Infiltration and redistribution of LNAPL into unsaturated layered porous media

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

Contamination of groundwater by harmful substances poses significant risks to both drinking water sources and aquatic ecosystems, making it a critical environmental concern. Most on-land spill events release organic molecules known as light non-aqueous phase liquids (LNAPLs), which then seep into the ground. Due to their low density and organic composition, they tend to float as they reach the water table. LNAPLs encompass a wide range of non-aqueous phase liquids, including various petroleum products, and can, over time, develop carcinogenic chemicals in water. However, due to frequent changes in hydraulic head, the confinement may fail to contain them, causing them to extend outward. When it contaminates water wells, people cannot reliably consume the water. The removal of dangerous contaminants from groundwater aquifers is made more challenging by LNAPLs. It is imperative to analyze the mechanisms governing LNAPL migration. As a response to this need and the associated dispersion of contaminants into adjacent aquifers, we have conducted a comprehensive qualitative literature review encompassing the years 2000-2022. Groundwater variability, soil structure, and precipitation have been identified as the three primary influential factors, ranked in the following order of significance. The rate of migration is shown to rise dramatically in response to changes in groundwater levels. Different saturation zones and confinement have a major effect on the lateral migration velocity. When the various saturation zones reach a balance, LNAPLs will stop moving. Although higher confinement slows the rate of lateral migration, it speeds up vertical migration. Beyond this, the lateral or vertical movement is also influenced by differences in the permeability of soil strata. Reduced mobility and tighter containment are the outcomes of migrating through fine-grained, low-porosity sand. The gaseous and liquid phases of LNAPLs move more quickly through coarse-grained soils. Due to the complexities and uncertainties associated with LNAPL behavior, accurately foreseeing the future spread of LNAPLs can be challenging. Although studies have utilized modeling techniques to simulate and predict LNAPL migration, the inherent complexities and uncertainties in the subsurface environment make it difficult to precisely predict the extent of LNAPL spread in the future. The granular soil structure considerably affects the porosity and pore pressure.

Spatial infiltration and redistribution of light crude oil in heterogeneous water-bearing soil layers under different hydrogeological processes

Two physical models were used to simulate the infiltration and redistribution process of light crude oil after leakage in a heterogeneous soil layer following water level variation and rainfall. Migration fronts and redistribution characteristics of oil during gravity seepage, water level variation, and rainfall were obtained using charge-coupled device (CCD) camera shooting and cyan-magenta-yellow‒black (CMYK)-based gray analysis, which were employed efficiently and at a low cost. Then, the influencing factors and migration mechanisms were examined. Finally, the soil water and oil contents were measured to verify the simulation results. The results are as follows: (1) the geologic lens and fine-coarse interface can intercept oil, resulting in a local highly contaminated area. (2) The crude oil infiltration path and velocity varied greatly with the different soil types and initial water contents. Within a certain range, the higher the initial water content is, the higher the lateral and vertical infiltration speeds. (3) The oil redistribution process was dominated by vertical infiltration under the condition of water level variation or rainfall, but oil-water displacement and the capillary pressure caused some oil to move horizontally near the geologic lens and fine-coarse interface. (4) Water level variation resulted in a synchronous rise or fall of the oil accumulation area, but rainfall caused it to move up. (5) Water level variation and rainfall imposed a certain influence on the periodic accumulation and release of crude oil in heterogeneous soil, especially in the presence of geologic lenses and lithologic interfaces.

Migration and redistribution of LNAPL in inclined stratified soil media

