uu.seUppsala University Publications
Change search
Link to record
Permanent link

Direct link
BETA
Niemi, Auli
Publications (10 of 85) Show all publications
Dessirier, B., Tsang, C.-F. & Niemi, A. (2018). A new scripting library for modeling flow and transport in fractured rock with channel networks. Computers & Geosciences, 111, 181-189
Open this publication in new window or tab >>A new scripting library for modeling flow and transport in fractured rock with channel networks
2018 (English)In: Computers & Geosciences, ISSN 0098-3004, E-ISSN 1873-7803, Vol. 111, p. 181-189Article in journal (Refereed) Published
Abstract [en]

Deep crystalline bedrock formations are targeted to host spent nuclear fuel owing to their overall low permeability. They are however highly heterogeneous and only a few preferential paths pertaining to a small set of dominant rock fractures usually carry most of the flow or mass fluxes, a behavior known as channeling that needs to be accounted for in the performance assessment of repositories. Channel network models have been developed and used to investigate the effect of channeling. They are usually simpler than discrete fracture networks based on rock fracture mappings and rely on idealized full or sparsely populated lattices of channels. This study reexamines the minimal requirements to describe a channel network in terms of groundwater flow and solute transport, leading to an extended description suitable for unstructured arbitrary networks of channels. An implementation of this formalism in a Python scripting library is presented and released along with this article. A new algebraic multigrid preconditioner delivers a significant speedup in the flow solution step compared to previous channel network codes. 3D visualization is readily available for verification and interpretation of the results by exporting the results to an open and free dedicated software. The new code is applied to three example cases to verify its results on full uncorrelated lattices of channels, sparsely populated percolation lattices and to exemplify the use of unstructured networks to accommodate knowledge on local rock fractures.

Keywords
channel network, groundwater flow, solute transport, graph theory
National Category
Oceanography, Hydrology and Water Resources
Research subject
Hydrology
Identifiers
urn:nbn:se:uu:diva-321572 (URN)10.1016/j.cageo.2017.11.013 (DOI)000423005900018 ()
Funder
Swedish Radiation Safety Authority, SSM2016-763
Available from: 2017-05-08 Created: 2017-05-08 Last updated: 2018-03-06Bibliographically approved
Dessirier, B., Tsang, C.-F. & Niemi, A. (2018). A new scripting library for modeling flow and transport in fractured rock with channel networks. Computers & Geosciences, 111, 181-189
Open this publication in new window or tab >>A new scripting library for modeling flow and transport in fractured rock with channel networks
2018 (English)In: Computers & Geosciences, ISSN 0098-3004, E-ISSN 1873-7803, Vol. 111, p. 181-189Article in journal (Refereed) Published
Abstract [en]

Deep crystalline bedrock formations are targeted to host spent nuclear fuel owing to their overall low permeability. They are however highly heterogeneous and only a few preferential paths pertaining to a small set of dominant rock fractures usually carry most of the flow or mass fluxes, a behavior known as channeling that needs to be accounted for in the performance assessment of repositories. Channel network models have been developed and used to investigate the effect of channeling. They are usually simpler than discrete fracture networks based on rock fracture mappings and rely on idealized full or sparsely populated lattices of channels. This study reexamines the fundamental parameter structure required to describe a channel network in terms of groundwater flow and solute transport, leading to an extended description suitable for unstructured arbitrary networks of channels. An implementation of this formalism in a Python scripting library is presented and released along with this article. A new algebraic multigrid preconditioner delivers a significant speedup in the flow solution step compared to previous channel network codes. 3D visualization is readily available for verification and interpretation of the results by exporting the results to an open and free dedicated software. The new code is applied to three example cases to verify its results on full uncorrelated lattices of channels, sparsely populated percolation lattices and to exemplify the use of unstructured networks to accommodate knowledge on local rock fractures.

