Pore-Network Modeling Dedicated to the Determination of the Petrophysical-Property Changes in the Presence of Reactive Fluid

SPE Journal ◽  
2010 ◽  
Vol 15 (03) ◽  
pp. 618-633 ◽  
Author(s):  
Lionnel Algive ◽  
Samir Bekri ◽  
Olga Vizika

Summary A pore-network model (PNM) is an efficient tool to account for phenomena occurring at the pore scale. Its explicit 3D network of pores interconnected by throats represents an easy way to consider the topology and geometry effects on upscaled and homogenized petrophysical parameters. In particular, this modeling approach is appropriate to study the rock/fluid interactions. It can provide quantitative information both on the effective transport property modifications caused by the reactions and on the structure evolution resulting from dissolution/precipitation mechanisms. The model developed is based on the resolution of the macroscopic reactive transport equation between the nodes of the network. By upscaling the results, we then determined the effective transport properties at the core scale. A sensitivity study on reactive and flow regimes has been conducted in the case of single-phase flow in the limit of long times. It has been observed that the mean reactive solute velocity and dispersion can vary up to one order of magnitude compared with the tracer values because of the concentration-profile heterogeneity at the pore scale resulting from the surface reactions. As for the reactive apparent coefficient, when the kinetics is limited by the mass transfer, it can decrease by several orders of magnitude with regard to that calculated by the usual perfect-mixing assumption. That is why scale factors should be added to the classical macroscopic transport equation implemented in reservoir simulators to predict accurately the reactive flow effects. For each study case, we also obtained the permeability variation vs. the porosity evolution in a physical way that accounts for reactive transport conditions. It appears that the wall-deformation pattern and its effect on petrophysical properties must be explained by considering both microscopic and macroscopic scales of the reactive transport, each one governed by a dimensionless number comparing reaction and transport characteristic times. This work helps improve the understanding of surface-reactions effects on reactive flow on the one hand and on permeability and porosity modifications on the other. Using the PNM approach, scale-factor parameters and permeability-vs.-porosity relations can be determined for various rock types and reactive flow regimes. Once integrated as inputs in a reservoir simulator, these relations form a powerful and convenient means of enhancing the modeling accuracy of the change in petrophysical properties during injection of a reactive fluid, such as brine rich in carbon dioxide (CO2).

2022 ◽  
Vol 3 ◽  
Author(s):  
Vitalii Starchenko

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.


2020 ◽  
Author(s):  
Einat Aharonov ◽  
Roi Roded ◽  
Ran Holtzman ◽  
Piotr Szymczak

<p>Dissolution by reactive fluid flow is a fundamental process in geological systems. It controls diagenesis and karst evolution and has broad implications for groundwater hydrology. Specifically, reactive flow controls the evolution of the void-space structure via the feedback between the reaction and transport. In some instances, advective transport rate is high compared to that of geochemical reactions (low Damkӧhler number, Da), such that the reactive fluid penetrates the system before its reactivity is exhausted, resulting in a relatively spatially-uniform dissolution. Despite the importance of low Da conditions, the emerging transformations in the medium structure, flow field, and its bulk properties are not well understood. Likewise, our ability to decipher diagenetic history and preexisting structure is lacking.</p><p>Here, using a network model, we investigate the evolution of heterogeneous and anisotropic medium during dissolution at low Da conditions. The numerical simulations show that the medium progressively becomes more homogeneous as well as isotropic, which consequently makes the flow field more uniform. Homogenization is particularly notable for anisotropic media, in which the transverse channels are wide relative to the channels parallel to the main flow direction. In this case, flow is initially focused within a few highly tortuous pathways, hence emphasizing the effect of dissolution on flow heterogeneity and tortuosity. The homogenization process is further enhanced when the surface reaction is transport-controlled—that is, when diffusion of dissolved ions away from the mineral surface to the bulk fluid is slow, reducing the reactivity adjacent to the surface: At first, since diffusive transport is more effective in narrow channels, they undergo faster dissolution, which selectively enlarges them leading to an initial steep rise in permeability. Later, however, as dissolution proceeds and the channels broaden, the overall dissolution rate drops, diminishing the growth rate of permeability. Our findings provide fundamental insights into reactive transport and hydrogeological processes in fractured and porous media.</p>


2016 ◽  
Vol 369 ◽  
pp. 125-130 ◽  
Author(s):  
Qing Rong Xiong ◽  
Andrey P. Jivkov

Adsorption of solutes in porous media is typically represented as an equilibrium process. The adsorption process is determined by the concentration of the solute in the fluid next to the solid grain in individual pores. The concentration near the solid surface is different from the average concentration within the pore due to the development of concentration gradients within the pore space. Simulations with pore network models are usually based on the average concentration in individual pore throats. To advance the realism of such models, we develop a relationship between pore-scale adsorption coefficients and corresponding coarse-grained adsorption parameters. A numerical scheme is proposed: firstly the solute transport within a single pore is simulated undergoing equilibrium adsorption at the pore wall, and secondly flux-averaged concentration breakthrough curves are obtained. By fitting the breakthrough curves, the coarse-grained adsorption parameters are determined from the pore-scale hydraulic and adsorption parameters.


