## Publications

### REFEREED PUBLICATIONS IN JOURNALS

23. | Cunningham, A. B., H. Class, A. Ebigbo, R. Gerlach, A. J. Phillips, and J. Hommel Field-scale modeling of microbially induced calcite precipitation, Computational Geosciences, 23/2, pp. 399-414, 2019. Abstract The biogeochemical process known as microbially induced calcite precipitation (MICP) is being investigated for engineering and material science applications. To model MICP process behavior in porous media, computational simulators must couple flow, transport, and relevant biogeochemical reactions. Changes in media porosity and permeability due to biomass growth and calcite precipitation, as well as their effects on one another must be considered. A comprehensive Darcy-scale model has been developed by Ebigbo et al. (Water Resour. Res. 48(7), W07519, 2012) and Hommel et al. (Water Resour. Res. 51, 3695–3715, 2015) and validated at different scales of observation using laboratory experimental systems at the Center for Biofilm Engineering (CBE), Montana State University (MSU). This investigation clearly demonstrates that a close synergy between laboratory experimentation at different scales and corresponding simulation model development is necessary to advance MICP application to the field scale. Ultimately, model predictions of MICP sealing of a fractured sandstone formation, located 340.8 m below ground surface, were made and compared with corresponding field observations. Modeling MICP at the field scale poses special challenges, including choosing a reasonable model-domain size, initial and boundary conditions, and determining the initial distribution of porosity and permeability. In the presented study, model predictions of deposited calcite volume agree favorably with corresponding field observations of increased injection pressure during the MICP fracture sealing test in the field. Results indicate that the current status of our MICP model now allows its use for further subsurface engineering applications, including well-bore cement sealing and certain fracture-related applications in unconventional oil and gas production. / Download |

22. | Hobé, A., D. Vogler, M.P. Seybold, A. Ebigbo, R.R. Settgast, and M.O. Saar Estimating Fluid Flow Rates through Fracture Networks using Combinatorial Optimization, Advances in Water Resources, 122, pp. 85-97, 2018. Abstract To enable fast uncertainty quantification of fluid flow in a discrete fracture network (DFN), we present two approaches to quickly compute fluid flow in DFNs using combinatorial optimization algorithms. Specifically, the presented Hanan Shortest Path Maxflow (HSPM) and Intersection Shortest Path Maxflow (ISPM) methods translate DFN geometries and properties to a graph on which a max flow algorithm computes a combinatorial flow, from which an overall fluid flow rate is estimated using a shortest path decomposition of this flow. The two approaches are assessed by comparing their predictions with results from explicit numerical simulations of simple test cases as well as stochastic DFN realizations covering a range of fracture densities. Both methods have a high accuracy and very low computational cost, which can facilitate much-needed in-depth analyses of the propagation of uncertainty in fracture and fracture-network properties to fluid flow rates. / Download |

21. | Niederau, J., A. Ebigbo, G. Marquart, J. Arnold, and C. Clauser On the impact of spatially heterogenous permeability on free convection in the Perth Basin, Australia, Geothermics, 66, pp. 119-133, 2017. Abstract We study the impact of spatially heterogeneous permeability on the formation and shape of hydrothermal porous flow convection in the Yarragadee Aquifer by modelling three simulation scenarios, each with differing permeability distributions. In all scenarios, the southern part of the model is characterised by convection rolls, while the north is dominated by a stable region of decreased temperatures at depth due to hydraulic interaction with shallower aquifers. This suggests that reservoir structure is a first-order controlling factor for the formation of the free con- vective system. The convective system adjusts to the spatially heterogeneous permeability distribution, yielding locally different convection patterns. / Download |

20. | Büsing, H., C. Vogt, A. Ebigbo, and N. Kitzsch Numerical study on CO2 leakage detection using electrical streaming potential data, Water Resour. Res, 53, pp. 1-15, 2017. Abstract We study the feasibility of detecting carbon dioxide (CO2) movement in the overburden of a storage reservoir due to CO2 leakage through an abandoned well by self-potential (SP) measurements at the surface. This is achieved with three-dimensional numerical (SP) modeling of two-phase fluid flow and electrokinetic coupling between flow and streaming potential. We find that, in typical leakage scenarios, for leaky and/or injection wells with conductive metal casing, self-potential signals originating from injection can be identified at the surface. As the injection signal is also observed at the leaky well with metal casing, SP monitoring can be applied for detecting abandoned wells. However, leakage signals are much smaller than the injection signal and thus masked by the latter. We present three alternatives to overcome this problem: (i) simulate the streaming potential of the nonleaky scenario and subtract the result from the measured streaming potential data; (ii) exploit the symmetry of the injection signal by analyzing the potential difference of dipoles with the dipole center at the injection well; or (iii) measure SP during periods where the injection is interrupted. In our judgement, the most promising approach for detecting a real-world CO2 leakage is by combining methods (i) and (ii), because this would give the highest signal from the leakage and omit signals originating from the injection well. Consequently, we recommend SP as monitoring method for subsurface CO2 storage, especially because a leakage can be detected shortly after the injection started even before CO2 arrives at the leaky well. / Download |

19. | Ebigbo, A., P.A. Lang, A. Paluszny, and R.W. Zimmerman Inclusion-based effective medium models for the permeability of a 3D fractured rock mass, Transport in Porous Media, 113/1, pp. 137-158, 2016. Abstract Effective permeability is an essential parameter for describing fluid flow through fractured rock masses. This study investigates the ability of classical inclusion-based effective medium models (following the work of Sævik et al. in Transp Porous Media 100(1):115–142, 2013. doi:10.1007/s11242-013-0208-0) to predict this permeability, which depends on several geometric properties of the fractures/networks. This is achieved by comparison of various effective medium models, such as the symmetric and asymmetric self-consistent schemes, the differential scheme, and Maxwell’s method, with the results of explicit numerical simulations of mono- and poly-disperse isotropic fracture networks embedded in a permeable rock matrix. Comparisons are also made with the Hashin–Shtrikman bounds, Snow’s model, and Mourzenko’s heuristic model (Mourzenko et al. in Phys Rev E 84:036–307, 2011. doi:10.1103/PhysRevE.84.036307). This problem is characterised by two small parameters, the aspect ratio of the spheroidal fractures, α, and the ratio between matrix and fracture permeability, κ. Two different regimes can be identified, corresponding to α/κ<1 and α/κ>1. The lower the value of α/κ, the more significant is flow through the matrix. Due to differing flow patterns, the dependence of effective permeability on fracture density differs in the two regimes. When α/κ≫1, a distinct percolation threshold is observed, whereas for α/κ≪1, the matrix is sufficiently transmissive that such a transition is not observed. The self-consistent effective medium methods show good accuracy for both mono- and polydisperse isotropic fracture networks. Mourzenko’s equation is very accurate, particularly for monodisperse networks. Finally, it is shown that Snow’s model essentially coincides with the Hashin–Shtrikman upper bound. / Download |

18. | Ebigbo, A., J. Niederau, G. Marquart, I. Dini, M. Thorwart, W. Rabbel, R. Pechnig, R. Bertani, and C. Clauser Influence of depth, temperature, and structure of a crustal heat source on the geothermal reservoirs of Tuscany: numerical modelling and sensitivity study, Geothermal Energy, 4/5, 2016. Abstract Granitoid intrusions are the primary heat source of many deep geothermal reservoirs in Tuscany. The depth and shape of these plutons, characterised in this study by a prominent seismic reflector (the K horizon), may vary significantly within the spatial scale of interest. In an exploration field, simulations reveal the mechanisms by which such a heat source influences temperature distribution. A simple analysis quantifies the sensitivity of potentially measurable indicators (i.e. vertical temperature profiles and surface heat flow) to variations in depth, temperature, and shape of the heat source within given ranges of uncertainty. / Download |

17. | Seidler, R., K. Padalkina, H.M. Buecker, A. Ebigbo, M. Herty, G. Marquart, and J. Niederau Optimal experimental design for reservoir property estimates in geothermal exploration, Computational Geomechanics, 20/375, 2016. Abstract During geothermal reservoir development, drilling deep boreholes turns out to be extremely expensive and risky. Thus, it is of great importance to work out the details of suitable borehole locations in advance. Here, given a set of existing boreholes, we demonstrate how a sophisticated numerical technique called optimal experimental design helps to find a location of an additional exploratory borehole that reduces risk and, ultimately, saves cost. More precisely, the approach minimizes the uncertainty when deducing the effective permeability of a buried reservoir layer from a temperature profile measured in this exploratory borehole. In this paper, we (1) outline the mathematical formulation in terms of an optimization problem, (2) describe the numerical implementation involving various software components, and (3) apply the method to a 3D numerical simulation model representing a real geothermal reservoir in northern Italy. Our results show that optimal experimental design is conceptually and computationally feasible for industrial-scale applications. For the particular reservoir and the estimation of permeability from temperature, the optimal location of the additional borehole coincides with regions of high flow rates and large deviations from the mean temperature of the reservoir layer in question. Finally, the presentation shows that, methodologically, the optimization method can be generalized from estimating permeability to finding any other reservoir properties. / Download |

