Spatial Moment Analysis Research Articles

Microbial community structure and function in a sandy soil column under dynamic hydrologic regimes

The capillary fringe (CF) exhibits steep chemical and redox gradients, representing a critical reactive interface in the subsurface. The spatial extent of the CF is temporally dynamic because of water table and soil moisture fluctuations. To simulate the spatio-temporal responses of the CF’s microbial community structure and activity, a 91 cm long sandy soil column experiment was conducted by imposing three different sequential hydrologic regimes: a stable water table separating the unsaturated and saturated zones (static regime), periodic pulse infiltration events with a constant water table position (infiltration regime), and water table fluctuations under three drainage-imbibition cycles of 6, 12, and 18 days (fluctuating regime). The investigation focused on the microbially-mediated turnover of nitrogen (N) and carbon (C). Geochemical profiles in the soil column were monitored continuously, while microbial DNA was extracted for 16S rDNA sequencing at the start, middle and end of the experiment. The observed microbial community compositions corresponded to five categories: (1) initial soil, (2) near-surface and (3) CF and saturated soil under the infiltration regime, (4) unsaturated and (5) saturated soil in the fluctuating regime. Based on co-occurrence network analysis, microbial physiologies matched the predominant geochemical conditions, and relatively stable community structures prevailed throughout the experiment. NO3- and NH4+ were more available to the microbes in the fluctuating regime’s alternating oxic/anoxic zone, leading to a more flexible network. The prevalent microbial metabolic pathways predicted using FAPROTAX were mainly associated with N and C metabolisms and significantly correlated with the taxonomic compositions. However, the dominant microbial groups were shared among all the soil samples. The limited spatial differences in microbial community structure under the fluctuating regime were likely due to a more pronounced role of microbial activity than microbial diversity, and similar drying-wetting conditions between the soil column and the filed site. Internally consistent geochemical and microbial depth stratifications were identified by (1) the redox potential and nitrogen species distributions, (2) spatial moment analysis of NO3- and NH4+ concentrations, (3) microbial community composition and (4) function. The five imposed hydrologic regimes generated distinct and coupled geochemical and microbial signatures. These signatures resulted from variations in the soil microbial community’s activity, rather than diversity or abundance.

Spatial Moment Analysis of Convective Mixing in Three-Dimensional Porous Media Using X-ray CT Images

Spatial Moment Analysis of Convective Mixing in Three-Dimensional Porous Media Using X-ray CT Images

Numerical analysis of spatial moment for colloid-facilitated contaminant transport through porous media

ABSTRACT Transport of contaminants in the presence of colloids is a three-phase problem in the groundwater porous media and colloids can also enhance the transport of contaminants by reducing the retardation factor. In this study, the numerical solution of the equilibrium model is obtained using a fully implicit Crank-Nicolson finite-difference method and validated with ten experimental results. Sensitivity and spatial moment analysis have been performed to investigate the mobility and spreading behavior of contaminants in the presence of colloids. Results suggest that enhancement of contaminant transport in the presence of colloids depends on the colloid concentration and the physical and chemical interaction of contaminants with the stationary solid matrix and suspended colloids. As the colloidal concentration increases, the mobility of contaminants augments, resulting in the early arrival of contaminants at the column outlet. The affinity of contaminants toward the stationary solid matrix of porous media can lead to the reduced transport of contaminants even under the presence of a significant amount of colloid concentration. In contrast, enhanced transport of contaminants is observed if the affinity of contaminants is more toward the suspended colloids.

Spatial Moment Analysis of Convective Mixing in Three-Dimensional Porous Media Using X-ray CT Images

Spatial Moment Analysis of Convective Mixing in Three-Dimensional Porous Media Using X-ray CT Images

Statistical analysis for biogeochemical processes in a sandy column with dynamic hydrologic regimes using spectral induced polarization (SIP) and self-potential (SP)

