Towards a digital twin for characterising natural source zone depletion: A feasibility study based on the Bemidji site

Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) may be a valid long-term management option at petroleum impacted sites. However, its future long-term reliability needs to be established. NSZD includes partitioning, biotic and abiotic degradation of LNAPL components plus multiphase fluid dynamics in the subsurface. Over time, LNAPL components are depleted and those partitioning to various phases change, as do those available for biodegradation. To accommodate these processes and predict trends and NSZD over decades to centuries, for the first time, we incorporated a multi-phase multi-component multi-microbe non-isothermal approach to representatively simulate NSZD at field scale. To validate the approach we successfully mimic data from the LNAPL release at the Bemidji site. We simulate the entire depth of saturated and unsaturated zones over the 27 years of post-release measurements. The study progresses the idea of creating a generic digital twin of NSZD processes and future trends. Outcomes show the feasibility and affordability of such detailed computational approaches to improve decision-making for site management and restoration strategies. The study provided a basis to progress a computational digital twin for complex subsurface systems

Towards characterizing LNAPL remediation endpoints

Remediating sites contaminated with light non-aqueous phase liquids (LNAPLs) is a demanding and often prolonged task. It is vital to determine when it is appropriate to cease engineered remedial efforts based on the long-term effectiveness of remediation technology options. For the first time, the long term effectiveness of a range of LNAPL remediation approaches including skimming and vacuum-enhanced skimming each with and without water table drawdown was simulated through a multi-phase and multi-component approach. LNAPL components of gasoline were simulated to show how component changes affect the LNAPL\u27s multi-phase behaviour and to inform the risk profile of the LNAPL. The four remediation approaches along with five types of soils, two states of the LNAPL specific mass and finite and infinite LNAPL plumes resulted in 80 simulation scenarios. Effective conservative mass removal endpoints for all the simulations were determined. As a key driver of risk, the persistence and mass removal of benzene was investigated across the scenarios. The time to effectively achieve a technology endpoint varied from 2 to 6 years. The recovered LNAPL in the liquid phase varied from 5% to 53% of the initial mass. The recovered LNAPL mass as extracted vapour was also quantified. Additional mass loss through induced biodegradation was not determined. Across numerous field conditions and release incidents, graphical outcomes provide conservative (i.e. more prolonged or greater mass recovery potential) LNAPL remediation endpoints for use in discussing the halting or continuance of engineered remedial efforts

On quantifying global carbon emission from oil contaminated lands over centuries

Petroleum releases into the subsurface contribute to global soil carbon emissions. Quantifying releases and changes in releases of carbon from soils over the lifetime of a spill is complex. Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) embodies all key mechanisms for transformation to carbon gases and their release from soils including partitioning, transport and degradation of petroleum components. Quantification of the interconnected behaviours of the soil microbiome, fluid flow, multi-component transport, partitioning, and biodegradation is crucial for understanding NSZD. Volatilization from LNAPL, aerobic biodegradation, methanogenesis, and heat production all lead to release of greenhouse gases to the atmosphere. To estimate carbon emissions, using a validated computational platform, we modelled the long term NSZD of four petroleum hydrocarbon types; crude oil, diesel, jet fuel and gasoline, to span the major products used globally. For two soil types, we estimated 150 years of carbon emissions from annual minor and 25 mostly major petroleum hydrocarbon land release incidents since 1950 – with an estimated released mass of ~9 million tonnes across the circumstances considered. Up to 2100 the mass of carbon emitted to the atmosphere is estimated to range from 4 to 6 Teragrams, with nearly 60 % currently released. Nomographs generated help predict the fate of LNAPL plumes and carbon emissions due to NSZD, which is crucially important to management of soil and groundwater contamination. The method provides a basis to include additionally identified and future petroleum releases. It is noted that the petroleum mixture composition, degradation rates, volatilization, and subsurface characteristics all can influence carbon emission estimations

Quantifying the benefits of in-time and in-place responses to remediate acute LNAPL release incidents

Acute large volume spills from storage tanks of petroleum hydrocarbons as light non aqueous phase liquids (LNAPLs) can contaminate soil and groundwater and may have the potential to pose explosive and other risks. In consideration of an acute LNAPL release scenario, we explore the value of a rapid remediation response, and the value of installing remediation infrastructure in close proximity to the spill location, in effecting greater recovery of LNAPL mass from the subsurface. For the first time, a verified three-dimensional multi-phase numerical framework and supercomputing resources was applied to explore the significance of in-time and in-place remediation actions. A sand aquifer, two release volumes and a low viscosity LNAPL were considered in key scenarios. The time of commencement of LNAPL remediation activities and the location of recovery wells were assessed requiring asymmetric computational considerations. The volume of LNAPL released considerably affected the depth of LNAPL penetration below the groundwater table, the radius of the plume over time and the recoverable LNAPL mass. The remediation efficiency was almost linearly correlated with the commencement time, but was a non-linear function of the distance of an extraction well from the spill release point. The ratio of the recovered LNAPL in a well located at the centre of the spill/release compared to a well located 5 m away was more than 3.5, for recovery starting only 7 days after the release. Early commencement of remediation with a recovery well located at the centre of the plume was estimated to recover 190 times more LNAPL mass than a one-month delayed commencement through a well 15 m away from the centre of the LNAPL plume. Optimally, nearly 40% of the initially released LNAPL could be recovered within two months of commencing LNAPL recovery actions

Advective and diffusive gas phase transport in vadose zones: Importance for defining vapour risks and natural source zone depletion of petroleum hydrocarbons