Leakage of light non-aqueous phase liquid (LNAPL) into soil can cause serious environmental issues. In this study, a two-dimensional device with adjustable dip angles was designed to investigate the migration and redistribution of LNAPL in natural inclined stratified soil media by the light transmission visualization (LTV) technology. The captured experimental images were processed to obtain the diesel distribution based on gray value which could represent the LNAPL saturation distribution. LNAPL may not be able to penetrate through the fine-coarse interface due to the capillary barrier effects. In this case, the vertical and horizontal migration distances (V and H), contaminated area (S), as well as deviation angle (γ) of centroid increased with the dip angle. Increasing the leakage amount to more than 30 mL would result in LNAPL breakthrough at the 10°-inclined interface, leading to much larger V, H, S, and γ than those at 10 mL, while 20-mL LNAPL failed to break through. In the latter case, a lower leakage rate than 10 mL/min would cause larger H and γ but similar V or S in the long term. This study could enrich the understanding of LNAPL contamination in vadose zone, providing reference for the prediction and treatment in realistic inclined contaminated sites.

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

Two-dimensional sand tank experiments were designed to investigate the retention process of the capillary zone during the migration of light non-aqueous phase liquid (LNAPL) pollutants. The fine sand and silt media experiments simulated the LNAPL migration process given a shallow point source leakage scenario. The results indicate that the LNAPL was retained in the capillary zone. A retention factor, based on the ratio of the change in the vertical migration velocity of the LNAPL front with time, was proposed to quantitatively characterize the retention effect. The retention factor and time satisfied the function of σ=A×exp(-kt). And the retention factor increased gradually with time, indicating the enhanced retention effect of capillary zone on the vertical migration of LNAPL. The concentration change rate was then used to investigate the LANPL redistribution process, which had a relationship with time of ν_c=B×ln(t)+C. The capillary zone could be divided from top to bottom into a weak retention zone (B > 0, v_c < 0), a strong retention zone (B < 0), and a barrier zone (B > 0, v_c > 0). The retention effect of capillary zone on LNAPL migration gradually strengthened during the vertical migration of LNAPL. In addition, the coefficient B had a relationship with the environmental factors (i.e., EC, pH, and ORP) of B=a×sin(b×α×β×γ)^c and the fitting coefficient R² of the function was above 0.913 for both media, indicating a strong correlation between the LNAPL redistribution process and the key environmental factors.

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

An analytical algorithm coupling free-phase migration, precipitation, and natural attenuation through volatilization and biodegradation (FPVB) was developed to calculate the flux of light nonaqueous-phase liquid (LNAPL) leaking from unsaturated zone to groundwater. Sandbox and soil column experiments were performed to identify the LNAPL migration characteristics and states to provide data to establish and verify FPVB algorithm. For free-phase migration, the Kinematic Oily Pollutant Transport (KOPT) model was used to determine LNAPL movement velocity and leakage time. The correlations of water saturation, residual LNAPL saturation and the cumulative dissolution ratio of residual LNAPL were described using an empirical formula for the precipitation leaching process. Equations for diesel volatilization kinetics and first order degradation were used to describe the natural attenuation processes. Coupling the algorithms for the different stages gave the final FPVB algorithm. The FPVB algorithm was used to describe the pollution situation at a real site, and the results were consistent with the actual situation. The FPVB algorithm could be used to quickly assess the scale and degree of pollution with little information on the parameters for the actual LNAPL leakage event.

Effects of ecohydrological interfaces on migrations and transformations of pollutants: A critical review

With the rapid development of society, the soil and water environments in many countries are suffering from severe pollution. Pollutants in different phases will eventually gather into the soil and water environments, and a series of migrations and transformations will take place at ecohydrological interfaces with water flow. However, it is still not clear how ecohydrological interfaces affect the migration and the transformation of pollutants. Therefore, this paper summarizes the physical, ecological, and biogeochemical characteristics of ecohydrological interfaces on the basis of introducing the development history of ecohydrology and the concept of ecohydrological interfaces. The effects of ecohydrological interfaces on the migration and transformation of heavy metals, organic pollutants, and carbon‑nitrogen‑phosphorus (C-N-P) pollutants are emphasized. Lastly, the prospects of applying ecohydrological interfaces for the removal of pollutants from the soil and water environment are put forward, including strengthening the ability to monitor and simulate ecohydrological systems at micro and macro scales, enhancing interdisciplinary research, and identifying main influencing factors that can provide theoretical basis and technical support.

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

