TY - JOUR
T1 - Numerical modelling and analysis of reactive flow and wormhole formation in fractured carbonate rocks
AU - Liu, Piyang
AU - Yao, Jun
AU - Couples, Gary Douglas
AU - Huang, Zhaoqin
AU - Sun, Hai
AU - Ma, Jingsheng
PY - 2017/11
Y1 - 2017/11
N2 - Combining the two-scale continuum model and the discrete fracture network model, a continuum-based model is developed that calculates the reactive flow of acid in carbonate rock with a complex fracture network. The locations of fractures in the model are explicitly defined and the method can capture complex geometric relationships. The governing equations are discretized by the finite-volume method, where the diffusion term and convection term are discretized using the two-point flux approximation (TPFA) scheme and the upwind scheme, respectively. The physical domain is discretized by Delaunay triangulation. To keep the robustness and efficiency of the TPFA scheme, the optimization algorithm is used to move the centroid node of the control volume to improve the orthogonality of the grids. Numerical simulations of reactive flow in 2D fractured porous media, in cases with simple and complex fracture arrays, under linear and radial flow conditions, are presented. In particular, a sensitivity analysis of the dissolution process with respect to the presence of fractures, fracture aperture, fracture distribution, and acid injection rate, is conducted.
AB - Combining the two-scale continuum model and the discrete fracture network model, a continuum-based model is developed that calculates the reactive flow of acid in carbonate rock with a complex fracture network. The locations of fractures in the model are explicitly defined and the method can capture complex geometric relationships. The governing equations are discretized by the finite-volume method, where the diffusion term and convection term are discretized using the two-point flux approximation (TPFA) scheme and the upwind scheme, respectively. The physical domain is discretized by Delaunay triangulation. To keep the robustness and efficiency of the TPFA scheme, the optimization algorithm is used to move the centroid node of the control volume to improve the orthogonality of the grids. Numerical simulations of reactive flow in 2D fractured porous media, in cases with simple and complex fracture arrays, under linear and radial flow conditions, are presented. In particular, a sensitivity analysis of the dissolution process with respect to the presence of fractures, fracture aperture, fracture distribution, and acid injection rate, is conducted.
U2 - 10.1016/j.ces.2017.06.027
DO - 10.1016/j.ces.2017.06.027
M3 - Article
SN - 0009-2509
VL - 172
SP - 143
EP - 157
JO - Chemical Engineering Science
JF - Chemical Engineering Science
ER -