Abstract:The finite-difference method, which is an important numerical tool for solving seismic wave equations, is widely applied in simulation, wave-equation-based migration and inversion. As the seismic wave phase plays a critical role in forward simulation and inversion, it should be preserved during wavefield simulation. In this paper, we propose a type of phase-preserving stereomodelling method, which is simultaneously symplectic and low numerical dispersive. First, we propose three new time-marching schemes for solving wave equations that are optimal symplectic partitioned Runge-Kutta schemes with minimized phase errors. Relevant simulations on a harmonic oscillator show that even after 200 000 temporal iterations, our schemes can still avoid the phase drifting issue that appears in other symplectic schemes. We use these symplectic schemes as time integrators, and a numerically low dispersive operator called the stereomodelling discrete operator as a spatial discretization approach to solve seismic wave equations. Theoretical analysis on the stability conditions shows that the new methods are more stable than previous methods. We also investigate the numerical dispersion relations of the methods proposed in this study. To further investigate phase accuracy, we compare the numerical solutions generated by the proposed methods with analytic solutions. Several numerical experiments indicate that our proposed methods are efficient for various models and perform well with perfectly matched layer boundary conditions.

Time discretization of an abstract problem applying to the linearized equations of coupled sound and heat flow

Time discretization of an initial value problem for a simultaneous abstract evolution equation applying to parabolic-hyperbolic phase-field systems

Energy-conserving time-discretisation of abstract dynamic problems with applications in continuum mechanics of solids

A stable and linear time discretization for a thermodynamically consistent model for two-phase incompressible flow

Coupled systems of fractional equations related to sound propagation: analysis and discussion

A Low-Dispersive Symplectic Partitioned Runge-Kutta Method for Solving Seismic-Wave Equations: I. Scheme and Theoretical Analysis

A novel method to study time fractional coupled systems of shallow water equations arising in ocean engineering

Time Discrete Approximation of Weak Solutions for Stochastic Equations of Geophysical Fluid Dynamics and Applications

On properties of a semi-explicit in time fourth-order vector compact scheme for the multidimensional acoustic wave equation

Compressible Euler Equations with Second Sound : Asymptotics of Discontinuous Solutions

A Discretization of Caputo Derivatives with Application to Time Fractional SDEs and Gradient Flows

Implicit Euler time discretization and FDM with Newton method in nonlinear heat transfer modeling

Power Series Method for Time-Integration in Computational Acoustics

Numerical solutions to time-fractional stochastic partial differential equations

First order decoupled method of the primitive equations of the ocean I: Time discretization

The Conservative Time High-Order AVF Compact Finite Difference Schemes for Two-Dimensional Variable Coefficient Acoustic Wave Equations

Initial-boundary value problems for discrete evolution equations: discrete linear Schrodinger and integrable discrete nonlinear Schrodinger equations

Numerical Analysis for a Hyperbolic PDE-Constrained Optimization Problem in Acoustic Full Waveform Inversion

Representation of Linear Waves in Inhomogeneous Media

A mimetic discretization of Westervelt's equation

On solvability of a time-fractional semilinear heat equation, and its quantitative approach to the classical counterpart