16. | Qin, C.-Z., S.M. Hassanizadeh, and A. Ebigbo Pore-scale network modeling of microbially induced calcium carbonate precipitation: Insight into scale dependence of biogeochemical reaction rates, Water Resources Research/52, 2016. Abstract The engineering of microbially induced calcium carbonate precipitation (MICP) has attracted much attention in a number of applications, such as sealing of CO2 leakage pathways, soil stabilization, and subsurface remediation of radionuclides and toxic metals. The goal of this work is to gain insight into pore-scale processes of MICP and scale dependence of biogeochemical reaction rates. This will help us develop efficient field-scale MICP models. In this work, we have developed a comprehensive pore-network model for MICP, with geochemical speciation calculated by the open-source PHREEQC module. A numerical pseudo-3-D micromodel as the computational domain was generated by a novel pore-network generation method. We modeled a three-stage process in the engineering of MICP including the growth of biofilm, the injection of calcium-rich medium, and the precipitation of calcium carbonate. A number of test cases were conducted to illustrate how calcite precipitation was influenced by different operating conditions. In addition, we studied the possibility of reducing the computational effort by simplifying geochemical calculations. Finally, the effect of mass transfer limitation of possible carbonate ions in a pore element on calcite precipitation was explored. / Download |

15. | Hommel, J., E. Lauchnor, R. Gerlach, A.B. Cunningham, A. Ebigbo, R. Helmig, and H. Class Investigating the influence of the initial biomass distribution and injection strategies on biofilm- mediated calcite precipitation in porous media, Transport in Porous Media, 2015. Abstract Attachment of bacteria in porous media is a complex mixture of processes resulting in the transfer and immobilization of suspended cells onto a solid surface within the porous medium. Quantifying the rate of attachment is difficult due to the many simultaneous processes possibly involved in attachment, including straining, sorption, and sedimentation, and the difficulties in measuring metabolically active cells attached to porous media. Preliminary experiments confirmed the difficulty associated with measuring active Sporosarcina pasteurii cells attached to porous media. However, attachment is a key process in applications of biofilm-mediated reactions in the subsurface such as microbially induced calcite precipitation. Independent of the exact processes involved, attachment determines both the distribution and the initial amount of attached biomass and as such the initial reaction rate. As direct experimental investigations are difficult, this study is limited to a numerical investigation of the effect of various initial biomass distributions and initial amounts of attached biomass. This is performed for various injection strategies, changing the injection rate as well as alternating between continuous and pulsed injections. The results of this study indicate that, for the selected scenarios, both the initial amount and the distribution of attached biomass have minor influence on the Ca2+ precipitation efficiency as well as the distribution of the precipitates compared to the influence of the injection strategy. The influence of the initial biomass distribution on the resulting final distribution of the precipitated calcite is limited, except for the continuous injection at intermediate injection rate. But even for this injection strategy, the Ca2+ precipitation efficiency shows no significant dependence on the initial biomass distribution. / Download |

14. | Hommel, J., E. Lauchnor, A. Phillips, R. Gerlach, A.B. Cunningham, R. Helmig, A. Ebigbo, and H. Class A revised model for microbially induced calcite precipitation: Improvements and new insights based on recent experiments, Water Resources Research, 51/5, pp. 3695-3715, 2015. Abstract The model for microbially induced calcite precipitation (MICP) published by Ebigbo et al. (2012) has been improved based on new insights obtained from experiments and model calibration. The challenge in constructing a predictive model for permeability reduction in the underground with MICP is the quantification of the complex interaction between flow, transport, biofilm growth, and reaction kinetics. New data from Lauchnor et al. (2015) on whole-cell ureolysis kinetics from batch experiments were incorporated into the model, which has allowed for a more precise quantification of the relevant parameters as well as a simplification of the reaction kinetics in the equations of the model. Further, the model has been calibrated objectively by inverse modeling using quasi-1D column experiments and a radial flow experiment. From the postprocessing of the inverse modeling, a comprehensive sensitivity analysis has been performed with focus on the model input parameters that were fitted in the course of the model calibration. It reveals that calcite precipitation and concentrations of math formula and math formula are particularly sensitive to parameters associated with the ureolysis rate and the attachment behavior of biomass. Based on the determined sensitivities and the ranges of values for the estimated parameters in the inversion, it is possible to identify focal areas where further research can have a high impact toward improving the understanding and engineering of MICP. / Download |

13. | Ebigbo, A., F. Golfier, and M. Quintard A coupled, pore-scale model for methanogenic microbial activity in underground hydrogen storage, Advances in Water Resources, 61, pp. 74-85, 2013. Abstract Underground hydrogen storage (UHS) as a means of energy storage is an efficient way of compensating for seasonal fluctuations in the availability of energy. One important factor which influences this technology is the activity of methanogenic microorganisms capable of utilising hydrogen and carbon dioxide for metabolism and leading to a change in the stored gas composition. A coupled, pore-scale model is presented which aids in the investigation of the mechanisms that govern the conversion of hydrogen to methane, i.e. advective hydrogen flow, its diffusion into microbial biofilms of multiple species, and its consumption within these biofilms. The model assumes that spherical grains are coated by a film of residual water and treats the biofilm development within each film in a quasi one-dimensional manner. A sample simulation using the presented model illustrates the biofilm growth process in these films as well as the competition between three different microbial species: methanogens, acetogens, and acetotrophs. / Download |

12. | Lange, T., M. Sauter, M. Heitfeld, K. Schetelig, K. Brosig, W. Jahnke, A. Kissinger, R. Helmig, A. Ebigbo, and H. Class Hydraulic fracturing in unconventional gas reservoirs: risks in the geological system, part 1, Environmental Earth Sciences, 70/8, pp. 3839-3853, 2013. Abstract Hydraulic fracturing of unconventional gas reservoirs rapidly developed especially in the USA to an industrial scale during the last decade. Potential adverse effects such as the deterioration of the quality of exploitable groundwater resources, areal footprints, or even the climate impact were not assessed. Because hydraulic fracturing has already been practised for a long time also in conventional reservoirs, the expansion into the unconventional domain was considered to be just a minor but not a technological step, with potential environmental risks. Thus, safety and environmental protection regulations were not critically developed or refined. Consequently, virtually no baseline conditions were documented before on-site applications as proof of evidence for the net effect of environmental impacts. Not only growing concerns in the general public, but also in the administrations in Germany promoted the commissioning of several expert opinions, evaluating safety, potential risks, and footprints of the technology in focus. The first two publications of the workgroup “Risks in the Geological System” of the independent “Information and Dialogue process on hydraulic fracturing” (commissioned by ExxonMobil Production Deutschland GmbH) comprises the strategy and approaches to identify and assess the potential risks of groundwater contamination of the exploitable groundwater system in the context of hydraulic fracturing operations in the Münsterland cretaceous basin and the Lower Saxony Basin, Germany. While being specific with respect to local geology and the estimation of effective hydraulic parameters, generalized concepts for the contamination risk assessment were developed. The work focuses on barrier effectiveness of different units of the overburden with respect to the migration of fracking fluids and methane, and considers fault zones as potential fluid pathway structures. / Download |

11. | Kissinger, A., R. Helmig, A. Ebigbo, H. Class, T. Lange, M. Sauter, M. Heitfeld, J. Klünker, and W. Jahnke Hydraulic fracturing in unconventional gas reservoirs: risks in the geological system, part 2, Environmental Earth Sciences, 70/8, pp. 3855-3873, 2013. Abstract Hydraulic fracturing is a method used for the production of unconventional gas resources. Huge amounts of so-called fracturing fluid (10,000–20,000 m3) are injected into a gas reservoir to create fractures in solid rock formations, upon which mobilised methane fills the pore space and the fracturing fluid is withdrawn. Hydraulic fracturing may pose a threat to groundwater resources if fracturing fluid or brine can migrate through fault zones into shallow aquifers. Diffuse methane emissions from the gas reservoir may not only contaminate shallow groundwater aquifers, but also escape into the atmosphere where methane acts as a greenhouse gas. The working group “Risks in the Geological System” as part of ExxonMobil’s hydrofracking dialogue and information dissemination processes was tasked with the assessment of possible hazards posed by migrating fluids as a result of hydraulic fracturing activities. In this work, several flow paths for fracturing fluid, brine and methane are identified and scenarios are set up to qualitatively estimate under what circumstances these fluids would leak into shallower layers. The parametrisation for potential hydraulic fracturing sites in North Rhine-Westphalia and Lower Saxony (both in Germany) is derived from literature using upper and lower bounds of hydraulic parameters. The results show that a significant fluid migration is only possible if a combination of several conservative assumptions is met by a scenario. / Download |

