首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Experiments designed to elucidate the pore-scale mechanisms of the dissolution of a residual non-aqueous phase liquid (NAPL), trapped in the form of ganglia within a porous medium, are discussed. These experiments were conducted using transparent glass micromodels with controlled pore geometry, so that the evolution of the size and shape of individual NAPL ganglia and, hence, the pore-scale mass transfer rates and mass transfer coefficients could be determined by image analysis. The micromodel design permitted reasonably accurate control of the pore water velocity, so that the mass transfer coefficients could be correlated in terms of a local (pore-scale) Peclet number. A simple mathematical model, incorporating convection and diffusion in a slit geometry was developed and used successfully to predict the observed mass transfer rates. For the case of non-wetting NAPL ganglia, water flow through the corners in the pore walls was seen to control the rate of NAPL dissolution, as recently postulated by Dillard and Blunt [Water Resour. Res. 36 (2000) 439–454]. Break-up of doublet non-wetting phase ganglia into singlet ganglia by snap-off in pore throats was also observed, confirming the interplay between capillarity and mass transfer. Additionally, the effect of wettability on dissolution mass transfer was demonstrated. Under conditions of preferential NAPL wettability, mass transfer from NAPL films covering the solid surfaces was seen to control the dissolution process. Supply of NAPL from the trapped ganglia to these films by capillary flow along pore corners was observed to result in a sequence of pore drainage events that increase the interfacial area for mass transfer. These observations provide new experimental evidence for the role of capillarity, wettability and corner flow on NAPL ganglia dissolution.  相似文献   

2.
Wettability profoundly affects not only the initial distribution of residual NAPL contaminants in natural soils, but also their subsequent dissolution in a flowing aqueous phase. Under conditions of preferential NAPL wettability, the residual NAPL phase is found within the smaller pores and in the form of continuous corner filaments and thick films on pore walls. Such films expose a much greater interfacial area for mass transfer than would be exposed by the same amount of non-wetting NAPL. Importantly, capillary and hydraulic continuity of NAPL filaments and thick films is essential for sustaining NAPL–water counterflow during the course of NAPL dissolution in flowing groundwater—a mechanism which maintains and even increases the interfacial area for mass transfer. Continued dissolution results in gradual thinning of the NAPL films, which may become unstable and rupture causing disconnection of the residual NAPL in the form of clusters. Using a pore network simulator, we demonstrate that NAPL film instability drastically modifies the microscopic configuration of residual NAPL, and hence the local hydrodynamic conditions and interfacial area for mass transfer, with concomitant effects on macroscopically observable quantities, such as the aqueous effluent concentration and the fractional NAPL recovery with time. These results strongly suggest that the disjoining pressure of NAPL films may exert an important, and hitherto unaccounted, control on the dissolution behaviour of a residual NAPL phase in oil wet systems.  相似文献   

3.
We investigated the dissolution of non-aqueous phase liquids (NAPLs) in a three-dimensional random sphere-pack medium using a pore-scale modeling approach to advance fundamental understanding and connect rigorously to microscale processes. Residual NAPL distributions were generated using a morphological approach and the entrapped non-wetting phase was quantitatively characterized by calculating volume, orientation, interfacial area, and shape of isolated NAPL regions. With a detailed aqueous-phase flow field obtained by a multiple-relaxation time lattice Boltzmann approach, we solved the advective–diffusive equation in the pore space using a high-resolution, adaptive-stencil finite-volume scheme and an operator-splitting algorithm. We show good agreement between the mass transfer rates predicted in the computational approach and previously published experimental observations. The pore-scale simulations presented in this work provide the first three-dimensional comparison to the considerable experimental work that has been performed to derive constitutive relations to quantify mass transfer from a residual NAPL to a flowing aqueous phase.  相似文献   

4.
Partitioning interwell tracer tests (PITTs) are a relatively new technique for measuring the amount of nonaqueous phase liquid (NAPL) within saturated porous media. In this work we examined the influence of mass transfer limitations on the accuracy of measured NAPL from PITTs. Two mathematical models were used along with laboratory column experiments to explore the influence of tracer partition coefficient, tracer detection limit, and injected tracer mass on NAPL measurements. When dimensionless mass transfer coefficients were small, NAPL measurement errors decreased with decreasing tracer partition coefficient, decreasing tracer detection limit, and increasing injected tracer mass. Extrapolating breakthrough curves exponentially reduced but did not eliminate systematic errors in NAPL measurement. Although transport in a single stream tube was used in the mathematical models and laboratory experiments, the results from this simplified domain were supported by data taken from a three-dimensional computational experiment, where the NAPL resided as large pool. Based on these results, we suggest guidelines for interpreting tracer breakthrough data to ascertain the importance of mass transfer limitations on NAPL measurements.  相似文献   