Quantifying the interlinked behaviour of the soil microbiome, fluid flow, multi-component transport and partitioning, and biodegradation is key to characterising vapour risks and natural source zone depletion (NSZD) of light non-aqueous phase liquid (LNAPL) petroleum hydrocarbons. Critical to vapour transport and NSZD is transport of gases through the vadose zone (oxygen from the atmosphere, volatile organic compounds (VOCs), methane and carbon dioxide from the zone of LNAPL biodegradation). Volatilisation of VOCs from LNAPL, aerobic biodegradation, methanogenesis and heat production all generate gas pressure changes that may lead to enhanced gas fluxes apart from diffusion. Despite the importance of the gaseous phase dynamics in the vadose zone processes, the relative pressure changes and consequent scales of advective (buoyancy and pressure driven) / diffusive transport is less studied. We use a validated multi-phase multi-component non-isothermal modelling framework to differentiate gas transport mechanisms. We simulate a multicomponent unweathered gasoline LNAPL with high VOC content to maximise the potential for pressure changes due to volatilisation and to enable the joint effects of methanogenesis and shallower aerobic biodegradation of vapours to be assessed, along with heat production. Considering a uniform fine sand profile with LNAPL resident in the water table capillary zone, results suggest that biodegradation plays the key role in gas phase formation and consequent pressure build-up. Results suggest that advection is the main transport mechanism over a thin zone inside the LNAPL/capillary region, where the effective gaseous diffusion is very low. In the bulk of the vadose zone above the LNAPL region, the pressure change is minimal, and gaseous diffusion is dominant. Even for high biodegradation rate cases, pressure build-up due to heat generation (inducing buoyancy effects) is smaller than the contribution of gas formation due to biodegradation. The findings are critical to support broader assumptions of diffusive transport being dominant in vapour transport and NSZD assessments

Evaluating an analytical model to predict subsurface LNAPL distributions and transmissivity from current and historic fluid levels in groundwater wells: Comparing results to numerical simulations

A recent analytical model predicts free, entrapped, and residual LNAPL saturations and the LNAPL transmissivity in the subsurface from current and historic fluid levels in groundwater wells. As such, the model accounts for effects of fluid level fluctuations in a well. The model was developed to predict LNAPL specific volumes and transmissivities from current fluid level measurements in wells and either recorded historic fluid level fluctuations in wells or estimates. An assumption is made in the model that the predictions are not dependent on whether the historic highest or lowest fluid level elevations in a well occur first. To test the assumption, we conduct two simulations with a modified multiphase flow numerical code TMVOC that incorporates relative permeability-saturation-capillary head relations employed in the model. In one simulation, the initial condition is for fluid levels in a well at the historic highest elevations. In the other simulation, the initial condition is for fluid levels in a well at the historic lowest elevations. We change the boundary conditions so both historical conditions occur followed by generating the current condition. Results from the numerical simulations are compared to model predictions and show the assumption in the analytical model is reasonable. The analytical model can be used to develop/refine conceptual site models and for assessing potential LNAPL recovery endpoints, especially on sites with fluctuating fluid levels in wells

Natural source zone depletion of LNAPL: A critical review supporting modelling approaches

Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) includes partitioning, transport and degradation of LNAPL components. NSZD is being considered as a site closure option during later stages of active remediation of LNAPL contaminated sites, and where LNAPL mass removal is limiting. To ensure NSZD meets compliance criteria and to design enhanced NSZD actions if required, residual risks posed by LNAPL and its long term behaviour require estimation. Prediction of long-term NSZD trends requires linking physicochemical partitioning and transport processes with bioprocesses at multiple scales within a modelling framework. Here we expand and build on the knowledge base of a recent review of NSZD, to establish the key processes and understanding required to model NSZD long term. We describe key challenges to our understanding, inclusive of the dominance of methanogenic or aerobic biodegradation processes, the potentially changeability of rates due to the weathering profile of LNAPL product types and ages, and linkages to underlying bioprocesses. We critically discuss different scales in subsurface simulation and modelling of NSZD. Focusing on processes at Darcy scale, 36 models addressing processes of importance to NSZD are investigated. We investigate the capabilities of models to accommodate more than 20 subsurface transport and transformation phenomena and present comparisons in several tables. We discuss the applicability of each group of models for specific site conditions

Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder

A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.Comment: 20 pages, 12 figures. submitted to Proc. SPIE, Astronomical Telescopes + Instrumentation (2014

Calibrating CHIME, A New Radio Interferometer to Probe Dark Energy

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements

MLPerf Inference Benchmark

Machine-learning (ML) hardware and software system demand is burgeoning. Driven by ML applications, the number of different ML inference systems has exploded. Over 100 organizations are building ML inference chips, and the systems that incorporate existing models span at least three orders of magnitude in power consumption and five orders of magnitude in performance; they range from embedded devices to data-center solutions. Fueling the hardware are a dozen or more software frameworks and libraries. The myriad combinations of ML hardware and ML software make assessing ML-system performance in an architecture-neutral, representative, and reproducible manner challenging. There is a clear need for industry-wide standard ML benchmarking and evaluation criteria. MLPerf Inference answers that call. In this paper, we present our benchmarking method for evaluating ML inference systems. Driven by more than 30 organizations as well as more than 200 ML engineers and practitioners, MLPerf prescribes a set of rules and best practices to ensure comparability across systems with wildly differing architectures. The first call for submissions garnered more than 600 reproducible inference-performance measurements from 14 organizations, representing over 30 systems that showcase a wide range of capabilities. The submissions attest to the benchmark's flexibility and adaptability.Comment: ISCA 202