Register      Login
Exploration Geophysics Exploration Geophysics Society
Journal of the Australian Society of Exploration Geophysicists
RESEARCH ARTICLE

Acoustic wave propagation simulation by reduced order modelling

Hadi Mahdavi Basir 1 Abdolrahim Javaherian 1 2 6 Zaher Hossein Shomali 2 3 Roohollah Dehghani Firouz-Abadi 4 Shaban Ali Gholamy 5
+ Author Affiliations
- Author Affiliations

1 Department of Petroleum Engineering, Amirkabir University of Technology, Tehran 15875-4413, Iran.

2 Institute of Geophysics, University of Tehran, Tehran 14155-6466, Iran.

3 Department of Earth Sciences, Uppsala University, Uppsala 75236, Sweden.

4 Department of Aerospace Engineering, Sharif University of Technology, Tehran 11365-11155, Iran.

5 Department of Geophysics, Exploration Directorate of National Iranian Oil Company, Tehran 19948-14695, Iran.

6 Corresponding author. Email: javaherian@aut.ac.ir

Exploration Geophysics 49(3) 386-397 https://doi.org/10.1071/EG16144
Submitted: 23 November 2016  Accepted: 3 May 2017   Published: 13 June 2017

Abstract

Wave propagation simulation, as an essential part of many algorithms in seismic exploration, is associated with high computational cost. Reduced order modelling (ROM) is a known technique in many different applications that can reduce the computational cost of simulation by employing an approximation of the model parameters. ROM can be carried out using different algorithms. The method proposed in this work is based on using the most important mode shapes of the model, which can be computed by an efficient numerical method. The numerical accuracy and computational performance of the proposed method were investigated over a simple synthetic velocity model and a portion of the SEG/EAGE velocity model. Different boundary conditions were discussed, and among them the random boundary condition had higher performance for applications like reverse time migration (RTM). The capability of the proposed method for RTM was evaluated and confirmed by the synthetic velocity model of SEG/EAGE. The results showed that the proposed ROM method, compared with the conventional finite element method (FEM), can decrease the computational cost of wave propagation modelling for applications with many simulations like the reverse time migration. Depending on the number of simulations, the proposed method can increase the computational efficiency by several orders of magnitudes.

Key words: acoustic wave propagation simulation, finite element method (FEM), reduced order modelling (ROM), seismic modelling.


References

Almuhaidib, A. M., and Toksöz, M. N., 2015, Finite difference elastic wave modeling with an irregular free surface using ADER scheme: Journal of Geophysics and Engineering, 12, 435–447
Finite difference elastic wave modeling with an irregular free surface using ADER scheme:Crossref | GoogleScholarGoogle Scholar |

Bai, Z., and Su, Y., 2005, Dimension reduction of large scale second-order dynamical systems via a second order Arnoldi method: SIAM Journal on Scientific Computing, 26, 1692–1709
Dimension reduction of large scale second-order dynamical systems via a second order Arnoldi method:Crossref | GoogleScholarGoogle Scholar |

Benner, P., 2004, Solving large-scale control problems: IEEE Control Systems Magazine, 24, 44–59
Solving large-scale control problems:Crossref | GoogleScholarGoogle Scholar |

Carcione, J. M., 1994, The wave equation in generalized coordinates: Geophysics, 59, 1911–1919
The wave equation in generalized coordinates:Crossref | GoogleScholarGoogle Scholar |

Carlberg, K., Cortial, J., Amsallem, D., Zahr, M., and Farhat, C. 2011, The GNAT non-linear model reduction method and its application to fluid dynamics problems: 6th AIAA Theoretical Fluid Mechanics Conference, Honolulu, Hawaii, 1–24.

Chaniotis, D., and Pai, M. A., 2005, Model reduction in power systems using Krylov subspace methods: IEEE Transactions on Power Systems, 20, 888–894

Clapp, R. G., 2009, Reverse-time migration with random boundaries: 79th Annual International Meeting, SEG, Expanded Abstracts, 2809–2813.

De Basabe, J., and Sen, M., 2009, New developments in the finite-element method for seismic modeling: The Leading Edge, 28, 562–567
New developments in the finite-element method for seismic modeling:Crossref | GoogleScholarGoogle Scholar |

Firouz-Abadi, R. D., Haddadpour, H., and Ghasemi, M., 2009, Reduced order modeling of liquid sloshing in 3D tanks using boundary element method: Engineering Analysis with Boundary Elements, 33, 750–761
Reduced order modeling of liquid sloshing in 3D tanks using boundary element method:Crossref | GoogleScholarGoogle Scholar |

Firouz-Abadi, R. D., Ghasemi, M., and Haddadpour, H., 2011, A modal approach to second order analysis of sloshing using boundary element method: Ocean Engineering, 38, 11–21
A modal approach to second order analysis of sloshing using boundary element method:Crossref | GoogleScholarGoogle Scholar |

Glover, K., 1984, All optimal Hankel-norm approximations of linear multivariable systems and their L, ∞-error bounds: International Journal of Control, 39, 1115–1193
All optimal Hankel-norm approximations of linear multivariable systems and their L, -error bounds:Crossref | GoogleScholarGoogle Scholar |

Grimme, E. J., 1997, Krylov projection methods for model reduction: Ph.D. thesis, University of Illinois, Urbana-Champaign.

Haddadpour, H., Behbahani-Nejad, M., and Firouz-Abadi, R. D., 2007, A reduced order aerodynamic model for aeroelastic analysis of complex configurations in incompressible flow: Journal of Aircraft, 44, 1015–1019
A reduced order aerodynamic model for aeroelastic analysis of complex configurations in incompressible flow:Crossref | GoogleScholarGoogle Scholar |

