共查询到20条相似文献，搜索用时 15 毫秒
1.
A common way to simulate fluid flow in porous media is to use Lattice Boltzmann (LB) methods. Permeability predictions from such flow simulations are controlled by parameters whose settings must be calibrated in order to produce realistic modelling results. Herein we focus on the simplest and most commonly used implementation of the LB method: the single-relaxation-time BGK model. A key parameter in the BGK model is the relaxation time τ which controls flow velocity and has a substantial influence on the permeability calculation. Currently there is no rigorous scheme to calibrate its value for models of real media. We show that the standard method of calibration, by matching the flow profile of the analytic Hagen-Poiseuille pipe-flow model, results in a BGK-LB model that is unable to accurately predict permeability even in simple realistic porous media (herein, Fontainebleau sandstone). In order to reconcile the differences between predicted permeability and experimental data, we propose a method to calibrate τ using an enhanced Transitional Markov Chain Monte Carlo method, which is suitable for parallel computer architectures. We also propose a porosity-dependent τ calibration that provides an excellent fit to experimental data and which creates an empirical model that can be used to choose τ for new samples of known porosity. Our Bayesian framework thus provides robust predictions of permeability of realistic porous media, herein demonstrated on the BGK-LB model, and should therefore replace the standard pipe-flow based methods of calibration for more complex media. The calibration methodology can also be extended to more advanced LB methods. 相似文献
2.
In this Comment on the recent work (Zhu and Ma, 2013) [11] by Zhu and Ma (ZM) we first show that all three local gray Lattice Boltzmann (GLB) schemes in the form (Zhu and Ma, 2013) [11]: GS (Chen and Zhu, 2008; Gao and Sharma, 1994) [1,4], WBS (Walsh et al., 2009) [12] and ZM, fail to get constant Darcy’s velocity in series of porous blocks. This inconsistency is because of their incorrect definition of the macroscopic velocity in the presence of the heterogeneous momentum exchange, while the original WBS model (Walsh et al., 2009) [12] does this properly. We improve the GS and ZM schemes for this and other related deficiencies. Second, we show that the “discontinuous velocity” they recover on the stratified interfaces with their WBS scheme is inherent, in different degrees, to all LBE Brinkman schemes, including ZM scheme. None of them guarantees the stress and the velocity continuity by their implicit interface conditions, even in the frame of the two-relaxation-times (TRT) collision operator where these two properties are assured in stratified Stokes flow, Ginzburg (2007) [5]. Third, the GLB schemes are presented in work (Zhu and Ma, 2013) [11] as the alternative ones to direct, Brinkman-force based (BF) schemes (Freed, 1998; Nie and Martys, 2007) [3,8]. Yet, we show that the BF-TRT scheme (Ginzburg, 2008) [6] gets the solutions of any of the improved GLB schemes for specific, viscosity-dependent choice of its one or two local relaxation rates. This provides the principal difference between the GLB and BF: while the BF may respect the linearity of the Stokes–Brinkman equation rigorously, the GLB-TRT cannot, unless it reduces to the BF via the inverse transform of the relaxation rates. Furthermore, we show that, in limited parameter space, “gray” schemes may run one another. From the practical point of view, permeability values obtained with the GLB are viscosity-dependent, unlike with the BF. Finally, the GLB shares with the BF a so-called anisotropy (Ginzburg, 2008; Nie and Martys, 2007) [6,8], that is, flow-direction-dependency in their effective viscosity corrections, related to the discretized spatial variation of the resistance forcing. 相似文献
3.
《Advances in water resources》1998,22(3):273-284
A series of multi-step outflow experiments was carried out to identify the unsaturated hydraulic properties of two homogeneous coarse-textured porous media (glass beads and sand). Because of the measured sharp fronts of water content decrease during these experiments the hydraulic functions are assumed to be represented by the complete van Genuchten–Mualem closed-form expressions with variable coefficients α, n, m and θr. The values of θs and Ks were measured directly. A sensitivity analysis with respect to α, n, and m shows that conditions of local identifiability are satisfied if measurements of water content at some inner points inside the column are considered. The inverse modelling technique consists of two steps: first, computation of objective function values based on water content data responses to obtain initial parameter estimations, and second, a more detailed parameter determination using a Levenberg–Marquardt scheme. In both steps a numerical model incorporating the hydraulic functions is utilized to simulate theoretical pressure head and water content distributions along the column. For both porous media unique solutions of the inverse problem could be obtained, and afterwards, the corresponding hydraulic functions were verified from additional drainage experiments.©1998 Elsevier Science Limited. All rights reserved 相似文献
4.
High resolution images acquired from X-ray μ-CT are able to map the internal structure of porous media on which multiphase flow properties can be computed. While the resolution of a few micrometers may be sufficient for capturing the pore space of many sandstones, most carbonates exhibit a large amount of microporosity; pores which are below the image resolution and are not resolved at specific resolution. Neglecting the effect of micropores on fluid flow and transport properties of these rocks can cause erroneous results in particular at partial saturations. Current image-based pore scale models typically only consider macropores for simulating fluid flow. In this paper, we quantify the effect of microporosity on the effective permeability of the wetting phase for heterogeneous model structures with varying amount of micro-to-macro porosity. A multi-scale numerical approach is proposed to couple an average effect of micropores with an explicit representation of macropores. The Brinkman equation is solved using a lattice Boltzmann formulation to facilitate the coupling of Darcy and Stokes equations in micropores and macropores, respectively. The results show good agreement between the fine scale solution and the results of the upscaled models in which microporous regions are homogenised. The paper analyses in particular the choice of the momentum sink parameter at low wetting phase saturations. It is shown that this parameter can be found using either a flux-based calculation of permeability of microporous regions or chosen purely on the basis of the effective permeability of these regions. 相似文献
5.
6.
《Advances in water resources》1998,22(3):257-272
This work presents a rigorous numerical validation of analytical stochastic models of steady state unsaturated flow in heterogeneous porous media. It also provides a crucial link between stochastic theory based on simplifying assumptions and empirical field and simulation evidence of variably saturated flow in actual or realistic hypothetical heterogeneous porous media. Statistical properties of unsaturated hydraulic conductivity, soil water tension, and soil water flux in heterogeneous soils are investigated through high resolution Monte Carlo simulations of a wide range of steady state flow problems in a quasi-unbounded domain. In agreement with assumptions in analytical stochastic models of unsaturated flow, hydraulic conductivity and soil water tension are found to be lognormally and normally distributed, respectively. In contrast, simulations indicate that in moderate to strong variable conductivity fields, longitudinal flux is highly skewed. Transverse flux distributions are leptokurtic. the moments of the probability distributions obtained from Monte Carlo simulations are compared to modified first-order analytical models. Under moderate to strong heterogeneous soil flux conditions (σ2y≥1), analytical solutions overestimate variability in soil water tension by up to 40% as soil heterogeneity increases, and underestimate variability of both flux components by up to a factor 5. Theoretically predicted model (cross-)covariance agree well with the numerical sample (cross-)covarianaces. Statistical moments are shown to be consistent with observed physical characteristics of unsaturated flow in heterogeneous soils.©1998 Elsevier Science Limited. All rights reserved 相似文献
7.
In the current paper, the effect of the geometrical characteristics of 2-D porous media on the relative permeability in immiscible two-phase flows is studied. The generation of the different artificial porous media is performed using a Boolean model based on a random distribution of overlapping circles/ellipses, the size and shape of which are chosen to satisfy the specific Minkowski functionals (i.e. volume fraction, solid line contour length, connectivity). The study aims to identify how each different Minkowski functional affects the relative permeability of each phase at various saturations of the non-wetting phase. A 2-D multi-relaxation time (MRT) lattice Boltzmann model (LBM) that can handle high density ratios is employed in the simulation. The relationship between the driving forces G and the relative permeabilities of the two phases for every artificial structure is quantified. It is found that for high non-wetting phase saturations (fully connected flow), a non-linear relationship exists between the non-wetting phase flow rate and the driving force, whilst this relationship becomes linear at higher magnitudes of the latter. The force magnitude required to approach the linear region is highly influenced by the pore size distribution and the connectivity of the solid phase. For lower non-wetting phase saturation values, its relative permeability in the linear regime decreases as the fraction of small pores in the structure increases and the non-wetting phase flow becomes disconnected. A strong influence of the solid phase connectivity is also observed. 相似文献
8.
In this paper we present a stochastic model reduction method for efficiently solving nonlinear unconfined flow problems in heterogeneous random porous media. The input random fields of flow model are parameterized in a stochastic space for simulation. This often results in high stochastic dimensionality due to small correlation length of the covariance functions of the input fields. To efficiently treat the high-dimensional stochastic problem, we extend a recently proposed hybrid high-dimensional model representation (HDMR) technique to high-dimensional problems with multiple random input fields and integrate it with a sparse grid stochastic collocation method (SGSCM). Hybrid HDMR can decompose the high-dimensional model into a moderate M-dimensional model and a few one-dimensional models. The moderate dimensional model only depends on the most M important random dimensions, which are identified from the full stochastic space by sensitivity analysis. To extend the hybrid HDMR, we consider two different criteria for sensitivity test. Each of the derived low-dimensional stochastic models is solved by the SGSCM. This leads to a set of uncoupled deterministic problems at the collocation points, which can be solved by a deterministic solver. To demonstrate the efficiency and accuracy of the proposed method, a few numerical experiments are carried out for the unconfined flow problems in heterogeneous porous media with different correlation lengths. The results show that a good trade-off between computational complexity and approximation accuracy can be achieved for stochastic unconfined flow problems by selecting a suitable number of the most important dimensions in the M-dimensional model of hybrid HDMR. 相似文献
9.
Liquid distributions in unsaturated porous media under different gravitational accelerations and corresponding macroscopic gaseous diffusion coefficients were investigated to enhance understanding of plant growth conditions in microgravity. We used a single-component, multiphase lattice Boltzmann code to simulate liquid configurations in two-dimensional porous media at varying water contents for different gravity conditions and measured gas diffusion through the media using a multicomponent lattice Boltzmann code. The relative diffusion coefficients (D rel) for simulations with and without gravity as functions of air-filled porosity were in good agreement with measured data and established models. We found significant differences in liquid configuration in porous media, leading to reductions in D rel of up to 25% under zero gravity. The study highlights potential applications of the lattice Boltzmann method for rapid and cost-effective evaluation of alternative plant growth media designs under variable gravity. 相似文献
10.
With the rapid growth of nanotechnology industry, nanomaterials as an emerging pollutant are gradually released into subsurface environments and become great concerns. Simulating the transport of nanomaterials in groundwater is an important approach to investigate and predict the impact of nanomaterials on subsurface environments. Currently, a number of transport models are used to simulate this process, and the outputs of these models could be inconsistent with each other due to conceptual model uncertainty. However, the performances of different models on simulating nanoparticles transport in groundwater are rarely assessed in Bayesian framework in previous researches, and these will be the primary objective of this study. A porous media column experiment is conducted to observe the transport of Titanium Dioxide Nanoparticles (nano-TiO2). Ten typical transport models which consider different chemical reaction processes are used to simulate the transport of nano-TiO2, and the observed nano-TiO2 breakthrough curves data are used to calibrate these models. For each transport model, the parameter uncertainty is evaluated using Markov Chain Monte Carlo, and the DREAM(ZS) algorithm is used to sample parameter probability space. Moreover, the Bayesian model averaging (BMA) method is used to incorporate the conceptual model uncertainty arising from different chemical reaction based transport models. The results indicate that both two-sites and nonequilibrium sorption models can well reproduce the retention of nano-TiO2 transport in porous media. The linear equilibrium sorption isotherm, first-order degradation, and mobile-immobile models fail to describe the nano-TiO2 retention and transport. The BMA method could instead provide more reliable estimations of the predictive uncertainty compared to that using a single model. 相似文献
11.
The attenuation characteristics of the Kinnaur area of the North West Himalayas were studied using local earthquakes that occurred during 2008–2009. Most of the analyzed events are from the vicinity of the Panjal Thrust (PT) and South Tibetan Detachment Thrust, which are well-defined tectonic discontinuities in the Himalayas. The frequency-dependent attenuation of P and S waves was estimated using the extended coda normalization method. Data from 64 local earthquakes recorded at 10 broadband stations were used. The coda normalization of the spectral amplitudes of P and S waves was done at central frequencies of 1.5, 3, 6, 9, and 12 Hz. Q p increases from about 58 at 1.5 Hz to 706 at 12 Hz, and Q s increases from 105 at 1.5 Hz to 1,207 at 12 Hz. The results show that the quality factors for both P and S waves (Q p and Q s) increase as a function of frequency according to the relation Q?=?Q o f n , where Q o is the corresponding Q value at 1 Hz frequency and “n” is the frequency relation parameter. We obtained Q p?=?(47?±?2)f (1.04±0.04) and Q s?=?(86?±?4)f (0.96±0.03) by fitting power law dependency model for the estimated values of the entire study region. The Q 0 and n values show that the region is seismically very active and the crust is highly heterogeneous. There was no systematic variation of values of Q p and Q s at different frequencies from one tectonic unit to another. As a consequence, average values of these parameters were obtained for each frequency for the entire region, and these were used for interpretation and for comparison with worldwide data. Q p values lie within the range of values observed for some tectonically active regions of the world, whereas Q s values were the lowest among the values compared for different parts of the world. Q s/Q p values were >1 for the entire range of frequencies studied. All these factors indicate that the crust is highly heterogeneous in the study region. The high Q s/Q p values also indicate that the region is partially saturated with fluids. 相似文献
12.
《Advances in water resources》2005,28(8):835-844
Permeability of porous media in subsurface environments is subject to potentially large uncertainties due to the heterogeneity of natural systems. In this study, a first-order reliability method (FORM) is combined with a lattice Boltzmann method (LBM) to estimate the permeability of randomly generated porous media. The proposed procedure provides an increased ease of addressing complex pore structures by employing LBM to model fluid flow, while inheriting the computational efficiency from FORM. Macroscale-equivalent permeability can thus be estimated with significantly reduced computational efforts, while maintaining a connection to the complex microscale fluid dynamics within a pore structure environment. Implemented on several randomly generated porous media domains, the proposed method provides 13–120 times the efficiency compared to Monte Carlo methods. 相似文献
13.
《Advances in water resources》2002,25(6):601-609
A mirror-image method is proposed in this paper to solve the boundary conditions in the lattice Boltzmann model proposed by Zhang et al. [Adv. Water Resour. 25 (2002) 1] for the advection and anisotropic dispersion of solute transport in porous media. Three types of boundary are considered: prescribed concentration boundary, prescribed flux boundary and prescribed concentration-gradient boundary. The accuracy of the proposed method is verified against benchmark problems and finite difference method. 相似文献
14.
Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law–based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions. 相似文献
15.
We study two asymptotic regimes of unstable miscible displacements in porous media, in the two limits, where a permeability-modified aspect ratio, RL=L/H(kv/kh)1/2, becomes large or small, respectively. The first limit is known as transverse (or vertical) equilibrium, the second leads to the problem of non-communicating layers (the Dykstra–Parsons problem). In either case, the problem reduces to the solution of a single integro-differential equation. Although at opposite limits of the parameter RL, the two regimes coincide in the case of equal viscosities, M=1. By comparison with high-resolution simulation we investigate the validity of these two approximations. The evolution of transverse averages, particularly under viscous fingering conditions, depends on RL. We investigate the development of a model to describe viscous fingering in weakly heterogeneous porous media under transverse equilibrium conditions, and compare with the various existing empirical models (such as the Koval, Todd–Longstaff and Fayers models). 相似文献
16.
This paper presents a stochastic model for multicomponent competitive monovalent cation exchange in hierarchical porous media. Reactive transport in porous media is highly sensitive to heterogeneities in physical and chemical properties, such as hydraulic conductivity (K), and cation exchange capacity (CEC). We use a conceptual model for multimodal reactive mineral facies and develop a Eulerian-based stochastic theory to analyze the transport of multiple cations in heterogeneous media with a hierarchical organization of reactive minerals. Numerical examples investigate the retardation factors and dispersivities in a chemical system made of three monovalent cations (Na+, K+, and Cs+). The results demonstrate how heterogeneity influences the transport of competitive monovalent cations, and highlight the importance of correlations between K and CEC. Further sensitivity analyses are presented investigating how the dispersion and retardation of each cation are affected by the means, variances, and integral scales of K and CEC. The volume fraction of organic matter is shown to be another important parameter. The Eulerian stochastic framework presented in this work clarifies the importance of each system parameters on the migration of cation plumes in formations with hierarchical organization of facies types. Our stochastic approach could be used as an alternative to numerical simulations for 3D reactive transport in hierarchical porous media, which become prohibitively expensive for the multicomponent applications considered in this work. 相似文献
17.
《Advances in water resources》2004,27(10):1017-1032
This paper presents a numerical solution for the effective conductivity of a periodic binary medium with cuboid inclusions located on an octahedral lattice. The problem is defined by five dimensionless geometric parameters and one dimensionless conductivity contrast parameter. The effective conductivity is determined by considering the flow through the “elementary flow domain” (EFD), which is an octant of the unitary domain of the periodic media. We derive practical bounds of interest for the six-dimensional parameter space of the EFD and numerically compute solutions at regular intervals throughout the entire bounded parameter space. A continuous solution of the effective conductivity within the limits of the simulated parameter space is then obtained via interpolation of the numerical results. Comparison to effective conductivities derived for random heterogeneous media demonstrate similarities and differences in the behavior of the effective conductivity in regular periodic (low entropy) vs. random (high entropy) media. The results define the low entropy bounds of effective conductivity in natural media, which is neither completely random nor completely periodic, over a large range of structural geometries. For aniso-probable inclusion spacing, the absolute bounds of Keff for isotropic inclusions are the Wiener bounds, not the Hashin-Shtrikman bounds. For isotropic inclusion and isoprobable conditions well below the percolation threshold, the results are in agreement with the self-consistent approach. For anisotropic cuboid inclusions, or at relatively close spacing in at least one direction (p > 0.2) (aniso-probable conditions), the effective conductivity of the periodic media is significantly different from that found in anisotropic random binary or Gaussian media. 相似文献
18.
In this paper the coupled equations governing the dynamic behavior of unsaturated soils are derived based on the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) [Gatmiri B. Analysis of fully coupled behavior of unsaturated porous medium under stress, suction and temperature gradient. Final report of CERMES-EDF, 1997] and Gatmiri et al. (1998) [Gatmiri B, Delage P, Cerrolaza M, UDAM: a powerful finite element software for the analysis of unsaturated porous media. Adv Eng Software 1998; 29(1): 29–43]. In this formulation, the solid skeleton displacements, water pressure and air pressure are presumed to be independent variables. The Boundary Integral formulations as well as fundamental solutions for such a dynamic u–pw–pa theory are presented in this paper for the first time. The boundary integral equations are derived via the use of the weighted residuals method in a way that permits an easy discretization and implementation in a Boundary Element code. Also, the associated two dimensional (2D) fundamental solutions for such deformable porous medium with linear elastic behavior are derived in Laplace transform domain using the method of Hörmander. Finally, some numerical results are presented to show the accuracy of the proposed solutions. The derived results are verified analytically by comparison with the previously introduced corresponding fundamental solutions in elastodynamic limiting case. 相似文献
19.
Mingjie Chen Zhiming Lu George A. Zyvoloski 《Stochastic Environmental Research and Risk Assessment (SERRA)》2008,22(4):587-596
This study is an extension of the stochastic analysis of transient two-phase flow in randomly heterogeneous porous media (Chen
et al. in Water Resour Res 42:W03425, 2006), by incorporating direct measurements of the random soil properties. The log-transformed intrinsic permeability, soil pore
size distribution parameter, and van Genuchten fitting parameter are treated as stochastic variables that are normally distributed
with a separable exponential covariance model. These three random variables conditioned on given measurements are decomposed
via Karhunen–Loève decomposition. Combined with the conditional eigenvalues and eigenfunctions of random variables, we conduct
a series of numerical simulations using stochastic transient water–oil flow model (Chen et al. in Water Resour Res 42:W03425,
2006) based on the KLME approach to investigate how the number and location of measurement points, different random soil properties,
as well as the correlation length of the random soil properties, affect the stochastic behavior of water and oil flow in heterogeneous
porous media. 相似文献
20.
Two-dimensional numerical simulations of two-phase (DNAPL-water) flow in spatially correlated random fields demonstrate the influence of nonwetting phase (NWP) relative permeability–saturation (kr,N–SW) relationships correlated to porous media intrinsic permeability (k). Both the volume of porous media invaded by the NWP and the length of time during which the NWP is migrating are under predicted if kr,N–k correlation is not accounted for in the model formulation. Not accounting for the kr,N–k correlation resulted in under predicting the volume of porous media invaded by up to approximately 10%, which is likely not significant for many practical applications. However, not accounting for the kr,N–k correlation resulted in under predicting field scale migration times by up to a factor of 4, which is likely significant in that the migration times are on the order of years to several decades for the DNAPL (1,2-DCE) considered in this study. The under prediction of migration times was greater for lower permeability aquifers. 相似文献