Converted-wave reflection traveltime inversion with free-surface multiples for ocean-bottom-node data

Abstract

Advancements in ocean-bottom-node (OBN) technologies have enabled the recording of high-quality multicomponent seismic data, which can provide crucial information about subsurface properties through elastic seismic imaging. However, imaging multicomponent data is challenging due to the requirement for P- and S-wave velocity models. Conventional processing relies on estimating the S-wave velocity model using a primary converted wave (C wave) because the pure S-wave primaries are not available if using standard air-gun sources. When a free surface is present, the C waves recorded with geophones can be contaminated by source-side P-wave multiples, especially water-layer reverberations, which are difficult to remove even using state-of-the-art wavefield decomposition methods. These C-wave multiples can cause significant mismatches between synthetic and observed data during reflection traveltime or waveform inversion if only primary P-to-S (PS) conversions are simulated. To address these issues and reliably build an S-wave velocity model for PS imaging, we develop a C-wave reflection traveltime inversion method that takes source-side free-surface multiples into account. The traveltime misfit of C-wave data is extracted using dynamic image warping to ensure a better match between the synthetic and observed data. The functional gradient of the objective function is derived using the adjoint state method. Based on sensitivity kernel analyses, a convenient gradient preconditioning strategy is developed to robustly separate the contribution of primary and multiple C waves during backpropagation. This strategy can effectively mitigate the high-wavenumber crosstalk associated with nonphysical wavepaths in the gradient and thus allow more accurate updating for the S-wave velocity model. The synthetic data and shallow-water OBN data from the East China Sea indicate that this approach can reliably recover the S-wave macrovelocity structures for PS imaging.