共查询到20条相似文献，搜索用时 31 毫秒
1.
This study formulates and analyzes continuous time random walk (CTRW) models in radial flow geometries for the quantification of non-local solute transport induced by heterogeneous flow distributions and by mobile–immobile mass transfer processes. To this end we derive a general CTRW framework in radial coordinates starting from the random walk equations for radial particle positions and times. The particle density, or solute concentration is governed by a non-local radial advection–dispersion equation (ADE). Unlike in CTRWs for uniform flow scenarios, particle transition times here depend on the radial particle position, which renders the CTRW non-stationary. As a consequence, the memory kernel characterizing the non-local ADE, is radially dependent. Based on this general formulation, we derive radial CTRW implementations that (i) emulate non-local radial transport due to heterogeneous advection, (ii) model multirate mass transfer (MRMT) between mobile and immobile continua, and (iii) quantify both heterogeneous advection in a mobile region and mass transfer between mobile and immobile regions. The expected solute breakthrough behavior is studied using numerical random walk particle tracking simulations. This behavior is analyzed by explicit analytical expressions for the asymptotic solute breakthrough curves. We observe clear power-law tails of the solute breakthrough for broad (power-law) distributions of particle transit times (heterogeneous advection) and particle trapping times (MRMT model). The combined model displays two distinct time regimes. An intermediate regime, in which the solute breakthrough is dominated by the particle transit times in the mobile zones, and a late time regime that is governed by the distribution of particle trapping times in immobile zones. These radial CTRW formulations allow for the identification of heterogeneous advection and mobile-immobile processes as drivers of anomalous transport, under conditions relevant for field tracer tests. 相似文献
2.
3.
We determine the relevance of Multi-Rate Mass Tansfer (MRMT) models (Haggerty and Gorelick, 1995) to general diffusive porosity structures. To this end, we introduce Structured INteracting Continua (SINC) models as the combination of a finite number of diffusion-dominated interconnected immobile zones exchanging with an advection-dominated mobile domain. It directly extends Multiple INteracting Continua framework (Pruess and Narasimhan, 1985) by introducing a structure in the immobile domain, coming for example from the dead-ends of fracture clusters or poorly-connected dissolution patterns. We demonstrate that, whatever their structure, SINC models can be made equivalent in terms of concentration in the mobile zone to a unique MRMT model. We develop effective shape-free numerical methods to identify its few dominant rates, that comply with any distribution of rates and porosities. We show that differences in terms of macrodispersion are not larger than 50% for approximate MRMT models with only one rate (double porosity models), and drop down to less than 0.1% for five rates MRMT models. Low-dimensional MRMT models accurately approach transport in structured diffusive porosities at intermediate and long times and only miss early responses. 相似文献
4.
The late tail of the breakthrough curve (BTC) of a conservative tracer in a regional-scale alluvial system is explored using Monte Carlo simulations. The ensemble numerical BTC, for an instantaneous point source injected into the mobile domain, has a heavy late tail transforming from power law to exponential due to a maximum thickness of clayey material. Haggerty et al.'s (2000) multiple-rate mass transfer (MRMT) method is used to predict the numerical late-time BTCs for solutes in the mobile phase. We use a simple analysis of the thicknesses of fine-grained units noted in boring logs to construct the memory function that describes the slow decline of concentrations at very late time. The good fit between the predictions and the numerical results indicates that the late-time BTC can be approximated by a summation of a small number of exponential functions, and its shape depends primarily on the thicknesses and the associated volume fractions of immobile water in "blocks" of fine-grained material. The prediction of the late-time BTC using the MRMT method relies on an estimate of the average advective residence time, t(ad). The predictions are not sensitive to estimation errors in t(ad), which can be approximated by L/v , where v is the arithmetic mean ground water velocity and L is the transport distance. This is the first example of deriving an analytical MRMT model from measured hydrofacies properties to predict the late-time BTC. The parsimonious model directly and quantitatively relates the observable subsurface heterogeneity to nonlocal transport parameters. 相似文献
5.
《Advances in water resources》2005,28(1):33-42
A numerical solution that is significantly more general than other semi-analytical solutions is presented for governing equations describing advective–dispersive transport with multirate mass transfer between mobile and immobile domains. The new solution approach is general in the sense that it does not impose any restrictive assumption on the spatial or temporal variability of advective and dispersive processes in the mobile domain. A single integro-differential equation (IDE) is developed for the concentration in the mobile domain by separating the concentration in the immobile domain from the set of two partial differential equations. The solution to the IDE requires the evaluation of a temporal integral of the concentration in the mobile domain, which is a function of the Laplace transform of the distribution of the mass transfer rate coefficient. The Laplace transform is not limited to flow fields with known constant velocities. The solutions for one- and two-dimensional examples obtained using the new approach agree with those obtained by existing semi-analytical and numerical approaches. 相似文献
6.
This paper deals with the lower order (first four) nonstationary statistical moments of the response of linear systems with random stiffness and random damping properties subject to random nonstationary excitation modeled as white noise multiplied by an envelope function. The method of analysis is based on a Markov approach using stochastic differential equations (SDE). The linear SDE with random coefficients subject to random excitation with deterministic initial conditions are transformed to an equivalent nonlinear SDE with deterministic coefficients and random initial conditions subject to random excitation. In this procedure, new SDE with random initial conditions, deterministic coefficients and zero forcing functions are introduced to represent the random variables. The joint statistical moments of the response are determined by considering an augmented dynamic system with state variables made up of the displacement and velocity vectors and the random variables of the structural system. The zero time-lag joint statistical moment equations for the augmented state vector are derived from the Itô differential formula. The statistical moment equations are ordinary nonlinear differential equations where hierarchy of moments appear. The hierarchy is closed by the cumulant neglect closure method applied at the fourth order statistical moment level. General formulation is given for multi-degree-of-freedom (MDOF) systems and the performance of the method in problems with nonstationary excitations and large variabilities is illustrated for a single-degree-of-freedom (SDOF) oscillator. 相似文献
7.
Diffusive mass exchange of non‐reactive substances in dual‐porosity porous systems – column experiments under saturated conditions 
下载免费PDF全文
下载免费PDF全文 Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
8.
The assumptions and analysis of Part 1 are extended to the three model soil-water profiles having different immobile-water contents up to 0.08 m3 m?3. Concentration profiles in the mobile zone obtained with a constant dispersion coefficient are compared with those calculated with the same dispersion coefficient and no immobile water. The average concentration profiles were calculated from the concentration of the solution in both the mobile and immobile phases. Increasing the immobile-water content decreased and shifted the average concentration profiles further into the soil. The values of the dispersion coefficient calculated with the immobile-water fraction ignored, were found to be water content dependent and many times greater than that of the input dispersion coefficient. 相似文献
9.
In modeling solute transport with mobile‐immobile mass transfer (MIMT), it is common to use an advection‐dispersion equation (ADE) with a retardation factor, or retarded ADE. This is commonly referred to as making the local equilibrium assumption (LEA). Assuming local equilibrium, Eulerian textbook treatments derive the retarded ADE, ostensibly exactly. However, other authors have presented rigorous mathematical derivations of the dispersive effect of MIMT, applicable even in the case of arbitrarily fast mass transfer. We resolve the apparent contradiction between these seemingly exact derivations by adopting a Lagrangian point of view. We show that local equilibrium constrains the expected time immobile, whereas the retarded ADE actually embeds a stronger, nonphysical, constraint: that all particles spend the same amount of every time increment immobile. Eulerian derivations of the retarded ADE thus silently commit the gambler's fallacy, leading them to ignore dispersion due to mass transfer that is correctly modeled by other approaches. We then present a particle tracking simulation illustrating how poor an approximation the retarded ADE may be, even when mobile and immobile plumes are continually near local equilibrium. We note that classic “LEA” (actually, retarded ADE validity) criteria test for insignificance of MIMT‐driven dispersion relative to hydrodynamic dispersion, rather than for local equilibrium. 相似文献
10.
A Practical Modeling Framework for Non‐Fickian Transport and Multi‐Species Sequential First‐Order Reaction 下载免费PDF全文
下载免费PDF全文 Daniel K. Burnell Jie Xu Scott K. Hansen Lawrence S. Sims Charles R. Faust 《Ground water》2018,56(4):524-540
Many studies indicate that small‐scale heterogeneity and/or mobile–immobile mass exchange produce transient non‐Fickian plume behavior that is not well captured by the use of the standard, deterministic advection‐dispersion equation (ADE). An extended ADE modeling framework is presented here that is based on continuous time random walk theory. It can be used to characterize non‐Fickian transport coupled with simultaneous sequential first‐order reactions (e.g., biodegradation or radioactive decay) for multiple degrading contaminants such as chlorinated solvents, royal demolition explosive, pesticides, and radionuclides. To demonstrate this modeling framework, new transient analytical solutions are derived and are inverted in Laplace space. Closed‐form, steady‐state, multi‐species analytical solutions are also derived for non‐Fickian transport in highly heterogeneous aquifers with linear sorption–desorption and matrix diffusion for use in spreadsheets. The solutions are general enough to allow different degradation rates for the mobile and immobile zones. The transient solutions for multi‐species transport are applied to examine the effects of source remediation on the natural attenuation of downgradient plumes of both parent and degradation products in highly heterogeneous aquifers. Results for representative settings show that the use of the standard, deterministic ADE can over‐estimate cleanup rates and under‐predict the cleanup timeframe in comparison to the extended ADE analytical model. The modeling framework and calculations introduced here are also applied for a 30 year groundwater cleanup program at a site in Palm Bay, Florida. The simulated plume concentrations using the extended ADE exhibited agreement with observed long concentration tails of trichloroethene, cis 1,2 DCE, and VC that remained above cleanup goals. 相似文献
11.
Many problems in hydraulics and hydrology are described by linear, time dependent partial differential equations, linearity being, of course, an assumption based on necessity.Solutions to such equations have been obtained in the past based purely on deterministic consideration. The derivation of such a solution requires that the initial conditions, the boundary conditions, and the parameters contained within the equations be stipulated in exact terms. It is obvious that the solution so derived is a function of these specified, values.There are at least four ways in which randomness enters the problem. i) the random initial value problem; ii) the random boundary value problem; iii) the random forcing problem when the non-homogeneous part becomes random and iv) the random parameter problem.Such randomness is inherent in the environment surrounding the system, the environment being endowed with a large number of degrees of freedom.This paper considers the problem of groundwater flow in a phreatic aquifer fed by rainfall. The goveming equations are linear second order partial differential equations. Explicit form solutions to this randomly forced equation have been derived in well defined regular boundaries. The paper also provides a derivation of low order moment equations. It contains a discussion on the parameter estimation problem for stochastic partial differential equations. 相似文献
12.
The paper presents Direct Numerical Simulations of an open channel flow laden with spherical particles at a bulk Reynolds number of 2941. The transport of thousands of mobile particles is simulated propagating over a rough bed which consists of immobile particles of the same size in hexagonal ordering. An Immersed Boundary Method is used for the numerical representation of the particles. With 22 points per diameter even the viscous scales of the flow are resolved at this Reynolds number. The reference run contains just as many fixed as mobile particles with a relative density slightly above the nominal threshold of incipient motion. Further runs were conducted with decreased mass loading and decreased Shields number together with a simulation containing only immobile particles. The variation of the parameters defining the mobile sediment yields a strong modification of particle–fluid as well as particle–particle interactions yielding different structures in space and time. This is assessed by means of appropriate statistical quantities addressing the continuous and the disperse phase. The results are in qualitative agreement with experimental observations at higher Reynolds number. 相似文献
13.
We consider a Bayesian model for inversion of observed amplitude variation with offset data into lithology/fluid classes, and study in particular how the choice of prior distribution for the lithology/fluid classes influences the inversion results. Two distinct prior distributions are considered, a simple manually specified Markov random field prior with a first-order neighbourhood and a Markov mesh model with a much larger neighbourhood estimated from a training image. They are chosen to model both horizontal connectivity and vertical thickness distribution of the lithology/fluid classes, and are compared on an offshore clastic oil reservoir in the North Sea. We combine both priors with the same linearized Gaussian likelihood function based on a convolved linearized Zoeppritz relation and estimate properties of the resulting two posterior distributions by simulating from these distributions with the Metropolis–Hastings algorithm. The influence of the prior on the marginal posterior probabilities for the lithology/fluid classes is clearly observable, but modest. The importance of the prior on the connectivity properties in the posterior realizations, however, is much stronger. The larger neighbourhood of the Markov mesh prior enables it to identify and model connectivity and curvature much better than what can be done by the first-order neighbourhood Markov random field prior. As a result, we conclude that the posterior realizations based on the Markov mesh prior appear with much higher lateral connectivity, which is geologically plausible. 相似文献
14.
Flow of nonvolatile nonaqueous phase liquid (NAPL) and aqueous phases that account for mobile, entrapped, and residual NAPL in variably saturated water-wet porous media is modeled and compared against results from detailed laboratory experiments. Residual saturation formation in the vadose zone is a process that is often ignored in multifluid flow simulators, which might cause an overestimation of the volume of NAPL that reaches the ground water. Mobile NAPL is defined as being continuous in the pore space and flows under a pressure gradient or gravitational body force. Entrapped NAPL is defined as being occluded by the aqueous phase, occurring as immobile ganglia surrounded by aqueous phase in the pore space and formed when NAPL is replaced by the aqueous phase. Residual NAPL is defined as immobile, nonwater entrapped NAPL that does not drain from the pore spaces and is conceptualized as being either continuous or discontinuous. Free NAPL comprises mobile and residual NAPL. The numerical model is formulated on mass conservation equations for oil and water, transported via NAPL and aqueous phases through variably saturated porous media. To account for phase transitions, a primary variable switching scheme is implemented for the oil-mass conservation equation over three phase conditions: (1) aqueous or aqueous-gas with dissolved oil, (2) aqueous or aqueous-gas with entrapped NAPL, and (3) aqueous or aqueous gas with free NAPL. Two laboratory-scale column experiments are modeled to verify the numerical model. Comparisons between the numerical simulations and experiments demonstrate the necessity to include the residual NAPL formation process in multifluid flow simulators. 相似文献
15.
The U.S. Department of Energy is currently studying Yucca Mountain, Nevada, as a potential site for a geological high-level waste repository. In the current conceptual models of radionuclide transport at Yucca Mountain, part of the transport path to pumping locations would be through an alluvial aquifer. Interactions with minerals in the alluvium are expected to retard the downstream migration of radionuclides, thereby delaying arrival times and reducing ground water concentrations. We evaluate the effectiveness of the alluvial aquifer as a transport barrier using the stochastic Lagrangian framework. A transport model is developed to account for physical and chemical heterogeneities and rate-limited mass transfer between mobile and immobile zones. The latter process is caused by small-scale heterogeneity and is thought to control the macroscopic-scale retardation in some field experiments. A geostatistical model for the spatially varying sorption parameters is developed from a site-specific database created from hydrochemical measurements and a calibrated modeling approach (Turner and Pabalan 1999). Transport of neptunium is considered as an example. The results are sensitive to the rate of transfer between mobile and immobile zones, and to spatial variability in the hydraulic conductivity. Chemical heterogeneity has only a small effect, as does correlation between hydraulic conductivity and the neptunium distribution coefficient. These results illustrate how general sensitivities can be explored with modest effort within the Lagrangian framework. Such studies complement and guide the application of more detailed numerical simulations. 相似文献
16.
研究了二维(2-D)介质粗糙面下方三维(3-D)金属目标的复合电磁散射问题.将表面积分方程(PMCHW)方程应用到介质粗糙面表面,电场积分方程(EFIE)应用于金属目标表面.基于矩量法,使用三角分域基函数(RWG)和伽略金法将表面积分方程离散为矩阵方程,并采用稳定的双共轭梯度迭代(BICGSTAB)算法对矩阵方程进行求解.针对矩量法(MOM)的高存储量和迭代过程中存在的矩阵向量积耗时的瓶颈,采用基于秩的多层矩阵分解法(MLUV),对矩阵元素进行压缩存储,以节省对计算机内存的需求,并加速迭代过程中的矩阵向量积运算.计算了高斯粗糙面下方球体的双站雷达散射截面积(RCS),并与最陡下降快速多级子算法(SDFMM)结果比较以验证该数值方法的正确性.最后分析了不同粗糙度、目标尺寸和目标位置对双站RCS的影响. 相似文献
17.
This study addresses estimation of net irrigation requirement over a growing season under climate uncertainty. An ecohydrological model, building upon the stochastic differential equation of soil moisture dynamics, is employed as a basis to derive new analytical expressions for estimating seasonal net irrigation requirement probabilistically. Two distinct irrigation technologies are considered. For micro irrigation technology, probability density function of seasonal net irrigation depth (SNID) is derived assessing transient behavior of a stochastic process which is time integral of dichotomous Markov process. Probability mass function of SNID which is a discrete random variable for traditional irrigation technology is also presented using a marked renewal process with quasi-exponentially-distributed time intervals. Comparing the results obtained from the presented models with those resulted from a Monte Carlo approach verified the significance of the probabilistic expressions derived and assumptions made. 相似文献
18.
The application of a modified version of dynamic TOPMODEL for two subcatchments at Plynlimon, Wales is described. Conservative chemical mixing within mobile and immobile stores has been added to the hydrological model in an attempt to simulate observed stream chloride concentrations. The model was not fully able to simulate the observed behaviour, in particular the short‐ to medium‐term dynamics. One of the primary problems highlighted by the study was the representation of dry deposition and cloud‐droplet‐deposited chloride, which formed a significant part of the long‐term chloride mass budget. Equifinality of parameter sets inhibited the ability to determine the effective catchment mixing volumes and coefficients or the most likely partition between occult mass inputs and chloride mass inputs determined by catchment immobile‐store antecedent conditions. Some success was achieved, in as much as some aspects of the dynamic behaviour of the signal were satisfactorily simulated, although spectral analysis showed that the model could not fully reproduce the 1/f power spectra of observed stream chloride concentrations with its implications of a wide distribution of residence times for water in the catchment. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
19.
Memory of a stochastic process implies its predictability, understood as a possibility to gain information on the future above
the random guess level. Here we search for memory in the mining-induced seismic process (MIS), that is, a process induced
or triggered by mining operations. Long memory is investigated by means of the Hurst rescaled range analysis, and the autocorrelation
function estimate is used to test for short memory. Both methods are complemented with result uncertainty analyses based on
different resampling techniques. The analyzed data comprise event series from Rudna copper mine in Poland. The studies show
that the interevent time and interevent distance processes have both long and short memory. MIS occurrences and locations
are internally interrelated. Internal relations among the sizes of MIS events are apparently weaker than those of other two
studied parameterizations and are limited to long term interactions. 相似文献
20.
Changes in the water table level result in variable water saturation and variable hydrological fluxes at the interface between the unsaturated and saturated zone. This may influence the transport and fate of contaminants in the subsurface. The objective of this study was to examine the impact of a decreasing and an increasing water table on solute transport. We conducted tracer experiments at downward flow conditions in laboratory columns filled with two different uniform porous media under static and transient flow conditions either increasing or decreasing the water table. Tracer breakthrough curves were simulated using a mobile–immobile transport model. The resulting transport parameters were compared to identify dominant transport processes. Changes in the water table level affected dispersivities and mobile water fractions depending on the direction of water table movement and the grain size of the porous media. In fine glass beads, the water flow velocity was similar to the decline rate of the water table, and the mobile water fraction was decreased compared with steady‐state saturated conditions. However, immobile water was negligible. In coarse glass beads, water flow was faster because of fingered flow in the unsaturated part, and the mobile water fraction was smaller than in the fine material. Here, a rising water table led to an even smaller mobile water fraction and increased solute spreading because of diffusive interaction with immobile water. We conclude that changes of the water table need to be considered to correctly simulate transport in the subsurface at the transition of the unsaturated–saturated zone. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献