10. | Cunningham, A.B., E. Lauchnor, J. Eldring, E. Esposito, A.C. Mitchell, R. Gerlach, A.J. Phillips, A. Ebigbo, and L.H. Spangler Abandoned well CO2 leakage mitigation using biologically induced mineralization: current progress and future directions, Greenhouse Gases: Science & Technology, 3, pp. 40-49, 2013. Abstract Methods of mitigating leakage or re-plugging abandoned wells before exposure to CO2are of high potential interest to prevent leakage of CO2 injected for geologic carbon sequestration in depleted oil and gas reservoirs where large numbers of abandoned wells are often present. While CO2resistant cements and ultrafine cements are being developed, technologies that can be delivered via low viscosity fluids could have significant advantages including the ability to plug small aperture leaks such as fractures or delamination interfaces. Additionally there is the potential to plug rock formation pore space around the wellbore in particularly problematic situations. We are carrying out research on the use of microbial biofilms capable of inducing the precipitation of crystalline calcium carbonate using the process of ureolysis. This method has the potential to reduce well bore permeability, coat cement to reduce CO2–related corrosion, and lower the risk of unwanted upward CO2 migration. In this spotlight, we highlight research currently underway at the Center for Biofilm Engineering (CBE) at Montana State University (MSU) in the area of ureolytic biomineralization sealing for reducing CO2 leakage risk. This research program combines two novel core testing systems and a 3-dimensional simulation model to investigate biomineralization under both radial and axial flow conditions and at temperatures and pressures which permit CO2 to exist in the supercritical state. This combination of modelling and experimentation is ultimately aimed at developing and verifying biomineralization sealing technologies and strategies which can successfully be applied at the field scale for carbon capture and geological storage (CCGS) projects. © 2013 Society of Chemical Industry and John Wiley & Sons, Ltd / Download |

9. | Helmig, R., B. Flemisch, M. Wolff, A. Ebigbo, and H. Class Model coupling for multiphase flow in porous media, Advances in Water Resources, 51/7, pp. 52-66, 2013. Abstract Numerical models for flow and transport in porous media are valid for a particular set of processes, scales, levels of simplification and abstraction, grids etc. The coupling of two or more specialised models is a method of increasing the overall range of validity while keeping the computational costs relatively low. Several coupling concepts are reviewed in this article with a focus on the authors’ work in this field. The concepts are divided into temporal and spatial coupling concepts, of which the latter is subdivided into multi-process, multi-scale, multi-dimensional, and multi-compartment coupling strategies. Examples of applications for which these concepts can be relevant include groundwater protection and remediation, carbon dioxide storage, nuclear-waste disposal, soil dry-out and evaporation processes as well as fuel cells and technical filters. / Download |

8. | Ebigbo, A., A. Phillips, R. Gerlach, R. Helmig, A.B. Cunningham, and H. Class Darcy-scale modeling of microbially induced carbonate mineral precipitation in sand columns, Water Resources Research, 48/7, W07519, 2012. Abstract [1] This investigation focuses on the use of microbially induced calcium carbonate precipitation (MICP) to set up subsurface hydraulic barriers to potentially increase storage security near wellbores of CO2 storage sites. A numerical model is developed, capable of accounting for carbonate precipitation due to ureolytic bacterial activity as well as the flow of two fluid phases in the subsurface. The model is compared to experiments involving saturated flow through sand-packed columns to understand and optimize the processes involved as well as to validate the numerical model. It is then used to predict the effect of dense-phase CO2 and CO2-saturated water on carbonate precipitates in a porous medium. / Download |

7. | Skjaelaaen, I., A. Ebigbo, M. Espedal, and R. Helmig A model for transport of hydrogen sulfide in oil- and water-saturated porous media, Computing and Visualization in Science, 13/6, pp. 265-273, 2010. Abstract In several oilfields, reservoir souring by generation of hydrogen sulfide (H2S) occurs in secondary recovery during which seawater is injected into originally sweet reservoirs. At the production site, high concentrations of H2S can cause severe damage to both equipment and human personnel. Proper modeling of H2S concentration in produced fluids can be useful for decision-making during field development design. We present a model for the transport of H2S in an oil- and water-saturated, water-wet porous medium. The different retardation mechanisms for the H2S are described. For the adsorption of H2S to rock, we include two distinct phases of adsorption. In addition, we introduce a functional relationship between adsorption capacity and permeability. As H2S mixes with oil, fractions become immobile as part of the residual oil. Communicated by Gabriel wittum. This article is dedicated to the memory of our dear colleague, friend and mentor, Magne Espedal, who passed away during the preparation of this manuscript. / Download |

6. | Ebigbo, A., R. Helmig, A.B. Cunningham, H. Class, and R. Gerlach Modelling biofilm growth in the presence of carbon dioxide and water flow in the subsurface, Advances in Water Resources, 33/7, pp. 762-781, 2010. Abstract The concentration of greenhouse gases – particularly carbon dioxide (CO2) – in the atmosphere has been on the rise in the past decades. One of the methods which have been proposed to help reduce anthropogenic CO2 emissions is the capture of CO2from large, stationary point sources and storage in deep geological formations. The caprock is an impermeable geological layer which prevents the leakage of stored CO2, and its integrity is of utmost importance for storage security. Due to the high pressure build-up during injection, the caprock in the vicinity of the well is particularly at risk of fracturing. Biofilms could be used as biobarriers which help prevent the leakage of CO2 through the caprock in injection well vicinity by blocking leakage pathways. The biofilm could also protect well cement from corrosion by CO2-rich brine. The goal of this paper is to develop and test a numerical model which is capable of simulating the development of a biofilm in a CO2 storage reservoir. This involves the description of the growth of the biofilm, flow and transport in the geological formation, and the interaction between the biofilm and the flow processes. Important processes which are accounted for in the model include the effect of biofilm growth on the permeability of the formation, the hazardous effect of supercritical CO2 on suspended and attached bacteria, attachment and detachment of biomass, and two-phase fluid flow processes. The model is tested by comparing simulation results to experimental data. / Download |

5. | van Noorden, T.L., I.S. Pop, A. Ebigbo, and R. Helmig An upscaled model for biofilm growth in a thin strip, Water Resources Research, 46/6, W06505, 2010. Abstract [1] The focus of this paper is the derivation of an effective model for biofilm growth in a porous medium and its effect on fluid flow. The starting point is a pore-scale model in which the local geometry of the pore is represented as a thin strip. The model accounts for changes in pore volume due to biomass accumulation. As the ratio of the width of the strip to its length approaches zero, we apply a formal limiting argument to derive a one-dimensional upscaled (effective) model. For a better understanding of the terms and parameters involved in the equations derived here, we compare these equations to a well-known core-scale model from the literature. / Download |

4. | Class, H., A. Ebigbo, R Helmig, H.K. Dahle, J.M. Nordbotten, M.A. Celia, P. Audigane, M. Darcis, J. Ennis-King, and Y. Fan A benchmark study on problems related to CO2 storage in geologic formations , Computational Geosciences, 13/4, Sp. Iss. SI, pp. 409-434, 2009. Abstract This paper summarises the results of a benchmark study that compares a number of mathematical and numerical models applied to specific problems in the context of carbon dioxide (CO2) storage in geologic formations. The processes modelled comprise advective multi-phase flow, compositional effects due to dissolution of CO2 into the ambient brine and non-isothermal effects due to temperature gradients and the Joule–Thompson effect. The problems deal with leakage through a leaky well, methane recovery enhanced by CO2 injection and a reservoir-scale injection scenario into a heterogeneous formation. We give a description of the benchmark problems then briefly introduce the participating codes and finally present and discuss the results of the benchmark study. / Download |

3. | Kopp, A., A. Ebigbo, A. Bielinski, H. Class, and R. Helmig Numerical simulation of temperature changes caused by CO2 injection in geological reservoirs, AAPG Studies in Geology, 59/26, pp. 439-456, 2009. Abstract Injection of CO 2 into the subsurface for geological storage has an effect on the temperature of the storage formation and the CO 2 itself. Numerical investigations are an essential tool in describing the relevant processes that determine such changes and the impact they may have on the migration and the storage mechanisms of CO 2 in the subsurface. This chapter focuses on the numerical simulation of such thermal effects and their consequences. Simulating the temperature changes in a storage site can be of interest for temperature-based monitoring. Determining whether or how such thermal effects change the transport of CO 2 in the formation is important for the success of a CO 2 storage effort. In particular, this chapter examines a leakage scenario and how temperature changes could affect the leakage flow. The second part of the chapter presents results of a complex reservoir-scale simulation. The target formation forms an anticlinal structure at a depth of about 570-900 m (1870-2953 ft). Strong temperature effects can be expected because of the possible tran Numerical simulation of temperature changes caused by CO2 injection in geological reservoirs. Available from: https://www.researchgate.net/publication/230838241_Numerical_simulation_of_temperature_changes_caused_by_CO2_injection_in_geological_reservoirs [accessed May 3, 2017]. / Download |

2. | Ebigbo, A., H. Class, and R. Helmig CO2 leakage through an abandoned well: problem- oriented benchmarks, Computational Geosciences, 11/2, pp. 103-115, 2007. Abstract The efficiency and sustainability of carbon dioxide (CO2) storage in deep geological formations crucially depends on the integrity of the overlying cap-rocks. Existing oil and gas wells, which penetrate the formations, are potential leakage pathways. This problem has been discussed in the literature, and a number of investigations using semi-analytical mathematical approaches have been carried out by other authors to quantify leakage rates. The semi-analytical results are based on a number of simplifying assumptions. Thus, it is of great interest to assess the influence of these assumptions. We use a numerical model to compare the results with those of the semi-analytical model. Then we ease the simplifying restrictions and include more complex thermodynamic processes including sub- and supercritical fluid properties of CO2 and non-isothermal as well as compositional effects. The aim is to set up problem-oriented benchmark examples that allow a comparison of different modeling approaches to the problem of CO2 leakage. / Download |

