← Back to Sessions

Quadrangular Adaptive Mesh for Elastic Wave Simulation in Smooth Anisotropic Media

Description: 

The management of smooth anisotropic media holds significant importance for seismic wave simulation and imaging. In such media, adaptive meshes constitute, in principle, an attractive representation to maximize simulation accuracy while minimizing computational cost. However, creating such meshes can be difficult, because the optimal local size of the elements is not clearly defined in the case of smooth anisotropic media.

In this work, we present a two-step algorithm to efficiently mesh 2D these media for spectral element method (SEM) simulations. Our algorithm yields quadrangular meshes which adapt the size of the elements to the local effective and directional S-wave velocity. It first relies on the octree decomposition introduced by Maréchal (2009) to divide the mesh until the size of each element edge is adapted to the local minimum wavelength that will be propagated. Then, a smoothing is applied to further optimize the size and shape of the elements, which increases the global time-step and, consequently, makes SEM simulations faster while keeping a good accuracy or even improving it in some cases. The application of our method on a section of the homogenized Groningen area shows that simulation time can be reduced by a factor up to 3.

Session: Numerical Modeling in Seismology: Developments and Applications - III

Type: Oral

Date: 5/1/2024

Presentation Time: 03:00 PM (local time)

Presenting Author: Marius

Student Presenter: Yes

Invited Presentation: 

Authors