This study focused on the processes of free infiltration, precipitation displacement, and natural attenuation of the LNAPL under the condition of near-surface leakage. Sandbox experiments were performed to explore the migration characteristics of LNAPL in the vadose zone with two media structures and the influences of the soil interface on the migration of LNAPL. The results indicate that the vertical migration velocity of the LNAPL infiltration front in medium and coarse sand was 1 order of magnitude higher than that in fine sand and that the LNAPL accumulated at the coarse–fine interface, which acted as the capillary barrier. Displacement of precipitation for LNAPL had little relationship with rainfall intensity and was obviously affected by medium particle size, where coarse sand (40.78%) > medium sand (20.5%) > fine sand (10%). The natural attenuation rate of the LNAPL in the vadose zone was related to the water content of the media; the natural attenuation rate of fine sand was higher. This study simulated the process of the LNAPL leakage from the near surface into the layered heterogeneous stratum, improved the understanding of the migration of the LNAPL under different stratum conditions, and can provide support for the treatment of LNAPL leakage events in the actual site.

Optimizing Laboratory Investigations of Saline Intrusion by Incorporating Machine Learning Techniques

Deriving saltwater concentrations from the light intensity values of dyed saline solutions is a long-established image processing practice in laboratory scale investigations of saline intrusion. The current paper presents a novel methodology that employs the predictive ability of machine learning algorithms in order to determine saltwater concentration fields. The proposed approach consists of three distinct parts, image pre-processing, porous medium classification (glass bead structure recognition) and saltwater field generation (regression). It minimizes the need for aquifer-specific calibrations, significantly shortening the experimental procedure by up to 50% of the time required. A series of typical saline intrusion experiments were conducted in homogeneous and heterogeneous aquifers, consisting of glass beads of varying sizes, to recreate the necessary laboratory data. An innovative method of distinguishing and filtering out the common experimental error introduced by both backlighting and the optical irregularities of the glass bead medium was formulated. This enabled the acquisition of quality predictions by classical, easy-to-use machine learning techniques, such as feedforward Artificial Neural Networks, using a limited amount of training data, proving the applicability of the procedure. The new process was benchmarked against a traditional regression algorithm. A series of variables were utilized to quantify the variance between the results generated by the two procedures. No compromise was found to the quality of the derived concentration fields and it was established that the proposed image processing technique is robust when applied to homogeneous and heterogeneous domains alike, outperforming the classical approach in all test cases. Moreover, the method minimized the impact of experimental errors introduced by small movements of the camera and the presence air bubbles trapped in the porous medium.

Experimental study on the transport characteristics of buried pipeline leakage and the performance of groundwater remediation system

Enhanced understanding of light non-aqueous phase liquid (LNAPL) infiltration into sandy porous medium is significant to the effective design of remediation strategies. A system for buried pipeline leakage in 2-D sandbox was conducted to investigate the migration of diesel through a sandy porous medium, and the system could also be conducted to investigate groundwater remediation. Two groups of experiments were carried out. The first experiment consisted of diesel infiltration into a fine sand matrix. We could notice that diesel spilled in dry sand layer at a constant speed and the diesel front kept longitudinal movement due to the gravity before it arrived at the edge of the capillary zone. The diesel front broadened as a whole because of the capillary force jacking after it reached the capillary zone. Finally, the bulk of the diesel was contained on top of the capillary zone. To protect groundwater, the second experiment consisted of remediating soils and groundwater. The results indicated that the voltage of electrocoagulation apparatus had a great influence on the treatment effect, and the removal rate of diesel was found to be more than 90% with a constant voltage of 20 V. The efficiency of groundwater remediation was influenced by the flow velocity, and it took 11 h when the flow velocity was 2.089 L/min. To summarize, the research was conducive to the study on diesel pollution control and pollution prediction.

