Modeling and analysis of the acidizing process in carbonate rocks using a two-phase thermal-hydrologic-chemical coupled model

Liu, Piyang; Yan, Xiangqiao; Ye, Jun; Sun, Shuyu

doi:10.1016/j.ces.2019.06.017

Cited by 46 publications

(34 citation statements)

References 74 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…(21)) to calculate the acid concentration. The operator splitting method is adopted here for solving the reactive-transport coupling problem as used in Maheshwari et al [27] and Liu et al [43]. In a time step, the diffustion-convection operator is solved firstly to calculate the acid concentration at time t*, at which the transport process has completed, but the reaction process has not yet started.…”

Section: Dimensionless Modelmentioning

confidence: 99%

“…The simulations presented in this work are performed with the number of grid fixed at 100 Â 40 Â 40 (100 cells in x direction, 40 cells in y and z directions). The grid size is established by comparing the simulation accuracy as described in Liu et al [43]. To put it briefly, the grid is refined until the simulation result is insensitive to any grid size changes.…”

Section: Dimensionless Modelmentioning

confidence: 99%

“…Based on previous works, Panga et al [26] presented a continuum model that can capture both the kinetic and mass transfer controlled regime simultaneously. After that, the continuum model has been widely used to investigate the effect of rock or fluid properties, including flow geometry (linear and radial) [30][31][32], completion methods [33], the type of injected acid (hydrochloric acid and self-diverting acid) [34][35][36][37][38], medium heterogeneity (the presence of vugs and fractures) [20,23,39,40], and reservoir temperature [41][42][43] on wormhole propagation and optimal injection volume.…”

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Three-dimensional simulation of acidizing process in carbonate rocks using the Darcy–Forchheimer framework

Liu

Ren

Kong

et al. 2020

Oil Gas Sci. Technol. – Rev. IFP Energies nouvelles

Self Cite

View full text Add to dashboard Cite

Acidizing is an economical and effective practice to remove the near wellbore damage, which is performed by injecting acid into the formation through the wellbore. The injected acid dissolves the rock, by which the permeability nearby the wellbore can be improved. For a carbonate reservoir, the injected acid dissolves some of the minerals and some narrow and long channels, named wormholes, are formed then. These wormholes can bypass the damaged zone and hence improve the productivity of the well. The process for acid dissolving rocks involves complex physicochemical change, including the chemical reactions at the pore scale and the fluid flow at Darcy scale. In this paper, a 3-D reactive flow model with non-Darcy framework is developed based on the two-scale continuum model, and is solved by using the finite volume method. Five types of dissolution patterns, named face dissolution, conical wormhole, wormhole, ramified wormhole, and uniform dissolution, are obtained as the injection velocity increases. The effect of non-Darcy flow on dissolution pattern and breakthrough volume is analyzed. It is found that there is no effect of non-Darcy on dissolution structure and breakthrough volume when the injection velocity is very low. However, when the injection velocity is very high, the generated wormhole has more branches when using the Forchheimer equation than using the Darcy equation. Moreover, the optimal injection velocity is found to be the same whether considering the non-Darcy flow or not.

show abstract

Section: Dimensionless Modelmentioning

confidence: 99%

Section: Dimensionless Modelmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Three-dimensional simulation of acidizing process in carbonate rocks using the Darcy–Forchheimer framework

Liu

Ren

Kong

et al. 2020

Oil Gas Sci. Technol. – Rev. IFP Energies nouvelles

Self Cite

View full text Add to dashboard Cite

show abstract

“…Note that advection terms in the Cahn-Hilliard and Navier-Stokes equations are discretized by the MINMOD scheme (Li et al 2015;Moukalled et al 2016;Liu et al 2019), which is a combination of the central difference, first-order and second-order upwind schemes. The MINMOD scheme not only achieves the second-order accuracy in space, but also preserves the physical properties of convection.…”

Section: Dimensionless Governing Equation and Numerical Schemementioning

confidence: 99%

Thermodynamically consistent modelling of two-phase flows with moving contact line and soluble surfactants