1. | Class, H., A. Bielinski, R. Helmig, A. Kopp, and A. Ebigbo Numerical simulation of CO2 storage in geological formations, Chemie Ingenieur Technik, 78/4, pp. 445-452, 2006. Abstract Die Speicherung von Kohlendioxid in geologischen Formationen wird derzeit als ein möglicher Beitrag zur Reduktion von Treibhausgaskonzentrationen in der Atmosphäre diskutiert und untersucht. Ein wichtiges Werkzeug begleitend zu experimentellen Untersuchungsmethoden ist die numerische Simulation. Hier wird ein Einblick in die physikalischen bzw. thermodynamischen Vorgänge im Untergrund während und nach einer Injektion von Kohlendioxid gegeben sowie deren modellkonzeptionelle Beschreibung durch nichtisotherme Mehrphasensysteme dargestellt. Die elementaren Schritte der physikalisch/mathematischen Modellbildung und numerischen Lösung der daraus entstehenden Gleichungssysteme werden erklärt. Anhand von Simulationsbeispielen werden dominierende Prozesse während und nach einer Injektion diskutiert. / Download |

### PROCEEDINGS REFEREED

2. | Hommel, J., A. Ebigbo, R. Gerlach, A. B. Cunningham, R. Helmig, and H. Class Finding a Balance between Accuracy and Effort For Modeling Biomineralization, Energy Procedia, 97, pp. 379-386, 2016. / Download |

1. | Hommel, J., A. B. Cunningham, R. Helmig, A. Ebigbo, and H. Class Numerical Investigation of Microbially Induced Calcite Precipitation as a Leakage Mitigation Technology, Energy Procedia, 40, pp. 392-397, 2013. / Download |

### THESES

2. | Ebigbo, A. Modelling of biofilm growth and its influence on CO2 and water (two-phase) flow in porous media, Dissertation University of Stuttgart, 131 pp., 2009. Abstract Bacterial biofilms are groups of microbial cells attached to surfaces and to each other. Cells in a biofilm are protected from adverse external conditions. In natural environments, this attached mode of growth is more successful than the suspended mode, and a major portion of microbial activity takes place at surfaces. In porous media, biofilms are used as bioreactors (e.g, in wastewater treatment) and as biobarriers (e.g., in enhanced oil recovery). They are also used in the containment and degradation of contaminants in groundwater aquifers. It has been proposed that biofilms be used as biobarriers for the mitigation of carbon dioxide (CO2) leakage from a geological storage reservoir. The concentration of greenhouse gases — particularly carbon dioxide (CO2) — in the atmosphere has been on the rise in the past decades. One of the methods which have been proposed to help reduce anthropogenic CO2 emissions is the capture of CO2 from large, stationary point sources and storage in deep geological formations. The caprock is an impermeable geological layer which prevents the leakage of stored CO2, and its integrity is of utmost importance for storage security. As mentioned above, biofilms could be used as biobarriers which help prevent the leakage of CO2 through the caprock in injection well vicinity. Due to the high pressure build-up during injection, the caprock in the vicinity of the well is particularly at risk of fracturing. The biofilm could also protect well cement from corrosion by CO2-rich brine. The goal of this work is to develop and test a numerical model which is capable of simulating the development of a biofilm in a CO2 storage reservoir. This involves the description of the growth of the biofilm, flow and transport in the geological formation, and the interaction between the biofilm and the flow processes. Important processes which are accounted for in the model include the effect of biofilm growth on the permeability of the formation, the hazardous effect of supercritical CO2 on suspended and attached bacteria, attachment and detachment of biomass, and two-phase fluid flow processes. The partial differential equations which describe the system are discretised in space with a vertex-centered finite volume method, and an implicit Euler scheme is used for time discretisation. The model is tested by comparing simulation results to experimental data. In a test case simulation, the model predicts the extent of biomass accumulation near an injection well and its effect on the permeability of the formation. The simulations show that the biobarrier is only effective for a limited amount of time. Regular injection of nutrients would be necessary to sustain the biofilm. In future work, the model could be extended to account for the active precipitation of minerals by the biofilm which would lead to a more enduring barrier. The model also needs to be extended to account for more than one growth-limiting factor. This would allow for the simulation of injection strategies which aim at growing a biofilm at some distance from the injection well. Biofilme, die in einem porösen Medium wachsen, blockieren Poren und verändern dabei die Eigenschaften des porösen Mediums. Diese veränderten Eigenschaften werden bei der biologischen Filtration (z. B. bei der Abwasserbehandlung), bei der biologischen Altlastensanierung (z. B. für die Erstellung hydraulischer Barrieren) und bei anderen Fragestellungen auf diesem Gebiet genutzt. Eine hydraulische Barriere biologischen Ursprungs könnte z. B. auch in einer geologischen Kohlendioxid-Lagerstätte eingesetzt werden, um das Entweichen von CO2 zu verhindern. CO2 ist das derzeit für am Wichtigsten erachtete anthropogene Treibhausgas. Die globale Erderwärmung wird demnach sehr stark durch die in den letzten Jahrzehnten stattfindende Anreicherung von anthropogenen Treibhausgasen in der Atmosphäre mitverursacht. Die Freisetzung von CO2 kann mit Hilfe effizienterer Technologien und alternativer Energiequellen reduziert werden. CO2-Emissionen können aber auch reduziert werden, indem man CO2 aus Kraftwerksabgasen abscheidet und in tiefen geologischen Formationen speichert. Bei den physikalischen Bedingungen, die in diesen unterirdischen Lagerstätten herrschen, liegt CO2 im überkritischen Zustand vor, gekennzeichnet durch eine hohe Dichte und geringe Viskosität. Diese Lagerstätten enthalten oft salzhaltiges Wasser, das dichter ist als CO2. Eine möglichst undurchlässige geologische Deckschicht verhindert das Aufsteigen des leichteren CO2 an die Erdoberfläche. Jedoch müssen, z. B. im Rahmen von Risikostudien, mögliche Störungen oder Risse in dieser Deckschicht betrachtet werden, die zu einem Entweichen des CO2 führen könnten. Die Deckschicht in der Nähe eines CO2-Injektionsbrunnens ist besonders gefährdet. Der hohe Druckanstieg während der ersten Injektionsphase, Zementkorrosion am Brunnen aufgrund des CO2-reichen Formationswassers und eventuelle Beschädigungen der Deckschicht während der Erstellung des Bohrlochs sind als mögliche Ursachen für gestörte Deckschichten zu nennen. Biobarrieren könnten verwendet werden, um solche Risiken zu minimieren, z. B. indem sie Risse in der Deckschicht abdichten oder den Bohrlochzement vor Korrosion schützen. Eine Biobarriere kann aus einem Biofilm selbst bestehen, aber auch aus vom Biofilm begüngstigten mineralischen Ablagerungen. Die vorliegende Arbeit behandelt im Wesentlichen die Entwicklung eines numerischen Modells, um die Anreicherung von mikrobieller Biomasse im Untergrund simulieren zu können. Das entwickelte Modell soll in der Lage sein, das Abdichten der beschädigten geologischen Deckschicht einer unterirdischen Kohlendioxid-Lagerstätte mit Hilfe von Biofilmen zu simulieren. Dafür müssen einerseits Strömungsprozesse und andererseits auch die mikrobielle Aktivität sowie die Interaktion dieser Vorgänge in porösen Medien richtig beschrieben werden. Die Anreicherung von Bakterien in einem porösen Medium beeinflusst die hydraulischen Eigenschaften des Mediums und als Folge davon auch die darin stattfindende Strömung. Im Gegenzug bestimmt die Strömung den Transport der mikrobiellen Nährstoffe und damit auch die Verteilung mikrobieller Wachstumsraten. Dementsprechend ist die richtige Beschreibung der Wechselwirkung zwischen Strömung und mikrobiellen Prozessen eine wesentliche Herausforderung in der Modellbildung. / Download |

1. | Ebigbo, A. Thermal Effects of Carbon Dioxide Sequestration in the Subsurface. Diplomarbeit, Institut für Wasserbau, MSc Thesis University of Stuttgart, 57 pp., 2005. Abstract In order to secure long-term storage of CO 2 in the subsurface, one has to be able to model the mass transport of CO2. The physical properties of CO2 have a very strong influence on the mass transport. These physical properties, in turn are dependent on temperature (and pressure). Hence, the modelling of heat transport in the subsurface during CO2 injection should be part of mass transport. In the simple set-up examined in this thesis, the Joule-Thompson cooling (as a result of the pressure drop) plays an important role near the injection point. The amount of cooling at the injection point depends on the temperature of the injected CO2 (and of course on the pressure difference that brought about the cooling). The temperature at that region achieves a stable value with time. The impermeable layer should be deeper than the depth at which the CO2 becomes gaseous as it is easier to trap less mobile supercritical or liquid CO2. / Download |

