Transport of dissolved methoxy-nonafluorobutane through a saturated column

Solute transport experiments were conducted in a one-dimensional saturated column using dissolved methoxy-nonafluorobutane (HFE-7100), a Novec engineered fluid developed by the 3M Corporation, as the solute. Novec engineered fluids are considered dense non-aqueous phase liquids (DNAPLs) because they are immiscible with water and have a specific gravity greater than one. The HFE-7100 fluid is safer and environmentally friendlier than common DNAPL contaminants such as tetrachloroethylene (PCE) or trichloroethylene (TCE); thus, it is an ideal substitute DNAPL for laboratory groundwater contamination research. Three sets of solute transport experiments were conducted. The first set of experiments was conducted in a glass-bead-packed column using dissolved HFE-7100 as the solute. The second set of experiments was conducted in a sand-packed column using dissolved HFE-7100 as the solute. The third set of experiments was conducted in a sand-packed column using dissolved PCE as the solute. The dissolved HFE-7100 column breakthrough concentrations were compared with dissolved PCE breakthrough concentrations. Results show that the one-dimensional solute transport equation was successful in describing the transport behavior of dissolved HFE-7100. This study demonstrates that the HFE-7100 fluid can be used as a safer substitute DNAPL for groundwater contaminant dissolution and transport research.     

An analytical solution to multi-component NAPL dissolution equations

This note presents a novel method for determining the changing composition of a multi-component NAPL body dissolving into moving groundwater, and the consequent changes in the aqueous phase solute concentrations in the surrounding pore water. A canonical system of coupled non-linear governing equations is derived which is suitable for representation of both pooled and residual configurations, and this is solved. Whereas previous authors have handled such problems numerically, it is shown that these governing equations succumb to analytical solution. By a suitable substitution, the equations become decoupled, and the problem collapses to a single first-order equation. The final result is expressed implicitly, with time as a function of the number of moles of the least soluble component, m₁. The number of moles of each other component is expressed explicitly in terms of m₁. It is shown that the time-m₁ relationship has a well behaved inverse. An example is given in which the analytic solution is verified against traditional finite difference analysis, and its computational efficiency is shown.     

Spreading of a NAPL lens in a double-porosity medium

Spreading of a non-aqueous phase liquid (NAPL) denser than water (DNAPL) lens (mound) in the unsaturated zone of double-porosity aquifer above an impervious plane boundary is investigated. The double-porosity aquifer is conceptualized as a fracture network surrounding pervious blocks. Vertical gravity equilibrium is assumed to prevail in each one of the two media, fractures and blocks. Through vertical integration, two coupled partial differential equations for the DNAPL content in each medium are obtained. The mass exchange rate between high- and low-permeability media is considered as a function of NAPL content. The dominant effect is gravity, whereas capillary forces are negligible. Analytical solutions for one-dimensional and axisymmetric problems of mound spreading are obtained.     

Effect of NAPL exposure on the wettability and two‐phase flow in a single rock fracture

Surface‐wetting properties are an important cause of changing the groundwater and two‐phase fluid flows. Various factors affecting the surface wettability were investigated in a parallel‐walled glass fracture with non‐aqueous phase liquid (NAPL) (gasoline, diesel, trichloroethylene, and creosote) wetted surfaces. First, the effect of the duration of NAPL exposure on wettability change was considered at pre‐wet fracture surfaces using the various NAPL species, and the result showed that the surface became hydrophobic after the exposure time of NAPL exceeded 2000 min. Second, the initial wetting state of the surface affected the timing when the wettability change begins as well as the extent of the wettability change in an NAPL‐wetted rock fractures. Under the dry condition, the wettability change was completed within a very short time of exposure to NAPL (~5 min), and then it finally reached the intermediate and weakly NAPL wetting (contact angle of 118°). Under the pre‐wet condition, a relatively long time of exposure (~5000 min) was needed to observe the obvious change of the surface wettability, which was changed up to strongly NAPL wetting (contact angle of 142°). Third, the wettability changed by NAPL exposure was stable and maintained for a long time, regardless of water flushing rate and temperature. Finally, the wettability change by the exposure of NAPL on parallel fracture surfaces was evaluated at various groundwater flow velocities. Result showed that groundwater flow velocity has an important impact upon measured contact angle. Although fracture surfaces were exposed to NAPL at the low groundwater flow velocity, the wettability was not changed from hydrophilic to hydrophobic when the contact time between NAPL and mineral surfaces was not sufficient owing to the pulse‐type movement of NAPL. This implies that the variation of exposure pattern due to groundwater flow on the wettability change can be an important factor affecting the wettability change of fracture surface and migration behaviour at natural fractured rock aquifers in case of NAPL spill. Copyright © 2015 John Wiley & Sons, Ltd.     

治理地下水有机污染抽出处理技术影响因素分析

抽出处理技术是治理地下水有机污染简便易行的方法,文中的简要介绍抽出处理方法去除污染机制之后,着重分析了讨论了导致该方法去除地下水有机污染耗时及运行费用高昂的原因,在此基础上,对地下水有机污染抽出处理技术作出全面详尽的评述。     

Numerical modeling of three-dimensional contaminant migration from dissolution of multicomponent NAPL pools in saturated porous media

