Comprehensive Two-Dimensional Gas Chromatography with Peak Tracking for Screening of Constituent Biodegradation in Petroleum UVCB Substances.

Petroleum substances, as archetypical UVCBs (substances of unknown or variable composition, complex reaction products, or biological substances), pose a challenge for chemical risk assessment as they contain hundreds to thousands of individual constituents. It is particularly challenging to determine the biodegradability of petroleum substances since each constituent behaves differently. Testing the whole substance provides an average biodegradation, but it would be effectively impossible to obtain all constituents and test them individually. To overcome this challenge, comprehensive two-dimensional gas chromatography (GC × GC) in combination with advanced data-handling algorithms was applied to track and calculate degradation half-times (DT50s) of individual constituents in two dispersed middle distillate gas oils in seawater. By tracking >1000 peaks (representing ∼53-54% of the total mass across the entire chromatographic area), known biodegradation patterns of oil constituents were confirmed and extended to include many hundreds not currently investigated by traditional one-dimensional GC methods. Approximately 95% of the total tracked peak mass biodegraded after 64 days. By tracking the microbial community evolution, a correlation between the presence of functional microbial communities and the observed progression of DT50s between chemical classes was demonstrated. This approach could be used to screen the persistence of GC × GC-amenable constituents of petroleum substance UVCBs.

Relative sensitivity of hydrodynamic, thermodynamic, and chemical processes for simulating the buoyant multiphase plume and surfacing flows of an oil and gas blowout.

Deepwater hydrocarbon releases experience complex chemical and physical processes. To assess simplifications of these processes on model predictions, we present a sensitivity analysis using simulations for the Deepwater Horizon oil spill. We compare the buoyant multiphase plume metrics (trap height, rise time etc), the hydrocarbon mass flowrates at the near-field plume termination and their mass fractions dissolved in the water column and reaching the water surface. The baseline simulation utilizes a 19-component hydrocarbon model, live-fluid state equations, hydrate dynamics, and heat and mass transfer. Other simulations turn-off each of these processes, with the simplest one using inert oil and methane gas. Plume metrics are the least sensitive to the modeled processes and can be matched by adjusting the release buoyancy flux. The mass flowrate metrics are more sensitive. Both liquid- and gas-phase mass transfer should be modeled for accurate tracking of soluble components (e.g. C1 - C7 hydrocarbons) in the environment.

An evaluation of models that estimate droplet size from subsurface oil releases.

Droplet size substantially affects the fate of oil released from deep subsea leaks. A baseline dataset of volume-median droplet diameters (d50), culled from ~250 laboratory observations, is used to validate seven droplet-size models. Four models compare reasonably well, having 95% confidence limits in d50 of ~±50%. Simulations with a near-field fate model (TAMOC) reveals that the four best-performing models, with d50 of 1.3-2.2 mm, agree similarly with observed fractionation of petroleum compounds in the water column during June 4-July 15, 2010. Model results suggest that, had a higher dose of dispersant been applied at the wellhead during Deepwater Horizon oil spill (DWH), the d50 would have dropped by an order of magnitude, reducing surfacing C1-C9 volatiles by 3.5×. Model uncertainty is found to be substantial for DWH-like blowouts treated with chemical dispersants, suggesting the need for further droplet-size model improvement.

Dynamics of Live Oil Droplets and Natural Gas Bubbles in Deep Water.

Explaining the dynamics of gas-saturated live petroleum in deep water remains a challenge. Recently, Pesch et al. [ Environ. Eng. Sci. 2018, 35 (4), 289-299] reported laboratory experiments on methane-saturated oil droplets under emulated deep-water conditions, providing an opportunity to elucidate the underlying dynamical processes. We explain these observations with the Texas A&M Oil spill/Outfall Calculator (TAMOC), which models the pressure-, temperature-, and composition-dependent interactions between oil-gas phase transfer; aqueous dissolution; and densities and volumes of liquid oil droplets, gas bubbles, and two-phase droplet-bubble pairs. TAMOC reveals that aqueous dissolution removed >95% of the methane from ∼3.5 mm live oil droplets within 14.5 min, prior to gas bubble formation, during the experiments of Pesch et al. Additional simulations indicate that aqueous dissolution, fluid density changes, and gas-oil phase transitions (ebullition, condensation) may all contribute to the fates of live oil and gas in deep water, depending on the release conditions. Illustrative model scenarios suggest that 5 mm diameter gas bubbles released at a <470 m water depth can transport methane, ethane, and propane to the water surface. Ethane and propane can reach the water surface from much deeper releases of 5 mm diameter live oil droplets, during which ebullition occurs at water depths of <70 m.

The treatment of biodegradation in models of sub-surface oil spills: A review and sensitivity study.

