High-resolution exploration of the submarine crust in the North Sea





A major breakthrough now makes it possible to explore the upper crust with unprecedented precision. Researchers from ISTerre have applied a high-resolution, multi-parameter visco-acoustic full waveform inversion (FWI) to a dense ocean-bottom sensor network deployed in the North Sea.



Figure 5.Cross-sections through the 3-D vertical velocity (a) initial and (b) final FWI models. The stratigraphic table shows the changes in lithology derived from the drilling in the central Graben structure (Munns 1985 ; Kristiansen et al. 2000). The reservoir R is located below the gas-cloud GC in the oil-bearing upper-cretaceous Tor and Hod formations T&H.


This study presents a 3D high-resolution characterization of the upper crust using FWI applied to dense seismic data. The subsurface is parameterized by P-wave velocity, density, quality factor, and VTI anisotropy. Starting from an initial velocity model obtained by reflection tomography and a density derived from Gardner’s law, the inversion uses the reciprocity principle and random source subsampling to reduce computational costs. The computations were carried out on the ADASTRA HPC platform, operated by CINES within the GENCIframework

The resulting models accurately reproduce P-wave velocity and density down to 2–3 km depth, validated by well log measurements, and reveal fine geological structures such as paleo-channels, low-velocity hydrocarbon-rich zones, and networks of fluid-filled faults. Although the study currently uses a visco-acoustic approximation, future extensions toward fully visco-elastic FWI with multicomponent data will allow better characterization of S-waves and fluid discrimination. These results demonstrate the potential of multiparameter FWI for detailed quantitative crustal characterization, bridging the gap between seismic imaging and reservoir analysis.




References :

Aubert, C., and al., 2025.
The MACIV multiscale seismic experiments in the French Massif Central (2023-2027) : deployment, data quality and availability
Annales Geophysicae, 68.

Bozdağ, E., Trampert, J., Tromp, J., 2011.
Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements
Geophysical Journal International, 185(2), 845–870.

Cao, J., Brossier, R., Górszczyk, A., Métivier, L., Virieux, J., 2022b.
3D multi-parameter full-waveform inversion for ocean-bottom seismic data using an efficient fluid-solid coupled spectral-element solver
Geophysical Journal International, 229(1), 671–703.

Espin, I., Salaun, N., Jiang, H., Reinier, M., 2023.
From FWI to ultra-high-resolution imaging
Leading Edge, 42(1), 16–23.

Kamath, N., Brossier, R., Métivier, L., Pladys, A., Yang, P., 2021.
Multiparameter full-waveform inversion of 3D ocean-bottom cable data from the Valhall field
Geophysics, 86(1), B15–B35.

Scientific Contacts :

  • L. Métivier – Researcher, Univ. Grenoble Alpes, ISTerre
  • R. Brossier – Researcher, Univ. Grenoble Alpes, ISTerre
  • A. Górszczyk – Researcher, Univ. Grenoble Alpes, ISTerre ; Institute of Geophysics, Polish Academy of Sciences, Warsaw, Poland
  • G. Provenzano – Researcher, Univ. Grenoble Alpes, ISTerre
  • S. Sambolian – Associate Professor at the Earth and Environment Institute of Strasbourg
  • A. Tarayoun – Researcher, Univ. Grenoble Alpes, ISTerre



Acknowledgements

The authors thank Espen Raknes, Ross Milne, and Einar Jarle Kjos from Aker BP for sharing the data and their assistance. This study was partially funded by the SEISCOPE consortium and benefited from access to HPC resources at GRICAD, Cray, and CINES/IDRIS/TGCC via GENCI.