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.

Modeling the Propagation of Diffusive-Viscous Waves Using Flux Corrected Transport-Finite Difference Method

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

Finite-difference Modelling in Two-Dimensional Anisotropic Media Using a Flux-Corrected Transport Technique

Modeling Of Rayleigh-Wave Propagation In Three-Dimensional Media With Topography Using Staggered-Grid Finite Difference Scheme

Seismic Wave Propagation in Viscoelastic and Saturated Rock

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

Modeling seismic wave propagation within complex structures

Viscoelastic wave finite-difference modeling in the presence of topography with adaptive free-surface boundary condition

Seismic wave modeling of fluid-saturated fractured porous rock: including fluid pressure diffusion effects of discretely distributed large-scale fractures

Frequency-dependent Q simulation and viscoacoustic reverse time migration based on the fractional Zener model

Finite difference frequency domain method with QR-decomposition-based complex-valued adaptive coefficients for 3D diffusive viscous wave modelling

3D frequency-domain modeling of marine controlled source electromagnetic responses with topography using finite volume method

Two-dimensional Time-Domain Finite-Difference Modeling for Viscoelastic Seismic Wave Propagation

Numerical Investigation of Rayleigh-wave Propagation on Topography Surface

Numerical simulation of cross-hole seismic wave field

A Wave Propagation Model with the Biot and the Fractional Viscoelastic Mechanisms

Finite Difference Modelling for Vertical Seismic Profile in Elastic Media

Modeling seismic wave propagation and amplification in 1D/2D/3D linear and nonlinear unbounded media

Acoustic wave propagation and scattering in visco-acoustic medium: Integral equation representation and De Wolf approximation

Wavefield simulation and analysis with finite difference method in partially saturated double-porosity media

Research on the forward-modeling method of viscoelastic Chapman extension in partially saturated rocks