Komatitsch, D., and Vilotte, J. P., 1998, The spectral element method: an efficient tool to simulate the seismic response of 2D and 3D geological structures: Bulletin of the Seismological Society of America, 88, 368–392

Liang, K., Ruess, M., and Abdalla, M., 2016, An eigenanalysis-based bifurcation indicator proposed in the framework of a reduced-order modeling technique for non-linear structural analysis: International Journal of Non-Linear Mechanics, 81, 129–138
An eigenanalysis-based bifurcation indicator proposed in the framework of a reduced-order modeling technique for non-linear structural analysis:Crossref | GoogleScholarGoogle Scholar |

Liu, S., Li, X., Wang, W., and Liu, Y., 2014, A mixed-grid finite element method with PML absorbing boundary conditions for seismic wave modelling: Journal of Geophysics and Engineering, 11, 055009
A mixed-grid finite element method with PML absorbing boundary conditions for seismic wave modelling:Crossref | GoogleScholarGoogle Scholar |

Lucia, D. J., Beran, P. S., and Silva, W. A., 2004, Reduced-order modeling: new approaches for computational physics: Progress in Aerospace Sciences, 40, 51–117
Reduced-order modeling: new approaches for computational physics:Crossref | GoogleScholarGoogle Scholar |

Mahdavi Basir, H., Javaherian, A., Shomali, Z. H., Dehghani Firouzabadi, R., and Rahimi Dalkhani, A., 2015, Using reduced order modeling algorithm for reverse time migration: Second EAGE Workshop on High Performance Computing for Upstream, HPC08b.

Moczo, P., Robertsson, J. O. A., and Eisner, L., 2007, The finite-difference time–domain method for modeling of seismic wave propagation: Advances in Geophysics, 48, 421–516
The finite-difference time–domain method for modeling of seismic wave propagation:Crossref | GoogleScholarGoogle Scholar |

Moore, B., 1981, Principal component analysis in linear systems: controllability, observability, and model reduction: IEEE Transactions on Automatic Control, 26, 17–32
Principal component analysis in linear systems: controllability, observability, and model reduction:Crossref | GoogleScholarGoogle Scholar |

Noorian, M. A., Firouz-Abadi, R. D., and Haddadpour, H., 2012, A reduced order model for liquid sloshing in tanks with flexible baffles using boundary element method: International Journal for Numerical Methods in Engineering, 89, 1652–1664
A reduced order model for liquid sloshing in tanks with flexible baffles using boundary element method:Crossref | GoogleScholarGoogle Scholar |

Oh, J. W., and Min, D. J., 2016, Multi-parameter full waveform inversion using Poisson’s ratio for elastic media: Exploration Geophysics, ,
Multi-parameter full waveform inversion using Poisson’s ratio for elastic media:Crossref | GoogleScholarGoogle Scholar |

Pereyra, V., 2013, Wave equation simulation in two-dimensions using a compressed modeler: American Journal of Computational Mathematics, 3, 231–241
Wave equation simulation in two-dimensions using a compressed modeler:Crossref | GoogleScholarGoogle Scholar |

Pereyra, V., 2016, Model order reduction with oblique projections for large scale wave propagation: Journal of Computational and Applied Mathematics, 295, 103–114
Model order reduction with oblique projections for large scale wave propagation:Crossref | GoogleScholarGoogle Scholar |

Qu, Z. Q., 2004, Model order reduction techniques: with applications in finite element analysis: Springer.

Robertsson, J. O., Bednar, B., Blanch, J., Kostov, C., and Manen, D., 2007, Introduction to the supplement on seismic modeling with applications to acquisition, processing, and interpretation: Geophysics, 72, SM1–SM4
Introduction to the supplement on seismic modeling with applications to acquisition, processing, and interpretation:Crossref | GoogleScholarGoogle Scholar |

Rochus, V., Rixen, D. J., and Golinval, J. C., 2005, Electrostatic coupling of mems structures: transient simulations and dynamic pull-in: Nonlinear Analysis, 63, e1619–e1633
Electrostatic coupling of mems structures: transient simulations and dynamic pull-in:Crossref | GoogleScholarGoogle Scholar |

Sarma, G. S., Mallick, K., and Gadhinglajkar, V. R., 1998, Nonreflecting boundary condition in finite-element formulation for an elastic wave equation: Geophysics, 63, 1006–1016
Nonreflecting boundary condition in finite-element formulation for an elastic wave equation:Crossref | GoogleScholarGoogle Scholar |

Shen, X., and Clapp, R. G., 2011, Random boundary condition for low-frequency wave propagation: SEG Technical Program Expanded Abstracts, 2962–2965.

Sirovich, L., 1987, Turbulence and dynamics of coherent structures: I: Quarterly of Applied Mathematics, 45, 561
Turbulence and dynamics of coherent structures: I:Crossref | GoogleScholarGoogle Scholar |

Symes, W., 2007, Reverse time migration with optimal checkpointing: Geophysics, 72, SM213–SM221
Reverse time migration with optimal checkpointing:Crossref | GoogleScholarGoogle Scholar |

Virieux, J., and Operto, S., 2009, An overview of full-waveform inversion in exploration geophysics: Geophysics, 74, WCC1–WCC26
An overview of full-waveform inversion in exploration geophysics:Crossref | GoogleScholarGoogle Scholar |

Wu, C., Bevc, D., and Pereyra, V., 2013, Model order reduction for efficient seismic modeling: 83rd Annual International Meeting, SEG, Expanded Abstracts, 3360–3364.