2014 ◽  
Vol 135 ◽  
pp. 321-335 ◽  
Author(s):  
Victoria Boyd ◽  
Hongkyu Yoon ◽  
Changyong Zhang ◽  
Mart Oostrom ◽  
Nancy Hess ◽  
...  

2021 ◽  
pp. 130917
Author(s):  
Xiangqian Wei ◽  
Wenzhi Li ◽  
Qiying Liu ◽  
Weitao Sun ◽  
Siwei Liu ◽  
...  

SPE Journal ◽  
2016 ◽  
Vol 21 (01) ◽  
pp. 280-292 ◽  
Author(s):  
John Lyons ◽  
Hadi Nasrabadi ◽  
Hisham A. Nasr-El-Din

Summary Fracture acidizing is a well-stimulation technique used to improve the productivity of low-permeability reservoirs and to bypass deep formation damage. The reaction of injected acid with the rock matrix forms etched channels through which oil and gas can then flow upon production. The properties of these etched channels depend on the acid-injection rate, temperature, reaction chemistry, mass-transport properties, and formation mineralogy. As the acid enters the formation, it increases in temperature by heat exchange with the formation and the heat generated by acid reaction with the rock. Thus, the reaction rate, viscosity, and mass transfer of acid inside the fracture also increase. In this study, a new thermal-fracture-acidizing model is presented that uses the lattice Boltzmann method to simulate reactive transport. This method incorporates both accurate hydrodynamics and reaction kinetics at the solid/liquid interface. The temperature update is performed by use of a finite-difference technique. Furthermore, heterogeneity in rock properties (e.g., porosity, permeability, and reaction rate) is included. The result is a model that can accurately simulate realistic fracture geometries and rock properties at the pore scale and that can predict the geometry of the fracture after acidizing. Three thermal-fracture-acidizing simulations are presented here, involving injection of 15 and 28 wt% of hydrochloric acid into a calcite fracture. The results clearly show an increase in the overall fracture dissolution because of the addition of temperature effects (increasing the acid-reaction and mass-transfer rates). It has also been found that by introducing mineral heterogeneity, preferential dissolution leads to the creation of uneven etching across the fracture surfaces, indicating channel formation.


2021 ◽  
Author(s):  
Hanbang Zou ◽  
Pelle Ohlsson ◽  
Edith Hammer

<p>Carbon sequestration has been a popular research topic in recent years as the rapid elevation of carbon emission has significantly impacted our climate. Apart from carbon capture and storage in e.g. oil reservoirs, soil carbon sequestration offers a long term and safe solution for the environment and human beings. The net soil carbon budget is determined by the balance between terrestrial ecosystem sink and sources of respiration to atmospheric carbon dioxide. Carbon can be long term stored as organic matters in the soil whereas it can be released from the decomposition of organic matter. The complex pore networks in the soil are believed to be able to "protect" microbial-derived organic matter from decomposition. Therefore, it is important to understand how soil structure impacts organic matter cycling at the pore scale. However, there are limited experimental studies on understanding the mechanism of physical stabilization of organic matter. Hence, my project plan is to create a heterogeneous microfluidic porous microenvironment to mimic the complex soil pore network which allows us to investigate the ability of organisms to access spaces starting from an initial ecophysiological precondition to changes of spatial accessibility mediated by interactions with the microbial community.</p><p>Microfluidics is a powerful tool that enables studies of fundamental physics, rapid measurements and real-time visualisation in a complex spatial microstructure that can be designed and controlled. Many complex processes can now be visualized enabled by the development of microfluidics and photolithography, such as microbial dynamics in pore-scale soil systems and pore network modification mimicking different soil environments – earlier considered impossible to achieve experimentally. The microfluidic channel used in this project contains a random distribution of cylindrical pillars of different sizes so as to mimic the variations found in real soil. The randomness in the design creates various spatial availability for microbes (preferential flow paths with dead-end or continuous flow) as an invasion of liquids proceeds into the pore with the lowest capillary entry pressure. In order to study the impact of different porosity in isolation of varying heterogeneity of the porous medium, different pore size chips that use the same randomly generated pore network is created. Those chips have the same location of the pillars, but the relative size of each pillar is scaled. The experiments will be carried out using sterile cultures of fluorescent bacteria, fungi and protists, synthetic communities of combinations of these, or a whole soil community inoculum. We will quantify the consumption of organic matter from the different areas via fluorescent substrates, and the bio-/necromass produced. We hypothesise that lower porosity will reduce the net decomposition of organic matter as the narrower pore throat limits the access, and that net decomposition rate at the main preferential path will be higher than inside branches</p>


Sign in / Sign up

Export Citation Format

Share Document