A methodology for optimizing modeling configuration in the numerical modeling of oil concentrations in underwater blowouts: a North Sea case study

A METHODOLOGY FOR OPTIMIZING MODELING CONFIGURATION IN THE NUMERICAL MODELING OF OIL CONCENTRATIONS IN UNDERWATER BLOWOUTS: A NORTH SEA CASE STUDY

Andrés Martínez^a,*, Ana J. Abascal^a, Andrés García^a, Beatriz Pérez-Díaz^a, Germán Aragón^a, Raúl Medina^a

^aIHCantabria - Instituto de Hidráulica Ambiental de la Universidad de Cantabria, Avda. Isabel Torres, 15, 39011 Santander, Spain

^* Corresponding author: martinezga@unican.es

Underwater oil and gas blowouts are not easy to repair. It may take months before the well is finally capped, releasing large amounts of oil into the marine environment. In addition, persistent oils (crude oil, fuel oil, etc.) break up and dissipate slowly, so they often reach the shore before the cleanup is completed, affecting vasts extension of seas-oceans, just as posing a major threat to marine organisms.

On account of the above, numerical modeling of underwater blowouts demands great computing power. High-resolution, long-term data bases of wind-ocean currents are needed to be able to properly model the trajectory of the spill at both regional (open sea) and local level (coastline), just as to account for temporal variability. Moreover, a large number of particles, just as a high-resolution grid, are unavoidable in order to ensure accurate modeling of oil concentrations, of utmost importance in risk assessment, so that threshold concentrations can be established (threshold concentrations tell you what level of exposure to a compound could harm marine organisms).

In this study, an innovative methodology has been accomplished for the purpose of optimizing modeling configuration: number of particles and grid resolution, in the modeling of an underwater blowout, with a view to accurately represent oil concentrations, especially when threshold concentrations are considered. In doing so, statistical analyses (dimensionality reduction and clustering techniques), just as numerical modeling, have been applied.

It is composed of the following partial steps: (i) classification of i representative clusters of forcing patterns (based on PCA and K-means algorithms) from long-term wind-ocean current hindcast data bases, so that forcing variability in the study area is accounted for; (ii) definition of j modeling scenarios, based on key blowout parameters (oil type, flow rate, etc.) and modeling configuration (number of particles and grid resolution); (iii) Lagrangian trajectory modeling of the combination of the i clusters of forcing patterns and the j modeling scenarios; (iv) sensitivity analysis of the Lagrangian trajectory model output: oil concentrations,  to modeling configuration; (v) finally, as a result, the optimal modeling configuration, given a certain underwater blowout (its key parameters), is provided.

It has been applied to a hypothetical underwater blowout in the North Sea, one of the world’s most active seas in terms of offshore oil and gas exploration and production. A 5,000 cubic meter per day-flow rate oil spill, flowing from the well over a 15-day period, has been modeled (assuming a 31-day period of subsequent drift for a 46-day modeling). Moreover, threshold concentrations of 0.1, 0.25, 1 and 10 grams per square meter have been applied in the sensitivity analysis. The findings of this study stress the importance of modeling configuration in accurate modeling of oil concentrations, in particular if lower threshold concentrations are considered.