1887

Abstract

Summary

For seismic modelling, imaging and inversion, finite-difference methods are still the workhorse of the industry despite their inability to meet the increasing demand for improved accuracy in subsurface imaging. Finite-element methods offer better accuracy but at a higher computational cost. A stress-velocity formulation with linear elements and an iterative method, defect correction, for inverting the mass matrix offers fourth-order super-convergence but is susceptible to numerical noise if waves in the wrong part of the dispersion curve are excited. We propose an improved source term that reduces that noise and investigate the accuracy of the method on structured triangular meshes as well as on unstructured rotated meshes. With an optimised source function, it is seen that the dispersive wavelengths can be avoided, giving the defect-correction approach a better performance than the mass-lumped formulation with only a marginal increase in compute effort.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201700769
2017-06-12
2024-04-26
Loading full text...

Full text loading...

References

  1. Ainsworth, M.
    [2014] Dispersive behaviour of high order finite element schemes for the one-way wave equation. Journal of Computational Physics, 259, 1–10.
    [Google Scholar]
  2. Hicks, G.J.
    [2002] Arbitrary source and receiver positioning in finite-difference schemes using Kaiser windowed sinc function. Geophysics, 67(1), 156–165.
    [Google Scholar]
  3. Komatitsch, D., Vilotte, J.P., Vai, R., Castillo-Covarrubias, J.M. and Sánchez-Sesma, F.J.
    [1999] The spectral element method for elastic wave equations—application to 2-D and 3-D seismic problems. International Journal for Numerical Methods in Engineering, 45 (9), 1139–1164.
    [Google Scholar]
  4. Mulder, W.A.
    [1996] A comparison between higher-order finite elements and finite differences for solving the wave equation. In: Désidéri, J.A., LeTallec, P., Oftate, E., Périaux, J. and Stein, E. (Eds.) Proceedings of the Second ECCOMAS Conference on Numerical Methods in Engineering, Paris, Sept. 9–13, 1996. John Wiley & Sons, Chichester, 344–350.
    [Google Scholar]
  5. [2013] New triangular mass-lumped finite elements of degree six for wave propagation. Progress In Electromagnetics Research, 141, 671–692.
    [Google Scholar]
  6. Mulder, W.A. and Shamasundar, R.
    [2016] Performance of continuous mass-lumped tetrahedral elements for elastic wave propagation with and without global assembly. Geophysical Journal International, 207(1), 414–421.
    [Google Scholar]
  7. Shamasundar, R. and Mulder, W.A.
    [2016a] Improving the accuracy of mass-lumped finite-elements in the first-order formulation of the wave equation by defect correction. Journal of Computational Physics, 322, 689–707.
    [Google Scholar]
  8. [2016b] Should we use the first- or second-order formulation with spectral elements for seismic modelling?Extended Abstract We P4 16, 78th EAGE Conference & Exhibition, Vienna, Austria, 2016.
    [Google Scholar]
  9. Zhebel, E., Minisini, S., Kononov, A. and Mulder, W.A.
    [2014] A comparison of continuous mass-lumped finite elements with finite differences for 3-D wave propagation. Geophysical Prospecting, 62(5), 1111–1125.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201700769
Loading
/content/papers/10.3997/2214-4609.201700769
Loading

Data & Media 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