Multiscale modelling of dual-porosity porous media; A computational pore-scale study for flow and solute transport.

Abstract Many environmental and agricultural applications involve the transport of water and dissolved constituents through aggregated soil profiles, or porous media that are structured, fractured or macroporous in other ways. During the past several decades, various process-based macroscopic models have been used to simulate contaminant transport in such media. Many of these models consider advective-dispersive transport through relatively large inter-aggregate pore domains, while exchange with the smaller intra-aggregate pores is assumed to be controlled by diffusion. Exchange of solute between the two domains is often represented using a first-order mass transfer coefficient, which is commonly obtained by fitting to observed data. This study aims to understand and quantify the solute exchange term by applying a dual-porosity pore-scale network model to relatively large domains, and analysing the pore-scale results in terms of the classical dual-porosity (mobile-immobile) transport formulation. We examined the effects of key parameters (notably aggregate porosity and aggregate permeability) on the main dual-porosity model parameters, i.e., the mobile water fraction ( ϕ m ) and the mass transfer coefficient ( α ). Results were obtained for a wide range of aggregate porosities (between 0.082 and 0.700). The effect of aggregate permeability was explored by varying pore throat sizes within the aggregates. Solute breakthrough curves (BTCs) obtained with the pore-scale network model at several locations along the domain were analysed using analytical solutions of the dual-porosity model to obtain estimates of ϕ m and α . An increase in aggregate porosity was found to decrease ϕ m and increase α , leading to considerable tailing in the BTCs. Changes in the aggregate pore throat size affected the relative flow velocity between the intra- and inter-aggregate domains. Higher flow velocities within the aggregates caused a change in the transport regime from diffusion dominated to more advection dominated. This change increased the exchange rate of solutes between the mobile and immobile domains, with a related increase in the value of the mass transfer coefficient and less tailing in the BTCs.

[1]  Rainer Horn,et al.  Three-dimensional quantification of intra-aggregate pore-space features using synchrotron-radiation-based microtomography , 2008 .

[2]  D. Wildenschild,et al.  X-ray imaging and analysis techniques for quantifying pore-scale structure and processes in subsurface porous medium systems , 2013 .

[3]  Kenneth Stuart Sorbie,et al.  Representation of multiscale heterogeneity via multiscale pore networks , 2013 .

[4]  Veerle Cnudde,et al.  Multi-scale, micro-computed tomography-based pore network models to simulate drainage in heterogeneous rocks , 2015 .

[5]  E. Perfect,et al.  Effects of organic and inorganic fertilization on soil aggregation in an Ultisol as characterized by synchrotron based X-ray micro-computed tomography , 2013 .

[6]  Amir Raoof,et al.  A New Method for Generating Pore-Network Models of Porous Media , 2010 .

[7]  Simulation of chloride transport in an aggregated soil using three conceptual models , 2014, Arabian Journal of Geosciences.

[8]  Veerle Cnudde,et al.  Imaging and image-based fluid transport modeling at the pore scale in geological materials : a practical introduction to the current state-of-the-art , 2016 .

[9]  N. Ferro,et al.  Dual-scale micro-CT assessment of soil structure in a long-term fertilization experiment , 2013 .

[10]  M. T. van Genuchten,et al.  Some Exact Solutions for Solute Transport Through Soils Containing Large Cylindrical Macropores , 1984 .

[11]  M. Maraqa,et al.  Prediction of mass-transfer coefficient for solute transport in porous media. , 2001, Journal of contaminant hydrology.

[12]  Amir Raoof,et al.  Upscaling Transport of Adsorbing Solutes in Porous Media: Pore‐Network Modeling , 2010 .

[13]  M. V. Genuchten,et al.  Models for simulating salt movement in aggregated field soils , 1986 .

[14]  Hannes Flühler,et al.  SUSCEPTIBILITY OF SOILS TO PREFERENTIAL FLOW OF WATER : A FIELD STUDY , 1994 .

[15]  S. H. Mansouri,et al.  Insights into the role of wettability in cathode catalyst layer of proton exchange membrane fuel cell; pore scale immiscible flow and transport processes , 2017 .

[16]  M. Mariño,et al.  Interlayer Diffusive Transfer and Transport of Contaminants in Stratified Formation. I: Theory , 1998 .

[17]  S. M. Hassanizadeh,et al.  Saturation‐dependent solute dispersivity in porous media: Pore‐scale processes , 2013 .

[18]  M. Shao,et al.  Evaluation of a simple technique for estimating two-domain transport parameters , 1999 .

[19]  Jan M. Nordbotten,et al.  Effect of Mean Network Coordination Number on Dispersivity Characteristics , 2012, Transport in Porous Media.

[20]  Faisal Khan,et al.  3D simulation of the permeability tensor in a soil aggregate on basis of nanotomographic imaging and LBE solver , 2011, Journal of Soils and Sediments.

[21]  S. Sjögersten,et al.  The effect of soil aggregate size on pore structure and its consequence on emission of greenhouse gases , 2013 .

[22]  J. Barker,et al.  Block-geometry functions characterizing transport in densely fissured media , 1985 .

[23]  Emil O. Frind,et al.  Contaminant transport in fractured porous media: Analytical solutions for a system of parallel fractures , 1982 .

[24]  E. A. Sudicky,et al.  Contaminant transport in fractured porous media: Analytical solution for a single fracture , 1981 .

[25]  S. Hassanizadeh,et al.  Retention and remobilization of colloids during steady‐state and transient two‐phase flow , 2013 .

[26]  John Maximilian Köhne,et al.  A review of model applications for structured soils: a) Water flow and tracer transport. , 2009, Journal of contaminant hydrology.

