Time-domain multidimensional deconvolution: A physically reliable and stable preconditioned implementation

David Vargas, Ivan Vasconcelos, Matteo Ravasi, Nick Luiken

Research output: Contribution to journalArticlepeer-review

13 Scopus citations


Multidimensional deconvolution constitutes an essential operation in a variety of geophysical scenarios at different scales ranging from reservoir to crustal, as it appears in applications such as surface multiple elimination, target-oriented redatuming, and interferometric body-wave retrieval just to name a few. Depending on the use case, active, microseismic, or teleseismic signals are used to reconstruct the broadband response that would have been recorded between two observation points as if one were a virtual source. Reconstructing such a response relies on the the solution of an ill-conditioned linear inverse problem sensitive to noise and artifacts due to incomplete acquisition, limited sources, and band-limited data. Typically, this inversion is performed in the Fourier domain where the inverse problem is solved per frequency via direct or iterative solvers. While this inversion is in theory meant to remove spurious events from cross-correlation gathers and to correct amplitudes, difficulties arise in the estimation of optimal regularization parameters, which are worsened by the fact they must be estimated at each frequency independently. Here we show the benefits of formulating the problem in the time domain and introduce a number of physical constraints that naturally drive the inversion towards a reduced set of stable, meaningful solutions. By exploiting reciprocity, time causality, and frequency-wavenumber locality a set of preconditioners are included at minimal additional cost as a way to alleviate the dependency on an optimal damping parameter to stabilize the inversion. With an interferometric redatuming example, we demonstrate how our time domain implementation successfully reconstructs the overburden-free reflection response beneath a complex salt body from noise-contaminated up-and down-going transmission responses at the target level.
Original languageEnglish (US)
Pages (from-to)3683
JournalRemote Sensing
Issue number18
StatePublished - Sep 15 2021

Bibliographical note

KAUST Repository Item: Exported on 2021-09-27
Acknowledgements: This research was funded by The sponsors of the Utrecht Consortium for Subsurface Imaging (UCSI).

ASJC Scopus subject areas

  • General Earth and Planetary Sciences


Dive into the research topics of 'Time-domain multidimensional deconvolution: A physically reliable and stable preconditioned implementation'. Together they form a unique fingerprint.

Cite this