SUMMARYThe capillary fringe (CF) is characterized by transient and steep redox gradients and is thought to be a hot spot for biogeochemical processes. Understanding chemical fate and transport in the CF is significant, however, biogeochemical dynamics at the CF are poorly understood because of the difficulty to measure representatively with high spatio-temporal resolution at depths under dynamic hydrologic regimes. Hydrogeophysics is a developing field that uses minimally intrusive and quick response methods to monitor hydrological properties. Two geoelectrical methods [spectral induced polarization (SIP) and self-potential (SP)], which are sensitive to the solid–liquid interfaces (SIP) and biogeochemical processes (SP) can address the above difficulty. The challenge lies on linking the geoelectrical responses with biogeochemical processes, where many different processes contribute to the signals. We conducted a soil column experiment under five hydrologic regimes focusing on nitrogen transformations with SIP and SP measurements: (1) a static regime with a stable water level; (2) an infiltration regime with periodic pulse infiltration events with a constant water level and (3) fluctuating regimes with water level fluctuations under three drying-wetting frequencies (6/12/18-day-cycle). This is the first large lab-scale work in a well-controlled and highly instrumented soil column. The dynamic hydrologic conditions stimulated complex biogeochemical processes at the CF, and therefore the SIP and SP signals result from many physical and biogeochemical processes. Therefore, we relied on statistical analysis in this study for a novel interpretation. Spearman correlation analysis supported water content played the most important role in real conductivity (σ′) dynamics in the vadose zone, whereas fluid conductivity dominated σ′ in the saturated zone. Both correlation analysis and spatial moment analysis implicated that water content was the driving factor for both σ′ and imaginary conductivity (σ″). A multiple linear regression model indicated the gradient of redox potential, the gradient of soil matric potential and water content were the three main influencing factors for the SP signals. We proposed that the water level fluctuation can efficiently facilitate microbial electron transfer through ions transport between the different redox zones, and aggregate redox processes to create SP signal gradients. Depth zonation analysis, using six environmental indexes (Eh and nitrogen species; water content; real conductivity; imaginary conductivity; SP signal; microbial community composition), suggested that water content induced by soil hydrology was the most dominant factor, captured by all the indexes. In turn, it led to indirect inference on biogeochemical processes and resultant geoelectrical signals. Applying geoelectrical methods to such biogeochemical processes will not only lead to a better understanding of the mechanistic meanings of the geoelectrical signals, but also build relationships between geoelectrical signals and biogeochemical parameters to facilitate a novel way to monitor biogeochemical processes.

Truncation effect on estimation of transport parameters for slug-injection tracer tests

For slug-injection tracer tests, tracer concentrations below the detection limit of the measurement instrument can cause truncation of the observed data. This study investigated the truncation effect on the estimation error of parameters based on analytical solutions and the results of laboratory-scale experiments. Spatial moment analysis was performed to estimate the measured total mass and transport parameters, including the pore velocity and the longitudinal and transverse dispersivities. Increasing the travel distance and detection limit caused the measured mass and dispersivities to be underestimated regardless of the dimensionality because hydrodynamic dispersion occurs with increasing travel distance, which smoothens the concentration front. The one- and two-dimensional cases showed that the truncation effect on the measured mass and longitudinal dispersivity depended on dimensionality. In contrast, the pore velocity showed no such dependence; the center of mass did not change as the unmeasured portion due to truncation was increased because the plume, which exhibited a Gaussian distribution, was truncated symmetrically. In the experiment, the measured mass and dispersivities likewise depended on the travel distance and detection limit, but there were large differences in the detection limit at which the dimensionless parameter reached a value of zero between the experimental results and analytical solution. This is because the initial plume in the experiment was of a finite size. Thus, experimental design factors such as the scale, device, and dimensionality should be considered to minimize the estimation error of transport parameters, excluding the pore velocity.

Confronting heterogeneous sorption and hydrodynamic dispersion on solute transport in a fracture-skin-matrix system using spatial moment analysis.

Understanding the intricacies of inter-dependency of fluid flow and solute transfer at the scale of a single fracture is limited by various simplifying assumptions employed for computational purposes. In the present study, the fracture-rock matrix interface is assumed to be consisting of a skin layer with sufficient mass transfer properties where non-linear adsorption is considered to be the limiting reaction among the various interfaces. A numerical model has been developed using implicit finite difference method with varying grids at the fracture-skin interface to capture the mass transfer during solute transport in the presence of non-linear Sips adsorption. The model was used to identify the critically influencing parameters on the temporal profiles of fluid velocity, macro-dispersion coefficient and dispersivity using the method of spatial moments. The results indicate that the presence of the skin has enhanced the mixing phenomenon as well as the sorptive mass transfer rates. Summary of the sensitivity analysis provides the critical factors to be considered while employing such comprehensive models for elucidating details at a small scale. The presence of fracture-skin evades the sensitive role played by (a) the fracture adsorption coefficient (with a reduced rate constant for adsorption and an increased rate constant for desorption) at early times, while resulting in an enhanced mixing characteristics at later times, and (b) maximum sorption capacity of fracture as a function of solute velocity; and in turn, it provides an improved control over the transportation of solutes through the fracture.