show/hide list of publications

### REFEREED PUBLICATIONS IN JOURNALS

23. | Cunningham, A. B., H. Class, A. Ebigbo, R. Gerlach, A. J. Phillips, and J. Hommel Field-scale modeling of microbially induced calcite precipitation, Computational Geosciences, 23/2, pp. 399-414, 2019. Abstract The biogeochemical process known as microbially induced calcite precipitation (MICP) is being investigated for engineering and material science applications. To model MICP process behavior in porous media, computational simulators must couple flow, transport, and relevant biogeochemical reactions. Changes in media porosity and permeability due to biomass growth and calcite precipitation, as well as their effects on one another must be considered. A comprehensive Darcy-scale model has been developed by Ebigbo et al. (Water Resour. Res. 48(7), W07519, 2012) and Hommel et al. (Water Resour. Res. 51, 3695–3715, 2015) and validated at different scales of observation using laboratory experimental systems at the Center for Biofilm Engineering (CBE), Montana State University (MSU). This investigation clearly demonstrates that a close synergy between laboratory experimentation at different scales and corresponding simulation model development is necessary to advance MICP application to the field scale. Ultimately, model predictions of MICP sealing of a fractured sandstone formation, located 340.8 m below ground surface, were made and compared with corresponding field observations. Modeling MICP at the field scale poses special challenges, including choosing a reasonable model-domain size, initial and boundary conditions, and determining the initial distribution of porosity and permeability. In the presented study, model predictions of deposited calcite volume agree favorably with corresponding field observations of increased injection pressure during the MICP fracture sealing test in the field. Results indicate that the current status of our MICP model now allows its use for further subsurface engineering applications, including well-bore cement sealing and certain fracture-related applications in unconventional oil and gas production. / Download |

22. | Hobé, A., D. Vogler, M.P. Seybold, A. Ebigbo, R.R. Settgast, and M.O. Saar Estimating Fluid Flow Rates through Fracture Networks using Combinatorial Optimization, Advances in Water Resources, 122, pp. 85-97, 2018. Abstract To enable fast uncertainty quantification of fluid flow in a discrete fracture network (DFN), we present two approaches to quickly compute fluid flow in DFNs using combinatorial optimization algorithms. Specifically, the presented Hanan Shortest Path Maxflow (HSPM) and Intersection Shortest Path Maxflow (ISPM) methods translate DFN geometries and properties to a graph on which a max flow algorithm computes a combinatorial flow, from which an overall fluid flow rate is estimated using a shortest path decomposition of this flow. The two approaches are assessed by comparing their predictions with results from explicit numerical simulations of simple test cases as well as stochastic DFN realizations covering a range of fracture densities. Both methods have a high accuracy and very low computational cost, which can facilitate much-needed in-depth analyses of the propagation of uncertainty in fracture and fracture-network properties to fluid flow rates. / Download |

21. | Niederau, J., A. Ebigbo, G. Marquart, J. Arnold, and C. Clauser On the impact of spatially heterogenous permeability on free convection in the Perth Basin, Australia, Geothermics, 66, pp. 119-133, 2017. Abstract We study the impact of spatially heterogeneous permeability on the formation and shape of hydrothermal porous flow convection in the Yarragadee Aquifer by modelling three simulation scenarios, each with differing permeability distributions. In all scenarios, the southern part of the model is characterised by convection rolls, while the north is dominated by a stable region of decreased temperatures at depth due to hydraulic interaction with shallower aquifers. This suggests that reservoir structure is a first-order controlling factor for the formation of the free con- vective system. The convective system adjusts to the spatially heterogeneous permeability distribution, yielding locally different convection patterns. / Download |

20. | Büsing, H., C. Vogt, A. Ebigbo, and N. Kitzsch Numerical study on CO2 leakage detection using electrical streaming potential data, Water Resour. Res, 53, pp. 1-15, 2017. Abstract We study the feasibility of detecting carbon dioxide (CO2) movement in the overburden of a storage reservoir due to CO2 leakage through an abandoned well by self-potential (SP) measurements at the surface. This is achieved with three-dimensional numerical (SP) modeling of two-phase fluid flow and electrokinetic coupling between flow and streaming potential. We find that, in typical leakage scenarios, for leaky and/or injection wells with conductive metal casing, self-potential signals originating from injection can be identified at the surface. As the injection signal is also observed at the leaky well with metal casing, SP monitoring can be applied for detecting abandoned wells. However, leakage signals are much smaller than the injection signal and thus masked by the latter. We present three alternatives to overcome this problem: (i) simulate the streaming potential of the nonleaky scenario and subtract the result from the measured streaming potential data; (ii) exploit the symmetry of the injection signal by analyzing the potential difference of dipoles with the dipole center at the injection well; or (iii) measure SP during periods where the injection is interrupted. In our judgement, the most promising approach for detecting a real-world CO2 leakage is by combining methods (i) and (ii), because this would give the highest signal from the leakage and omit signals originating from the injection well. Consequently, we recommend SP as monitoring method for subsurface CO2 storage, especially because a leakage can be detected shortly after the injection started even before CO2 arrives at the leaky well. / Download |

19. | Ebigbo, A., P.A. Lang, A. Paluszny, and R.W. Zimmerman Inclusion-based effective medium models for the permeability of a 3D fractured rock mass, Transport in Porous Media, 113/1, pp. 137-158, 2016. Abstract Effective permeability is an essential parameter for describing fluid flow through fractured rock masses. This study investigates the ability of classical inclusion-based effective medium models (following the work of Sævik et al. in Transp Porous Media 100(1):115–142, 2013. doi:10.1007/s11242-013-0208-0) to predict this permeability, which depends on several geometric properties of the fractures/networks. This is achieved by comparison of various effective medium models, such as the symmetric and asymmetric self-consistent schemes, the differential scheme, and Maxwell’s method, with the results of explicit numerical simulations of mono- and poly-disperse isotropic fracture networks embedded in a permeable rock matrix. Comparisons are also made with the Hashin–Shtrikman bounds, Snow’s model, and Mourzenko’s heuristic model (Mourzenko et al. in Phys Rev E 84:036–307, 2011. doi:10.1103/PhysRevE.84.036307). This problem is characterised by two small parameters, the aspect ratio of the spheroidal fractures, α, and the ratio between matrix and fracture permeability, κ. Two different regimes can be identified, corresponding to α/κ<1 and α/κ>1. The lower the value of α/κ, the more significant is flow through the matrix. Due to differing flow patterns, the dependence of effective permeability on fracture density differs in the two regimes. When α/κ≫1, a distinct percolation threshold is observed, whereas for α/κ≪1, the matrix is sufficiently transmissive that such a transition is not observed. The self-consistent effective medium methods show good accuracy for both mono- and polydisperse isotropic fracture networks. Mourzenko’s equation is very accurate, particularly for monodisperse networks. Finally, it is shown that Snow’s model essentially coincides with the Hashin–Shtrikman upper bound. / Download |

18. | Ebigbo, A., J. Niederau, G. Marquart, I. Dini, M. Thorwart, W. Rabbel, R. Pechnig, R. Bertani, and C. Clauser Influence of depth, temperature, and structure of a crustal heat source on the geothermal reservoirs of Tuscany: numerical modelling and sensitivity study, Geothermal Energy, 4/5, 2016. Abstract Granitoid intrusions are the primary heat source of many deep geothermal reservoirs in Tuscany. The depth and shape of these plutons, characterised in this study by a prominent seismic reflector (the K horizon), may vary significantly within the spatial scale of interest. In an exploration field, simulations reveal the mechanisms by which such a heat source influences temperature distribution. A simple analysis quantifies the sensitivity of potentially measurable indicators (i.e. vertical temperature profiles and surface heat flow) to variations in depth, temperature, and shape of the heat source within given ranges of uncertainty. / Download |

17. | Seidler, R., K. Padalkina, H.M. Buecker, A. Ebigbo, M. Herty, G. Marquart, and J. Niederau Optimal experimental design for reservoir property estimates in geothermal exploration, Computational Geomechanics, 20/375, 2016. Abstract During geothermal reservoir development, drilling deep boreholes turns out to be extremely expensive and risky. Thus, it is of great importance to work out the details of suitable borehole locations in advance. Here, given a set of existing boreholes, we demonstrate how a sophisticated numerical technique called optimal experimental design helps to find a location of an additional exploratory borehole that reduces risk and, ultimately, saves cost. More precisely, the approach minimizes the uncertainty when deducing the effective permeability of a buried reservoir layer from a temperature profile measured in this exploratory borehole. In this paper, we (1) outline the mathematical formulation in terms of an optimization problem, (2) describe the numerical implementation involving various software components, and (3) apply the method to a 3D numerical simulation model representing a real geothermal reservoir in northern Italy. Our results show that optimal experimental design is conceptually and computationally feasible for industrial-scale applications. For the particular reservoir and the estimation of permeability from temperature, the optimal location of the additional borehole coincides with regions of high flow rates and large deviations from the mean temperature of the reservoir layer in question. Finally, the presentation shows that, methodologically, the optimization method can be generalized from estimating permeability to finding any other reservoir properties. / Download |