5.
Nonaqueous phase liquid (NAPL) is a long-term source of ground water contamination as the pollutant slowly partitions into the air and water phases. The objective of this work was to study the efficacy of aqueous surfactant solution to enhance the dissolution of a residual NAPL below the capillary fringe, hence reducing the time needed for aquifer restoration. An analytical technique was developed to measure the concentration of NAPL in a nonionic surfactant. Soil column experiments simulated conditions in the saturated soil where a NAPL may become trapped as a discontinuous immobile phase. Experimental results indicate that dissolution was a rate-limited process, approaching equilibrium concentrations after 24 hours. The relative permeability of the aqueous phase initially decreased as surfactant was injected, but increased over time as the saturation of residual NAPL was reduced through mass transfer into the surfactant-enhanced aqueous phase. These findings suggest that enhancing the aqueous phase with a nonionic surfactant may significantly enhance the in situ recovery or residual NAPL.  相似文献   

6.
Interphase mass transfer in porous media is commonly modeled using Sherwood number expressions that are developed in terms of fluid and porous medium properties averaged over some representative elementary volume (REV). In this work the influence of sub-grid scale properties on interphase mass transfer was investigated using a two-dimensional pore network model. The focus was on assessing the impact of (i) NAPL saturation, (ii) interfacial area (iii) NAPL spatial distribution at the pore scale, (iv) grain size heterogeneity, (v) REV or domain size and (vi) pore scale heterogeneity of the porous media on interphase mass transfer. Variability of both the mass transfer coefficient that explicitly accounts for the interfacial area and the mass transfer coefficient that lumps the interfacial area was examined. It was shown that pore scale NAPL distribution and its orientation relative to the flow direction have significant impact on flow bypassing and the interphase mass transfer coefficient. This results in a complex non-linear relationship between interfacial area and the REV-based interphase mass transfer rate. Hence, explicitly accounting for the interfacial area does not eliminate the uncertainty of the mass transfer coefficient. It was also shown that, even for explicitly defined flow patterns, changing the domain size over which the mass transfer process is defined influences the extent of NAPL bypassing and dilution and, consequently, the interphase mass transfer. It was also demonstrated that the spatial variability of pore scale parameters such as pore throat diameters may result in different rates of interphase mass transfer even for the same pore size distribution index.  相似文献   

7.
The generation of vapor‐phase contaminant plumes within the vadose zone is of interest for contaminated site management. Therefore, it is important to understand vapor sources such as non‐aqueous‐phase liquids (NAPLs) and processes that govern their volatilization. The distribution of NAPL, gas, and water phases within a source zone is expected to influence the rate of volatilization. However, the effect of this distribution morphology on volatilization has not been thoroughly quantified. Because field quantification of NAPL volatilization is often infeasible, a controlled laboratory experiment was conducted in a two‐dimensional tank (28 cm × 15.5 cm × 2.5 cm) with water‐wet sandy media and an emplaced trichloroethylene (TCE) source. The source was emplaced in two configurations to represent morphologies encountered in field settings: (1) NAPL pools directly exposed to the air phase and (2) NAPLs trapped in water‐saturated zones that were occluded from the air phase. Airflow was passed through the tank and effluent concentrations of TCE were quantified. Models were used to analyze results, which indicated that mass transfer from directly exposed NAPL was fast and controlled by advective‐dispersive‐diffusive transport in the gas phase. However, sources occluded by pore water showed strong rate limitations and slower effective mass transfer. This difference is explained by diffusional resistance within the aqueous phase. Results demonstrate that vapor generation rates from a NAPL source will be influenced by the soil water content distribution within the source. The implications of the NAPL morphology on volatilization in the context of a dynamic water table or climate are discussed.  相似文献   

8.
In this work, we have applied a group-contribution activity-coefficient model, UNIFAC, and the solubility of alcohols in water to estimate partition coefficients for alcohol tracers between water and nonaqueous-phase liquids (NAPLs). The effects of temperature and mutual solubility between NAPL and aqueous phases on the estimation of partition coefficients were also investigated. By comparing the estimated results with experimental partition coefficients for 30 alcohol tracers between 10 NAPLs and water, we found that: i) the UNIFAC-solubility method, in which the UNIFAC model in its infinite-dilution form is applied to the NAPL phase and the solubility of tracers in water is used for estimation of the activity coefficient in the aqueous phase, works better than the UNIFAC model; ii) a linear relation between the logarithm of partition coefficients and the logarithm of tracer solubility in water is observed for those tracers having a similar chemical structure (i.e. the same number of branched methyl groups). This can serve as a useful tool for quick selection of the tracers that exhibit the desired partition coefficients; iii) the effect of mutual solubility between NAPL and aqueous phases can be neglected because such miscibility is very small, usually of the order of 10−3 mole/mole unit; and iv) temperature variation between 15° and 25°C does not significantly affect partition coefficients.  相似文献   