Exploring the Impacts of Source Condition Uncertainties on Far‐Field Brine Leakage Plume Predictions in Geologic Storage of CO2: Integrating Intermediate‐Scale Laboratory Testing With Numerical Modeling

AbstractNatural fissures/faults or pressure‐induced fractures in the caprock confining injected CO2 have been identified as a potential leakage pathways of far‐field native brine contaminating underground sources of drinking water. Developing models to simulate brine propagation through the overlaying formations and aquifers is essential to conduct reliable pre‐ and post‐risk assessments for site selection and operation, respectively. One of the primary challenges of performing such simulations is lack of adequate information about source conditions, such as hydro‐structural properties of caprock fracture/fault zone and the permeability field of the storage formation. This research investigates the impact of source condition uncertainties on the accuracy of leaking brine plume predictions. Prediction models should be able to simulate brine leakage and transport in complex multilayered geologic systems with interacting regional natural and leakage flows. As field datasets are not readily available for model testing and validation, three comprehensive intermediate‐scale laboratory experiments were used to generate high‐resolution spatiotemporal data on brine plume development under different leakage scenarios. Experimental data were used to validate a flow and transport model developed using existing code FEFLOW to simulate brine plume under varying source conditions. Spatial moment analysis was conducted to evaluate how uncertainty in source conditions impacts brine migration predictions. Results showed that inaccurately prescribing the permeability field of storage formation and caprock fractures in models can cause errors in leakage pathway and spread predictions up to ∼19% and ∼100%, respectively. These findings will help in selecting and characterizing storage sites by factoring in potential risks to shallow groundwater resources.

Analysis of soil wetting patterns in subsurface drip irrigation systems – Indoor alfalfa experiments

When using subsurface drip irrigation (SDI) to grow field crops, concerns may exist for water availability to plants, because emitters and plants cannot be lined up perfectly. The spatial distribution of water around SDI emitters is largely determined by soil hydraulic properties and characteristics of root water uptake (RWU). Experiments were performed in an indoor facility to quantify the effects of alfalfa RWU on characteristics of wetting patterns in three soil textures. In all soils, the majority of RWU took place in a layer below the emitters, where soil moisture was most dynamic. Unique wetting patterns developed in the three soils in response to irrigation events. However, when irrigation applications were aimed to replace potential evapotranspiration, these unique patterns disappeared within 24 hr. When irrigation was supplied in excess for a 21-d period, an effective wetted zone developed over time. Both the extent of and moisture distribution within the wetted zone varied between the three soil textures. An analysis of spatial moments showed that the moisture plume characteristics stabilized within seven days of irrigation and RWU. Zones of effective drying were observed in all three soils during the 21-d period. The alfalfa crop was not negatively affected by the presence of these zones, because (compensatory) RWU was concentrated around the SDI emitters. The results of these experiments imply that the design of SDI line spacing can be less focused on soil type when root system size is considered. This does require the presence of sufficient early season moisture for the crop to develop a complete root system.

Evidence of Preferential Flow Activation in the Vadose Zone via Geophysical Monitoring

Preferential pathways allow rapid and non-uniform water movement in the subsurface due to strong heterogeneity of texture, composition, and hydraulic properties. Understanding the importance of preferential pathways is crucial, because they have strong impact on flow and transport hydrodynamics in the unsaturated zone. Particularly, improving knowledge of the water dynamics is essential for estimating travel time through soil to quantify hazards for groundwater, assess aquifer recharge rates, improve agricultural water management, and prevent surface stormflow and flooding hazards. Small scale field heterogeneities cannot be always captured by the limited number of point scale measurements collected. In order to overcome these limitations, noninvasive geophysical techniques have been widely used in the last decade to predict hydrodynamic processes, due to their capability to spatialize hydrogeophysical properties with high resolution. In the test site located in Bari, Southern Italy, the geophysical approach, based on electrical resistivity tomography (ERT) monitoring, has been implemented to detect preferential pathways triggered by an artificial rainfall event. ERT-derived soil moisture estimations were obtained in order to quantitatively predict the water storage (m3m−3), water velocity (ms−1), and spread (m2) through preferential pathways by using spatial moments analysis.

