共查询到20条相似文献，搜索用时 203 毫秒
1.
The fractional Brownian motion (fBm) and fractional Lévy motion (fLm) can easily describe the geometry and the statistical structure of hydraulic conductivity (K) for real-world. However, the fBm and fLm models have not been systematically evaluated when building the K field for a low-permeability site. In this study, both the fBm and fLm are used to simulate the low-K field at NingCheGu (NCG), Tianjin, China. Groundwater flow and solute transport are then computed using MODFLOW and MT3DMS, respectively, and the influence of the fBm/fLm models for K on groundwater flow and solute transport is discussed. Results show that the fLm fits better the statistics of the low-K medium than fBm, and the random logarithmic K (LnK) field generated by fLm is more stable because the resultant LnK field captures more of the measured properties at the field site than that generated by fBm. In contrast, the LnK generated by fBm is more likely to form both high-K channels and low-K barriers. The fBm therefore predicts more extreme behaviours in flow and transport, including the preferential flow, low-concentration blocks and solute retention. The overall groundwater renewal period and solute travel time for the fLm simulation are slightly shorter than those for fBm. The impacts of the fLm and fBm models on the statistics of the resultant LnK fields and the dynamics of groundwater flow and solute transport revealed by this study shed light on the selection and evaluation of the fractional probability distribution models in capturing the K fields for low-K media. 相似文献
2.
Permissibility of fractal exponents and models of band-limited two-point functions for fGn and fBm random fields 总被引：1，自引：1，他引：0
The fractional Gaussian noise (fGn) and fractional Brownian motion (fBm) random field models have many applications in the environmental sciences. An issue of practical interest is the permissible range and the relations between different fractal exponents used to characterize these processes. Here we derive the bounds of the covariance exponent for fGn and the Hurst exponent for fBm based on the permissibility theorem by Bochner. We exploit the theoretical constraints on the spectral density to construct explicit two-point (covariance and structure) functions that are band-limited fractals with smooth cutoffs. Such functions are useful for modeling a gradual cutoff of power-law correlations. We also point out certain peculiarities of the relations between fractal exponents imposed by the mathematical bounds. Reliable estimation of the correlation and Hurst exponents typically requires measurements over a large range of scales (more than 3 orders of magnitude). For isotropic fractals and partially isotropic self-affine processes the dimensionality curse is partially lifted by estimating the exponent from measurements along fixed directions. We derive relations between the fractal exponents and the one-dimensional spectral density exponents, and we illustrate the relations using measurements of paper roughness.The author would like to acknowledge helpful comments from two anonymous referees. 相似文献
3.
Apparent multifractality and scale‐dependent distribution of data sampled from self‐affine processes
Shlomo P. Neuman 《水文研究》2011,25(11):1837-1840
It has been previously demonstrated theoretically and numerically by the author that square or absolute increments of data sampled from fractional Brownian/Lévy motion (fBm/fLm), or of incremental data sampled from fractional Gaussian/Lévy noise (fGn/fLn), exhibit apparent/spurious multifractality. Here, we generalize these previous development in a way that (a) rigorously subordinates (truncated) fLn to fGn or, in a statistically equivalent manner, (truncated) fLm to fBm; (b) extends the analysis to a wider class of subordinated self‐affine processes; (c) provides a simple way to generate such processes and (d) explains why the distribution of corresponding increments tends to evolve from heavy tailed at small lags (separation distances or scales) to Gaussian at larger lags. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
4.
《Advances in water resources》1999,23(1):49-57
Travel-time statistics for non-reacting tracers in fractal and multifractal media are addressed through numerical simulations. The logarithm of hydraulic conductivity is modeled using fractional Brownian motion (fBm) and more recently developed multifractal model based on bounded fractional Levy motion (bfLm). These models have been shown previously to accurately reproduce statistical properties of large conductivity datasets. The ensemble-mean travel time increases nearly linearly with travel distance and the variance in the travel time increases nearly parabolically with travel distance. This is consistent with near-field analytical approximations developed for non-fractal media and suggests that these analytical results may have some degree of robustness to non-ideal features in the random-field models. The magnitudes of the travel-time moments are dependent on the system size. For fBm media, this size dependence can be explained using an effective variance that increases with increasing size of the flow system. However, the magnitudes of the travel-time moments are also sensitive to other non-ideal effects such as deviations from Gaussian behavior. This sensitivity illustrates the need for careful aquifer characterization and conditional numerical simulation in practical situations requiring accurate estimates of uncertainty in the plume position. 相似文献
5.
Kaili Wang Guanhua Huang 《Advances in water resources》2011,34(6):671-683
The effect of aquifer heterogeneity on flow and solute transport in two-dimensional isotropic porous media was analyzed using the Monte Carlo method. The two-dimensional logarithmic permeability (ln K) was assumed to be a non-stationary random field with its increments being a truncated fractional Lévy motion (fLm). The permeability fields were generated using the modified successive random additions (SRA) algorithm code SRA3DC [http://www.iamg.org/CGEditor/index.htm]. The velocity and concentration fields were computed respectively for two-dimensional flow and transport with a pulse input using the finite difference codes of MODFLOW 2000 and MT3DMS. Two fLm control parameters, namely the width parameter (C) and the Lévy index (α), were varied systematically to examine their effect on the resulting permeability, flow velocity and concentration fields. We also computed the first- and second-spatial moments, the dilution index, as well as the breakthrough curves at different control planes with the corresponding concentration fields. In addition, the derived breakthrough curves were fitted using the continuous time random walk (CTRW) and the traditional advection-dispersion equation (ADE). Results indicated that larger C and smaller α both led to more heterogeneous permeability and velocity fields. The Lévy-stable distribution of increments in ln K resulted in a Lévy-stable distribution of increments in logarithm of the velocity (ln v). Both larger C and smaller α created sharper leading edges and wider tailing edges of solute plumes. Furthermore, a relatively larger amount of solute still remained in the domain after a relatively longer time transport for smaller α values. The dilution indices were smaller than unity and increased as C increased and α decreased. The solute plume and its second-spatial moments increased as C increased and α decreased, while the first-spatial moments of the solute plume were independent of C and α values. The longitudinal macrodispersivity was scale-dependent and increased as a power law function of time. Increasing C and decreasing α both resulted in an increase in longitudinal macrodispersivity. The transport in such highly heterogeneous media was slightly non-Gaussian with its derived breakthrough curves being slightly better fitted by the CTRW than the ADE, especially in the early arrivals and late-time tails. 相似文献
6.
Shlomo P. Neuman 《水文研究》2010,24(15):2056-2067
Many earth and environmental variables appear to be self‐affine (monofractal) or multifractal with spatial (or temporal) increments having exceedance probability tails that decay as powers of − α where 1 < α ≤ 2. The literature considers self‐affine and multifractal modes of scaling to be fundamentally different, the first arising from additive and the second from multiplicative random fields or processes. We demonstrate theoretically that data having finite support, sampled across a finite domain from one or several realizations of an additive Gaussian field constituting fractional Brownian motion (fBm) characterized by α = 2, give rise to positive square (or absolute) increments which behave as if the field was multifractal when in fact it is monofractal. Sampling such data from additive fractional Lévy motions (fLm) with 1 < α < 2 causes them to exhibit spurious multifractality. Deviations from apparent multifractal behaviour at small and large lags are due to nonzero data support and finite domain size, unrelated to noise or undersampling (the causes cited for such deviations in the literature). Our analysis is based on a formal decomposition of anisotropic fLm (fBm when α = 2) into a continuous hierarchy of statistically independent and homogeneous random fields, or modes, which captures the above behaviour in terms of only E + 3 parameters where E is Euclidean dimension. Although the decomposition is consistent with a hydrologic rationale proposed by Neuman (2003), its mathematical validity is independent of such a rationale. Our results suggest that it may be worth checking how closely would variables considered in the literature to be multifractal (e.g. experimental and simulated turbulent velocities, some simulated porous flow velocities, landscape elevations, rain intensities, river network area and width functions, river flow series, soil water storage and physical properties) fit the simpler monofractal model considered in this paper (such an effort would require paying close attention to the support and sampling window scales of the data). Parsimony would suggest associating variables found to fit both models equally well with the latter. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
7.
This paper proposes a numerical method to simulate oil spill trajectories, which are affected by the combination of advection, turbulent diffusion and mechanical spreading process, based on a particle tracking algorithm. Recent studies have shown that the trajectories of drifters on the ocean surface have a fractal structure that is far from being described using ordinary Brownian motion. Thus, in modeling the diffusion process, a discrete method has been employed for the generation of fractional Brownian motion (fBm) to illustrate superdiffusive transport. The algorithm is implemented to predict oil slick trajectories following the “Arteaga” oil spill accident that occurred near the Dalian coastal region in 2005. When compared with the observed data and the results of traditional diffusion modeling, the numerical results based on the fBm model are encouraging. 相似文献
8.
Intrinsic random fields of order k, defined as random fields whose high-order increments (generalized increments of order k) are second-order stationary, are used in spatial statistics to model regionalized variables exhibiting spatial trends, a feature that is common in earth and environmental sciences applications. A continuous spectral algorithm is proposed to simulate such random fields in a d-dimensional Euclidean space, with given generalized covariance structure and with Gaussian generalized increments of order k. The only condition needed to run the algorithm is to know the spectral measure associated with the generalized covariance function (case of a scalar random field) or with the matrix of generalized direct and cross-covariances (case of a vector random field). The algorithm is applied to synthetic examples to simulate intrinsic random fields with power generalized direct and cross-covariances, as well as an intrinsic random field with power and spline generalized direct covariances and Matérn generalized cross-covariance. 相似文献
9.
On the reconstruction of structural and functional properties in random heterogeneous media 总被引：1，自引：0，他引：1
Flow and transport in porous media is determined by its structure. Beside spatial correlation, especially the connectivity of heterogeneous conductivities is acknowledged to be a key factor. This has been demonstrated for well defined random fields having different topological properties. Yet, it remains an open question which morphological measures carry sufficient information to actually predict flow and transport in porous media. We analyze flow and transport in classical, two-dimensional random fields showing different topology and we determine a selection of structural characteristics including classical two-point statistics, chord-length distribution and Minkowski functions (four-point statistics) including the Euler number as a topological measure. Using the approach of simulated annealing for global optimization we generate analog random fields that are forced to reproduce one or several of theses structural characteristics. Finally we evaluate in how far the generated analogons reproduce the original flow and transport behavior as well as some more elaborate structural characteristics including percolation probabilities and the pair connectivity function. The results confirm that two-point statistics is insufficient to capture functional properties since it is not sensitive to connectivity. In contrast, the combination of Minkowski functions and chord length distributions carries sufficient information to reproduce the breakthrough curve of a conservative solute. Hence, global topology provided by the Euler number together with local clustering provided by the chord length distribution seems to be a powerful condensation of structural complexity with respect to functional properties. 相似文献
10.
11.
Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
12.
Hakan Sirin 《Stochastic Environmental Research and Risk Assessment (SERRA)》2006,20(6):381-390
Solute plume subjected to field scale hydraulic conductivity heterogeneity shows a large dispersion/macrodispersion, which is the manifestation of existing fields scale heterogeneity on the solute plume. On the other hand, due to the scarcity of hydraulic conductivity measurements at field scale, hydraulic conductivity heterogeneity can only be defined statistically, which makes the hydraulic conductivity a random variable/function. Random hydraulic conductivity as a parameter in flow equation makes the pore flow velocity also random and the ground water solute transport equation is a stochastic differential equation now. In this study, the ensemble average of stochastic ground water solute transport equation is taken by the cumulant expansion method in order to upscale the laboratory scale transport equation to field scale by assuming pore flow velocity is a non stationary, non divergence-free and unsteady random function of space and time. Besides the stochastic explanation of macrodispersion and the velocity correction term obtained by Kavvas and Karakas (J Hydrol 179:321–351, 1996) before a new velocity correction term, which is a function of mean pore flow velocity divergence, is obtained in this study due to strict second order cumulant expansion (without omitting any term after the expansion) performed. The significance of the new velocity correction term is investigated on a one dimensional transport problem driven by a density dependent flow field. 相似文献
13.
The results of a field study of the small-scale spatial structure of the depth of shallow seasonal snowcovers in prairie and arctic environments are presented. It is shown that the spatial distribution of snow depth is fractal at small scales, becoming random at scales beyond some limiting length. This is due to the autocorrelation of depth at small sampling distances. The transition of fractal to random behaviour is indexed by a ‘cutoff length’, which is defined by the intersection of the ‘fractal’ slope and horizontal tangent of a logarithmic plot of the standard deviation of depth versus sampling distance. The magnitude of the cutoff length is related to the degree of macroscopic variability of the underlying topography. An increase in length due to the effects of macroscopic topographic variability on snowcover accumulation is confirmed by de-trending field measurements. The de-trended data shown a cutoff length for wheat stubble and fallow surfaces of approximately 30 m, which is consistent with the distance determined from measurements on ‘flat’ fields. The implications of the transition of snow depth from fractal to random structure on the scales of snow sampling and modelling are presented. The cutoff length may provide a statistic for stratifying shallow snowcovers, by linking snowcover properties to the underlying topography. 相似文献
14.
Alberto Bellin Yoram Rubin 《Stochastic Environmental Research and Risk Assessment (SERRA)》1996,10(4):253-278
This paper describes a new method for generating spatially-correlated random fields. Such fields are often encountered in hydrology and hydrogeology and in the earth sciences. The method is based on two observations: (i) spatially distributed attributes usually display a stationary correlation structure, and (ii) the screening effect of measurements leads to the sufficiency of a small search neighborhood when it comes to projecting measurements and data in space. The algorithm which was developed based on these principles is called HYDRO_GEN, and its features and properties are discussed in depth. HYDRO_GEN is found to be accurate and extremely fast. It is also versatile: it can simulate fields of different nature, starting from weakly stationary fields with a prescribed covariance and ending with fractal fields. The simulated fields can display statistical isotropy or anisotropy. 相似文献
15.
Paul S ADDISION 《国际泥沙研究》2009,24(4):439-454
An accurate particle tracking method using FBMINC （new fractional Brownian motion） is outlined. It generates non-Fickian diffusion rather than Fickian diffusion as traditional particle tracking model does. The FBMINC model is based on fractional Brownian motion （fl3m） which is generalization of regular Brownian motion. The two models of fBms （FBM model and FBMINC model） were explored and the differences of the two models are compared from the three aspects： the standard deviation of each step, the small memory and the effect of the number of particles in the cloud. The results show the FBMINC model is a better model as it produces more accurate statistics. The effect of simple shear dispersion for both Brownian and fBm was investigated. The power law scaling of fBm shear dispersion was correctly identified. In addition, a scaling coettieient was found numerically. The FBMINC model is then used for producing both superdiffusive and subdiffusive particle paths, therefore, the non-Fickian diffusion of soil particle clouds can be modelled. The particle clouds represent soil contaminant are released in an idealised coastal bay and the fBm particle tracking method is used for simulation the particle cloud spreading in the bay. There is a noticeable increase in the spreading rate of the cloud. In addition, owning to the spreading rate of the cloud, a noticeable part of it has escaped the bay area and transported downstream. The variation of the Hurst exponent can lead to an area of the flow being affected by a contaminant cloud which is not picked up by the regular Brownian motion models. The purpose of this paper is to bring soil transport engineers a new angle on soil particle transport research in fluids. Using FBM1NC particle tracking model allows more flexibility in simulation of diffusion in soil contaminant spread in coastal bay or ocean surface. 相似文献
16.
《Advances in water resources》2004,27(3):207-222
An Eulerian perturbation approach was applied to develop a method of moment for solute transport in a nonstationary, fractured medium. The conceptualized fractured medium is described through a dual-porosity model. Stochastic governing equations for mean concentration and concentration covariance were analytically derived to the first-order accuracy of log-conductivity variance and solved with a numerical method––a finite difference method. The developed method is called a numerical Eulerian method of moment (NEMM). This method was compared with the stationary transport theory [Water Resour. Res. 36(7) (2000) 1665] for predicting mean concentration and its spatial moments. The comparison indicated that the two methods matched very well in predicting first and second spatial moments. NEMM solutions were also compared with Monte Carlo simulations for solute transport in stationary fractured media. The results of the two methods were consistent for calculating small log conductivity variance. The theory was then used to study effects of various parameters and nonstationarity of the medium on flow and transport processes. Results indicated that medium nonstationarity would significantly influence the solute transport process. The nonstationary transport theory relaxes many assumptions adopted in stationary theories and paves the way for applying the NEMM to many environmental projects, especially in analyzing uncertainty of solute transport. 相似文献
17.
Subsurface fluid flow and solute transport take place in a multiscale heterogeneous environment. Neither these phenomena nor their host environment can be observed or described with certainty at all scales and locations of relevance. The resulting ambiguity has led to alternative conceptualizations of flow and transport and multiple ways of addressing their scale and space–time dependencies. We focus our attention on four approaches that give rise to nonlocal representations of advective and dispersive transport of nonreactive tracers in randomly heterogeneous porous or fractured continua. We compare these approaches theoretically on the basis of their underlying premises and the mathematical forms of the corresponding nonlocal advective–dispersive terms. One of the four approaches describes transport at some reference support scale by a classical (Fickian) advection–dispersion equation (ADE) in which velocity is a spatially (and possibly temporally) correlated random field. The randomness of the velocity, which is given by Darcy’s law, stems from random fluctuations in hydraulic conductivity (and advective porosity though this is often disregarded). Averaging the stochastic ADE over an ensemble of velocity fields results in a space–time-nonlocal representation of mean advective–dispersive flux, an approach we designate as stnADE. A closely related space–time-nonlocal representation of ensemble mean transport is obtained upon averaging the motion of solute particles through a random velocity field within a Lagrangian framework, an approach we designate stnL. The concept of continuous time random walk (CTRW) yields a representation of advective–dispersive flux that is nonlocal in time but local in space. Closely related to the latter are forms of ADE entailing fractional derivatives (fADE) which leads to representations of advective–dispersive flux that are nonlocal in space but local in time; nonlocality in time arises in the context of multirate mass transfer models, which we exclude from consideration in this paper. We describe briefly each of these four nonlocal approaches and offer a perspective on their differences, commonalities, and relative merits as analytical and predictive tools. 相似文献
18.
Yong Zhang Dongbao Zhou Maosheng Yin HongGuang Sun Wei Wei Shiyin Li Chunmiao Zheng 《水文研究》2020,34(25):5104-5122
Modelling pollutant transport in water is one of the core tasks of computational hydrology, and various physical models including especially the widely used nonlocal transport models have been developed and applied in the last three decades. No studies, however, have been conducted to systematically assess the applicability, limitations and improvement of these nonlocal transport models. To fill this knowledge gap, this study reviewed, tested and improved the state-of-the-art nonlocal transport models, including their physical background, mathematical formula and especially the capability to quantify conservative tracers moving in one-dimensional sand columns, which represents perhaps the simplest real-world application. Applications showed that, surprisingly, neither the popular time-nonlocal transport models (including the multi-rate mass transfer model, the continuous time random walk framework and the time fractional advection-dispersion equation), nor the spatiotemporally nonlocal transport model (ST-fADE) can accurately fit passive tracers moving through a 15-m-long heterogeneous sand column documented in literature, if a constant dispersion coefficient or dispersivity is used. This is because pollutant transport in heterogeneous media can be scale-dependent (represented by a dispersion coefficient or dispersivity increasing with spatiotemporal scales), non-Fickian (where plume variance increases nonlinearly in time) and/or pre-asymptotic (with transition between non-Fickian and Fickian transport). These different properties cannot be simultaneously and accurately modelled by any of the transport models reviewed by this study. To bypass this limitation, five possible corrections were proposed, and two of them were tested successfully, including a time fractional and space Hausdorff fractal model which minimizes the scale-dependency of the dispersion coefficient in the non-Euclidean space, and a two-region time fractional advection-dispersion equation which accounts for the spatial mixing of solute particles from different mobile domains. Therefore, more efforts are still needed to accurately model transport in non-ideal porous media, and the five model corrections proposed by this study may shed light on these indispensable modelling efforts. 相似文献
19.
We investigate effective solute transport in a chemically heterogeneous medium subject to temporal fluctuations of the flow conditions. Focusing on spatial variations in the equilibrium adsorption properties, the corresponding fluctuating retardation factor is modeled as a stationary random space function. The temporal variability of the flow is represented by a stationary temporal random process. Solute spreading is quantified by effective dispersion coefficients, which are derived from the ensemble average of the second centered moments of the normalized solute distribution in a single disorder realization. Using first-order expansions in the variances of the respective random fields, we derive explicit compact expressions for the time behavior of the disorder induced contributions to the effective dispersion coefficients. Focusing on the contributions due to chemical heterogeneity and temporal fluctuations, we find enhanced transverse spreading characterized by a transverse effective dispersion coefficient that, in contrast to transport in steady flow fields, evolves to a disorder-induced macroscopic value (i.e., independent of local dispersion). At the same time, the asymptotic longitudinal dispersion coefficient can decrease. Under certain conditions the contribution to the longitudinal effective dispersion coefficient shows superdiffusive behavior, similar to that observed for transport in s stratified porous medium, before it decreases to its asymptotic value. The presented compact and easy to use expressions for the longitudinal and transverse effective dispersion coefficients can be used for the quantification of effective spreading and mixing in the context of the groundwater remediation based on hydraulic manipulation and for the effective modeling of reactive transport in heterogeneous media in general. 相似文献
20.
Sergei A. Fomin Vladimir A. ChugunovToshiyuki Hashida 《Advances in water resources》2011,34(2):205-214
The paper provides an introduction to fundamental concepts of mathematical modeling of mass transport in fractured porous heterogeneous rocks. Keeping aside many important factors that can affect mass transport in subsurface, our main concern is the multi-scale character of the rock formation, which is constituted by porous domains dissected by the network of fractures. Taking into account the well-documented fact that porous rocks can be considered as a fractal medium and assuming that sizes of pores vary significantly (i.e. have different characteristic scales), the fractional-order differential equations that model the anomalous diffusive mass transport in such type of domains are derived and justified analytically. Analytical solutions of some particular problems of anomalous diffusion in the fractal media of various geometries are obtained. Extending this approach to more complex situation when diffusion is accompanied by advection, solute transport in a fractured porous medium is modeled by the advection-dispersion equation with fractional time derivative. In the case of confined fractured porous aquifer, accounting for anomalous non-Fickian diffusion in the surrounding rock mass, the adopted approach leads to introduction of an additional fractional time derivative in the equation for solute transport. The closed-form solutions for concentrations in the aquifer and surrounding rocks are obtained for the arbitrary time-dependent source of contamination located in the inlet of the aquifer. Based on these solutions, different regimes of contamination of the aquifers with different physical properties can be readily modeled and analyzed. 相似文献