首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 406 毫秒
1.
Hydraulic conductivity distribution and plume initial source condition are two important factors affecting solute transport in heterogeneous media. Since hydraulic conductivity can only be measured at limited locations in a field, its spatial distribution in a complex heterogeneous medium is generally uncertain. In many groundwater contamination sites, transport initial conditions are generally unknown, as plume distributions are available only after the contaminations occurred. In this study, a data assimilation method is developed for calibrating a hydraulic conductivity field and improving solute transport prediction with unknown initial solute source condition. Ensemble Kalman filter (EnKF) is used to update the model parameter (i.e., hydraulic conductivity) and state variables (hydraulic head and solute concentration), when data are available. Two-dimensional numerical experiments are designed to assess the performance of the EnKF method on data assimilation for solute transport prediction. The study results indicate that the EnKF method can significantly improve the estimation of the hydraulic conductivity distribution and solute transport prediction by assimilating hydraulic head measurements with a known solute initial condition. When solute source is unknown, solute prediction by assimilating continuous measurements of solute concentration at a few points in the plume well captures the plume evolution downstream of the measurement points.  相似文献   

2.
Large peatland complexes dominate the landscape of the James Bay Lowland in subarctic Ontario, Canada. However, there is not a thorough understanding of the hydrological processes occurring in these important systems, particularly how ladder fens connect large domed bogs to the aquatic ecosystems that drain the peatland complex. Ladder fens consist of a pool‐rib topography where flow downgradient is controlled by the peat ribs. Within the ribs, low‐lying preferential flow paths typically enhance the transmission of water, whereas the elevated ridge microforms impede water flow to downgradient aquatic ecosystems. To assess the hydrological connectivity, we study the role of the water table, peat transmissivity, and microtopography of a small ladder fen for 3 summers (2013–2015) in the James Bay Lowland. The system was manipulated with a sustained hydrological forcing (water addition) to the upslope boundary of the fen during 2014 (38 m3/day) and 2015 (30 m3/day). There was an exponential increase in transmissivity towards the peat surface due to extremely high‐hydraulic conductivities within the upper few centimeters of the peat deposit. At the maximum water table, the saturated hydraulic conductivity of the 0.1 m layer of peat below the water table varied depending on peat microtopography (preferential flow paths = 42–598 m/day and ridges = 16–52 m/day), resulting in high‐hydrological connectivity periods. Furthermore, during 2015, there was an abnormally large amount of precipitation (300 mm vs. long‐term average ~ 100 mm) that resulted in complete surface water connectivity of the site. This caused rapid movement of water from the head of system to the outlet (~15 hr) and runoff ratios >1, compared to low‐water table periods (runoff ratio ~ 0.05). This study highlights the profound importance of the transmissivity–water table feedback mechanism in ladder fens, on controlling the water retention and drainage of large peatland complexes.  相似文献   

3.
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.  相似文献   

4.
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.  相似文献   

5.
Peat soils are heterogeneous, anisotropic porous media. Compared to mineral soils, there is still limited understanding of physical and solute transport properties of fen peat soils. In this study, we aimed to explore the effect of soil anisotropy on solute transport in degraded fen peat. Undisturbed soil cores, taken in vertical and horizontal direction, were collected from one drained and one restored fen peatland both in a comparable state of soil degradation. Saturated hydraulic conductivity (K s) and chemical properties of peat were determined for all soil cores. Miscible displacement experiments were conducted under saturated steady state conditions using potassium bromide as a conservative tracer. The results showed that (1) the K s in vertical direction (K sv) was significantly higher than that in horizontal direction (Ksh), indicating that K s of degraded fen peat behaves anisotropically; (2) pronounced preferential flow occurred in vertical direction with a higher immobile water fraction and a higher pore water velocity; (3) the 5% arrival time (a proxy for the strength of preferential flow) was affected by soil anisotropy as well as study site. A strong correlation was found between 5% arrival time and dispersivity, K s and mobile water fraction; (4) phosphate release was observed from drained peat only. The impact of soil heterogeneity on phosphate leaching was more pronounced than soil anisotropy. The soil core with the strongest preferential flow released the highest amount of phosphate. We conclude that soil anisotropy is crucial in peatland hydrology but additional research is required to fully understand anisotropy effects on solute transport.  相似文献   

