Abstract:The amount, location, and form of NAPL in contaminated vadose zones are controlled by the spatial distribution of water saturation and soil permeability, the NAPL spill scenario, water infiltration events, and vapor transport. To evaluate the effects of these processes, we used the three-phase flow simulator STOMP, which includes a new permeability-liquid saturation-capillary pressure (k-S-P) constitutive model. This new constitutive model considers three NAPL forms: free, residual, and trapped. A 2-D vertical cross-section with five stratigraphic layers was assumed, and simulations were performed for seven cases. The conceptual model of the soil heterogeneity was based upon the stratigraphy at the Hanford carbon tetrachloride (CT) spill site. Some cases considered co-disposal of NAPL with large volumes of wastewater, as also occurred at the Hanford CT site. In these cases, the form and location of NAPL were most strongly influenced by high water discharge rates and NAPL evaporation to the atmosphere. In order to investigate the impact of heterogeneity, the hydraulic conductivity within the lower permeability layer was modeled as a realization of a random field having three different classes. For six extreme cases of 100 realizations, the CT mass that reached the water table varied by a factor of two, and was primarily controlled by the degree of lateral connectivity of the low conductivity class within the lowest permeability layer. The grid size at the top boundary had a dramatic impact on NAPL diffusive flux just after the spill event when the NAPL was present near the ground surface. NAPL evaporation with a fine grid spacing at the top boundary decreased CT mass that reached the water table by 74%, compared to the case with a coarse grid spacing, while barometric pumping had a marginal effect for the case of a continuous NAPL spill scenario considered in this work. For low water infiltration rate scenarios, the distribution of water content prior to a NAPL spill event decreased CT mass that reached the water table by 98% and had a significant impact on the formation of trapped NAPL. For all cases simulated, use of the new constitutive model that allows the formation of residual NAPL increased the amount of NAPL retained in the vadose zone. Density-driven advective gas flow from the ground surface controlled vapor migration in strongly anisotropic layers, causing NAPL mass flux to the lower layer to be reduced. These simulations indicate that consideration of the formation of residual and trapped NAPLs and dynamic boundary conditions (e.g., areas, rates, and periods of different NAPL and water discharge and fluctuations of atmospheric pressure) in the context of full three-phase flow are needed, especially for NAPL spill events at the ground surface. In addition, NAPL evaporation, density-driven gas advection, and NAPL vertical movement enhanced by water flow must be considered in order to predict NAPL distribution and migration in the vadose zone.

Analysis of the LNAPL Migration Process in the Vadose Zone under Two Different Media Conditions

Study on Influence of Geological Heterogeneity on Migration of LNAPL in Contaminated Site Through Numerical Analysis

Migration behaviour of LNAPL in fractures filled with porous media: Laboratory experiments and numerical simulations

Retention effect and mode of capillary zone on the migration process of LNAPL pollutants based on experimental exploration

Experimental study on non-aqueous phase liquid multiphase flow characteristics and controlling factors in heterogeneous porous media

Experiment study on multiphase flow of LNAPL polluted soil

The Study of Variations in Groundwater Occurrence Influencing the Migration Characteristics of LNAPL

Numerical Investigation on Light Non–aqueous Phase Liquid Flow in the Vadose Zone Considering Porosity Effects on Soil Hydraulic Properties

Numerical Modeling of Leakage, Transport and Remediation of Mixed DNAPL and LNAPL in the Unsaturated Clayey Soil Underlain by Saturated Sandy Soil

Numerical Analysis and Centrifugal Modeling of LNAPLs Transport in Subsurface System

Assessing the impact of water infiltration on LNAPL mobilization in sand column using simplified image analysis method

Study on the Behavior of Transportation of LNAPL in Soil and Groundwater

Experimental study on migration characteristics of LNAPL in the aquitard under pumping conditions

Evaluation of factors causing lateral migration of light non-aqueous phase liquids (LNAPLs) in onshore oil spill accidents

Influence of Flow Velocity and Spatial Heterogeneity on DNAPL Migration in Porous Media: Insights from Laboratory Experiments and Numerical Modelling

Evaluation for the Behavior of LNAPL Released in Vadose Zone As Source of Groundwater Contamination

Construction, application and validation of a new algorithm for determining light nonaqueous-phase liquid fluxes in unsaturated zones

Laboratory testing and numerical simulation of multiphase flow in unsaturated soils

Pattern of LNAPL's migration and distribution in 2-D auriferous media

Effect of soil moisture dynamics on dense nonaqueous phase liquid (DNAPL) spill zone architecture in heterogeneous porous media

Effect of Soil Fabric on Transport of a LNAPL Through Unsaturated Fine-Grained Soils: A Centrifugal Model Study