A meshless method to simulate solute transport in heterogeneous porous media

We derive a meshless numerical method based on smoothed particle hydrodynamics (SPH) for the simulation of conservative solute transport in heterogeneous geological formations. We demonstrate that the new proposed scheme is stable, accurate, and conserves global mass. We evaluate the performance of the proposed method versus other popular numerical methods for the simulation of one- and two-dimensional dispersion and two-dimensional advective–dispersive solute transport in heterogeneous porous media under different Pèclet numbers. The results of those benchmarks demonstrate that the proposed scheme has important advantages over other standard methods because of its natural ability to control numerical dispersion and other numerical artifacts. More importantly, while the numerical dispersion affecting traditional numerical methods creates artificial mixing and dilution, the new scheme provides numerical solutions that are “physically correct”, greatly reducing these artifacts.     

Variability of volume strain in bounded heterogeneous media

The deformation of the solid matrix affects the fluid pore pressure and flow by altering the pore volume. Such interaction in turn affects the storage of groundwater in the void space. Obviously, this subject is of interest in groundwater hydrology. This paper describes an investigation of the effect of aquifer heterogeneity on the variability of the fluid pressure head and solid's volume strain, where the assumption of a constant vertical total stress leads to a relatively simple relationship between changes in solid's volume strain and fluid pressure head. To solve the problem analytically, focus is placed on the one‐dimensional models. It is found from our closed‐form solutions that the variance and correlation length of the log hydraulic conductivity are important in increasing the variability of pressure head and solid's volume strain. It is hoped that our findings will provide a basic framework for understanding and quantifying field‐scale volume strain processes and be useful in stimulating further research in this area. Copyright © 2011 John Wiley & Sons, Ltd.     

Density-dependent dispersion in heterogeneous porous media Part I: A numerical study

In this paper, we describe carefully conducted numerical experiments, in which a dense salt solution vertically displaces fresh water in a stable manner. The two-dimensional porous media are weakly heterogeneous at a small scale. The purpose of these simulations, conducted for a range of density differences, is to obtain accurate concentration profiles that can be used to validate nonlinear models for high-concentration-gradient dispersion. In this part we focus on convergence of the computations, in numerical and statistical sense, to ensure that the uncertainty in the results is small enough.Concentration variances are computed, which give estimates of the uncertainty in local concentration values. These local variations decrease with increasing density contrast. For tracer transport, obtained longitudinal dispersivities are in accordance with analytical findings. In the case of high-density contrasts, stabilizing gravity forces counteract the growth of dispersive fingers, decreasing the effective width of the transition zone. For small log-permeability variances, the decrease of the apparent dispersivity that is found is in agreement with laboratory results for homogeneous columns.     

峨眉山超级地幔柱对四川盆地烃源岩热演化的影响

四川盆地的构造、热演化与峨眉山超级地幔柱有密切关系.峨眉山超级地幔柱对四川盆地中二叠统之下的烃源岩热演化有着十分重要的影响.在四川盆地热历史恢复的基础上,研究了峨眉山超级地幔柱对盆地内烃源岩,特别是中二叠统之下的古生界烃源岩热演化的影响.结果表明,中二叠统及下伏烃源岩的热演化受中晚二叠世发生在盆地西南方向的峨眉山超级地幔柱的影响巨大,且具地区差异性.即在靠近峨眉山地幔柱中心的地区,有机质迅速成熟并达到其成熟度的最高值(以H1井为代表),古生界烃源岩迅速进入过成熟,此后未有二次生烃;而远离峨眉山地幔柱的盆地大部分地区,古生界烃源岩在二叠纪以来具有多次生烃过程.中生界烃源岩热演化,主要和前陆盆地阶段的构造过程包括前陆沉积和断裂的逆冲推覆等相关. 在烃源岩有机质成熟度演化史的基础上,从盆地热史和烃源岩热演化的角度指出了下二叠统及之下烃源层在四川盆地不同地区油气勘探中的不同意义.     

Stochastic modeling of solute transport in 3-D heterogeneous porous media with random source condition