A three-dimensional model for contaminant transport resulting from the dissolution of multicomponent nonaqueous phase liquid (NAPL) pools in three-dimensional saturated subsurface formations is developed. The solution is obtained numerically by a finite-difference scheme, and it is suitable for homogeneous porous media with unidirectional interstitial velocity. Each dissolved component may undergo first-order decay and may sorb under local equilibrium conditions. It is also assumed that the dissolution process is mass transfer limited. The nonaqueous phase activity coefficients of the NAPL pool components are evaluated at each time step. The model behavior is illustrated through a synthetic example with a NAPL pool consisting of a mixture of TCA (1,1,2-trichloroethane) and TCE (trichloroethylene). The numerical solution presented in this work is in good agreement with a recently developed analytical solution for the special case of a single component NAPL pool. The results indicate the importance of accounting for the necessary changes in the organic phase activity which significantly affects the equilibrium aqueous solubility.Notation C liquid phase solute concentration (solute mass/liquid volume) (M L^–3) - C _s single component aqueous saturation concentration (solubility) (M L^–3) - C ^w equilibrium aqueous solubility (M L^–3) - D molecular diffusion coefficient (L² t ^–1) - D _e effective molecular diffusion coefficient (L² t ^–1) - D _x longitudinal hydrodynamic dispersion coefficient (L² t ^–1) - D _y lateral hydrodynamic dispersion coefficient (L² t ^–1) - D _z hydrodynamic dispersion coefficient in the vertical direction (L² t ^–1) - I() integer mode arithmetic operator - k local mass transfer coefficient (Lt ^–1) - k ^* average mass transfer coefficient (Lt ^–1) - L length - l _x ,l _y pool dimensions inx andy directions (L) - ll _x ,l _y x andy Cartesian coordinates of the pool origin (L) - M number of moles remaining in a pool (moles) - M initial number of moles (moles) - n finite-difference scheme time level - R retardation factor (dimensionless) - t time (t) - U _x average interstitial velocity (Lt ^–1) - x, y, z spatial Cartesian coordinates (L) - X dimensionless mole fraction - dimensionless activity coefficient - _w viscosity of water (=0.8904 cp at 25°C) - decay coefficient (t ^–1) - ^* tortuosity ( 1) - i,j, k finite-difference scheme grid indicators - p component number indicator - P total number of components - s pure single component - o nonaqueous phase - w aqueous phase     

A multiphase flow and multispecies reactive transport model for DNAPL-involved Compound Specific Isotope Analysis

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.     

Slow gas expansion in saturated natural porous media by gas injection and partitioning with non-aqueous phase liquids

The partitioning of volatile non-aqueous phase liquid (NAPL) compounds to a discontinuous gas phase can result in the expansion of that gas phase, and the resulting gas flow can significantly affect the mass transfer from NAPL source zones. This recently reported gas flow generated by the spontaneous expansion of a discontinuous gas phase has not been extensively characterized in the literature. This study measured the expansion rate of a single gas cluster in a 1.1 mm sand above a pool of trans-1,2-dichloroethene (tDCE) in small-scale flow cell experiments. To characterize the gas flow, gas injection experiments in three sizes of sand were conducted at very slow injection rates typical of gas flow rates produced by gas expansion due to NAPL partitioning. Gas cluster spontaneous expansion rates above a tDCE pool were found to be 0.34 ± 0.02 and 0.29 ± 0.01 mL/day in duplicate experiments, which is sufficiently slow to result in discontinuous gas flow in porous media with a grain size diameter greater than 0.02 mm. Measured capillary pressures during gas injection showed patterns consistent with discontinuous gas flow, and identified multiple fragmentation events and expansion by coalescence with trapped clusters. The combination of pressure data and light transmission images were used to identify fragmentation and obtain direct measurements of the critical cluster length (i.e. the length at which withdrawal of the gas phase from a pore space occurs) in quasi-two-dimensional porous media for the first time. The measured critical cluster lengths were 1.4–3.6, 3.2–6.0 and 2.8–6.5 cm in 1.1, 0.7 and 0.5 mm sands, respectively. These values agreed well with estimates of the critical cluster length made using previously reported equations, and parameters derived from the medium’s capillary pressure-saturation relationship.     

A novel experimental procedure to investigate the biodegradation of NAPL under unsaturated conditions

Soils need to be thoroughly investigated regarding their potential for the natural attenuation of non-aqueous phase liquids (NAPL). Laboratory investigations truly representative of degradation processes in field conditions are difficult to implement for porous media partially saturated with water, NAPL and air. We propose an innovative protocol to investigate degradation processes under steady-state vadose zone conditions. Experiments are carried out in glass columns filled with a sand and, as bacteria source, a soil from a diesel-fuel-polluted site. Water and NAPL (n-hexadecane diluted in heptamethylnonane (HMN)) are added to the porous medium in a two-step procedure using ceramic membranes placed at the bottom of the column. This procedure results, for appropriate experimental conditions, in a uniform distribution of the two fluids (water and NAPL) throughout the column. In a biodegradation experiment non-biodegradable HMN is used to provide NAPL mass, while keeping biodegradable n-hexadecane small enough to monitor its rapid degradation. Biodegradation is followed as a function of time by measuring oxygen consumption, using a respirometer. Degradative activity is controlled by diffusive transfers in the porous network, of oxygen from the gas phase to the water phase and of n-hexadecane from the NAPL phase to the water phase.     

环境地球物理学及其现状与进展

简述了环境地球物理学的发展过程及其研究内容和方法。归纳总结出目前环境地球物理学研究的5大趋势:①勘探技术的进步和数据处理方法的不断提高;②非水相流体的研究是环境地球物理学研究的重要领域;③地质灾害预报和环境污染监测、治理仍是环境地球物理学研究的主要内容;④特殊环境下的地球物理技术有待发展;⑤生态环境研究是环境地球物理学今后研究的新热点。在分析我国环境地球物理学研究现状的基础上,提出了发展我国环境地球物理学的若干对策与建议。