Electrical monitoring of saline tracers to reveal subsurface flow pathways in a flat ditch-drained field

Time-lapse electrical resistivity imaging (ERI) combined with salt tracers can enhance our understanding of subsurface flow and transport in ditch-drained agroecosystems. In a nearly level field adjacent to a drainage ditch, we installed 192 electrodes in a 72-m2 plot along with five hydrometric stations to enable time-lapse ERI monitoring during storms. Immediately prior to a wintertime rainstorm, we applied 5 kg of KBr to a shallow trench located roughly 11 m upgradient from the ditch. Antecedent conditions prior to the storm were wet, assuring a hydrologically connected state that favored tracer dissolution and movement in the subsurface. We monitored the evolution of the tracer plume with ERI over a seven-day period. Four-dimensional (4D) delineation of precipitation and tracer transport using time-lapse ERI showed that the tracer plume moved preferentially downward through a discontinuity in the argillic horizon within the first 57 hrs. Over the following 79 hrs, the tracer plume exhibited steady, lateral movement towards the drainage ditch within a more permeable sandy horizon situated just below the argillic. A spatial moment analysis of logarithmic changes in conductivity in the sand layer revealed seepage velocities ranging from 0.0027 to 0.0482 m d−1; these values were positively correlated with hydraulic gradients, suggesting laminar flow processes typical of sandy aquifers. Increasing groundwater velocities with hydraulic gradients suggest that solute transfers from groundwater to ditch are likely to be faster in storms, especially during wet periods. Even so, predominantly slow lateral transport within the sandy horizon, even during storms, indicates that critical source areas of pollutant loss are likely within a few meters of the ditch.

Development of quantitative metrics of plume migration at geologic CO2 storage sites

AbstractWe develop plume migration metrics based on spatial moment analysis methods that quantify the spatio‐temporal evolution of plumes at geologic CO2 storage sites. The metrics are generalized to handle any 3‐D scalar attribute field values. Within the geologic CO2 storage context, these can be parameters such as CO2 saturation, effective pressure, overpressure, dissolved CO2 concentration, total dissolved solids, pH, and other attributes that are critical for assessing risks. The metrics are comprehensive in that they can effectively handle and account for complex continuous and discontinuous plumes and intra‐plume migration. We demonstrate the metrics on simulated CO2 plumes injected into flat and tilted reservoirs with homogeneous and heterogeneous permeability fields. Using these idealized reservoir scenarios, we demonstrate the information that the metrics extract, showing that the metrics elucidate nuances in plume migration not apparent by standard approaches to the scalar fields values. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Feasibility Assessment of Long-Term Electrical Resistivity Monitoring of a Nitrate Plume.

Long-term monitoring solutions at contaminated sites are necessary to track plume migration and evaluate the performance of remediation efforts. Electrical resistivity imaging (ERI) can potentially provide information about plume dynamics; however, the feasibility and likelihood of success are seldom evaluated before conducting a field study. Coupling flow and transport models with geoelectrical models provide a powerful way to assess the potential effectiveness of an actual ERI field campaign. We present a coupled approach for evaluating the feasibility of monitoring nitrate migration and remediation using 4D time-lapse ERI at a legacy nuclear waste facility. This kilometer-scale study focuses on depths below the water table (∼70 m). A flow and transport model is developed to perform simulations of nitrate migration and removal via a hypothetical pump-and-treat system. A tracer injection is also simulated at the leading edge of the nitrate plume to enhance the conductivity contrast between the native subsurface and the groundwater fluids. Images of absolute bulk conductivity provide limited information concerning plume migration while time-lapse difference images, which remove the static effects of geology, provide more useful information concerning plume dynamics over time. A spatial moment analysis performed on flow and transport and ERI models matches well during the tracer injection; however, inversion regularization smoothing otherwise limits the value in terms of locating the center of mass. We find that the addition of a tracer enables ERI to characterize plume dynamics during pump-and-treat operations, and late-time ERI monitoring provides a conservative estimate of nitrate plume boundaries in this synthetic study.

Interaction of dissolution, sorption and biodegradation on transport of BTEX in a saturated groundwater system: Numerical modeling and spatial moment analysis