During probabilistic analysis of flow and transport in porous media, the uncertainty due to spatial heterogeneity of governing parameters are often taken into account. The randomness in the source conditions also play a major role on the stochastic behavior in distribution of the dependent variable. The present paper is focused on studying the effect of both uncertainty in the governing system parameters as well as the input source conditions. Under such circumstances, a method is proposed which combines with stochastic finite element method (SFEM) and is illustrated for probabilistic analysis of concentration distribution in a 3-D heterogeneous porous media under the influence of random source condition. In the first step SFEM used for probabilistic solution due to spatial heterogeneity of governing parameters for a unit source pulse. Further, the results from the unit source pulse case have been used for the analysis of multiple pulse case using the numerical convolution when the source condition is a random process. The source condition is modeled as a discrete release of random amount of masses at fixed intervals of time. The mean and standard deviation of concentration is compared for the deterministic and the stochastic system scenarios as well as for different values of system parameters. The effect of uncertainty of source condition is also demonstrated in terms of mean and standard deviation of concentration at various locations in the domain.     

Comparison of theory and experiment for solute transport in weakly heterogeneous bimodal porous media

In this work, the influence of non-equilibrium effects on solute transport in a weakly heterogeneous medium is discussed. Three macro-scale models (upscaled via the volume averaging technique) are investigated: (i) the two-equation non-equilibrium model, (ii) the one-equation asymptotic model and (iii) the one-equation local equilibrium model. The relevance of each of these models to the experimental system conditions (duration of the pulse injection, dispersivity values…) is analyzed. The numerical results predicted by these macroscale models are compared directly with the experimental data (breakthrough curves). Our results suggest that the preasymptotic zone (for which a non-Fickian model is required) increases as the solute input pulse time decreases. Beyond this limit, the asymptotic regime is recovered. A comparison with the results issued from the stochastic theory for this regime is performed. Results predicted by both approaches (volume averaging method and stochastic analysis) are found to be consistent.     

Comparison of theory and experiments for dispersion in homogeneous porous media

Modeling dispersion in homogeneous porous media with the convection–dispersion equation commonly requires computing effective transport coefficients. In this work, we investigate longitudinal and transverse dispersion coefficients arising from the method of volume averaging, for a variety of periodic, homogeneous porous media over a range of particle Péclet (Pe_p) numbers. Our objective is to validate the upscaled transverse dispersion coefficients and concentration profiles by comparison to experimental data reported in the literature, and to compare the upscaling approach to the more common approach of inverse modeling, which relies on fitting the dispersion coefficients to measured data. This work is unique in that the exact microscale geometry is available; thus, no simplifying assumptions regarding the geometry are required to predict the effective dispersion coefficients directly from theory. Transport of both an inert tracer and non-chemotactic bacteria is investigated for an experimental system that was designed to promote transverse dispersion. We highlight the occurrence of transverse dispersion coefficients that (1) depart from power-law behavior at relatively low Pe_p values and (2) are greater than their longitudinal counterparts for a specific range of Pe_p values. The upscaling theory provides values for the transverse dispersion coefficient that are within the 98% confidence interval of the values obtained from inverse modeling. The mean absolute error between experimental and upscaled concentration profiles was very similar to that between the experiments and inverse modeling. In all cases the mean absolute error did not exceed 12%. Overall, this work suggests that volume averaging can potentially be used as an alternative to inverse modeling for dispersion in homogeneous porous media.     

Application of stress-pore pressure coupling theory for porous media to the Xinfengjiang reservoir earthquakes

Theory of the coupling of stress-pore pressure in the saturated, elastic porous media is used in the study of the formation mechanism of the Xinfengjiang reservoir-induced earthquakes. Based on the results, it is believed that compared with the mechanism of additional stress in the vicinity of the reservoir, the mechanism of the coupling of additional stress and pore pressure may be more well-founded for the occurrence of reservoir-induced earthquakes.     

An extended stability criterion for density-driven flows in homogeneous porous media

