The influence of wettability on NAPL dissolution fingering

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.     

Impact of NAPL architecture on interphase mass transfer: A pore network study

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.     

Effect of NAPL film stability on the dissolution of residual wetting NAPL in porous media: A pore-scale modeling study

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.     

Evaluation of an upscaled model for DNAPL dissolution kinetics in heterogeneous aquifers

Estimates of contaminant fluxes from DNAPL sources as a function of time and DNAPL mass reduction are important to assess the long-term sustainability and costs of monitored natural attenuation and to determine the benefits of partial source removal. We investigate the accuracy of the upscaled mass transfer function (MTF) proposed by Parker and Park [Parker JC, Park E. Modeling field-scale dense nonaqueous phase liquid dissolution kinetics in heterogeneous aquifers. WRR 2004;40:W05109] to describe field-scale dissolved phase fluxes from DNAPL sources for a range of scenarios generated using high-resolution 3-D numerical simulations of DNAPL infiltration and long-term dissolved phase transport. The results indicate the upscaled MTF is capable of accurately describing field-scale DNAPL dissolution rates as a function of time. For finger-dominated source regions, an empirical mass depletion exponent in the MTF takes on values greater than one which results in predicted mass flux rates that decrease continuously with diminishing DNAPL mass over time. Lens-dominated regions exhibit depletion exponents less than one, which results in more step-function like mass flux versus time behavior. Mass fluxes from DNAPL sources exhibiting both lens- and finger-dominated subregions were less accurately described by the simple MTF, but were well described by a dual-continuum model of the same form for each subregion. The practicality of calibrating a dual-continuum model will likely depend on the feasibility of obtaining spatially resolved field measurements of contaminant fluxes or concentrations associated with the subregions using multilevel sampling or some other means.     

文安地震前后的转换函数参数变化分析

为了提取地震前后的磁场异常变化,本文用转换函数方法分析了文安地震前后,距离震中分别为52km和298km的静海地磁台和泰安地磁台不同周期的转换函数参数随时间的变化过程,发现静海地磁台转换函数参数的数据异常点随着时间的向后推移从短周期向比较长的周期迁移,而泰安地磁台的转换函数参数没有出现类似现象。     

Implications of rate-limited mass transfer for aquifer storage and recovery

Pressure to decrease reliance on surface water storage has led to increased interest in aquifer storage and recovery (ASR) systems. Recovery efficiency, which is the ratio of the volume of recovered water that meets a predefined standard to total volume of injected fluid, is a common criterion of ASR viability. Recovery efficiency can be degraded by a number of physical and geochemical processes, including rate-limited mass transfer (RLMT), which describes the exchange of solutes between mobile and immobile pore fluids. RLMT may control transport behavior that cannot be explained by advection and dispersion. We present data from a pilot-scale ASR study in Charleston, South Carolina, and develop a three-dimensional finite-difference model to evaluate the impact of RLMT processes on ASR efficiency. The modeling shows that RLMT can explain a rebound in salinity during fresh water storage in a brackish aquifer. Multicycle model results show low efficiencies over one to three ASR cycles due to RLMT degrading water quality during storage; efficiencies can evolve and improve markedly, however, over multiple cycles, even exceeding efficiencies generated by advection-dispersion only models. For an idealized ASR model where RLMT is active, our simulations show a discrete range of diffusive length scales over which the viability of ASR schemes in brackish aquifers would be hindered.     

岩石压应力传递系数计算及影响机制

钻孔应变观测是地壳岩石应变相对测量的主要方法,对其附加应力场的研究,可以有效分析岩石压应力向探头的传递过程,可对应变仪有关力学参数的设计和观测数据的物理解释提供帮助。考虑到外加应力在岩石、水泥和仪器钢筒等3层介质中传递的实际情况,根据弹性力学中的厚壁圆筒方程,推导轴对称应力作用下岩石压应力的传递系数公式,分析影响机制,并给出有关变化规律曲线和定量指标。     

Geochemical modeling to aid experimental design for multiple isotope tracer studies of coupled dissolution and precipitation reaction kinetics