Keywords
Channel network, Groundwater flow, Solute transport, Graph theory
National Category
Earth and Related Environmental Sciences
Identifiers
urn:nbn:se:uu:diva-343655 (URN)10.1016/j.cageo.2017.11.013 (DOI)000423005900018 ()
Available from: 2018-05-09 Created: 2018-05-09 Last updated: 2018-05-09Bibliographically approved
Figueiredo, B., Tsang, C.-F., Rutqvist, J. & Niemi, A. (2018). Corrigendum to “The effects of nearby fractures on hydraulically induced fracture propagation and permeability changes”[Eng. Geol. 228 (2017) 197–213]. Engineering Geology, 239, 344-344
Open this publication in new window or tab >>Corrigendum to “The effects of nearby fractures on hydraulically induced fracture propagation and permeability changes”[Eng. Geol. 228 (2017) 197–213]
2018 (English)In: Engineering Geology, ISSN 0013-7952, E-ISSN 1872-6917, Vol. 239, p. 344-344Article in journal (Other academic) Published
National Category
Geology
Identifiers
urn:nbn:se:uu:diva-359373 (URN)10.1016/j.enggeo.2018.03.006 (DOI)000432769300030 ()
Funder
The Geological Survey of Sweden (SGU), 1724EU, Horizon 2020, 636811
Note

WoS title: The effects of nearby fractures on hydraulically induced fracture propagation and permeability changes (vol 228, pg 197, 2017)

Correction to: Engineering Geology, vol. 228, pages 197-213. DOI: 10.1016/j.enggeo.2017.08.011

Available from: 2018-09-04 Created: 2018-09-04 Last updated: 2018-09-04Bibliographically approved
Hedayati, M., Wigston, A., Wolf, J. L., Rebscher, D. & Niemi, A. (2018). Impacts of SO2 gas impurity within a CO2 stream on reservoir rock of a CCS pilot site: Experimental and modelling approach. International Journal of Greenhouse Gas Control, 70, 32-44
Open this publication in new window or tab >>Impacts of SO2 gas impurity within a CO2 stream on reservoir rock of a CCS pilot site: Experimental and modelling approach
Show others...
2018 (English)In: International Journal of Greenhouse Gas Control, ISSN 1750-5836, E-ISSN 1878-0148, Vol. 70, p. 32-44Article in journal (Refereed) Published
Abstract [en]

In order to evaluate chemical impacts of SO2 impurity on reservoir rock during CO2 capture and storage in deep saline aquifers, several batch reactor experiments were performed on laboratory scale using core rock samples from the pilot CO2 injection site in Heletz. In this experiment, the samples were exposed to pure N-2(g), pure CO2(g), and CO2(g) with an impurity of 1.5% SO2(g) under reservoir conditions for pressure and temperature (14.5 MPa, 60 degrees C). Based on the set-up and the obtained experimental results, a batch chemical model was established using the numerical simulation program TOUGHREACT V3.0-OMP. Comparing laboratory and simulation data provides a better understanding of the rock-brine-gas interactions. In addition, it offers an evaluation of the capability of the model to predict chemical interactions in the target injection reservoir during exposure to pure and impure CO2. The best match between the geochemical model and experimental data was achieved when the reactive surface area of minerals in the model was adjusted in order to calibrate the kinetic rates of minerals. The simulations indicated that SO2(g) tends to dissolve rather quickly and oxidizes under a kinetic control. Hence, it has a stronger effect on the acidity of the brine than pure CO2(g) and as a result, increased mineral dissolution and caused the precipitation of sulfate and sulfide minerals. Ankerite, dolomite, and siderite, the most abundant carbonates in the sandstone rock sample, are subject to stronger dissolution in the presence of SO2 gas. The performed simulations confirmed a slower dissolution rate for ankerite and siderite than for dolomite. The model reproduced the precipitation of pyrite and anhydrite as observed in the laboratory. The dissolution of dolomite observed in the batch reaction test with pure N-2 is assumed to be due to slight contamination with oxygen and modelling supported this. The inclusion of SO2 increased the porosity over that of the pure CO2 case, and is thus considered to increase the permeability and injectivity of the reservoir as well. Exposure to SO2 also increased the concentration of trace elements. The calibrated kinetic parameters determined in this study will be used to model the injection and long-term behavior of CO2 at the Heletz field site, and may be used for similar geologic reservoirs.