16. | Qin, C.-Z., S.M. Hassanizadeh, and A. Ebigbo Pore-scale network modeling of microbially induced calcium carbonate precipitation: Insight into scale dependence of biogeochemical reaction rates, Water Resources Research/52, 2016. Abstract The engineering of microbially induced calcium carbonate precipitation (MICP) has attracted much attention in a number of applications, such as sealing of CO2 leakage pathways, soil stabilization, and subsurface remediation of radionuclides and toxic metals. The goal of this work is to gain insight into pore-scale processes of MICP and scale dependence of biogeochemical reaction rates. This will help us develop efficient field-scale MICP models. In this work, we have developed a comprehensive pore-network model for MICP, with geochemical speciation calculated by the open-source PHREEQC module. A numerical pseudo-3-D micromodel as the computational domain was generated by a novel pore-network generation method. We modeled a three-stage process in the engineering of MICP including the growth of biofilm, the injection of calcium-rich medium, and the precipitation of calcium carbonate. A number of test cases were conducted to illustrate how calcite precipitation was influenced by different operating conditions. In addition, we studied the possibility of reducing the computational effort by simplifying geochemical calculations. Finally, the effect of mass transfer limitation of possible carbonate ions in a pore element on calcite precipitation was explored. / Download |

15. | Hommel, J., E. Lauchnor, R. Gerlach, A.B. Cunningham, A. Ebigbo, R. Helmig, and H. Class Investigating the influence of the initial biomass distribution and injection strategies on biofilm- mediated calcite precipitation in porous media, Transport in Porous Media, 2015. Abstract Attachment of bacteria in porous media is a complex mixture of processes resulting in the transfer and immobilization of suspended cells onto a solid surface within the porous medium. Quantifying the rate of attachment is difficult due to the many simultaneous processes possibly involved in attachment, including straining, sorption, and sedimentation, and the difficulties in measuring metabolically active cells attached to porous media. Preliminary experiments confirmed the difficulty associated with measuring active Sporosarcina pasteurii cells attached to porous media. However, attachment is a key process in applications of biofilm-mediated reactions in the subsurface such as microbially induced calcite precipitation. Independent of the exact processes involved, attachment determines both the distribution and the initial amount of attached biomass and as such the initial reaction rate. As direct experimental investigations are difficult, this study is limited to a numerical investigation of the effect of various initial biomass distributions and initial amounts of attached biomass. This is performed for various injection strategies, changing the injection rate as well as alternating between continuous and pulsed injections. The results of this study indicate that, for the selected scenarios, both the initial amount and the distribution of attached biomass have minor influence on the Ca2+ precipitation efficiency as well as the distribution of the precipitates compared to the influence of the injection strategy. The influence of the initial biomass distribution on the resulting final distribution of the precipitated calcite is limited, except for the continuous injection at intermediate injection rate. But even for this injection strategy, the Ca2+ precipitation efficiency shows no significant dependence on the initial biomass distribution. / Download |

14. | Hommel, J., E. Lauchnor, A. Phillips, R. Gerlach, A.B. Cunningham, R. Helmig, A. Ebigbo, and H. Class A revised model for microbially induced calcite precipitation: Improvements and new insights based on recent experiments, Water Resources Research, 51/5, pp. 3695-3715, 2015. Abstract The model for microbially induced calcite precipitation (MICP) published by Ebigbo et al. (2012) has been improved based on new insights obtained from experiments and model calibration. The challenge in constructing a predictive model for permeability reduction in the underground with MICP is the quantification of the complex interaction between flow, transport, biofilm growth, and reaction kinetics. New data from Lauchnor et al. (2015) on whole-cell ureolysis kinetics from batch experiments were incorporated into the model, which has allowed for a more precise quantification of the relevant parameters as well as a simplification of the reaction kinetics in the equations of the model. Further, the model has been calibrated objectively by inverse modeling using quasi-1D column experiments and a radial flow experiment. From the postprocessing of the inverse modeling, a comprehensive sensitivity analysis has been performed with focus on the model input parameters that were fitted in the course of the model calibration. It reveals that calcite precipitation and concentrations of math formula and math formula are particularly sensitive to parameters associated with the ureolysis rate and the attachment behavior of biomass. Based on the determined sensitivities and the ranges of values for the estimated parameters in the inversion, it is possible to identify focal areas where further research can have a high impact toward improving the understanding and engineering of MICP. / Download |

13. | Ebigbo, A., F. Golfier, and M. Quintard A coupled, pore-scale model for methanogenic microbial activity in underground hydrogen storage, Advances in Water Resources, 61, pp. 74-85, 2013. Abstract Underground hydrogen storage (UHS) as a means of energy storage is an efficient way of compensating for seasonal fluctuations in the availability of energy. One important factor which influences this technology is the activity of methanogenic microorganisms capable of utilising hydrogen and carbon dioxide for metabolism and leading to a change in the stored gas composition. A coupled, pore-scale model is presented which aids in the investigation of the mechanisms that govern the conversion of hydrogen to methane, i.e. advective hydrogen flow, its diffusion into microbial biofilms of multiple species, and its consumption within these biofilms. The model assumes that spherical grains are coated by a film of residual water and treats the biofilm development within each film in a quasi one-dimensional manner. A sample simulation using the presented model illustrates the biofilm growth process in these films as well as the competition between three different microbial species: methanogens, acetogens, and acetotrophs. / Download |

12. | Lange, T., M. Sauter, M. Heitfeld, K. Schetelig, K. Brosig, W. Jahnke, A. Kissinger, R. Helmig, A. Ebigbo, and H. Class Hydraulic fracturing in unconventional gas reservoirs: risks in the geological system, part 1, Environmental Earth Sciences, 70/8, pp. 3839-3853, 2013. Abstract Hydraulic fracturing of unconventional gas reservoirs rapidly developed especially in the USA to an industrial scale during the last decade. Potential adverse effects such as the deterioration of the quality of exploitable groundwater resources, areal footprints, or even the climate impact were not assessed. Because hydraulic fracturing has already been practised for a long time also in conventional reservoirs, the expansion into the unconventional domain was considered to be just a minor but not a technological step, with potential environmental risks. Thus, safety and environmental protection regulations were not critically developed or refined. Consequently, virtually no baseline conditions were documented before on-site applications as proof of evidence for the net effect of environmental impacts. Not only growing concerns in the general public, but also in the administrations in Germany promoted the commissioning of several expert opinions, evaluating safety, potential risks, and footprints of the technology in focus. The first two publications of the workgroup “Risks in the Geological System” of the independent “Information and Dialogue process on hydraulic fracturing” (commissioned by ExxonMobil Production Deutschland GmbH) comprises the strategy and approaches to identify and assess the potential risks of groundwater contamination of the exploitable groundwater system in the context of hydraulic fracturing operations in the Münsterland cretaceous basin and the Lower Saxony Basin, Germany. While being specific with respect to local geology and the estimation of effective hydraulic parameters, generalized concepts for the contamination risk assessment were developed. The work focuses on barrier effectiveness of different units of the overburden with respect to the migration of fracking fluids and methane, and considers fault zones as potential fluid pathway structures. / Download |

11. | Kissinger, A., R. Helmig, A. Ebigbo, H. Class, T. Lange, M. Sauter, M. Heitfeld, J. Klünker, and W. Jahnke Hydraulic fracturing in unconventional gas reservoirs: risks in the geological system, part 2, Environmental Earth Sciences, 70/8, pp. 3855-3873, 2013. Abstract Hydraulic fracturing is a method used for the production of unconventional gas resources. Huge amounts of so-called fracturing fluid (10,000–20,000 m3) are injected into a gas reservoir to create fractures in solid rock formations, upon which mobilised methane fills the pore space and the fracturing fluid is withdrawn. Hydraulic fracturing may pose a threat to groundwater resources if fracturing fluid or brine can migrate through fault zones into shallow aquifers. Diffuse methane emissions from the gas reservoir may not only contaminate shallow groundwater aquifers, but also escape into the atmosphere where methane acts as a greenhouse gas. The working group “Risks in the Geological System” as part of ExxonMobil’s hydrofracking dialogue and information dissemination processes was tasked with the assessment of possible hazards posed by migrating fluids as a result of hydraulic fracturing activities. In this work, several flow paths for fracturing fluid, brine and methane are identified and scenarios are set up to qualitatively estimate under what circumstances these fluids would leak into shallower layers. The parametrisation for potential hydraulic fracturing sites in North Rhine-Westphalia and Lower Saxony (both in Germany) is derived from literature using upper and lower bounds of hydraulic parameters. The results show that a significant fluid migration is only possible if a combination of several conservative assumptions is met by a scenario. / Download |