6.
Timms WA  Hendry MJ 《Ground water》2008,46(4):616-628
Characterizing and predicting reactive solute transport in low hydraulic conductivity ( K ) clay-rich media is challenging because the very long transport time for solutes renders conventional column tests impractical. In this study, a centrifugation technique was developed to assess the transport of a simple aqueous solution (NaCl) by accelerating flow by centrifugal force through low K (1.1 × 10−11 m/s) core samples. Duplicate cores (52-mm length × 33-mm diameter) were centrifuged at 330 × g for 90 d to model the migration of saline pore water (0.5 M NaCl) under in situ conditions through an approximately 17-m-thick clay prototype over approximately 24,000 years. A PHREEQC one-dimensional reactive solute transport code simulated effluent breakthrough of the NaCl during centrifugation, with best-fit cation exchange coefficients similar to batch tests. The calibrated code was used to predict solute profile development over the long term in the prototype or simulated field-scale conditions. Chromatographic separation of solutes due to ion exchange was evident over several meters in the simulated prototype and the field profile. The applicability of centrifugation methods to predict transport of more complex suites of reactive solutes over the long term is yet to be verified.  相似文献   

7.
Little is known about solute transport in peats, despite the obvious importance of solute transport on eco‐hydrological processes in both managed and natural peatlands. To address this lack of knowledge, we investigated solute transport processes in an agricultural fen peat using a conservative KBr tracer. The main aim of the study was to elucidate solute transport behaviour in general in this peat, with a more specific aim of investigating whether preferential or bypassing flow occurred. The tracer moved through the peat more rapidly than expected, and the pattern of movement showed clear evidence of plot‐scale bypassing flow. The data also provide evidence that bypassing flow occurs in pores at smaller scales. The implications of this study for management of wetland pastures in the Somerset Moors in south‐west England are discussed. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

8.
Numerical simulations of variable-density flow and solute transport have been conducted to investigate dense plume migration for various configurations of 2D fracture networks. For orthogonal fractures, simulations demonstrate that dispersive mixing in fractures with small aperture does not stabilize vertical plume migration in fractures with large aperture. Simulations in non-orthogonal 2D fracture networks indicate that convection cells form and that they overlap both the porous matrix and fractures. Thus, transport rates in convection cells depend on matrix and fracture flow properties. A series of simulations in statistically equivalent networks of fractures with irregular orientation show that the migration of a dense plume is highly sensitive to the geometry of the network. If fractures in a random network are connected equidistantly to the solute source, few equidistantly distributed fractures favor density-driven transport. On the other hand, numerous fractures have a stabilizing effect, especially if diffusive transport rates are high. A sensitivity analysis for a network with few equidistantly distributed fractures shows that low fracture aperture, low matrix permeability and high matrix porosity impede density-driven transport because these parameters reduce groundwater flow velocities in both the matrix and the fractures. Enhanced molecular diffusion slows down density-driven transport because it favors solute diffusion from the fractures into the low-permeability porous matrix where groundwater velocities are smaller. For the configurations tested, variable-density flow and solute transport are most sensitive to the permeability and porosity of the matrix, which are properties that can be determined more accurately than the geometry and hydraulic properties of the fracture network, which have a smaller impact on density-driven transport.  相似文献   

9.
Heterogeneity is prevalent in aquifers and has an enormous impact on contaminant transport in groundwater. Numerical simulations are an effective way to deal with heterogeneity directly by assigning different hydraulic property values to each numerical grid block. Because hydraulic properties vary on different scales, but they cannot be sampled exhaustively and the number of numerical grid blocks is limited by computational considerations, the dispersive effects of unmodeled heterogeneity need to be accounted for. Dispersion tensors can be used to model the dispersion caused by unmodeled heterogeneity. The concept of block-effective macrodispersion tensors for modeling the effects of small-scale variability on solute transport introduced by Rubin et al. [Rubin Y, Sun A, Maxwell R, Bellin A. The concept of block-effective macrodispersivity and a unified approach for grid-scale- and plume-scale-dependent transport. J Fluid Mech 1999;395:161–80] is extended in this paper for use with reactive solutes. The tensors are derived for reactive solutes with spatially variable retardation factors and for solutes experiencing spatially uniform rate-limited sorption. The longitudinal block-effective macrodispersion coefficient is largest for perfect negative correlation between the log-hydraulic conductivity and the retardation factor. Because dispersion tensors, as they are usually implemented in numerical simulations, produce symmetric spreading, the applicability of the concept depends on the portion of the plume asymmetry caused by small-scale variability. The presented results show that the concept is applicable for rate-limited sorption for block sizes of one and two integral scales.  相似文献   

