Finite-element modelling of seismic wave propagation on tetrahedra requires meshes that accurately follow interfaces between impedance contrasts or surface topography and have element sizes proportional to the local velocity. We explain a mesh generation approach by example. Starting from a finite-difference representation of the velocity model, triangulated surfaces are generated along impedance discontinuities. These define subdomains that are meshed independently and in parallel, honouring the local velocity values. The resulting volumetric meshes are merged into a single mesh. The approach is flexible, efficient, scalable and capable of producing quality meshes.


Article metrics loading...

Loading full text...

Full text loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error