It is a challenge to make thorough but efficient experimental designs for the coupled mineral dissolution and precipitation studies in a multi-mineral system, because it is difficult to speculate the best experimental duration, optimal sampling schedule, effects of different experimental conditions, and how to maximize the experimental outputs prior to the actual experiments. Geochemical modeling is an efficient and effective tool to assist the experimental design by virtually running all scenar...     

Prediction of railway ground vibrations: Accuracy of a coupled lumped mass model for representing the track/soil interaction

Recent advances in railway-induced ground vibrations showed that the track/soil interaction plays an important role in the low frequency range. This paper contributes to the numerical analysis of train/track/foundation dynamics by presenting the accuracy of a coupled lumped mass (CLM) model devoted to the railway foundations and to the track/soil coupling. Following a summary of the background and the advantages of the CLM model, the coupling strategy is quantified through two application cases. Firstly, the dynamic track deflection is calculated for different railway lines considering various degrees of complexities of foundations. Then, the foundation responses are compared depending on whether detailed coupling is introduced or not. The benefit of the proposed model is emphasized by presenting free-field ground vibration responses generated by a tram and a high-speed train, obtained by a revisited two-step prediction model developed by the authors.     

Instrumentation for measuring and recording mass and energy transfer from volcanoes to atmosphere

A new approach to the understanding of volcanic activity is attempted. It consists in establishing tentative correlations between all measurable physical, chemical and geophysical eruptive parameters. The first stage of this program is to devise and design gauges and recorders of different types to be used on different types of existing eruptive vents, whatever their shape, size, violence and accessibility. The methods adopted and the first results obtained on the physical parameters of eruptive gases are described here.     

Time scale of hydrothermal water-rock reactions in Yellowstone National Park based on radium isotopes and radon

We have measured ²²⁴Ra (3.4 d), ²²⁸Ra (5.7 yr), and ²²⁶Ra (1620 yr) and chloride in hot spring waters from the Norris-Mammoth Corridor, Yellowstone National Park. Two characteristic cold-water components mix with the primary hydrothermal water: one for the travertine-depositing waters related to the Mammoth Hot Springs and the other for the sinter-depositing Norris Geyser Basin springs. The Mammoth Hot Springs water is a mixture of the primary hydrothermal fluid with meteoric waters flowing through the Madison Limestone, as shown by the systematic decrease of the (²²⁸Ra/²²⁶Ra) activity ratio proceeding northward. The Norris Geyser Basin springs are mixtures of primary hydrothermal water with different amounts of cold meteoric water with no modification of the primary hydrothermal (²²⁸Ra/²²⁶Ra) activity ratio. Using a solution and recoil model for radium isotope supply to the primary hydrothermal water, a mean water-rock reaction time prior to expansion at 350°C and supply to the surface is 540 years assuming that 250 g of water are involved in the release of the radium from one gram of rock. The maximum reaction time allowed by our model is 1150 years.     

Simulation of flood hydrographs based on lumped and semi-distributed models for two tropical catchments in Kenya

Abstract

The design and operation of flood management systems require computation of flood hydrographs for both design floods and flood forecasting purposes, since observed data are usually inadequate for these tasks. This is particularly relevant for most developing countries, i.e. mainly for tropical catchments. One possible way of obtaining information about flood hydrographs is through the use of rainfall—runoff models. Two such models, namely the Bochum model and the Nash Cascade—Diskin Infiltration model, which are semi-distributed and lumped models, respectively, were used in the present study. These models were applied to two catchments in Kenya with drainage areas of 6.71 km² and 26.03 km². A set of 13 selected rainfall—runoff events was used to calibrate and validate the models. The physical parameters required by the models were derived from catchment characteristics using GIS and remote sensing data while the conceptual parameters were obtained by optimization. The flood hydrographs simulated using the parameters so derived indicated that it is possible to use the two models in this tropical environment.     

Effects of rate-limited mass transfer on water sampling with partially penetrating wells

