Level set simulation of coupled advection‐diffusion and pore structure evolution due to mineral precipitation in porous media

Water Resources Research - Tập 44 Số 12 - 2008
Xiaoyi Li1, Hai Huang1, Paul Meakin1
1Idaho National Laboratory, Idaho Falls, Idaho, USA

Tóm tắt

A pore‐scale simulation technique based on level set interface tracking has been developed for modeling coupled reactive flow and structure evolution in porous media and fracture apertures. Advection, diffusion, and mineral precipitation resulting in changes in pore geometries are treated simultaneously by solving fully coupled fluid flow and reactive solute transport equations. The reaction‐induced evolution of solid grain surfaces is captured using a level set method, and a subgrid representation of the interface based on the level set approach is used instead of a pixel representation of the interface often used in cellular automata and lattice‐Boltzmann simulations. Precipitation processes within a 2‐D porous medium represented by nonoverlapping discs were simulated under various flow conditions and reaction rates, and the resulting changes of pore geometry are discussed. The simulation results indicate that under reaction‐limited conditions, precipitation is nearly uniform over the grain surfaces. However, this is no longer true when reaction is relatively fast and diffusion is the dominant transport process. In such cases, precipitation occurs mostly near the throat inlet and results in rapid permeability reduction with only a small reduction of porosity. In the case of fast reaction with transport dominated by advection (which is mostly likely in engineered remediation applications), solute can be delivered deep into fracture apertures and precipitation occurs mostly along preferential flow paths. Quantitative relationships between permeability and porosity under various flow conditions and reaction rates are also reported.

Từ khóa


Tài liệu tham khảo

10.1006/jcph.2002.7092

10.1016/j.jcp.2003.09.020

10.1016/0009-2509(95)00121-K

10.1016/S0013-7952(97)00044-6

10.1029/95WR03615

10.1016/0167-6636(82)90010-2

10.1103/PhysRevE.71.016311

10.2113/gselements.2.6.351

10.1029/2003GL017396

10.1029/2001WR000246

Economides M. J., 2000, Reservoir Stimulation

10.1006/jcph.1999.6236

10.1029/1998JB900029

10.1029/2005WR004204

10.1029/2005GL023899

10.1103/PhysRevE.65.036318

10.1029/2003JB002504

10.1029/2004GL021107

Kim S. G., 2002, Phase field modeling of directional eutectic solidification and comparison with experiments, Int. J. Cast Metals Res., 15, 241, 10.1080/13640461.2003.11819490

10.1016/j.jcrysgro.2003.08.078

10.1021/es025985u

10.1063/1.1844471

Lichtner P. C., 1996, Continuum formulation of multicomponent‐multiphase reactive transport, Reactive Transp. Porous Media, 34, 1

Mandelbrot B. B., 1984, The Fractal Geometry of Nature

10.1016/j.jcp.2005.03.031

Meakin P., 1998, Fractals, Scaling and Growth Far From Equilibrium

10.1080/00986449608936530

10.1063/1.1702607

10.1063/1.1713333

10.1016/0017-9310(91)90221-Y

10.1006/jcph.2000.6636

Press W. H., 1992, Numerical Recipes in FORTRAN: The Art of Scientific Computing

10.1144/jgs.157.1.239

10.1002/aic.690350713

10.1016/0009-2509(93)80031-K

10.1017/S0022112001004025

10.1146/annurev.fluid.35.101101.161105

10.1103/PhysRevLett.94.018501

10.1016/j.epsl.2005.09.017

10.1006/jcph.1994.1155

10.1016/S0045-7930(97)00053-4

10.1016/j.jcp.2006.08.013

10.1029/2005WR004770

10.1016/S0017-9310(02)00194-1

10.1006/jcph.1999.6294

10.1006/jcph.2001.6916

10.1071/EN05023

10.1103/PhysRevLett.47.1400

10.1029/2002JB001979

10.1016/j.apgeochem.2003.11.003

10.1016/j.jcp.2005.04.024

10.1029/95JB00958

Thông tin
Thông tin xuất bản

Water Resources Research

Tập 44 Số 12

Thông tin tác giả