Abstract:ABSTRACT The curvilinear-grid finite-difference method (FDM), which uses curvilinear coordinates to discretize the nonplanar interface geometry, is extended to simulate acoustic and seismic-wave propagation across the fluid–solid interface at the sea bottom. The coupled acoustic velocity-pressure and elastic velocity-stress formulation that governs wave propagation in seawater and solid earth is expressed in curvilinear coordinates. The formulation is solved on a collocated grid by alternative applications of forward and backward MacCormack finite difference within a fourth-order Runge–Kutta temporal integral scheme. The shape of a fluid–solid interface is discretized by a curvilinear grid to enable a good fit with the topographic interface. This good fit can obtain a higher numerical accuracy than the staircase approximation in the conventional FDM. The challenge is to correctly implement the fluid–solid interface condition, which involves the continuity of tractions and the normal component of the particle velocity, and the discontinuity (slipping) of the tangent component of the particle velocity. The fluid–solid interface condition is derived for curvilinear coordinates and explicitly implemented by a domain-decomposition technique, which splits a grid point on the fluid–solid interface into one grid point for the fluid wavefield and another one for the solid wavefield. Although the conventional FDM that uses effective media parameters near the fluid–solid interface to implicitly approach the boundary condition conflicts with the fluid–solid interface condition. We verify the curvilinear-grid FDM by conducting numerical simulations on several different models and compare the proposed numerical solutions with independent solutions that are calculated by the Luco-Apsel-Chen generalized reflection/transmission method and spectral-element method. Besides, the effects of a nonplanar fluid–solid interface and fluid layer on wavefield propagation are also investigated in a realistic seafloor bottom model. The proposed algorithm is a promising tool for wavefield propagation in heterogeneous media with a nonplanar fluid–solid interface.

Finite difference method for acoustic-elastic coupled equations of seismic waves in the frequency domain

2D Frequency-Domain Finite-Difference Acoustic Wave Modeling Using Optimized Perfectly Matched Layers

Finite Difference Seismic Forward Modeling Method for Fluid–solid Coupled Media with Irregular Seabed Interface

A Multi-Block Finite Difference Method for Seismic Wave Equation in Auxiliary Coordinate System with Irregular Fluid–solid Interface

High-order finite-difference seismic forward modeling method for fluid-solid boundary coupling media

Numerical simulation of elastic wave based on the staggered grid finite difference method

Efficient and Stable Finite Difference Modelling of Acoustic Wave Propagation in Variable-density Media

Optimized staggered-grid finite-difference method in time-space domain based on exact time evolution schemes

An Efficient and high accuracy finite-difference scheme for the acoustic wave equation in 3D heterogeneous media

The Method Of Finite Difference Of High Precision Elastic Wave Equations In The Frequency Domain And Wave-Field Simulation

Finite-difference Time-Domain Scheme for Nonlinear Acoustic Waves in Solids

Multiscale finite-difference method for frequency-domain acoustic wave modeling

An Optimal Finite-Difference Method Based on the Elongated Stencil for 2D Frequency-Domain Acoustic-Wave Modeling

A finite-element algorithm with a perfectly matched layer boundary condition for seismic modelling in a diffusive-viscous medium

Numerical simulation of acoustic waves propagation by finite element method based on optimized matrices

Efficient Modeling of Fractional Laplacian Viscoacoustic Wave Equation with Fractional Finite Difference Method

A high order finite difference method for the elastic wave equation in bounded domains with nonconforming interfaces

An Improved 25‐point Finite‐difference Scheme for Frequency‐domain Elastic Wave Modelling

3D Seismic-Wave Modeling with a Topographic Fluid–Solid Interface at the Sea Bottom by the Curvilinear-Grid Finite-Difference Method

Forward modeling of seismic wave fields in complex anisotropic media using finite difference method

Complex-valued Adaptive-Coefficient Finite-Difference Frequency-Domain Method for Wavefield Modeling Based on the Diffusive-Viscous Wave Equation