Nonequilibrium concentration type curves are numerically developed and sensitivity analyses are performed to examine the relationships between effluent concentrations in partially penetrating monitoring/extraction wells, the vertical plume shape, and the mass transfer characteristics of the aquifer. The governing two-dimensional, axisymmetric nonequilibrium solute transport equation is solved in three stages using an operator-splitting approach. In the first two stages, the advection and dispersion terms are solved with the Eulerian-Lagrangian method, based on the backward method of characteristics for advection and the standard implicit Galerkin finite element method for dispersion. In the third step, the first-order, immobile-mobile domain mass transfer term is computed analytically for both two-site and lognormally distributed, multirate models. Effluent concentration variations with time and contour plots of the pore water concentration distribution in the aquifer are compared for a wide range of field- and laboratory-measured mass transfer rates, various plume shapes, and relevant physical/chemical parameter values, including pumping rate, vertical anisotropy ratio, retardation factor, and porosity. The simulation results show that rate-limited mass transfer can have a significant impact on sample and aquifer pore water concentrations during three-dimensional transport to a partially penetrating well. An alternative dimensionless form of the nonequilibrium solute transport equation is derived to illustrate the key parameter groupings that quantify rate-limited sorption effects and show the relative importance of individual parameters. A hypothetical field application example demonstrates the fitting of dimensional type curves to discrete-interval sampling data in order to evaluate the mass transfer characteristics of an aquifer and shows how type curve superposition can be used to model complex plume shapes.     

Influence of an endothermic phase transition on mass transfer between the upper and the lower mantle

Geochemical data indicate that two major reservoirs 1–2 Ga in age are present in the mantle. The upper mantle, feeding mid-ocean ridges, is depleted in chemical elements carried away into the continental crust. The lower mantle, feeding hotspot plumes, is close in composition to primordial matter. The 660-km depth of an endothermic phase transition in olivine has been considered over the last two decades as a possible boundary between the reservoirs. In this period, many models of mantle convection were constructed that used values of the phase transition parameters, which led to temporal (up to 1 Gyr long) convection layerings and periodic avalanche-induced mantle intermixing events. However, laboratory measurements with new improved instrumentation give other values of the phase transition parameters that require a revision of the majority of the existence of large-scale avalanches in the Earth’s history becomes disputable. The paper is devoted to comprehensive study of the phase transition effect on the structure of mantle flows with different values of phase transition parameters and Rayleigh numbers; in particular, the mass transfer through the phase boundary is calculated for different regimes of steady-state convection.     

Numerical analysis of simultaneous heat and mass transfer during infiltration into frozen ground

This paper concerns heat transfer and water infiltration into a homogeneous, unsaturated, frozen soil. A numerical analysis of the processes, which is based on the local volume averaging formulation of transport phenomena in porous media, is described. Simultaneous heat and mass transfers with phase changes are considered. The results of the simulations provide an insight into the mechanics of the infiltration process. They suggest that infiltration involves primarily two flow regimes, a transient regime and a quasi-steady-state regime. Once the quasi-steady-state regime is reached, the energy needed to increase soil temperature at depth is supplied by latent heat released from freezing of water in the upper layers of soil. The effects of surface saturation, initial soil saturation and initial soil temperature on infiltration are examined.     

Influence of mass transfer resistance on detection of nonaqueous phase liquids with partitioning tracer tests

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.     

Seismic response of bridge piers on pile groups for different soil damping models and lumped parameter representations of the foundation

This paper presents a wide parametric study aimed at elucidating the influence, on the computed seismic response of bridge piers, of two related aspects of the model: (1) the adoption of the classical hysteretic or the causal Biot's damping models for the soil and (2) the use of two different lumped parameter models of different complexity and accuracy to approximate the impedances of the pile foundation. A total of 2072 cases, including different superstructures, pile foundations, soil deposits, and seismic input signals, are studied. The results are presented so that the influence of the different parameters involved in the analysis can be assessed. From an engineering point of view, both lumped parameter models provide, in general, sufficiently low errors. The choice of the most adequate model for each case will depend not only on the configuration of the structure and the soil-foundation system but also on the assumed soil damping model, whose influence on the computed seismic responses is relevant in many cases. The nonphysical behaviour provided by the classical hysteretic damping model for the soil at zero frequency generates issues in the process of fitting the impedance functions. It is also found that larger deck displacements are predicted by Biot's model due to the higher damping at low frequencies provided by the classical hysteretic damping model.     

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.