Toward a Practical Appraisal for Waveform Tomography of Band- and Offset-Limited Marine Seismic Data

• Published in: IEEE transactions on geoscience and remote sensing 2022, Vol.60, p.1-14
• Main Authors: Gras, Claudia, Jimenez-Tejero, Clara Estela, Sallares, Valenti, Melendez, Adria, Ranero, Cesar R.
• Format: Article
• Language: English
• Subjects: Acoustic waves; Analytical models; Data models; Datasets; downward continuation (DC); full-waveform inversion; Geology; Geometry; High resolution; inverse theory; Inversions; marine seismic data; Model accuracy; Nonlinear systems; Nonlinearity; numerical modeling; Objective function; Ocean floor; P waves; Propagation; Records; Resolution; Seismic data; Seismic surveys; Seismic velocities; Seismograms; signal analysis; Streaming media; Tomography; Travel time; travel-time tomography (TTT); Velocity; velocity models; Wave velocity; waveform tomography; Waveforms; Wavelengths; Workflow
• ISSN: 0196-2892; 1558-0644
• DOI: 10.1109/TGRS.2021.3097966

We present a generalized workflow to retrieve high-resolution P-wave velocity ( V_{\mathrm{ P}} ) models of complex Earth's subsurface structures from traditional marine near-vertical seismic reflection experiments. These records have typically offsets too short to map refraction phases and lack low-frequency information. The workflow is composed of three steps: 1) downward continuation (DC) of seismic records to the seafloor to recover diving wave information; 2) travel-time tomography (TTT) of first arrivals obtained from DC data, to retrieve a kinematically correct model; and 3) full-waveform inversion (FWI) of the original streamer dataset, starting with the model obtained with TTT and sequentially introducing higher wavenumber details into the model. We show that the TTT allows overcoming the issues associated with the nonlinearity intrinsic to FWI. We also disentangle envelope and phase from the waveform to choose the objective function most suitable for FWI. We assess the accuracy of initial models and predict the quality of the FWI results by quantifying the early arrival cycle skipping between original and simulated data. The efficiency of the workflow is tested with a challenging synthetic target model, containing vertical boundaries with strong velocity contrasts and velocity inversions embedded in a checkerboard-like pattern. We show that workflow steps 1) and 2) provide a TTT model that is not cycle skipped at the frequencies available in most marine seismic experiments and thus allow step 3) FWI to obtain high-resolution V_{\mathrm{ P}} models of the subsurface using band- and offset-limited field datasets, traditionally collected in marine airgun and streamer acquisitions.