10. | Cunningham, A.B., E. Lauchnor, J. Eldring, E. Esposito, A.C. Mitchell, R. Gerlach, A.J. Phillips, A. Ebigbo, and L.H. Spangler Abandoned well CO2 leakage mitigation using biologically induced mineralization: current progress and future directions, Greenhouse Gases: Science & Technology, 3, pp. 40-49, 2013. Abstract Methods of mitigating leakage or re-plugging abandoned wells before exposure to CO2are of high potential interest to prevent leakage of CO2 injected for geologic carbon sequestration in depleted oil and gas reservoirs where large numbers of abandoned wells are often present. While CO2resistant cements and ultrafine cements are being developed, technologies that can be delivered via low viscosity fluids could have significant advantages including the ability to plug small aperture leaks such as fractures or delamination interfaces. Additionally there is the potential to plug rock formation pore space around the wellbore in particularly problematic situations. We are carrying out research on the use of microbial biofilms capable of inducing the precipitation of crystalline calcium carbonate using the process of ureolysis. This method has the potential to reduce well bore permeability, coat cement to reduce CO2–related corrosion, and lower the risk of unwanted upward CO2 migration. In this spotlight, we highlight research currently underway at the Center for Biofilm Engineering (CBE) at Montana State University (MSU) in the area of ureolytic biomineralization sealing for reducing CO2 leakage risk. This research program combines two novel core testing systems and a 3-dimensional simulation model to investigate biomineralization under both radial and axial flow conditions and at temperatures and pressures which permit CO2 to exist in the supercritical state. This combination of modelling and experimentation is ultimately aimed at developing and verifying biomineralization sealing technologies and strategies which can successfully be applied at the field scale for carbon capture and geological storage (CCGS) projects. © 2013 Society of Chemical Industry and John Wiley & Sons, Ltd / Download |

9. | Helmig, R., B. Flemisch, M. Wolff, A. Ebigbo, and H. Class Model coupling for multiphase flow in porous media, Advances in Water Resources, 51/7, pp. 52-66, 2013. Abstract Numerical models for flow and transport in porous media are valid for a particular set of processes, scales, levels of simplification and abstraction, grids etc. The coupling of two or more specialised models is a method of increasing the overall range of validity while keeping the computational costs relatively low. Several coupling concepts are reviewed in this article with a focus on the authors’ work in this field. The concepts are divided into temporal and spatial coupling concepts, of which the latter is subdivided into multi-process, multi-scale, multi-dimensional, and multi-compartment coupling strategies. Examples of applications for which these concepts can be relevant include groundwater protection and remediation, carbon dioxide storage, nuclear-waste disposal, soil dry-out and evaporation processes as well as fuel cells and technical filters. / Download |

8. | Ebigbo, A., A. Phillips, R. Gerlach, R. Helmig, A.B. Cunningham, and H. Class Darcy-scale modeling of microbially induced carbonate mineral precipitation in sand columns, Water Resources Research, 48/7, W07519, 2012. Abstract [1] This investigation focuses on the use of microbially induced calcium carbonate precipitation (MICP) to set up subsurface hydraulic barriers to potentially increase storage security near wellbores of CO2 storage sites. A numerical model is developed, capable of accounting for carbonate precipitation due to ureolytic bacterial activity as well as the flow of two fluid phases in the subsurface. The model is compared to experiments involving saturated flow through sand-packed columns to understand and optimize the processes involved as well as to validate the numerical model. It is then used to predict the effect of dense-phase CO2 and CO2-saturated water on carbonate precipitates in a porous medium. / Download |

7. | Skjaelaaen, I., A. Ebigbo, M. Espedal, and R. Helmig A model for transport of hydrogen sulfide in oil- and water-saturated porous media, Computing and Visualization in Science, 13/6, pp. 265-273, 2010. Abstract In several oilfields, reservoir souring by generation of hydrogen sulfide (H2S) occurs in secondary recovery during which seawater is injected into originally sweet reservoirs. At the production site, high concentrations of H2S can cause severe damage to both equipment and human personnel. Proper modeling of H2S concentration in produced fluids can be useful for decision-making during field development design. We present a model for the transport of H2S in an oil- and water-saturated, water-wet porous medium. The different retardation mechanisms for the H2S are described. For the adsorption of H2S to rock, we include two distinct phases of adsorption. In addition, we introduce a functional relationship between adsorption capacity and permeability. As H2S mixes with oil, fractions become immobile as part of the residual oil. Communicated by Gabriel wittum. This article is dedicated to the memory of our dear colleague, friend and mentor, Magne Espedal, who passed away during the preparation of this manuscript. / Download |

6. | Ebigbo, A., R. Helmig, A.B. Cunningham, H. Class, and R. Gerlach Modelling biofilm growth in the presence of carbon dioxide and water flow in the subsurface, Advances in Water Resources, 33/7, pp. 762-781, 2010. Abstract The concentration of greenhouse gases – particularly carbon dioxide (CO2) – in the atmosphere has been on the rise in the past decades. One of the methods which have been proposed to help reduce anthropogenic CO2 emissions is the capture of CO2from large, stationary point sources and storage in deep geological formations. The caprock is an impermeable geological layer which prevents the leakage of stored CO2, and its integrity is of utmost importance for storage security. Due to the high pressure build-up during injection, the caprock in the vicinity of the well is particularly at risk of fracturing. Biofilms could be used as biobarriers which help prevent the leakage of CO2 through the caprock in injection well vicinity by blocking leakage pathways. The biofilm could also protect well cement from corrosion by CO2-rich brine. The goal of this paper is to develop and test a numerical model which is capable of simulating the development of a biofilm in a CO2 storage reservoir. This involves the description of the growth of the biofilm, flow and transport in the geological formation, and the interaction between the biofilm and the flow processes. Important processes which are accounted for in the model include the effect of biofilm growth on the permeability of the formation, the hazardous effect of supercritical CO2 on suspended and attached bacteria, attachment and detachment of biomass, and two-phase fluid flow processes. The model is tested by comparing simulation results to experimental data. / Download |

5. | van Noorden, T.L., I.S. Pop, A. Ebigbo, and R. Helmig An upscaled model for biofilm growth in a thin strip, Water Resources Research, 46/6, W06505, 2010. Abstract [1] The focus of this paper is the derivation of an effective model for biofilm growth in a porous medium and its effect on fluid flow. The starting point is a pore-scale model in which the local geometry of the pore is represented as a thin strip. The model accounts for changes in pore volume due to biomass accumulation. As the ratio of the width of the strip to its length approaches zero, we apply a formal limiting argument to derive a one-dimensional upscaled (effective) model. For a better understanding of the terms and parameters involved in the equations derived here, we compare these equations to a well-known core-scale model from the literature. / Download |

4. | Class, H., A. Ebigbo, R Helmig, H.K. Dahle, J.M. Nordbotten, M.A. Celia, P. Audigane, M. Darcis, J. Ennis-King, and Y. Fan A benchmark study on problems related to CO2 storage in geologic formations , Computational Geosciences, 13/4, Sp. Iss. SI, pp. 409-434, 2009. Abstract This paper summarises the results of a benchmark study that compares a number of mathematical and numerical models applied to specific problems in the context of carbon dioxide (CO2) storage in geologic formations. The processes modelled comprise advective multi-phase flow, compositional effects due to dissolution of CO2 into the ambient brine and non-isothermal effects due to temperature gradients and the Joule–Thompson effect. The problems deal with leakage through a leaky well, methane recovery enhanced by CO2 injection and a reservoir-scale injection scenario into a heterogeneous formation. We give a description of the benchmark problems then briefly introduce the participating codes and finally present and discuss the results of the benchmark study. / Download |

3. | Kopp, A., A. Ebigbo, A. Bielinski, H. Class, and R. Helmig Numerical simulation of temperature changes caused by CO2 injection in geological reservoirs, AAPG Studies in Geology, 59/26, pp. 439-456, 2009. Abstract Injection of CO 2 into the subsurface for geological storage has an effect on the temperature of the storage formation and the CO 2 itself. Numerical investigations are an essential tool in describing the relevant processes that determine such changes and the impact they may have on the migration and the storage mechanisms of CO 2 in the subsurface. This chapter focuses on the numerical simulation of such thermal effects and their consequences. Simulating the temperature changes in a storage site can be of interest for temperature-based monitoring. Determining whether or how such thermal effects change the transport of CO 2 in the formation is important for the success of a CO 2 storage effort. In particular, this chapter examines a leakage scenario and how temperature changes could affect the leakage flow. The second part of the chapter presents results of a complex reservoir-scale simulation. The target formation forms an anticlinal structure at a depth of about 570-900 m (1870-2953 ft). Strong temperature effects can be expected because of the possible tran Numerical simulation of temperature changes caused by CO2 injection in geological reservoirs. Available from: https://www.researchgate.net/publication/230838241_Numerical_simulation_of_temperature_changes_caused_by_CO2_injection_in_geological_reservoirs [accessed May 3, 2017]. / Download |

2. | Ebigbo, A., H. Class, and R. Helmig CO2 leakage through an abandoned well: problem- oriented benchmarks, Computational Geosciences, 11/2, pp. 103-115, 2007. Abstract The efficiency and sustainability of carbon dioxide (CO2) storage in deep geological formations crucially depends on the integrity of the overlying cap-rocks. Existing oil and gas wells, which penetrate the formations, are potential leakage pathways. This problem has been discussed in the literature, and a number of investigations using semi-analytical mathematical approaches have been carried out by other authors to quantify leakage rates. The semi-analytical results are based on a number of simplifying assumptions. Thus, it is of great interest to assess the influence of these assumptions. We use a numerical model to compare the results with those of the semi-analytical model. Then we ease the simplifying restrictions and include more complex thermodynamic processes including sub- and supercritical fluid properties of CO2 and non-isothermal as well as compositional effects. The aim is to set up problem-oriented benchmark examples that allow a comparison of different modeling approaches to the problem of CO2 leakage. / Download |

