Full text loading...
Frequency-domain Full Waveform Inversion (FWI) is amenable to efficient processing of full-azimuth long-offset seabed acquisition carried out with sparsely-deployed ocean bottom nodes (OBNs) because the inversion can be performed with a few discrete frequencies. However, computing the solution of the forward problem efficiently in the frequency domain with linear algebra solvers remains a challenge for large computational domains (tens to hundreds of millions of parameters). We apply 3D frequency-domain FWI to the Gorgon OBN dataset in the NorthWestern shelf, Australia. We solve the forward problem with the multifrontal direct solver MUMPS, which includes four key features: A parallelism combining message-passing interface and multithreading, block low-rank compression, mixed precision arithmetic and efficient block processing of sparse sources. The dataset involves 650 OBNs and 400,000 sources. Mono-parameter visco-acoutic VTI FWI for vertical wavespeed is performed with a classical frequency continuation between 1.7 Hz and 13 Hz. The target covers an area ranging from 260 km2 (frequency > 8.5 Hz) to 705 km2 (frequency < 8.5 Hz) for a maximum depth of 8 km. FWI highlights the bounding faults of the Gorgon horst at reservoir depths as well as several intra-horst faults and several horizons of the Mungaroo formation down 7 km depth.