10.
Illman WA  Berg SJ  Yeh TC 《Ground water》2012,50(3):421-431
The main purpose of this paper was to compare three approaches for predicting solute transport. The approaches include: (1) an effective parameter/macrodispersion approach (Gelhar and Axness 1983); (2) a heterogeneous approach using ordinary kriging based on core samples; and (3) a heterogeneous approach based on hydraulic tomography. We conducted our comparison in a heterogeneous sandbox aquifer. The aquifer was first characterized by taking 48 core samples to obtain local-scale hydraulic conductivity (K). The spatial statistics of these K values were then used to calculate the effective parameters. These K values and their statistics were also used for kriging to obtain a heterogeneous K field. In parallel, we performed a hydraulic tomography survey using hydraulic tests conducted in a dipole fashion with the drawdown data analyzed using the sequential successive linear estimator code (Yeh and Liu 2000) to obtain a K distribution (or K tomogram). The effective parameters and the heterogeneous K fields from kriging and hydraulic tomography were used in forward simulations of a dipole conservative tracer test. The simulated and observed breakthrough curves and their temporal moments were compared. Results show an improvement in predictions of drawdown behavior and tracer transport when the K tomogram from hydraulic tomography was used. This suggests that the high-resolution prediction of solute transport is possible without collecting a large number of small-scale samples to estimate flow and transport properties that are costly to obtain at the field scale.  相似文献   

11.
Zheng C  Gorelick SM 《Ground water》2003,41(2):142-155
Several recent studies at the Macrodispersion Experiment (MADE) site in Columbus, Mississippi, have indicated that the relative preferential flowpaths and flow barriers resulting from decimeter-scale aquifer heterogeneities appear to have a dominant effect on plume-scale solute transport. Numerical experiments are thus conducted in this study to explore the key characteristics of solute transport in two-dimensional flow fields influenced by decimeter-scale preferential flowpaths. A hypothetical but geologically plausible network of 10 cm wide channels of high hydraulic conductivity is used to represent the relative preferential flowpaths embedded in an otherwise homogeneous aquifer. When the hydraulic conductivity in the channels is 100 times greater than that in the remaining portion of the aquifer, the calculated concentration distributions under three source configurations all exhibit highly asymmetrical, non-Gaussian patterns. These patterns, with peak concentrations close to the source and extensive spreading downgradient, resemble that observed at the MADE site tracer tests. When the contrast between the channel and nonchannel hydraulic conductivities is reduced to 30:1 from 100:1, the calculated mass distribution curve starts to approach a Gaussian one with the peak concentration near the central portion of the plume. Additional analysis based on a field-scale model demonstrates that the existence of decimeter-scale preferential flowpaths can have potentially far-reaching implications for ground water remediation. Failure to account for them in numerical simulation could lead to overestimation of the effectiveness of the remedial measure under consideration.  相似文献   

12.
Solute discharge moments (mean and variance) are computed using numerical modeling of flow and advective transport in two-dimensional heterogeneous aquifers and are compared to theoretical results. The solute discharge quantifies the temporal evolution of the total contaminant mass crossing a certain compliance boundary. In addition to analyzing the solute discharge moments within a classical absolute dispersion framework, we also analyze relative dispersion formulation, whereby plume meandering (deviation from mean flow path caused by velocity variations at scales larger than plume size) is removed. This study addresses some important issues related to the computation of solute discharge moments from random walk particle tracking experiments, and highlights some of the important differences between absolute and relative dispersion frameworks. Relative dispersion formulation produces maximum uncertainty that coincides with the peak mean discharge. Absolute dispersion, however, results in earlier arrival of the uncertainty peak as compared to the first moment peak. Simulations show that the standard deviation of solute discharge in a relative dispersion framework requires increasingly large temporal sampling windows to smooth out some of the large fluctuations in breakthrough curves associated with advective transport. Using smoothing techniques in particle tracking to distribute the particle mass over a volume rather than at a point significantly reduces the noise in the numerical simulations and removes the need to use large temporal windows. Same effect can be obtained by adding a local dispersion process to the particle tracking experiments used to model advective transport. The effect of the temporal sampling window bears some relevance and important consequences for evaluating risk-related parameters. The expected value of peak solute discharge and its standard deviation are very sensitive to this sampling window and so will be the risk distribution relying on such numerical models.  相似文献   