Biodegradation is important for the fate of oil spilled in marine environments, yet parameterization of biodegradation varies across oil spill models, which usually apply constant first-order decay rates to multiple pseudo-components describing an oil. To understand the influence of model parameterization on the fate of subsurface oil droplets, we reviewed existing algorithms and rates and conducted a model sensitivity study. Droplets were simulated from a blowout at 2000 m depth and were either treated with sub-surface dispersant injection (2% dispersant to oil ratio) or untreated. The most important factor affecting oil fate was the size of the droplets, with biodegradation contributing substantially to the fate of droplets ≤0.5 mm. Oil types, which were similar, had limited influence on simulated oil fate. Model results suggest that knowledge of droplet sizes and improved estimation of pseudo-component biodegradation rates and lag times would enhance prediction of the fate and transport of subsurface oil.

Oil spill modeling in deep waters: Estimation of pseudo-component properties for cubic equations of state from distillation data.

Deep-water oil spills represent a major, localized threat to marine ecosystems. Multi-purpose computer models have been developed to predict the fate of spilled oil. These models include databases of pseudo-components from distillation cut analysis for hundreds of oils, and have been used for guiding response action, damage assessment, and contingency planning for marine oil spills. However, these models are unable to simulate the details of deep-water, high-pressure chemistry. We present a new procedure to calculate the chemical properties necessary for such simulations that we validate with 614 oils from the ADIOS oil library. The calculated properties agree within 20.4% with average values obtained from data for measured compounds, for 90% of the chemical properties. This enables equation-of-state calculations of dead oil density, viscosity, and interfacial tension. This procedure enables development of comprehensive oil spill models to predict the behavior of petroleum fluids in the deep sea.

Simulating Gas-Liquid-Water Partitioning and Fluid Properties of Petroleum under Pressure: Implications for Deep-Sea Blowouts.

With the expansion of offshore petroleum extraction, validated models are needed to simulate the behaviors of petroleum compounds released in deep (>100 m) waters. We present a thermodynamic model of the densities, viscosities, and gas-liquid-water partitioning of petroleum mixtures with varying pressure, temperature, and composition based on the Peng-Robinson equation-of-state and the modified Henry's law (Krychevsky-Kasarnovsky equation). The model is applied to Macondo reservoir fluid released during the Deepwater Horizon disaster, represented with 279-280 pseudocomponents, including 131-132 individual compounds. We define >n-C8 pseudocomponents based on comprehensive two-dimensional gas chromatography (GC × GC) measurements, which enable the modeling of aqueous partitioning for n-C8 to n-C26 fractions not quantified individually. Thermodynamic model predictions are tested against available laboratory data on petroleum liquid densities, gas/liquid volume fractions, and liquid viscosities. We find that the emitted petroleum mixture was ∼29-44% gas and ∼56-71% liquid, after cooling to local conditions near the broken Macondo riser stub (∼153 atm and 4.3 °C). High pressure conditions dramatically favor the aqueous dissolution of C1-C4 hydrocarbons and also influence the buoyancies of bubbles and droplets. Additionally, the simulated densities of emitted petroleum fluids affect previous estimates of the volumetric flow rate of dead oil from the emission source.

Resolving biodegradation patterns of persistent saturated hydrocarbons in weathered oil samples from the Deepwater Horizon disaster.

Biodegradation plays a major role in the natural attenuation of oil spills. However, limited information is available about biodegradation of different saturated hydrocarbon classes in surface environments, despite that oils are composed mostly of saturates, due to the limited ability of conventional gas chromatography (GC) to resolve this compound group. We studied eight weathered oil samples collected from four Gulf of Mexico beaches 12-19 months after the Deepwater Horizon disaster. Using comprehensive two-dimensional gas chromatography (GC × GC), we successfully separated, identified, and quantified several distinct saturates classes in these samples. We find that saturated hydrocarbons eluting after n-C22 dominate the GC-amenable fraction of these weathered samples. This compound group represented 8-10%, or 38-68 thousand metric tons, of the oil originally released from Macondo well. Saturates in the n-C22 to n-C29 elution range were found to be partly biodegraded, but to different relative extents, with ease of biodegradation decreasing in the following order: n-alkanes > methylalkanes and alkylcyclopentanes+alkylcyclohexanes > cyclic and acyclic isoprenoids. We developed a new quantitative index designed to characterize biodegradation of >n-C22 saturates. These results shed new light onto the environmental fate of these persistent, hydrophobic, and mostly overlooked compounds in the unresolved complex mixtures (UCM) of weathered oils.

Oxygenated weathering products of Deepwater Horizon oil come from surprising precursors.

Following the release of crude oil from the Macondo well in 2010, a wide range of weathering processes acted on the spilled oil. A recent study revealed that samples from this spill were oxidized into oxygenated hydrocarbons (OxHC) comprising more than 50% of the extracted hydrocarbons. The precursors of these compounds were not identified despite using a wide range of analytical tools, including gas chromatography (GC). To search for these precursors, over 40 samples were analyzed by comprehensive two-dimensional gas chromatography (GC×GC), one of the largest studies of its kind to date. Partial least squares regression was employed to elucidate the GC×GC peaks that could be the precursors of OxHC in our samples. We found that the formation of OxHC correlated with the disappearance of saturated hydrocarbons, including alkylcyclopentanes, alkyl cyclohexanes, alkylated bicyclic saturated compounds, tricyclic terpanpoids, and alkylbenzenes. These results indicate a previously under-reported chemodynamic process in oil spill weathering.