National Category
Geochemistry
Identifiers
urn:nbn:se:uu:diva-340872 (URN)10.1016/j.ijggc.2018.01.003 (DOI)000428773100004 ()
Funder
EU, FP7, Seventh Framework Programme, 309102, 309067
Available from: 2018-02-05 Created: 2018-02-05 Last updated: 2018-06-01Bibliographically approved
Tsang, C.-F., Figueiredo, B. & Niemi, A. (2018). Importance of stress effects on inputs to fracture network models used for subsurface flow and transport studies. International Journal of Rock Mechanics And Mining Sciences, 101, 13-17
Open this publication in new window or tab >>Importance of stress effects on inputs to fracture network models used for subsurface flow and transport studies
2018 (English)In: International Journal of Rock Mechanics And Mining Sciences, ISSN 1365-1609, E-ISSN 1873-4545, Vol. 101, p. 13-17Article in journal (Refereed) Published
National Category
Geosciences, Multidisciplinary
Identifiers
urn:nbn:se:uu:diva-334138 (URN)10.1016/j.ijrmms.2017.11.012 (DOI)000418717200002 ()
Funder
Swedish Radiation Safety Authority, SSM2016-763The Geological Survey of Sweden (SGU), 1724
Available from: 2017-11-21 Created: 2017-11-21 Last updated: 2018-01-18Bibliographically approved
Rasmusson, K., Rasmusson, M., Tsang, Y., Benson, S., Hingerl, F., Fagerlund, F. & Niemi, A. (2018). Residual trapping of carbon dioxide during geological storage: insight gained through a pore-network modeling approach. International Journal of Greenhouse Gas Control, 74, 62-78
Open this publication in new window or tab >>Residual trapping of carbon dioxide during geological storage: insight gained through a pore-network modeling approach
Show others...
2018 (English)In: International Journal of Greenhouse Gas Control, ISSN 1750-5836, E-ISSN 1878-0148, Vol. 74, p. 62-78Article in journal (Refereed) Published
Abstract [en]

To reduce emissions of the greenhouse gas CO2 to the atmosphere, sequestration in deep saline aquifers is a viable strategy. Residual trapping is a key containment process important to the success of CO2 storage operations. While residual trapping affects CO2 migration over large scales, it is inherently a pore-scale process. Pore-network models (PNMs), capturing such processes, are useful for our understanding of residual trapping, and for upscaling trapping parameters for larger scale models. A PNM for simulation of quasi-static two-phase flow; CO2 intrusion (drainage) followed by water flooding (imbibition) was developed. It accounts for pore-scale displacement mechanisms, and was used to investigate residual CO2 trapping. The sensitivity of the residual CO2 saturation to several parameters was studied, to validate a trapping behavior in agreement with earlier studies. Then the PNM was calibrated to core sample data and used to simulate drainage-imbibition scenarios with different turning point saturations. From these the initial-residual saturation curves of CO2 in Heletz sandstone were estimated, essential for future macroscopic-scale simulations. Further, the occurrence of different pore-scale mechanisms were quantified and the size distribution of the residual clusters was shown to exhibit a bimodal appearance. The findings improve the understanding of residual trapping in Heletz sandstone.

National Category
Earth and Related Environmental Sciences
Identifiers
urn:nbn:se:uu:diva-327991 (URN)10.1016/j.ijggc.2018.04.021 (DOI)000434428100007 ()
Funder
EU, FP7, Seventh Framework Programme, 309067Swedish Energy Agency, 43526-1
Available from: 2017-08-15 Created: 2017-08-15 Last updated: 2018-08-30Bibliographically approved
Rasmusson, M., Rasmusson, K., Fagerlund, F., Tsang, Y. & Niemi, A. (2018). The impact of co-contaminant SO2, versus salinity and thermodynamic conditions, on residual CO2 trapping during geological storage. Greenhouse Gases: Science and Technology, 8(6), 1053-1065
Open this publication in new window or tab >>The impact of co-contaminant SO2, versus salinity and thermodynamic conditions, on residual CO2 trapping during geological storage
Show others...
2018 (English)In: Greenhouse Gases: Science and Technology, E-ISSN 2152-3878, Vol. 8, no 6, p. 1053-1065Article in journal (Refereed) Published
Abstract [en]

During geological storage in deep saline aquifers, immobilization of CO2 in reservoir rock determines both storage safety and capacity. Assessment of the sensitivity of residual trapping to different parameters (interfacial tension and contact angles) and the storage conditions affecting these is therefore of great importance. One aspect of concern is the presence of co-contaminants such as SO2 in the injected gas. Using experimentally measured values of interfacial tensions and contact angles, we apply pore-network modelling (which accounts for pore-scale mechanisms such as snap-off, cooperative pore body filling and piston-type displacement) to a generic sandstone network to quantify the impact of SO2 co-injection on residual CO2 trapping, and its relative importance as compared to the influences of thermodynamic conditions and salinity. We show that the presence of small amounts of SO2 in the injected CO2 has a notable positive effect on the amount of CO2 becoming residually trapped (similar to 3% increase at 1 wt% SO2). However, this effect is small compared to that of the brine salinity (similar to 20% decrease in residually trapped CO2 over the salinity range 0.2 to 5 M NaCl). Still, co-injection of SO2 could potentially favour the residual trapping of CO2 in reservoir rocks, especially at storage sites with inclined aquifers where the CO2 is set to migrate hydro-dynamically over long distances. The salinity of the resident brine is of primary importance during storage site selection. Furthermore, sensitivity analysis shows that the advancing contact angle strongly impacts residual CO2 trapping. 