1. | Class, H., A. Bielinski, R. Helmig, A. Kopp, and A. Ebigbo Numerical simulation of CO2 storage in geological formations, Chemie Ingenieur Technik, 78/4, pp. 445-452, 2006. Abstract Die Speicherung von Kohlendioxid in geologischen Formationen wird derzeit als ein möglicher Beitrag zur Reduktion von Treibhausgaskonzentrationen in der Atmosphäre diskutiert und untersucht. Ein wichtiges Werkzeug begleitend zu experimentellen Untersuchungsmethoden ist die numerische Simulation. Hier wird ein Einblick in die physikalischen bzw. thermodynamischen Vorgänge im Untergrund während und nach einer Injektion von Kohlendioxid gegeben sowie deren modellkonzeptionelle Beschreibung durch nichtisotherme Mehrphasensysteme dargestellt. Die elementaren Schritte der physikalisch/mathematischen Modellbildung und numerischen Lösung der daraus entstehenden Gleichungssysteme werden erklärt. Anhand von Simulationsbeispielen werden dominierende Prozesse während und nach einer Injektion diskutiert. / Download |

### PROCEEDINGS REFEREED

2. | Hommel, J., A. Ebigbo, R. Gerlach, A. B. Cunningham, R. Helmig, and H. Class Finding a Balance between Accuracy and Effort For Modeling Biomineralization, Energy Procedia, 97, pp. 379-386, 2016. / Download |

1. | Hommel, J., A. B. Cunningham, R. Helmig, A. Ebigbo, and H. Class Numerical Investigation of Microbially Induced Calcite Precipitation as a Leakage Mitigation Technology, Energy Procedia, 40, pp. 392-397, 2013. / Download |

### THESES

2. | Ebigbo, A. Modelling of biofilm growth and its influence on CO2 and water (two-phase) flow in porous media, Dissertation University of Stuttgart, 131 pp., 2009. Abstract Bacterial biofilms are groups of microbial cells attached to surfaces and to each other. Cells in a biofilm are protected from adverse external conditions. In natural environments, this attached mode of growth is more successful than the suspended mode, and a major portion of microbial activity takes place at surfaces. In porous media, biofilms are used as bioreactors (e.g, in wastewater treatment) and as biobarriers (e.g., in enhanced oil recovery). They are also used in the containment and degradation of contaminants in groundwater aquifers. It has been proposed that biofilms be used as biobarriers for the mitigation of carbon dioxide (CO2) leakage from a geological storage reservoir. The concentration of greenhouse gases — particularly carbon dioxide (CO2) — in the atmosphere has been on the rise in the past decades. One of the methods which have been proposed to help reduce anthropogenic CO2 emissions is the capture of CO2 from large, stationary point sources and storage in deep geological formations. The caprock is an impermeable geological layer which prevents the leakage of stored CO2, and its integrity is of utmost importance for storage security. As mentioned above, biofilms could be used as biobarriers which help prevent the leakage of CO2 through the caprock in injection well vicinity. Due to the high pressure build-up during injection, the caprock in the vicinity of the well is particularly at risk of fracturing. The biofilm could also protect well cement from corrosion by CO2-rich brine. The goal of this work is to develop and test a numerical model which is capable of simulating the development of a biofilm in a CO2 storage reservoir. This involves the description of the growth of the biofilm, flow and transport in the geological formation, and the interaction between the biofilm and the flow processes. Important processes which are accounted for in the model include the effect of biofilm growth on the permeability of the formation, the hazardous effect of supercritical CO2 on suspended and attached bacteria, attachment and detachment of biomass, and two-phase fluid flow processes. The partial differential equations which describe the system are discretised in space with a vertex-centered finite volume method, and an implicit Euler scheme is used for time discretisation. The model is tested by comparing simulation results to experimental data. In a test case simulation, the model predicts the extent of biomass accumulation near an injection well and its effect on the permeability of the formation. The simulations show that the biobarrier is only effective for a limited amount of time. Regular injection of nutrients would be necessary to sustain the biofilm. In future work, the model could be extended to account for the active precipitation of minerals by the biofilm which would lead to a more enduring barrier. The model also needs to be extended to account for more than one growth-limiting factor. This would allow for the simulation of injection strategies which aim at growing a biofilm at some distance from the injection well. Biofilme, die in einem porösen Medium wachsen, blockieren Poren und verändern dabei die Eigenschaften des porösen Mediums. Diese veränderten Eigenschaften werden bei der biologischen Filtration (z. B. bei der Abwasserbehandlung), bei der biologischen Altlastensanierung (z. B. für die Erstellung hydraulischer Barrieren) und bei anderen Fragestellungen auf diesem Gebiet genutzt. Eine hydraulische Barriere biologischen Ursprungs könnte z. B. auch in einer geologischen Kohlendioxid-Lagerstätte eingesetzt werden, um das Entweichen von CO2 zu verhindern. CO2 ist das derzeit für am Wichtigsten erachtete anthropogene Treibhausgas. Die globale Erderwärmung wird demnach sehr stark durch die in den letzten Jahrzehnten stattfindende Anreicherung von anthropogenen Treibhausgasen in der Atmosphäre mitverursacht. Die Freisetzung von CO2 kann mit Hilfe effizienterer Technologien und alternativer Energiequellen reduziert werden. CO2-Emissionen können aber auch reduziert werden, indem man CO2 aus Kraftwerksabgasen abscheidet und in tiefen geologischen Formationen speichert. Bei den physikalischen Bedingungen, die in diesen unterirdischen Lagerstätten herrschen, liegt CO2 im überkritischen Zustand vor, gekennzeichnet durch eine hohe Dichte und geringe Viskosität. Diese Lagerstätten enthalten oft salzhaltiges Wasser, das dichter ist als CO2. Eine möglichst undurchlässige geologische Deckschicht verhindert das Aufsteigen des leichteren CO2 an die Erdoberfläche. Jedoch müssen, z. B. im Rahmen von Risikostudien, mögliche Störungen oder Risse in dieser Deckschicht betrachtet werden, die zu einem Entweichen des CO2 führen könnten. Die Deckschicht in der Nähe eines CO2-Injektionsbrunnens ist besonders gefährdet. Der hohe Druckanstieg während der ersten Injektionsphase, Zementkorrosion am Brunnen aufgrund des CO2-reichen Formationswassers und eventuelle Beschädigungen der Deckschicht während der Erstellung des Bohrlochs sind als mögliche Ursachen für gestörte Deckschichten zu nennen. Biobarrieren könnten verwendet werden, um solche Risiken zu minimieren, z. B. indem sie Risse in der Deckschicht abdichten oder den Bohrlochzement vor Korrosion schützen. Eine Biobarriere kann aus einem Biofilm selbst bestehen, aber auch aus vom Biofilm begüngstigten mineralischen Ablagerungen. Die vorliegende Arbeit behandelt im Wesentlichen die Entwicklung eines numerischen Modells, um die Anreicherung von mikrobieller Biomasse im Untergrund simulieren zu können. Das entwickelte Modell soll in der Lage sein, das Abdichten der beschädigten geologischen Deckschicht einer unterirdischen Kohlendioxid-Lagerstätte mit Hilfe von Biofilmen zu simulieren. Dafür müssen einerseits Strömungsprozesse und andererseits auch die mikrobielle Aktivität sowie die Interaktion dieser Vorgänge in porösen Medien richtig beschrieben werden. Die Anreicherung von Bakterien in einem porösen Medium beeinflusst die hydraulischen Eigenschaften des Mediums und als Folge davon auch die darin stattfindende Strömung. Im Gegenzug bestimmt die Strömung den Transport der mikrobiellen Nährstoffe und damit auch die Verteilung mikrobieller Wachstumsraten. Dementsprechend ist die richtige Beschreibung der Wechselwirkung zwischen Strömung und mikrobiellen Prozessen eine wesentliche Herausforderung in der Modellbildung. / Download |

1. | Ebigbo, A. Thermal Effects of Carbon Dioxide Sequestration in the Subsurface. Diplomarbeit, Institut für Wasserbau, MSc Thesis University of Stuttgart, 57 pp., 2005. Abstract In order to secure long-term storage of CO 2 in the subsurface, one has to be able to model the mass transport of CO2. The physical properties of CO2 have a very strong influence on the mass transport. These physical properties, in turn are dependent on temperature (and pressure). Hence, the modelling of heat transport in the subsurface during CO2 injection should be part of mass transport. In the simple set-up examined in this thesis, the Joule-Thompson cooling (as a result of the pressure drop) plays an important role near the injection point. The amount of cooling at the injection point depends on the temperature of the injected CO2 (and of course on the pressure difference that brought about the cooling). The temperature at that region achieves a stable value with time. The impermeable layer should be deeper than the depth at which the CO2 becomes gaseous as it is easier to trap less mobile supercritical or liquid CO2. / Download |