Least-squares Fourier reconstruction is basically a discrete linear inverse problem that attempts to recover the Fourier spectrum of the seismic wave-field from irregularly sampled data along the spatial coordinates. The estimated Fourier coefficients are then used to reconstruct the data in a regular grid via a standard inverse Fourier transform (IDFT or IFFT). Unfortunately, this kind of inverse problem is usually under-determined and ill-conditioned. For this reason the LS Fourier reconstruction with minimum norm (FRMN) adopts a damped least-squares inversion to retrieve a unique and stable solution. In this work we show how the damping can introduce artefacts on the reconstructed 3D data. To quantitatively describe this issue, we introduce the concept of “extended” model resolution matrix (EMRM) and we formulate the reconstruction problem as an appraisal problem. Through the simultaneous analysis of the EMRM and of the noise term, we discuss the limits of the FRMN reconstruction and we assess the validity of the reconstructed data and the possible bias introduced by the inversion process. Also, we can guide the parametrization of the forward problem to minimize the occurrence of unwanted artefacts. A simple synthetic example and real data from a 3D marine common shot gather are used to discuss our approach and to show the results of FRMN reconstruction.

Appraisal problem in the 3D least squares Fourier seismic data reconstruction

Stucchi E.;
2015-01-01

Abstract

Least-squares Fourier reconstruction is basically a discrete linear inverse problem that attempts to recover the Fourier spectrum of the seismic wave-field from irregularly sampled data along the spatial coordinates. The estimated Fourier coefficients are then used to reconstruct the data in a regular grid via a standard inverse Fourier transform (IDFT or IFFT). Unfortunately, this kind of inverse problem is usually under-determined and ill-conditioned. For this reason the LS Fourier reconstruction with minimum norm (FRMN) adopts a damped least-squares inversion to retrieve a unique and stable solution. In this work we show how the damping can introduce artefacts on the reconstructed 3D data. To quantitatively describe this issue, we introduce the concept of “extended” model resolution matrix (EMRM) and we formulate the reconstruction problem as an appraisal problem. Through the simultaneous analysis of the EMRM and of the noise term, we discuss the limits of the FRMN reconstruction and we assess the validity of the reconstructed data and the possible bias introduced by the inversion process. Also, we can guide the parametrization of the forward problem to minimize the occurrence of unwanted artefacts. A simple synthetic example and real data from a 3D marine common shot gather are used to discuss our approach and to show the results of FRMN reconstruction.
Ciabarri, F.; Mazzotti, A.; Stucchi, E.; Bienati, N.
File in questo prodotto:
File Dimensione Formato  
Geoph Prosp 2014 3D Fourier Reconstruction.pdf

solo utenti autorizzati

Descrizione: articolo finale
Tipologia: Versione finale editoriale
Licenza: NON PUBBLICO - Accesso privato/ristretto
Dimensione 4.1 MB
Formato Adobe PDF
4.1 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
paper_040414_72.pdf

accesso aperto

Descrizione: Post-print
Tipologia: Documento in Post-print
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 6.69 MB
Formato Adobe PDF
6.69 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11568/611069
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 1
social impact