We present a procedure for the numerical computation of the seismicwaves in realistic 3-D geomodels. The procedure is based on a parallel implementation of the Fourier pseudo-spectral method (FPSM) on staggered grids. FPSM solves spatial derivatives with optimal accuracy in domains which are discretized in easy-to-handle, structured meshes. The application of the presented procedure is therefore advantageous, in terms of both accuracy and simplicity of use, whenever the geo-model can be adequately sampled by a structured grid. To be applicable to a wide range of geophysical problems, the procedure implements the most important features needed for modelling the seismic waves propagation, that is intrinsic attenuation, seismic velocity anisotropy and irregular topography. The intrinsic attenuation is implemented by means of the generalized zener body GZB mechanical model. Orthorhombic symmetry is considered for the description of the seismic anisotropy. Irregular topography is approximated by a staircase of cubic cells and treated using an original approach based on discrete Fourier transforms of arbitrary length. Sharp impedance contrasts in the geo-model are handled by suitable averaging of the material properties to reduce phase misalignments and staircase effects due to the structured grid. Absorbing boundaries based on the convolutional perfectly matching layer technique are adopted. Finally, the efficiency of the proposed procedure is illustrated by a number of tests and examples.

Numerical simulation of seismic wave propagation in realistic 3-D geo-models with a Fourier pseudo-spectral method

Klin P.;Priolo E.;
2010-01-01

Abstract

We present a procedure for the numerical computation of the seismicwaves in realistic 3-D geomodels. The procedure is based on a parallel implementation of the Fourier pseudo-spectral method (FPSM) on staggered grids. FPSM solves spatial derivatives with optimal accuracy in domains which are discretized in easy-to-handle, structured meshes. The application of the presented procedure is therefore advantageous, in terms of both accuracy and simplicity of use, whenever the geo-model can be adequately sampled by a structured grid. To be applicable to a wide range of geophysical problems, the procedure implements the most important features needed for modelling the seismic waves propagation, that is intrinsic attenuation, seismic velocity anisotropy and irregular topography. The intrinsic attenuation is implemented by means of the generalized zener body GZB mechanical model. Orthorhombic symmetry is considered for the description of the seismic anisotropy. Irregular topography is approximated by a staircase of cubic cells and treated using an original approach based on discrete Fourier transforms of arbitrary length. Sharp impedance contrasts in the geo-model are handled by suitable averaging of the material properties to reduce phase misalignments and staircase effects due to the structured grid. Absorbing boundaries based on the convolutional perfectly matching layer technique are adopted. Finally, the efficiency of the proposed procedure is illustrated by a number of tests and examples.
2010
Numerical solutions, Seismic anisotropy, Seismic attenuation, Computational seismology, Wave propagation.
File in questo prodotto:
File Dimensione Formato  
FPSM_GJI2010.pdf

non disponibili

Tipologia: Altro materiale allegato
Licenza: Non specificato
Dimensione 618.56 kB
Formato Adobe PDF
618.56 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

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/20.500.14083/3401
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 39
  • ???jsp.display-item.citation.isi??? 37
social impact