Interaction of various physical, chemical and biological transport processes plays an important role in deciding the fate and migration of contaminants in groundwater systems. In this study, a numerical investigation on the interaction of various transport processes of BTEX in a saturated groundwater system is carried out. In addition, the multi-component dissolution from a residual BTEX source under unsteady flow conditions is incorporated in the modeling framework. The model considers Benzene, Toluene, Ethyl Benzene and Xylene dissolving from the residual BTEX source zone to undergo sorption and aerobic biodegradation within the groundwater aquifer. Spatial concentration profiles of dissolved BTEX components under the interaction of various sorption and biodegradation conditions have been studied. Subsequently, a spatial moment analysis is carried out to analyze the effect of interaction of various transport processes on the total dissolved mass and the mobility of dissolved BTEX components. Results from the present numerical study suggest that the interaction of dissolution, sorption and biodegradation significantly influence the spatial distribution of dissolved BTEX components within the saturated groundwater system. Mobility of dissolved BTEX components is also found to be affected by the interaction of these transport processes.

Mathematical Modeling on Mobility and Spreading of BTEX in a Discretely Fractured Aquifer System Under the Coupled Effect of Dissolution, Sorption, and Biodegradation

This paper presents a numerical model to investigate the migration of BTEX (benzene, toluene, ethylbenzene, and xylene) within a fractured aquifer system at the scale of a single fracture under the coupled effect of various transport processes such as dissolution, sorption, biodegradation. The developed model also considers the influence of equilibrium- and kinetic-controlled sorption scenarios on BTEX transport. Further, the transport characteristics of dissolved BTEX within the fracture–matrix system (FMS) are obtained by carrying out the spatial moment analysis on the concentration profiles simulated. In this study, the spatial moment analysis is conducted to estimate the following transport characteristics of dissolved BTEX: (a) velocity within the fracture, (b) dispersion coefficient within the fracture, (c) dissolved mass within the matrix. In order to investigate the sensitivity of various input parameters, two sensitivity indices are computed based on the variation in the velocity of dissolved BTEX constituent within the fracture (SI_vel) and on the variation in the dissolved mass of BTEX constituents within the matrix (SI_mat). Results from the present simulation study suggest that the sorption and biodegradation reactions influence the concentration distribution of highly soluble BTEX constituents (benzene, toluene) within the FMS significantly. The influence of biodegradation on the migration of BTEX within the FMS is found to be more when it co-occurs with the sorption reaction. The effect of sorption and biodegradation reactions on the mobility and dispersion coefficient of dissolved BTEX constituents within the FMS is found to be significant during the early simulation period.

Micro-positron emission tomography for measuring sub-core scale single and multiphase transport parameters in porous media

Accurate descriptions of heterogeneity in porous media are important for understanding and modeling single phase (e.g. contaminant transport, saltwater intrusion) and multiphase (e.g. geologic carbon storage, enhanced oil recovery) transport problems. Application of medical imaging to experimentally quantify these processes has led to significant progress in material characterization and understanding fluid transport behavior at laboratory scales. While widely utilized in cancer diagnosis and management, cardiology, and neurology, positron emission tomography (PET) has had relatively limited applications in earth science. This study utilizes a small-bore micro-PET scanner to image and quantify the transport behavior of pulses of a conservative aqueous radiotracer injected during single and multiphase flow experiments in two heterogeneous Berea sandstone cores. The cores are discretized into axial-parallel streamtubes, and using the reconstructed micro-PET data, expressions are derived from spatial moment analysis for calculating sub-core tracer flux and pore water velocity. Using the flux and velocity measurements, it is possible to calculate porosity and saturation from volumetric flux balance, and calculate permeability and water relative permeability from Darcy’s law. Second spatial moment analysis enables measurement of sub-core solute dispersion during both single phase and multiphase experiments. A numerical simulation model is developed to verify the assumptions of the streamtube dimension reduction technique. A variation of the reactor ratio is presented as a diagnostic metric to efficiently determine the validity of the streamtube approximation in core and column-scale experiments. This study introduces a new method to quantify sub-core permeability, relative permeability, and dispersion. These experimental and analytical methods provide a foundation for future work on experimental measurements of differences in transport behavior across scales.

Spatial moment analysis of multispecies contaminant transport in porous media

Spatial moment analysis of multispecies contaminant transport in porous media

Predicting suitability of different scale-dependent dispersivities for reactive solute transport through stratified porous media