Place, publisher, year, edition, pages
WILEY PERIODICALS, INC, 2018
Keywords
capillary trapping, carbon dioxide, co-injection, impurities, pore-network model, sulfur dioxide
National Category
Oceanography, Hydrology and Water Resources
Identifiers
urn:nbn:se:uu:diva-371116 (URN)10.1002/ghg.1816 (DOI)000451042600006 ()
Funder
Swedish Research Council, 2010-3657EU, FP7, Seventh Framework Programme, 282900EU, FP7, Seventh Framework Programme, 309102
Available from: 2018-12-28 Created: 2018-12-28 Last updated: 2018-12-28Bibliographically approved
Wolf, J. L., Niemi, A., Bensabat, J., May, F., Ruetters, H. & Rebscher, D. (2017). 2D reactive transport simulations of CO2 streams containing impurities in a saline aquifer, Heletz, Israel. In: Dixon, T Laloui, L Twinning, S (Ed.), 13Th International Conference on Greenhouse Gas Control Technologies, Ghgt-13: . Paper presented at 13th International Conference on Greenhouse Gas Control Technologies (GHGT), NOV 14-18, 2016, Lausanne, SWITZERLAND (pp. 3537-3546). Elsevier
Open this publication in new window or tab >>2D reactive transport simulations of CO2 streams containing impurities in a saline aquifer, Heletz, Israel
Show others...
2017 (English)In: 13Th International Conference on Greenhouse Gas Control Technologies, Ghgt-13 / [ed] Dixon, T Laloui, L Twinning, S, Elsevier, 2017, p. 3537-3546Conference paper, Published paper (Refereed)
Abstract [en]

In order to evaluate the chemical impacts of CO2 stream impurities on reservoir rocks, 2D reactive transport simulations using the code TOUGHREACT V3.0 were performed. The underlying reservoir properties are based on in-situ data from the CO2 injection test site Heletz, Israel. Two different CO2 compositions (mole fractions 99 % CO2 + 1 % SO2 and 98.8 % CO2 + 1 % SO2 + 0.2 % NO2, respectively) were chosen to represent oxidising impurities. Different modelling approaches, namely trace gas transport (TGT) and additional brine injection (ABI), were applied to investigate the influence of these modelling approaches on qualitative and quantitative simulation results. The simulations using either approach show an accumulation of SO2 and NO2 close to the injection well due to the preferential dissolution of these acidic impurities compared to CO2. Both modelling approaches indicate the same general chemical impact and related mineral reactions. Within the affected rock volume a distinct ankerite to anhydrite conversion occurs, which slightly enhances porosity. While the same qualitative conclusions independently from the chosen modelling approach were obtained, the quantitative magnitude of mineral conversion and the spatial extent of impurity affected rock material depend on the chosen modelling approach and thus need further investigation with respect to e.g. validation by field test data. (C) 2017 The Authors. Published by Elsevier Ltd.

Place, publisher, year, edition, pages
Elsevier, 2017
Series
Energy Procedia, ISSN 1876-6102 ; 114
Keywords
Impurities, reactive transport, SO2, NO2, TOUGHREACT V3.0, CCS
National Category
Geosciences, Multidisciplinary
Identifiers
urn:nbn:se:uu:diva-349328 (URN)10.1016/j.egypro.2017.03.1483 (DOI)000419147303068 ()
Conference
13th International Conference on Greenhouse Gas Control Technologies (GHGT), NOV 14-18, 2016, Lausanne, SWITZERLAND
Funder
EU, FP7, Seventh Framework Programme, 309102
Available from: 2018-04-26 Created: 2018-04-26 Last updated: 2018-04-26Bibliographically approved
Tian, L., Wilkinson, R., Yang, Z., Power, H., Fagerlund, F. & Niemi, A. (2017). Gaussian Process Emulators for Quantifying Uncertainty in CO2 Spreading Predictions in Heterogeneous Media. Computers & Geosciences
Open this publication in new window or tab >>Gaussian Process Emulators for Quantifying Uncertainty in CO2 Spreading Predictions in Heterogeneous Media
Show others...
2017 (English)In: Computers & Geosciences, ISSN 0098-3004, E-ISSN 1873-7803Article in journal (Refereed) Published
Abstract [en]