Stability of density-driven flows is a challenging problem with current applications in major areas like energy exploration, water pollution, nuclear and oil industries. The mathematical model for such flows is a system of coupled non linear partial differential equations. To study the physical stability of the system, we consider steady-state flow and perturb the solution of the full system of equations (without Boussinesq approximation) and investigate how it evolves in time: if the solution does not grow indefinitely, the system is called stable. The perturbations are treated as being the result of sub-scale interactions between the velocity field and the solute mass. Making use of a two-scale expansion of the solution, we derived extended stability criteria that include the effects of density, viscosity and flow velocity in flow configurations aligned parallel as well as orthogonal to gravity forces. Numerical simulations with the numerical simulator d³f$d^{3}f$ are presented to test the theoretical stability criteria.     

Modeling the effects of nonlinear equilibrium sorption on the transport of solute plumes in saturated heterogeneous porous media

Transport of sorbing solutes in 2D steady and heterogeneous flow fields is modeled using a particle tracking random walk technique. The solute is injected as an instantaneous pulse over a finite area. Cases of linear and Freundlich sorption isotherms are considered. Local pore velocity and mechanical dispersion are used to describe the solute transport mechanisms at the local scale. This paper addresses the impact of the degree of heterogeneity and correlation lengths of the log-hydraulic conductivity field as well as negative correlation between the log-hydraulic conductivity field and the log-sorption affinity field on the behavior of the plume of a sorbing chemical. Behavior of the plume is quantified in terms of longitudinal spatial moments: center-of-mass displacement, variance, 95% range, and skewness. The range appears to be a better measure of the spread in the plumes with Freundlich sorption because of plume asymmetry. It has been found that the range varied linearly with the travelled distance, regardless of the sorption isotherm. This linear relationship is important for extrapolation of results to predict behavior beyond simulated times and distances. It was observed that the flow domain heterogeneity slightly enhanced the spreading of nonlinearly sorbing solutes in comparison to that which occurred for the homogeneous flow domain, whereas the spreading enhancement in the case of linear sorption was much more pronounced. In the case of Freundlich sorption, this enhancement led to further deceleration of the solute plume movement as a result of increased retardation coefficients produced by smaller concentrations. It was also observed that, except for plumes with linear sorption, correlation between the hydraulic conductivity and the sorption affinity fields had minimal effect on the spatial moments of solute plumes with nonlinear sorption.     

含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究

本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.     

Fluid dispersion effects on density-driven thermohaline flow and transport in porous media

This study introduces the dispersive fluid flux of total fluid mass to the density-driven flow equation to improve thermohaline modeling of salt and heat transports in porous media. The dispersive fluid flux in the flow equation is derived to account for an additional fluid flux driven by the density gradient and mechanical dispersion. The coupled flow, salt transport and heat transport governing equations are numerically solved by a fully implicit finite difference method to investigate solution changes due to the dispersive fluid flux. The numerical solutions are verified by the Henry problem and the thermal Elder problem under a moderate density effect and by the brine Elder problem under a strong density effect. It is found that increment of the maximum ratio of the dispersive fluid flux to the advective fluid flux results in increasing dispersivity for the Henry problem and the brine Elder problem. The effects of the dispersive fluid flux on salt and heat transports under high density differences and high dispersivities are more noticeable than under low density differences and low dispersivities. Values of quantitative indicators such as the Nusselt number, mass flux, salt mass stored and maximum penetration depth in the brine Elder problem show noticeable changes by the dispersive fluid flux. In the thermohaline Elder problem, the dispersive fluid flux shows a considerable effect on the shape and the number of developed fingers and makes either an upwelling or a downwelling flow in the center of the domain. In conclusion, for the general case that involves strong density-driven flow and transport modeling in porous media, the dispersive fluid flux should be considered in the flow equation.     

Density-dependent dispersion in heterogeneous porous media Part II: Comparison with nonlinear models

The results of a series of high-resolution numerical experiments are used to test and compare three nonlinear models for high-concentration-gradient dispersion. Gravity stable miscible displacement is considered. The first model, introduced by Hassanizadeh, is a modification of Fick’s law which involves a second-order term in the dispersive flux equation and an additional dispersion parameter β. The numerical experiments confirm the dependency of β on the flow rate. In addition, a dependency on travelled distance is observed. The model can successfully be applied to nearly homogeneous media (σ² = 0.1), but additional fitting is required for more heterogeneous media.The second and third models are based on homogenization of the local scale equations describing density-dependent transport. Egorov considers media that are heterogeneous on the Darcy scale, whereas Demidov starts at the pore-scale level. Both approaches result in a macroscopic balance equation in which the dispersion coefficient is a function of the dimensionless density gradient. In addition, an expression for the concentration variance is derived. For small σ², Egorov’s model predictions are in satisfactory agreement with the numerical experiments without the introduction of any new parameters. Demidov’s model involves an additional fitting parameter, but can be applied to more heterogeneous media as well.     