13.
Despite that discrete flow features (DFFs, e.g. fractures and faults) are common features in the subsurface, few studies have explored the influence of DFFs on solute plumes in otherwise permeable rocks (e.g. sandstone, limestone), compared to low-permeability rock settings (e.g. granite and basalt). DFFs can provide preferential flow pathways (i.e. ‘preferential flow features’; PFFs), or can act to impede flow (i.e. ‘barrier flow features’; BFFs). This research uses a simple analytical expression and numerical modelling to explore how a single DFF influences the steady-state distributions of solute plumes in permeable aquifers. The analysis quantifies the displacement and widening (or narrowing) of a steady-state solute plume as it crosses a DFF in idealised, 1 × 1 m moderately permeable rock aquifers. Previous research is extended by accounting for DFFs as 2D flow features, and including BFF situations. A range of matrix-DFF permeability ratios (0.01 to 100) and DFF apertures (0.25 mm to 2 cm), typical of sedimentary aquifers containing medium-to-large fractures, are considered. The results indicate that for the conceptual models considered here, PFFs typically have a more significant influence on plume distributions than BFFs, and the impact of DFFs on solute plumes generally increases with increasing aperture. For example, displacement of peak solute concentration caused by DFFs exceeds 20 cm in some PFF cases, compared to a maximum of 0.64 cm in BFF cases. PFFs widen plumes up to 9.7 times, compared to a maximum plume widening of 2.0 times in BFF cases. Plumes crossing a PFF are less symmetrical, and peak solute concentrations beneath PFFs are up to two orders of magnitude lower than plumes in BFF cases. This study extends current knowledge of the attenuating influence of DFFs in otherwise permeable rocks on solute plume characteristics, through evaluation of 2D flow effects in DFFs for a variety of DFF apertures, and by considering BFF situations.  相似文献   

14.
A number of experimental studies have tackled the issue of solute transport parameter assessments either in the laboratory or in the field. But yet, the behavior of a plume in the field under density driven forces, is not well known due to possible development of instabilities. Some field tracer tests on the fate of plumes denser than native groundwater such as those encountered under waste disposal facilities, have pointed out the processes of sinking and splitting at the early stage of migration. The process of dispersion was widely investigated, but the range of dispersivity values obtained from either experimental tests, or numerical and theoretical calculations is still very large, even for the same type of aquifers. These discrepancies were considered to be essentially caused by soil heterogeneities and scale effects. In the meantime, studies on the influence of sinking and fingering have remained more scarce. The objective of the work is to analyze how transport parameters such as dispersivities can be affected by unstable conditions, which lead to plume sinking and fingering. A series of tracer tests were carried out to study under natural conditions, the transport of a dense chloride solution injected in a shallow two-layered aquifer. Two types of experiments were performed: in the first type, source injection was such that the plume could travel downward from one layer to the other of higher pore velocity, and in the second one, the migration took place only in the faster layer. The results suggest some new insights in the processes occurring at the early stages of a dense plume migration moving in a stratified aquifer under groundwater fluctuations, which can be summarized through the following points: (i) Above a stability criterion threshold, a fingering process and a multi modal plume transport take place, but local dispersivities can be cautiously derived, using breakthrough curves matching. (ii) When water table is subject to some cycling or rising, the plume can be significantly distorted in the transverse direction, leading to unusual values of the ratio between longitudinal and transverse dispersivities. (iii) Under stable conditions, for example in the case of straightforward injection in the faster aquifer layer, longitudinal dispersivity is greater than the transverse component as usually encountered, and the obtained transport parameters are closed to macro dispersivity values, which reach their asymptotic limit at very short distances. (iv) The classical scale effect about the varying dispersivity at short distances could be a process mainly due to the distance required for a plume stabilization.  相似文献   

