Staggered-Grid Finite-Difference Wavefield Simulation for Two-Phase Media

In staggered-grid finite-difference wavefield simulation for two-phase media, we implemented staggered-grid programs with different orders of accuracy, employing 2nd, 4th, and higher-order spatial derivatives to enhance computational efficiency and numerical precision. The algorithm incorporates a novel approach to handle medium heterogeneity through convolutional perfectly matched layer (CPML) boundary conditions and adaptive material parameter grids. Key functions include velocity-stress staggered grid updating schemes with stress correction for fluid-solid coupling. We utilized specialized processing for non-uniform media that dynamically adjusts finite-difference coefficients based on local material properties, resulting in more accurate simulation outcomes. The simulation results underwent comprehensive analysis and interpretation, including waveform comparisons and dispersion error assessments, further advancing the understanding of wavefield simulation in two-phase porous media with complex heterogeneities.