3D simulation of the permeability tensor in a soil aggregate on basis of nanotomographic imaging and LBE solver
Tóm tắt
The purpose of this hydropedological work is to investigate the feasibility of a method to calculate permeability of soil peds on the matrix pore scale resolution. This paper focuses on imaging of the intra-aggregate microstructure and, based on the three-dimensional (3D) images, quantification of the pore network connectivity and permeability tensor. Finally, lattice Boltzmann equation (LBE) simulations of Navier–Stokes flow in the thus derived pore network allow to compute the heterogeneous 3D flow velocity field. Nanotomographic X-ray absorption mode imaging of a single soil ped has been performed at the TOMCAT beamline of the Swiss Light Source synchrotron facility with 0.74 μm spatial resolution. Segmentation of the 3D Nano-CT images into solid phase and pore space allowed to study the statistical properties of the connected pore network. The thus derived pore network data were used as direct input for the software package GeoDict integrated with a LBE algorithm to perform saturated water flow modeling. The soil ped features quantified from the tomographic images were pore and grain size distributions (PSD, GSD), porosity, percolation tensor in x, y, z-direction and percolation pathways (macropores). The PSD frequency has a peak in the range 3–8 μm (mesopores) contributing 50% of the total. In general, the matrix pore structure and, hence, saturated flow field velocity of our sample is highly anisotropic. LBE-simulated pore scale fluid flow is used to ultimately determine gross parameters (coefficients) of Darcy’s law such as the saturated hydraulic conductivity, K
s. The thus simulated K
s = 105 ± 24 cm day−1 of a single soil ped (mean of three ROIs and directions) is one order of magnitude higher than the value of K
s = 12 ± 2 cm day−1 predicted by a classical pedotranfer function approach. This K
s underestimation by PTF is typical for macropore flow. Integration of 3D image evaluation with the LBE approach as an essential step toward understanding the highly heterogeneous intra-aggregate microstructure is now possible on the submicron scale. Modeling the resulting anisotropic fluid advection field is, however, just a first step for 3D models of biogeochemical reactions at soil interfaces on this pore scale.
Tài liệu tham khảo
Bartoli F, Bird NRA, Gomendy V, Vivier H, Niquet S (1999) The relation between silty soil structures and their mercury porosimetry curve counterparts: fractals and percolation. Eur J Soil Sci 50:9–22
Baveye PC, Laba M, Otten W, Bouckaert L, Sterpaio PD, Goswami RR, Grinev D, Houston A, Hu Y, Liu J, Mooney S, Pajor R, Sleutel S, Tarquis A, Wang W, Wei Q, Sezgin M (2010) Observer-dependent variability of the thresholding step in the quantitative analysis of soil images and X-ray microtomography data. Geoderma 157:51–63
Bronnikov AV (2002) Theory of quantitative phase-contrast computed tomography. J Opt Soc Am A 19:472–480
Calboun PS, Duszyk BS, Heath DG, Carley JC, Fishman EK (1999) Three-dimensional volume rendering of spiral CT data: theory and method. RadioGraphics 19:745–764
Campbell GS (1985) Soil physics with BASIC—transport models for soil–plant systems. Developments in Soil Science, vol 14. Elsevier, Amsterdam, pp 150
Carminati A, Kaestner A, Flühler H, Lehmann P, Or D, Lehmann E, Stampanoni M (2007) Hydraulic contacts controlling water flow across grains. Phys Rev E 76:026311. doi:10.1103/PhysRevE.76.026311
Culligan KA, Wildenschild D, Christensen BSB, Gray WG, Rivers ML (2006) Pore-scale characteristics of multiphase flow in porous media: a synchrotron-based CMT comparison of air-water and oil-water experiments. Adv Water Resour 29:227–238
Fourie W, Said R, Young P, Barnes DL (2007) The simulation of pore scale fluid flow with real world geometries obtained from X-ray computed tomography. Proc COMSOL Conf 2007, Newton, MA, USA, pp 6 (http://www.comsol.com/papers/3039/download/Fourie.pdf)
Fox AL, Eisenhauer DE, Dosskey MG (2005) Modeling water and sediment trapping by vegetated filters using VFSMOD: comparing methods for estimating infiltration parameters. ASAE Pap No. 052118. http://abe.ufl.edu/carpena/files/pdf/software/vfsmod/ModelingWaterAndSedimentTrapping.pdf
Gallagher NC Jr, Wise GL (1981) A theoretical analysis of the properties of median filters. IEEE Trans ASSP 29:1136–1141
Ginzburg I, d’Humiéres D (2003) Multireflection boundary conditions for lattice Boltzmann models. Phys Rev E 68:066614. doi:10.1103/PhysRevE.68.066614
Ginzburg I, Steiner K (2003) Lattice Boltzmann model for free-surface and its applications to filling processes in casting. J Comput Phys 185:61–99
Graf von der Schulenburg DA, Pintelon TRR, Picioreanu C, Van Loosdrecht MCM, Johns ML (2009) Three-dimensional simulations of biofilm growth in porous media. AICHE J 55:494–504
Guber AK, Pachepsky YA (2010) Multimodeling with pedotransfer functions, documentation and user manual for PTF calculator (CalcPTF), version 2.0. Environmental Microbial and Food Safety Laboratory, Hydrology and Remote Sensing Laboratory, Beltsville Agricultural Research Center, USDA-ARS. http://www.ars.usda.gov/Research/docs.htm?docid=21211
Guber AK, Pachepsky YA, Shein E, Rawls WJ (2004) Soil aggregates and water retention. Dev Soil Sci 30:143–151
Humières DD, Ginzburg I, Krafczyk M, Lallemand P, Luo LS (2002) Multiple-relaxation-time lattice Boltzmann models in three dimensions. Philos Trans R Soc Lond A 360:437–451
Inamuro T, Yoshino M, Ogino F (1999) Lattice Boltzmann simulation of flows in a three-dimensional porous structure. Int J Numer Meth Fluids 29:737–748
Kaestner A, Lehmann E, Stampanoni M (2008) Imaging and image processing in porous media research. Adv Wat Resour 31:1174–1187
Knutson CE, Werth CJ, Valocchi AJ (2001) Pore-scale modeling of dissolution from variably distributed nonaqueous phase liquid blobs. Water Resour Res 37:2951–2963
Kutilek M, Nielsen DR (2007) Interdisciplinary of hydropedology. Geoderma 138:252–260
Marone F, Hintermüller C, McDonald S, Abela R, Mikuljan G, Stampanoni M (2009) X-ray tomographic microscopy at TOMCAT. J Phys Conf Ser 186:012042. doi:10.1088/1742-6596/186/1/012042
Marone F, Münch B, Stampanoni M (2010) Fast reconstruction algorithm dealing with tomography artifacts. In Stock SR (ed) Developments in X-ray tomography VII. SPIE Proc 7804. doi:10.1117/12.859703
McDonald SA, Marone F, Hintermüller C, Mikuljan G, David C, Stampanoni M (2011) Phase contrast X-ray tomographic microscopy for biological and materials science applications. Adv Eng Mater 13:116–121
Murshed MM, Klapp SA, Szeder T, Huthwelker T, Stampanoni M, Marone F, Hintermüller C, Bohrmann G, Kuhs WF, Kersten M (2008) Natural gas hydrate investigations by synchrotron radiation X-ray cryo-tomographic microscopy (SRXCTM). Geophys Res Lett 35:L23612. doi:10.1029/2008GL035460
Pachepsky Y, Van Genuchten M (2011) Pedotransfer functions. In: Glinski J, Horabik J, Lipiec J (eds) Encyclopedia of agrophysics. Springer, Berlin
Pal NR, Pal SK (1993) A review of image segmentation techniques. Pattern Recogn 29:1277–1294
Pan C, Luo LS, Miller CT (2006) An evaluation of lattice Boltzmann schemes for porous medium flow simulation. Comput Fluids 35:898–909
Peth S (2010) Chapter 3—Applications of microtomography in soils and sediments. Developments in Soil Science, vol. 34. Elsevier, Amsterdam, pp 73–101
Peth S, Horn R, Beckmann F, Donath T, Fischer J, Smucker AJM (2008) Three-dimensional quantification of intra-aggregate pore-space features using synchrotron-radiation-based microtomography. Soil Sci Soc Am J 72:897–907
Pfrang A, Schladitz K, Wiegmann A, Schimmel TH (2007) Calculation of the evolution of surface area and free volume during the infiltration of fiber felts. Chem Vap Deposition 13:705–715
Porter ML, Wildenschild D (2010) Image analysis algorithms for estimating porous media multiphase flow variables from computed microtomography data: a validation study. Comput Geosci 14:15–30
Rennert T, Totsche KU, Heister K, Kersten M, Thieme J (2011) Advanced spectroscopic, microscopic and tomographic characterization techniques to study biogeochemical interfaces in soil. J Soils Sediments. doi:10.1007/s11368-011-0417-5
San José Martínez F, Martín MA, Caniego FJ, Tuller M, Guber A, Pachepsky Y, García-Gutiérrez C (2010) Multifractal statistics of discretized X-ray CT visualizations for the characterization of soil macropore structures. Geoderma 156:32–42
Sheppard AP, Sok RM, Averdunk H (2004) Techniques for image enhancement and segmentation of tomographic images of porous materials. Phys A Stat Mech Appl 339:145–151
Silin D, Tomutsa L, Benson SM, Patzek TW (2011) Microtomography and pore-scale modeling of two-phase fluid distribution. Transp Porous Med 86:495–515
Stampanoni M, Grošo A, Isenegger A, Mikuljan G, Chen Q, Bertrand A, Henein S, Betemps R, Frommherz U, Böhler P, Meister D, Lange M, Abela R (2006) Developments in X-ray tomography V. Proc SPIE 638:63180M1-14
Succi S (2001) The lattice Boltzmann equation: for fluid dynamics and beyond. Clarendon, Gloucestershire
Thoemmes G, Becker J, Junk M, Vaikuntam AK, Kehrwald D, Klar A, Steiner K, Wiegmann A (2009) A lattice Boltzmann method for immiscible multiphase flow simulations using the level set method. J Comput Phys 228:1139–1156
Tietje O, Hennings V (1996) Accuracy of the saturated hydraulic conductivity prediction by pedo-transfer functions compared to the variability within FAO textural classes. Geoderma 69:71–84
Tippkötter R, Eickhorst T, Taubner H, Greder B, Rademaker G (2009) Detection of soil water in macropores of undistrurbed soil using microfocus X-ray tube computerized tomography (mCT). Soil Tillage Res 105:12–20
Wildenschild D, Hopmans JW, Rivers ML, Kent AJR (2005) Quantitative analysis of flow processes in a sand using synchrotron-based X-ray microtomography. Vadose Zone J 4:112–126
Zhang X, Bengough AG, Deeks LK, Crawford JW, Young IM (2002) A novel three-dimensional lattice Boltzmann model for solute transport in variably saturated porous media. Water Resour Res 38:1167. doi:10.1029/2001WR000982