[27]  Miguel A. Mariño,et al.  Interlayer Diffusive Transfer and Transport of Contaminants in Stratified Formation. II: Analytical Solutions , 1998 .

[28]  G. Matthews,et al.  Investigating the effects of wettability and pore size distribution on aggregate stability: the role of soil organic matter and the humic fraction , 2012 .

[29]  Peyman Mostaghimi,et al.  Insights into non-Fickian solute transport in carbonates , 2013, Water resources research.

[30]  Julia Frankfurter.  Conceptual Models Of Flow And Transport In The Fractured Vadose Zone , 2016 .

[31]  N. Jarvis.  A review of non‐equilibrium water flow and solute transport in soil macropores: principles, controlling factors and consequences for water quality , 2007 .

[32]  Pascale Royer,et al.  Contaminant transport in fractured porous media , 2003 .

[33]  S. H. Mansouri,et al.  Effects of Porosity and Water Saturation on the Effective Diffusivity of a Cathode Catalyst Layer , 2017 .

[34]  Veerle Cnudde,et al.  High-resolution X-ray computed tomography in geosciences: A review of the current technology and applications , 2013 .

[35]  Jirka Šimůnek,et al.  Modeling Nonequilibrium Flow and Transport Processes Using HYDRUS , 2008 .

[36]  Mark L. Brusseau,et al.  The influence of solute size, pore water velocity, and intraparticle porosity on solute dispersion and transport in soil , 1993 .

[37]  K. Abbaspour,et al.  Application of a two‐dimensional model to simulate flow and transport in a macroporous agricultural soil with tile drains , 2001 .

[38]  G. Bolt,et al.  Chapter 9: Movement of Solutes in Soil: Principles of Adsorption/Exchange Chromatography , 1979 .

[39]  Ayaz Mehmani,et al.  The effect of microporosity on transport properties in porous media , 2014 .

[40]  Adrian Sheppard,et al.  Imaged-based multiscale network modelling of microporosity in carbonates , 2014 .

[41]  S. Hassanizadeh,et al.  Oxidation of trichloroethylene, toluene, and ethanol vapors by a partially saturated permeable reactive barrier. , 2014, Journal of contaminant hydrology.

[42]  M. K. Shukla,et al.  Soil Physics: An Introduction , 2013 .

[43]  Tim Schmitz,et al.  Introduction To Environmental Soil Physics , 2016 .

[44]  M. V. Genuchten.  General approach for modeling solute transport in structured soils , 1985 .

[45]  Suzanne E. Allaire,et al.  Quantifying preferential flow in soils: A review of different techniques , 2009 .

[46]  Ľ. Lichner,et al.  Impact of long‐term mineral and organic fertilizer application on the water stability, wettability and porosity of aggregates obtained from two loamy soils , 2015 .

[47]  Horst H. Gerke,et al.  Preferential flow descriptions for structured soils , 2006 .

[48]  János Urai,et al.  Multi-scale characterization of porosity in Boom Clay (HADES-level, Mol, Belgium) using a combination of X-ray μ-CT, 2D BIB-SEM and FIB-SEM tomography , 2015 .

[49]  M. V. Genuchten,et al.  Mass transfer studies in sorbing porous media. I. Analytical solutions , 1976 .

[50]  K. H. Coats,et al.  Dead-End Pore Volume and Dispersion in Porous Media , 1964 .

[51]  A. Valocchi,et al.  Constraints on the Validity of Equilibrium and First-Order Kinetic Transport Models in Structured Soils , 1986 .

[52]  M. V. Genuchten,et al.  Review and comparison of models for describing non-equilibrium and preferential flow and transport in the vadose zone , 2003 .

[53]  Marco Bittelli.  Conceptual Models of Flow and Transport in the Fractured Vadose Zone. , 2002 .

[54]  Anders Rasmuson,et al.  Exact solution of a model for diffusion in particles and longitudinal dispersion in packed beds , 1980 .

[55]  A. Kravchenko,et al.  Soil aggregate geometry: Measurements and morphology , 2015 .

[56]  R. Yamamoto,et al.  Lattice-Boltzmann method combined with smoothed-profile method for particulate suspensions. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[57]  B. Wilson,et al.  Aggregate hierarchy and carbon mineralization in two Oxisols of New South Wales, Australia , 2015 .

[58]  Brent Clothier,et al.  In Situ Measurement of the Effective Transport Volume for Solute Moving Through Soil , 1992 .

[59]  J. Tisdall,et al.  Organic matter and water‐stable aggregates in soils , 1982 .

[60]  Muhammad Sahimi,et al.  Flow and Transport in Porous Media and Fractured Rock: From Classical Methods to Modern Approaches , 1995 .

[61]  M. Balhoff,et al.  A forward analysis on the applicability of tracer breakthrough profiles in revealing the pore structure of tight gas sandstone and carbonate rocks , 2015 .

[62]  Mark N. Goltz,et al.  Using the method of moments to analyze three-dimensional diffusion-limited solute transport from tem , 1987 .

[63]  J. Šimůnek,et al.  The effects of preferential flow and soil texture on risk assessments of a NORM waste disposal site. , 2010, Journal of hazardous materials.

[64]  S. M. Hassanizadeh,et al.  A new formulation for pore‐network modeling of two‐phase flow , 2012 .

[65]  Kazuo Hayashi,et al.  High Resolution X-Ray Computed Tomography. , 1993 .

[66]  T. Vogel,et al.  Using Oxygen‐18 to Study the Role of Preferential Flow in the Formation of Hillslope Runoff , 2010 .

[67]  Amir Raoof,et al.  PoreFlow: A complex pore-network model for simulation of reactive transport in variably saturated porous media , 2013, Comput. Geosci..