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 <i>Deepwater Horizon</i> disaster, represented with 279–280 pseudocomponents, including 131–132 individual compounds. We define ><i>n</i>-C₈ pseudocomponents based on comprehensive two-dimensional gas chromatography (GC × GC) measurements, which enable the modeling of aqueous partitioning for <i>n</i>-C₈ to <i>n</i>-C₂₆ 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 C₁−C₄ 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