Search

Light non-aqueous phase liquids (LNAPLs), which include various petroleum products, are a significant source of groundwater contamination globally. Once introduced into the subsurface, these contaminants tend to accumulate in the vadose zone, causing chronic soil and water pollution. The vadose zone often contains lens-shaped bodies with diverse properties that can significantly influence the migration and distribution of LNAPLs. Understanding the interaction between LNAPLs and these lens-shaped bodies is crucial for developing effective environmental management and remediation strategies. Prior research has primarily focused on LNAPL behavior in homogeneous media, with less emphasis on the impact of heterogeneous conditions introduced by lens-shaped bodies. To investigate the impact of lens-shaped structures on the migration of LNAPLs and to assess the specific effects of different types of lens-shaped structures on the distribution characteristics of LNAPL migration, this study simulates the LNAPL leakage process using an indoor two-dimensional sandbox. Three distinct test groups were conducted: one with no lens-shaped aquifer, one with a low-permeability lens, and one with a high-permeability lens. This study employs a combination of oil front curve mapping and high-density resistivity imaging techniques to systematically evaluate how the presence of lens-shaped structures affects the migration behavior, distribution patterns, and corresponding resistivity anomalies of LNAPLs. The results indicate that the migration rate and distribution characteristics of LNAPLs are influenced by the presence of a lens in the gas band of the envelope. The maximum vertical migration distances of the LNAPL are as follows: high-permeability lens (45 cm), no lens-shaped aquifer (40 cm), and low-permeability lens (35 cm). Horizontally, the maximum migration distances of the LNAPL to the upper part of the lens body decreases in the order of low-permeability lens, high-permeability lens, and no lens-shaped aquifer. The low-permeability lens impedes the vertical migration of the LNAPL, significantly affecting its migration path. It creates a flow around effect, hindering the downward migration of the LNAPL. In contrast, the high-permeability lens has a weaker retention effect and creates preferential flow paths, promoting the downward migration of the LNAPL. Under conditions with no lens-shaped aquifer and a high-permeability lens, the region of positive resistivity change rate is symmetrical around the axis where the injection point is located. Future research should explore the impact of various LNAPL types, lens geometries, and water table fluctuations on migration patterns. Incorporating numerical simulations could provide deeper insights into the mechanisms controlling LNAPL migration in heterogeneous subsurface environments.

Contamination from light nonaqueous phase fluids (LNAPLs) and their derivatives during mining, production, and transportation has become a concern. Scholars have extensively studied LNAPL contamination, but the role of water content variation on its migration process in the unsaturated zone has not been sufficiently researched. The specific issue addressed in this study is the impact of water content on the migration of light nonaqueous phase liquids (LNAPLs) in sandy soils, a critical yet under-researched aspect of subsurface contamination. To tackle this, we employed indoor simulated vertical, one-dimensional, multiphase flow soil column experiments, utilizing the orthogonal experimental method to systematically evaluate the effects of varying water contents on the occurrence state and migration rate of LNAPLs. The experimental results indicate the following: (1) The migration rate of LNAPL exhibits an L-shaped trend during subsurface imbibition and a nonlinear relationship with migration time. The migration rate and migration time of surface infiltration have a linear growth relationship. (2) The residual rate of LNAPL is negatively correlated with water content and positively correlated with oil content in the homogeneous non-saturated state. With the increase in the amount of leaked oil, 40% of the leaked LNAPL is sorbed within the soil. (3) When the water content of the test medium is below 14%, and the oil content is below 11%, LNAPL appears in the unsaturated zone in a solid phase. As the water content increases, the adsorption rate of the oil phase gradually decreases and eventually reaches the oil saturation point. (4) When the water content of the medium exceeds 8%, over time, LNAPL will be subject to oil–water interfacial tension, and the rate of LNAPL movement first decreases and then increases, displaying nonlinear growth. The innovation of this work lies in the comprehensive analysis of LNAPL migration under controlled laboratory conditions, providing results that enhance the understanding of LNAPL behavior in sandy soils. These quantitative insights are crucial for developing targeted remediation strategies for LNAPL-induced pollution in the unsaturated zone.

The spatial variability of soil properties has always been a significant research field in geoscience. The types of soil properties cover a wide range, but most studies have focused on the spatial variability of soil physicochemical properties over the past decades. Studies on soil hydraulic characteristics are limited, and most of them are limited to the farmland scale. However, the spatial variability of regional soil properties (soil texture and hydraulic properties) is valuable for the study of sedimentation processes and soil water transport. Therefore, here, the spatial variation of six soil properties (sand, silt, clay content, bulk density, saturated water content and saturated hydraulic conductivity) in the typical alluvial plain area of the lower Yellow River is quantitatively studied, by using classical statistics, geostatistics and single fractal and multifractal methods. This study mainly quantitatively analysed the spatial variability of different soil properties and compared four research methods. Although the coefficient of variation, nugget coefficient, single fractal dimension and multifractal spectral width can reflect spatial variability, diverse conclusions are drawn (on variability) if different methods are used, and the different soil properties show large disparities. These four methods show a different variation order of soil properties, but there are some common conclusions based on analysis and judgment. In general, the silt content in the study area is stable, mainly originating from loess transported by Yellow River erosion, which is also reflected in the Kriging interpolation maps under the geostatistical models. The variation in bulk density and saturated water content is weak, and the spatial variability of sand and clay content is moderate. In addition, the saturated hydraulic conductivity fluctuates violently. This may be related to the differences in local topography, human activity and the content of sand and clay, each of which significantly affects the saturated hydraulic conductivity. Classical statistics has a limitation because it fails to corelate with spatial location. Due to the small sample capacity and calculation error of lag distance, the accuracy of geostatistics and single fractal dimensions needs to be improved. Multifractal spectral analysis does not need to consider the normality of data and can quantitatively represent local characteristics; therefore, its results have high reliability.