9.
 A stochastic simulation is performed to study multiphase flow and contaminant transport in fractal porous media with evolving scales of heterogeneity. Numerical simulations of residual NAPL mass transfer and subsequent transport of dissolved and/or volatilized NAPL mass in variably saturated media are carried out in conjunction with Monte Carlo techniques. The impact of fractal dimension, plume scale and anisotropy (stratification) of fractal media on relative dispersivities is investigated and discussed. The results indicate the significance of evolving scale of porous media heterogeneity to the NAPL transport in the subsurface. In general, the fractal porous media enhance the dispersivities of NAPL mass plume transport in both the water phase and the gas phase while the influence on the water phase is more significant. The porous media with larger fractal dimension have larger relative dispersivities. The aqueous horizontal dispersivity exhibits a most significant increase against the plume scale.  相似文献   

10.
Laboratory experiments and numerical simulations in homogeneous porous media were used to investigate the influence of porous medium wettability on the formation and growth of preferential dissolution pathways, dissolution fingers, during nonaqueous phase liquid (NAPL) dissolution. As the porous medium became increasingly NAPL-wet, dissolution fingers grew wider and slower. This result was observed in physical experiments with 0% and 100% NAPL-wet conditions and confirmed with numerical simulations at these and intermediate wettabilities. A previously derived expression for an upscaled mass transfer rate coefficient that accounts for the growth of dissolution fingers was used to quantify the effect of fingering on overall NAPL removal rates. For the test cases evaluated, NAPL dissolution fingering controlled the overall rate of NAPL dissolution after the dissolution front moved 4 cm in 0% NAPL-wet conditions and 18 cm in 100% NAPL-wet conditions. Thus, even in completely NAPL-wet media dissolution fingering may control the overall rate of NAPL dissolution after relatively short travel distances. The importance of NAPL dissolution fingering in heterogeneous systems with spatially varying NAPL saturations, though, remains an important question for future work.  相似文献   

11.
This study presents a multiphase flow and multispecies reactive transport model for the simultaneous simulation of NAPL and groundwater flow, dissolution, and reactive transport with isotope fractionation, which can be used for better interpretation of NAPL-involved Compound Specific Isotope Analysis in 3D heterogeneous hydrogeologic systems. The model was verified for NAPL-aqueous phase equilibrium partitioning, aqueous phase multi-chain and multi-component reactive transport, and aqueous phase multi-component transport with isotope fractionation. Several illustrative examples are presented to investigate the effect of DNAPL spill rates, degradation rate constants, and enrichment factors on the temporal and spatial distribution of the isotope signatures of chlorinated aliphatic hydrocarbon groundwater plumes. The results clearly indicate that isotope signatures can be significantly different when considering multiphase flow within the source zone. A series of simulations indicate that degradation and isotope enrichment compete with dissolution to determine the isotope signatures in the source zone: isotope ratios remain the same as those of the source if dissolution dominates the reaction, while heavy isotopes are enriched in reactants along groundwater plume flow paths when degradation becomes dominant. It is also shown that NAPL composition can change from that of the injected source due to the partitioning of components between the aqueous and NAPL phases even when degradation is not allowed in NAPL phase. The three-dimensional simulation is presented to mechanistically illustrate the complexities in determining and interpreting the isotopic signatures with evolving DNAPL source architecture.  相似文献   

