Statistics – Computation
Scientific paper
Nov 2008
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2008geoji.175..541b&link_type=abstract
Geophysical Journal International, Volume 175, Issue 2, pp. 541-559.
Statistics
Computation
8
Numerical Solutions, Computational Seismology, Wave Propagation
Scientific paper
A new numerical technique for solving 2-D elastodynamic equations based on a finite-volume frequency-domain approach is proposed. This method has been developed as a tool to perform 2-D elastic frequency-domain full-waveform inversion. In this context, the system of linear equations that results from the discretization of the elastodynamic equations is solved with a direct solver, allowing efficient multiple-source simulations at the partial expense of the memory requirement. The discretization of the finite-volume approach is through triangles. Only fluxes with the required quantities are shared between the cells, relaxing the meshing conditions, as compared to finite-element methods. The free surface is described along the edges of the triangles, which can have different slopes. By applying a parsimonious strategy, the stress components are eliminated from the discrete equations and only the velocities are left as unknowns in the triangles. Together with the local support of the P0 finite-volume stencil, the parsimonious approach allows the minimizing of core memory requirements for the simulation. Efficient perfectly matched layer absorbing conditions have been designed for damping the waves around the grid. The numerical dispersion of this FV formulation is similar to that of O(Δx2) staggered-grid finite-difference (FD) formulations when considering structured triangular meshes. The validation has been performed with analytical solutions of several canonical problems and with numerical solutions computed with a well-established FD time-domain method in heterogeneous media. In the presence of a free surface, the finite-volume method requires 10 triangles per wavelength for a flat topography, and fifteen triangles per wavelength for more complex shapes, well below the criteria required by the staircase approximation of O(Δx2) FD methods. Comparisons between the frequency-domain finite-volume and the O(Δx2) rotated FD methods also show that the former is faster and less memory demanding for a given accuracy level, an attractive feature for frequency-domain seismic inversion. We have thus developed an efficient method for 2-D P-SV-wave modelling on structured triangular meshes as a tool for frequency-domain full-waveform inversion. Further work is required to improve the accuracy of the method on unstructured meshes.
Brossier Romain
Operto Stéphane
Virieux Jean
No associations
LandOfFree
Parsimonious finite-volume frequency-domain method for 2-D P-SV-wave modelling does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Parsimonious finite-volume frequency-domain method for 2-D P-SV-wave modelling, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Parsimonious finite-volume frequency-domain method for 2-D P-SV-wave modelling will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-746679