Numerical modelling and analysis of reactive flow and wormhole formation in fractured carbonate rocks

Reactive Flows in Porous Media: Challenges in Theoretical and Numerical Methods

We review theoretical and computational research, primarily from the past 10 years, addressing the flow of reactive fluids in porous media. The focus is on systems where chemical reactions at the solid-fluid interface cause dissolution of the surrounding porous matrix, creating nonlinear feedback mechanisms that can often lead to greatly enhanced permeability. We discuss insights into the evolution of geological forms that can be inferred from these feedback mechanisms, as well as some geotechnical applications such as enhanced oil recovery, hydraulic fracturing, and carbon sequestration. Until recently, most practical applications of reactive transport have been based on Darcy-scale modeling, where averaged equations for the flow and reactant transport are solved. We summarize the successes and limitations of volume averaging, which leads to Darcy-scale equations, as an introduction to pore-scale modeling. Pore-scale modeling is computationally intensive but offers new insights as well as tests of averaging theories and pore-network models. We include recent research devoted to validation of pore-scale simulations, particularly the use of visual observations from microfluidic experiments.

Reactive Transport Simulation of Cavern Formation along Fractures in Carbonate Rocks

Karst cavities and caves are often present along fractures in limestone reservoirs and are of significance for oil and gas exploration. Understanding the formation and evolution of caves in fractured carbonate rocks will enhance oil and gas exploration and development. Herein, a reactive transport model was established considering both the matrix and fractures. Different factors affecting the dissolution along fractures were considered in the simulation of matrix–fracture carbonate rocks, including the magnitude and characteristic length of the matrix porosity heterogeneity, intersecting fractures, and complex fracture network. The results show that a strong heterogeneity of the matrix porosity significantly affects the cave formation along the fracture and the existence of fractures increases the heterogeneity due to the high permeability as well as the dissolution area. The characteristic length of the matrix porosity heterogeneity affects the cave location and shape. The larger permeability of intersecting fractures or the matrix greatly increases the cave size, leading to the formation of large, connected cave areas. A complex fracture network leads to more developed karst dissolution caves. The topology of the fracture network and preferential flow dominate the distribution of caves and alleviate the effect of the matrix heterogeneity.

A Numerical Simulator for Modeling the Coupling Processes of Subsurface Fluid Flow and Reactive Transport Processes in Fractured Carbonate Rocks

Water–rock interactions can alter rock properties through chemical reactions during subsurface transport processes like geological CO2 sequestration (GCS), matrix acidizing, and waterflooding in carbonate formations. Dynamic changes in rock properties cause a failure of waterflooding and GCS and could also dramatically affect the efficiency of the acidizing. Efficient numerical simulations are thus essential to the optimized design of those subsurface processes. In this paper, we develop a three-dimensional (3D) numerical model for simulating the coupled processes of fluid flow and chemical reactions in fractured carbonate formations. In the proposed model, we employ the Stokes–Brinkman equation for momentum balance, which is a single-domain formulation for modeling fluid flow in fractured porous media. We then couple the Stokes–Brinkman equation with reactive-transport equations. The model can be formulated to describe linear as well as radial flow. We employ a decoupling procedure that sequentially solves the Stokes–Brinkman equation and the reactive transport equations. Numerical experiments show that the proposed method can model the coupled processes of fluid flow, solute transport, chemical reactions, and alterations of rock properties in both linear and radial flow scenarios. The rock heterogeneity and the mineral volume fractions are two important factors that significantly affect the structure of conductive channels.

Matrix acidizing: a fouling mitigation process in oil and gas wells

Fouling mitigation in underground reservoirs enhances the permeability and the flow capacity of production or injection wells and is carried out by reservoir stimulation methods such as matrix acidizing. This process is known as the most significant method used to improve the production or injection indices of oil and gas wells as well as water and steam wells. Here, different aspects of this process, its chemical advances and novel high-technologies are compared and discussed in order to reveal their advantages and determine under what conditions they are applicable. Knowledge for adapting the proper acid treatment with the well characteristics is another issue that has been considered in this paper. The final goal is to present the state-of-the-art fouling mitigation methods based on novel experiments, simulations and investigations in order to emphasize the engineering aspects of fouling mitigation in oil and gas wells by matrix acidizing.

3-D Modelling and Experimental Comparison of Reactive Flow in Carbonates under Radial Flow Conditions

We use a two-scale continuum model to simulate reactive flow and wormhole formation in carbonate rocks under 3-D radial flow conditions. More specifically, we present a new structure-property relationship based on the fractal geometry theory, to describe the evolution of local permeability, pore radius, and specific area with porosity variation. In the numerical calculation, to improve the convergence rate, the heterogeneous medium in question is extended by adding a thin layer of homogeneous porous medium to its inlet. We compare the simulation results with the available experimental observations and find that they are qualitatively consistent with each other. Additionally, sensitivity analysis of the dissolution process with respect to acid injection rate and rock heterogeneity, including heterogeneity magnitude and correlation length, is presented.