12.
Neat ethanol (75.7 L) was released into the upper capillary zone in a continuous-flow, sand-packed aquifer tank (8.2 m3) with an average seepage velocity of 0.75 m/day. This model aquifer system contained a residual nonaqueous phase liquid (NAPL) that extended from the capillary zone to 10 cm below the water table. Maximum aqueous concentrations of ethanol were 20% v/v in the capillary zone and 0.08% in the saturated zone at 25 and 30 cm downgradient from the emplaced NAPL source, respectively. A bench-scale release experiment was also conducted for a similar size spill (scaled to the plan area). The concentrations of ethanol in ground water for both the bench- and pilot-scale experiments were consistent with advective–dispersive limited mass transfer from the capillary to the saturated zone. Concentrations of monoaromatic hydrocarbons and isooctane increased in the pore water of the capillary zone as a result of both redistribution of residual NAPL (confirmed by visualization) and enhanced hydrocarbon dissolution due to the cosolvent effect exerted by ethanol. In the tank experiment, higher hydrocarbon concentrations in ground water were also attributed to decreased hydrocarbon biodegradation activity caused by preferential microbial utilization of ethanol and the resulting depletion of oxygen. These results infer that spills of highly concentrated ethanol will be largely confined to the capillary zone due to its buoyancy, and ethanol concentrations in near-source zone ground water will be controlled by mass transfer limitations and hydrologic conditions. Furthermore, highly concentrated ethanol releases onto pre-existing NAPL will likely exacerbate impacts to ground water, due to NAPL mobilization and dissolution, and decreased bioattenuation of hydrocarbons.  相似文献   

13.
A non-equilibrium, two-phase, three-component compositional model for the simulation of alcohol flooding has been developed and tested. Inter-phase mass transfer algorithms allow for transfer of all three components at high concentrations and high mass flux rates using a two-film model. The model has been used to simulate alcohol floods where the alcohol has an affinity for either the water-rich phase, or the organic-rich phase. Calibration, using experimental effluent data from an alcohol flood which used a 2-propanol (IPA)-water-tetrachlorethene (PCE) ternary system, indicates that inter-phase mass transfer parameters can be non-unique. Sensitivity studies, completed using the non-equilibrium model for the IPA-water-PCE system, indicate that experimentally derived organic-rich phase composition data should lead to better estimates of the non-wetting phase film thickness. For alcohol flooding experiments where the primary mechanism of non-aqueous phase liquid (NAPL) removal is enhanced dissolution, near-equilibrium conditions may be achieved with NAPL recovery similar for conditions of near-equilibrium and equilibrium. However, for systems where remobilization is the primary mechanism of NAPL recovery, it is expected that although local conditions may approach equilibrium, the resulting NAPL recovery can be significantly lower than would be attained if equilibrium conditions persisted.  相似文献   

14.
《Advances in water resources》2005,28(11):1143-1158
Biodegrading plumes in groundwater are often typified by relatively reactive zones around the fringes and less reactive zones in the core. A high degree of refinement is required at the fringes if a model is to be of use in improving the conceptual understanding of plumes. Two strategies for dealing with the potentially high computational demands are (i) parallel processing, where the workload is shared between multiple processors, and (ii) locally adaptive remeshing, where a refined area of the grid tracks the moving plume fringes through the domain. The partial differential equation toolbox, UG (Unstructured Grids) offers advanced numerical tools including adaptive remeshing, sparse matrix storage schemes, and multigrid solvers. It embraces many of these features within a parallel processing environment. This paper reports on a recent development of UG to simulate field scale reactive biogeochemistry including Monod kinetics, NAPL dissolution, mineral precipitation/dissolution and ion exchange. The non-linear multicomponent reactive transport system is solved with the fully coupled method. Test cases have been used for verification of the new capability. The paper illustrates an application to a 3D field site. It is demonstrated that both adaptive remeshing and parallel processing can improve efficiency and in turn facilitate the incorporation of a more complex set of species and reactions such that understanding of plume processes is enhanced.  相似文献   

15.
Innovative remediation studies were conducted between 1994 and 2004 at sites contaminated by nonaqueous phase liquids (NAPLs) at Hill and Dover AFB, and included technologies that mobilize, solubilize, and volatilize NAPL: air sparging (AS), surfactant flushing, cosolvent flooding, and flushing with a complexing-sugar solution. The experiments proved that aggressive remedial efforts tailored to the contaminant can remove more than 90% of the NAPL-phase contaminant mass. Site-characterization methods were tested as part of these field efforts, including partitioning tracer tests, biotracer tests, and mass-flux measurements. A significant reduction in the groundwater contaminant mass flux was achieved despite incomplete removal of the source. The effectiveness of soil, groundwater, and tracer based characterization methods may be site and technology specific. Employing multiple methods can improve characterization. The studies elucidated the importance of small-scale heterogeneities on remediation effectiveness, and fomented research on enhanced-delivery methods. Most contaminant removal occurs in hydraulically accessible zones, and complete removal is limited by contaminant mass stored in inaccessible zones. These studies illustrated the importance of understanding the fluid dynamics and interfacial behavior of injected fluids on remediation design and implementation. The importance of understanding the dynamics of NAPL-mixture dissolution and removal was highlighted. The results from these studies helped researchers better understand what processes and scales are most important to include in mathematical models used for design and data analysis. Finally, the work at these sites emphasized the importance and feasibility of recycling and reusing chemical agents, and enabled the implementation and success of follow-on full-scale efforts.  相似文献   