In this paper, the behavior of breakthrough curves (BTCs) for reactive solute transport through stratified porous media is investigated. A physical laboratory model for layered porous media was constructed, in which thin layer of gravel was sandwiched in between two thick layers of natural soil. Gravel layer and natural soil layers were hydraulically connected as single porous continuum. A constant source of tracer was connected through gravel layer and elucidated at different sampling points in the direction of flow. Flexible multiprocess non-equilibrium (MPNE) transport equation with scale-dependent dispersivity function was used to simulate experimental BTCs of reactive solute transport through layered porous media. The values of equilibrium sorption coefficient and other input parameters were obtained experimentally. The simulation of BTC was performed using MPNE model with scale-dependent dispersivity. The simulation of different scale-dependent dispersivities was then compared and it was found that for field scale of estimation of dispersivity, asymptotic and exponential dispersivity functions performed better. In continuation to the comparison of simulated BTCs obtained using different models, spatial moment analysis of each aforesaid scale-dependent dispersivity model was also done. Spatial moment analysis provides the information related to mean solute mass, rate of mass travel, and mean plume dispersion. Linear and constant dispersivities showed higher variance as compared to asymptotic and exponential dispersion functions. This supports the field applicability of asymptotic and exponential dispersivity functions. The BTCs were also found to elucidate a nonzero concentration with time, which was clearly affected by physical non-equilibrium. In natural condition, such information is required in effective aquifer remediation process.

Characterization of reactive transport by 3-D electrical resistivity tomography (ERT) under unsaturated conditions

The leaching of nitrate from intensively used arable soil is of major concern in many countries. In this study, we show how time lapse electrical resistivity tomography (ERT) can be used to characterize spatially heterogeneous processes of ion production, consumption, and transport in soils. A controlled release fertilizer was introduced into an undisturbed soil core in a laboratory lysimeter and subjected to infiltration events. The production of ions resulting from processes associated with nitrification and their transport through the soil core was observed by time lapse ERT and analysis of seepage water samples from a multicompartment sampler. ERT images show development and propagation of a high-conductivity plume from the fertilizer source zone. Molar amounts of nitrate produced in and exported from the soil core could be well reproduced by time lapse ERT using a spatial moment analysis. Furthermore, we observed that several shape measures of local breakthrough-curves (BTCs) of seepage water conductivity and nitrate derived by effluent analyses and BTCs of bulk conductivity derived by ERT are highly correlated, indicating the preservation of spatial differences of the plume breakthrough in the ERT data. Also differences between nitrate breakthrough and a conservative tracer breakthrough can be observed by ERT. However, the estimation of target ion concentrations by ERT is error bound and the smoothing algorithm of the inversion masks spatial conductivity differences. This results in difficulties reproducing spatial differences of ion source functions and variances of travel times. Despite the observed limitations, we conclude that time lapse ERT can be qualitatively and quantitatively informative with respect to processes affecting the fate of nitrate in arable soils.

Experimental analysis of spatial correlation effects on capillary trapping of supercritical CO2 at the intermediate laboratory scale in heterogeneous porous media

AbstractSeveral numerical studies have demonstrated that the heterogeneous nature of typical sedimentary formations can favorably dampen the accumulation of mobile CO2 phase underneath the caprock. Core flooding experiments have also shown that contrasts in capillary entry pressure can lead to buildup of nonwetting fluid phase (NWP) at interfaces between facies. Explicit representation of geological heterogeneity at the intermediate (cm‐to‐m) scale is a powerful approach to identify the key mechanisms that control multiphase flow dynamics in porous media. The ability to carefully control flow regime and permeability contrast at a scale that is relevant to CO2 plume dynamics in saline formations offers valuable information to understand immiscible displacement processes and provides a benchmark for mathematical models. To provide insight into the impact of capillary heterogeneity on flow dynamics and trapping efficiency of supercritical CO2 under successive drainage and imbibition conditions, we present an experimental investigation conducted in a synthetic sand reservoir. By mimicking the interplay of governing forces at reservoir conditions via application of surrogate fluids, we performed three immiscible displacement experiments to observe the entrapment of NWP in heterogeneous porous media. Capillary trapping performance is evaluated for each scenario through spatial and temporal variations of NWP saturation; for this reason we adopted X‐ray attenuation to precisely measure phase saturation throughout the flow domain and apply spatial moment analysis. The sweeping performance of two different permeability fields with comparable variance but distinct spatial correlation was compared against a homogeneous base case with equivalent mean permeability by means of spatial moment analysis.