et al. 2019

Self Cite

View full text Add to dashboard Cite

Droplet dynamics on a solid substrate is significantly influenced by surfactants. It remains a challenging task to model and simulate the moving contact line dynamics with soluble surfactants. In this work, we present a derivation of the phase-field moving contact line model with soluble surfactants through the first law of thermodynamics, associated thermodynamic relations and the Onsager variational principle. The derived thermodynamically consistent model consists of two Cahn–Hilliard type of equations governing the evolution of interface and surfactant concentration, the incompressible Navier–Stokes equations and the generalized Navier boundary condition for the moving contact line. With chemical potentials derived from the free energy functional, we analytically obtain certain equilibrium properties of surfactant adsorption, including equilibrium profiles for phase-field variables, the Langmuir isotherm and the equilibrium equation of state. A classical droplet spread case is used to numerically validate the moving contact line model and equilibrium properties of surfactant adsorption. The influence of surfactants on the contact line dynamics observed in our simulations is consistent with the results obtained using sharp interface models. Using the proposed model, we investigate the droplet dynamics with soluble surfactants on a chemically patterned surface. It is observed that droplets will form three typical flow states as a result of different surfactant bulk concentrations and defect strengths, specifically the coalescence mode, the non-coalescence mode and the detachment mode. In addition, a phase diagram for the three flow states is presented. Finally, we study the unbalanced Young stress acting on triple-phase contact points. The unbalanced Young stress could be a driving or resistance force, which is determined by the critical defect strength.

show abstract

“…Before, lots of studies have been conducted to study the carbonate matrix acidizing process. Both experimental (Bernadiner et al, 1992;Wang et al, 1993;Frick et al, 1994;Fredd and Fogler, 1999;Izgec et al, 2008;Furui et al, 2012a;Furui et al, 2012b;Dong, 2012;Aidagulov et al, 2018;Qiu et al, 2018) and numerical (Kanaka and Panga, 2003;Panga et al, 2005;Kalia and Balakotaiah, 2007;Kalia and Balakotaiah, 2009;Maheshwari et al, 2013;Maheshwari et al, 2016b;Liu et al, 2019) results show that an optimum acid injection rate widely exists, at which the minimum amount of acid solution is consumed to achieve a breakthrough (Jia et al, 2021). Meanwhile, different dissolution patterns (Hoefner and Fogler, 1989;Fredd and Fogler, 1999;Huang et al, 2000;Izgec et al, 2008;Aidagulov et al, 2018) can also be observed through visualization techniques at different acid injection velocities.…”

Section: Introductionmentioning

confidence: 98%

Numerical Studies and Analyses on the Acidizing Process in Vug Carbonate Rocks

et al. 2021

Self Cite

View full text Add to dashboard Cite

A vug porosity system, in addition to a matrix, is also the target rocks for the acidizing. In this work, the acidizing process in two typical core-scale separate-vug porosity systems is studied in detail. Numerous cases are conducted to discuss a parametric study on the acidizing process and hydraulic behavior. Results indicate that the presence of vug reduces the pore volume of acid solution consumed to achieve a breakthrough, which is consistent with experimental observations. Increasing the vug diameter and porosity decreases pore volume to breakthrough both for a vugular carbonate rock and isolated vug carbonate rock. In comparison, the acid mass does not change a lot. Typical dissolution patterns can also be observed in the acidizing process when a vug exists. Compared to matrix dissolution patterns, the presence of vug induces wormhole to pass through the vug region.

show abstract

Modeling and analysis of the acidizing process in carbonate rocks using a two-phase thermal-hydrologic-chemical coupled model

Cited by 46 publications

References 74 publications

Three-dimensional simulation of acidizing process in carbonate rocks using the Darcy–Forchheimer framework

Three-dimensional simulation of acidizing process in carbonate rocks using the Darcy–Forchheimer framework

Thermodynamically consistent modelling of two-phase flows with moving contact line and soluble surfactants

Numerical Studies and Analyses on the Acidizing Process in Vug Carbonate Rocks

Contact Info

Product

Resources

About