16.
Organic contaminants present as nonaqueous phase liquids (NAPLs) in the subsurface often pose a long-term risk to human health and the environment. Investigating the distribution of NAPLs in porous media remains a major challenge in risk assessment and management of contaminated sites. Conventional soil coring and monitoring wells have been widely used over past decades as the primary means of subsurface investigation to determine NAPL extent. Known limitations of conventional approaches have led us to explore an alternative or a complementary technique to provide high-quality information of NAPL source zone architecture. This work advances an imaging tool for a variety of organic NAPL contaminants in unconsolidated soils through magnetic resonance imaging (MRI) of frozen cores. Using trichloroethylene (TCE) and o-xylene as model species, we illustrate that discriminatory freezing of water, while keeping the NAPL in a liquid state, enables high-resolution qualitative delineation of NAPL distribution within porous media. This novel approach may help improve site conceptual models and consequentially lead to highly tailored, more efficient remedial measures.  相似文献   

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

18.
A macroscopic transport model is developed, following the Taylor shear dispersion analysis procedure, for a 2D laminar shear flow between parallel plates possessing a constant specified concentration. This idealized geometry models flow with contaminant dissolution at pore-scale in a contaminant source zone and flow in a rock fracture with dissolving walls. We upscale a macroscopic transient transport model with effective transport coefficients of mean velocity, macroscopic dispersion, and first-order mass transfer rate. To validate the macroscopic model the mean concentration, covariance, and wall concentration gradient are compared to the results of numerical simulations of the advection–diffusion equation and the Graetz solution. Results indicate that in the presence of local-scale variations and constant concentration boundaries, the upscaled mean velocity and macrodispersion coefficient differ from those of the Taylor–Aris dispersion, and the mass transfer flux described by the first-order mass transfer model is larger than the diffusive mass flux from the constant wall. In addition, the upscaled first-order mass transfer coefficient in the macroscopic model depends only on the plate gap and diffusion coefficient. Therefore, the upscaled first-order mass transfer coefficient is independent of the mean velocity and travel distance, leading to a constant pore-scale Sherwood number of 12. By contrast, the effective Sherwood number determined by the diffusive mass flux is a function of the Peclet number for small Peclet number, and approaches a constant of 10.3 for large Peclet number.  相似文献   

19.
A two-dimensional numerical transport model is developed to determine the effect of aquifer anisotropy and heterogeneity on mass transfer from a dense nonaqueous phase liquid (DNAPL) pool. The appropriate steady state groundwater flow equation is solved implicitly whereas the equation describing the transport of a sorbing contaminant in a confined aquifer is solved by the alternating direction implicit method. Statistical anisotropy in the aquifer is introduced by two-dimensional, random log-normal hydraulic conductivity field realizations with different directional correlation lengths. Model simulations indicate that DNAPL pool dissolution is enhanced by increasing the mean log-transformed hydraulic conductivity, groundwater flow velocity, and/or anisotropy ratio. The variance of the log-transformed hydraulic conductivity distribution is shown to be inversely proportional to the average mass transfer coefficient.  相似文献   

20.
Spills and releases of hydrocarbons may result in zones of nonaqueous phase liquids (NAPL) within soils and groundwater. The NAPL will change, or “weather” over time due to a range of physical, chemical, and biological processes. Hydrocarbon constituents in environmental samples collected from NAPL-impacted groundwater wells, sediments, or soils vary in composition over time due to weathering. The changing composition can be used to estimate mass depletion rates and trends for the bulk NAPL and for individual constituent chemicals relative to marker constituents which are less susceptible to weathering. Methods for the selection of marker constituents and for quantitatively and conservatively estimating NAPL depletion rates and trends over time are shown. Estimates are included for two sites with different NAPL mixtures present (crude oil and gasoline/diesel-range product), for which depletion of half the initial total NAPL is estimated at 13.6 ± 2.9 years and 7.3 ± 1.8 years respectively, and with no active NAPL remediation at either site. Similar methods for oil or NAPL depletion estimates have often relied on a prior-identified suite of presumed-conserved marker constituents. The method presented here includes steps which identify the best set of analyzed candidate marker constituents in a NAPL mixture. This can confirm prior-selected markers but is particularly useful for NAPL mixtures in which no prior-identified marker constituents are present.  相似文献   

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

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