15.
A comprehensive numerical study was undertaken to investigate transport of a variable-density, conservative solute plume in an unconfined coastal aquifer subject to high and low frequency oceanic forcing. The model combined variable-density saturated flow for groundwater and solute transport, and wave hydrodynamics from a 2D Navier–Stokes solver. A sinusoidal tidal signal was specified by implementing time-varying heads at the seaward boundary. The solute plume behavior was investigated under different oceanic forcing conditions: no forcing, waves, tide, and combined waves and tide. For each forcing condition, four different injected solute densities (freshwater, brackish water, seawater, brine) were used to investigate the effects of density on the transport of the injected plume beneath and across the beach face. The plume’s low-order spatial moments were computed, viz., mass, centroid, variance and aspect ratio. The results confirmed that both tide- and wave-forcing produce an upper saline plume beneath the beach face in addition to the classical saltwater wedge. For the no-forcing and tide-only cases (during rising tides), an additional small circulation cell below the beach face was observed. Oceanic forcing affects strongly the solute plume’s flow path, residence time and discharge rate across the beach face, as well as its spreading. For the same oceanic forcing, solute plumes with different densities follow different trajectories from the source to the discharge location (beach face). The residence time and plume spreading increased with plume density. It was concluded that simulations that neglect the effect of waves or tides cannot reproduce accurately solute plume dispersion and also, in the case of coasts with small waves or tides, the solute residence time in the aquifer.  相似文献   

16.
Study of the mobility of contaminants in an aquifer is an important issue for the proper remediation of contaminated groundwater. Determination of associated solute transport parameters therefore is essential for investigation of the extent to which groundwater can be contaminated. This study aimed at determining solute transport parameters for an unconfined sandy aquifer at a laboratory scale through various tracer tests using a conservative solute as a tracer. Tracer tests consisted of both well‐tracer tests (single and double wells) and an aquifer tracer test using a plume‐capturing device such as time domain reflectometry (TDR). The results showed that longitudinal dispersivities estimated from the single and double well‐tracer tests were 2·2 cm and 13·5 cm for a travel distance of 9·3 cm and 13·5 cm from the injection point respectively. These results agreed reasonably well with the results of the aquifer tracer test. The solute transport parameters obtained at multiple points in the aquifer through the aquifer tracer test revealed that the dispersivity length was proportional to the travel distance by a factor of 0·3, which was moderately higher than the value of 0·1 given in the literature. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

17.
  相似文献   

18.
We analyze the impact of a linear trend in the mean log-conductivity on the transport of a conservative tracer in a bounded domain. The effects of such a linear trend on solute transport were analyzed in depth for unbounded domains (Rubin and Seong, Water Resour Res 30(11):2901–2911, 1994; Indelman and Rubin, Water Resour Res 31(5):1257–1265, 1995; Water Resour Res 32(5):1257–1265, 1996), whereas studies concerning this special case of medium nonstationarity in finite domains usually focus on head or flow statistics (Guadagnini et al., Stoch Environ Res Risk Assess, 17:394–407, 2003). In this study both ensemble and effective plume moments are provided for an instantaneous release of a solute through a linear source normal to the mean flow direction, by taking into account different sizes of the source. The analysis involving a steady velocity field spatially nonstationary is developed by using the stochastic finite element method. Results show that ensemble moments are affected by increasing trends both parallel and normal to the mean flow direction, but the impact on effective plume moments is very different. A parallel trend does not seem to influence the effective second moments; while a normal trend, although modifies the transverse effective moment only weakly, strongly increases the longitudinal one, especially for large initial sizes of the source. Furthermore, the increase of the particle displacement variance produced by a parallel trend in the finite domain disagrees with the results obtained in an unbounded domain, due to the boundary conditions here considered making both head and velocity moments nonstationary and nonsymmetric.  相似文献   

19.
A 3D ERT study of solute transport in a large experimental tank   总被引：2，自引：0，他引：2  
A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m3. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes.  相似文献   

20.
Using a slab of Massillon Sandstone, laboratory-scale solute tracer experiments were carried out to test numerical simulations using the Advection–Dispersion Equation (ADE). While studies of a similar nature exist, our work differs in that we combine: (1) experimentation in naturally complex geologic media, (2) X-ray absorption imaging to visualize and quantify two-dimensional solute transport, (3) high resolution transport property characterization, with (4) numerical simulation. The simulations use permeability, porosity, and solute concentration measured to sub-centimeter resolution. While bulk breakthrough curve characteristics were adequately matched, large discrepancies exist between the experimental and simulated solute concentration fields. Investigation of potential experimental errors suggests that the failure to fit solute concentration fields may lie in loss of intricate connectivity within the cross-bedded sandstone occurring at scales finer than our property characterization measurements (i.e., sub-centimeter).  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号