Study on diesel vertical migration characteristics and mechanism in water-bearing sand stratum using an automated resistivity monitoring system

Oil spills frequently occur on both land and sea. Petroleum in mobile phase will cause serious pollution in the sediment and can form a secondary pollution source. Therefore, it is very important to study the migration of petroleum in sediments ideally in a rapid and simplified approach. The release of diesel was simulated using fine beach sand to construct a model aquifer, and dynamic monitoring was carried out using an automated monitoring system including a resistivity probe originally developed by our research group. The mobile phase migration fronts were determined accurately using wavelet analysis method combined with resistivity curve method. Then, a relationship between resistivity and the joint oil-water content was established. The main conclusions were as follows. The seepage velocity of the diesel with high mobility at the initial stage of infiltration was faster, followed by a period when gravity seepage was dominant, and finally a redistribution period at the later stage, which was mainly an oil-water displacement process. The resistivity trends for diesel infiltration in different water-saturated soil layers varied with depth. The resistivity in the vadose zone fluctuated significantly, increasing initially and later decreasing. The resistivity change in the capillary zone was relatively small and constant in the initial stage; then, it increased and subsequently decreased. The resistivity in the saturated zone was basically unchanged with depth, and the value became slightly larger than the background value over time. Overall, for a large volume of mobile phase diesel leakage, the arrival migration fronts can be detected by wavelet analysis combined with resistivity curves. The thickness of the oil slick in the capillary zone can be estimated by resistivity changes. The relationships between resistivity and both the moisture content and oil-water joint saturation are in agreement with the linear models. The research results provide basic data and a new data processing method for monitoring of contaminated sites following major oil spills using the resistivity method.

Improved light-transmission method for the study of LNAPL migration and distribution rule

In this study, diesel was selected as a pollutant to study the migration and distribution rule of light non-aqueous phase liquid (LNAPL) in the simulated vadose zone. Saturation was regarded as a critical parameter to reflect the LNAPL migration and distribution rule. To get LNAPL saturation distribution figures, an image processing method of saturation was established to improve light-transmission technology, which can deal with digital camera images. Results showed that the vadose zone contains three areas from top to bottom, named dry media, transition zone and capillary zone. The system has two interfaces at which moisture content increased significantly. The significant increase in moisture content induced two apparent horizontal LNAPL diffusions in the two interfaces. Furthermore, the highly saturated LNAPL was mostly distributed near the wet interface, which lay between the dry media and the transition zone. Moreover, the downstream expansion of LNAPL in the capillary zone was promoted by groundwater flow, yet cutting off LNAPL supply could stop the downstream expansion after a period of time. The accuracy of this image processing method of saturation was verified by mass balance theory and reported a relative error of 4.38%.

Two-dimensional laboratory simulation of LNAPL infiltration and redistribution in the vadose zone

A quantitative two-dimensional laboratory experiment was conducted to investigate the immiscible flow of a light non-aqueous phase liquid (LNAPL) in the vadose zone. An image analysis technique was used to determine the two-dimensional saturation distribution of LNAPL, water and air during LNAPL infiltration and redistribution. Vertical water saturation variations were also continuously monitored with miniature resistivity probes. LNAPL and water pressures were measured using hydrophobic and hydrophilic tensiometers. This study is limited to homogeneous geological conditions, but the unique experimental methods developed will be used to examine more complex systems. The pressure measurements and the quantification of the saturation distribution of all the fluids in the entire flow domain under transient conditions provide quantitative data essential for testing the predictive capability of numerical models. The data are used to examine the adequacy of the constitutive pressure-saturation relations that are used in multiphase flow models. The results indicate that refinement of these commonly used hydraulic relations is needed for accurate model prediction. It is noted in particular that, in three-fluid phase systems, models should account for the existence of a residual NAPL saturation occurring after NAPL drainage. This is of notable importance because residual NAPL can act as a non negligible persistent source of contamination.