We explore the use of Gaussian process emulators (GPE) in the numerical simulation of CO2 injection into a deep heterogeneous aquifer. The model domain is a two-dimensional, log-normally distributed stochastic permeability field. We first estimate the cumulative distribution functions (CDFs) of the CO2 breakthrough time and the total CO2 mass using a computationally expensive Monte Carlo (MC) simulation. We then show that we can accurately reproduce these CDF estimates with a GPE, using only a small fraction of the computational cost required by traditional MC simulation. In order to build a GPE that can predict the simulator output from a permeability field consisting of 1000s of values, we use a truncated Karhunen-Loève (K-L) expansion of the permeability field, which enables the application of the Bayesian functional regression approach. We perform a cross-validation exercise to give an insight of the optimization of the experiment design for selected scenarios: we find that it is sufficient to use 100s values for the size of the training set and that it is adequate to use as few as 15 K-L components. Our work demonstrates that GPE with truncated K-L expansion can be effectively applied to uncertainty analysis associated with modeling of multiphase flow and transport processes in heterogeneous media.

Keywords
CO2, Bayesian, Permeability, KL expansion, Monte Carlo, Cumulative distribution function, Uncertainty analysis
National Category
Geosciences, Multidisciplinary Computer Sciences
Identifiers
urn:nbn:se:uu:diva-298748 (URN)10.1016/j.cageo.2017.04.006 (DOI)000404697000011 ()
Funder
EU, FP7, Seventh Framework Programme, 227286EU, FP7, Seventh Framework Programme, 282900EU, FP7, Seventh Framework Programme, 309067
Available from: 2016-07-06 Created: 2016-07-06 Last updated: 2018-02-06Bibliographically approved
Basirat, F., Yang, Z. & Niemi, A. (2017). Pore-scale modeling of wettability effects on CO2–brine displacement during geological storage. Advances in Water Resources, 109, 181-195
Open this publication in new window or tab >>Pore-scale modeling of wettability effects on CO2–brine displacement during geological storage
2017 (English)In: Advances in Water Resources, ISSN 0309-1708, E-ISSN 1872-9657, Vol. 109, p. 181-195Article in journal (Refereed) Published
Abstract [en]

Wetting properties of reservoir rocks and caprocks can vary significantly, and they strongly influence geological storage of carbon dioxide in deep saline aquifers, during which CO2 is supposed to displace the resident brine and to become permanently trapped. Fundamental understanding of the effect of wettability on CO2-brine displacement is thus important for improving storage efficiency and security. In this study, we investigate the influence of wetting properties on two-phase flow of CO2 and brine at the pore scale. A numerical model based on the phase field method is implemented to simulate the two-phase flow of CO2-brine in a realistic pore geometry. Our focus is to study the pore-scale fluid-fluid displacement mechanisms under different wetting conditions and to quantify the effect of wettability on macroscopic parameters such as residual brine saturation, capillary pressure, relative permeability, and specific interfacial area. Our simulation results confirm that both the trapped wetting phase saturation and the normalized interfacial area increase with decreasing contact angle. However, the wetting condition does not appear to influence the CO2 breakthrough time and saturation. We also show that the macroscopic capillary pressures based on the pressure difference between inlet and outlet can differ significantly from the phase averaging capillary pressures for all contact angles when the capillary number is high ( log Ca > -5). This indicates that the inlet-outlet pressure difference may not be a good measure of the continuum-scale capillary pressure. In addition, the results show that the relative permeability of CO2 can be significantly lower in strongly water-wet conditions than in the intermediate-wet conditions.

Place, publisher, year, edition, pages
Elsevier, 2017
National Category
Oceanography, Hydrology and Water Resources
Research subject
Hydrology
Identifiers
urn:nbn:se:uu:diva-315304 (URN)10.1016/j.advwatres.2017.09.004 (DOI)000416037100014 ()
Funder
EU, FP7, Seventh Framework Programme, 309067Swedish Research Council, 637-2014-445Swedish Energy Agency, 43526-1
Available from: 2017-02-14 Created: 2017-02-14 Last updated: 2018-02-23Bibliographically approved
Organisations

Search in DiVA

Show all publications