Bayesian uncertainty quantification for flows in heterogeneous porous media using reversible jump Markov chain Monte Carlo methods

In this paper, we study the uncertainty quantification in inverse problems for flows in heterogeneous porous media. Reversible jump Markov chain Monte Carlo algorithms (MCMC) are used for hierarchical modeling of channelized permeability fields. Within each channel, the permeability is assumed to have a log-normal distribution. Uncertainty quantification in history matching is carried out hierarchically by constructing geologic facies boundaries as well as permeability fields within each facies using dynamic data such as production data. The search with Metropolis–Hastings algorithm results in very low acceptance rate, and consequently, the computations are CPU demanding. To speed-up the computations, we use a two-stage MCMC that utilizes upscaled models to screen the proposals. In our numerical results, we assume that the channels intersect the wells and the intersection locations are known. Our results show that the proposed algorithms are capable of capturing the channel boundaries and describe the permeability variations within the channels using dynamic production history at the wells.     

Generation of three dimensional flow fields for statistically anisotropic heterogeneous porous media

A methodology for generating three dimensional (3D) flow fields for statistically anisotropic heterogeneous porous media is presented and demonstrated. The simulated flow fields are shown to exhibit the input spatial correlation structure and observe mass continuity. Sample flow fields are presented in the form of cross sectional slices of the 3D formation. These cross sections demonstrate visually the characteristics of subsurface flow. The method was found to be faster than traditional techniques in terms of its computational requirements. Given this method, it is possible to generate the large number of realizations of a velocity field necessary to compute high order statistics in transport problems.     

An experimental study of dispersion of an interactive tracer in chemically heterogeneous medium at a column scale

This paper aims at examining the increase of phenol adsorption breakthrough curves spreading caused by the chemical heterogeneity of granular activated carbon fixed beds. The local and the thermodynamic equilibrium assumption, as well as the nonlinear adsorption obeying to Langmuir isotherm, are considered. This study particularly tempts to link the reduced variance of phenol breakthrough curves to a measurable quantifying parameter of the chemical heterogeneity. The investigated artificial heterogeneous media are prepared by alternating layers of two types of granular activated carbon, active and non-active ones, that have similar physical properties. On the one hand, the chemical heterogeneity is quantified by the active layer relative thickness of the column length, l₁/L. On the other hand, it is quantified by the mean value of the probability distribution γ. The latter also represents the mean active grains mass ratio of the total medium mass, hence the medium mean capacity. The obtained results show an increase in the reduced variance and thus the effective global dispersion with the heterogeneity; the increase is as important as the medium capacity decreases. However, the dispersion increase achieves a limit value, even when the heterogeneity increases. The results are statistically modelled using a regression equation function of the capacity variation in terms of γ and the chemical heterogeneity in terms of l₁/L. The relationship combining the medium capacity and the chemical heterogeneity is obtained. The relationship implicitly takes into account the effect of the column length.     

Permeability dependence of streaming potential coefficient in porous media

In theory, the streaming potential coefficient depends not only on the zeta potential but also on the permeability of the rocks that partially determines the surface conductivity of the rocks. However, in practice, it is hard to show the permeability dependence of streaming potential coefficients because of the variation of zeta potential from sample to sample. To study permeability dependence of streaming potential, including the effects of the variation of the zeta potential and surface conductance due to the difference in mineral compositions between samples, we perform measurements on 12 consolidated samples, including natural and artificial samples saturated with 7 different NaCl solutions to determine the streaming potential coefficients. The results have shown that the streaming potential coefficients strongly depend on the permeability of the samples for low fluid conductivity. When the fluid conductivity is larger than than 0.50 S/m for the natural samples or 0.25 S/m for the artificial ceramic samples, the streaming potential coefficient is independent of permeability. This behavior is quantitatively explained by a theoretical model.