Forward Modeling of Rayleigh Surface Waves in Two-Phase Media

This forward modeling simulation focuses on Rayleigh surface waves in two-phase media, providing valuable insights into their propagation characteristics. The implementation employs advanced algorithms including PML (Perfectly Matched Layer) boundary conditions to effectively absorb outgoing waves and minimize artificial reflections, creating more realistic subsurface scenarios. Key optimizations include: - Staggered-grid finite-difference method for wave equation discretization - Convolutional Perfectly Matched Layer (CPML) implementation for broadband absorption - Memory-optimized matrix operations for handling complex media parameters The simulation utilizes high-performance computing techniques to significantly accelerate computation, enabling rapid result generation and analysis. Through this modeling approach, researchers can systematically investigate Rayleigh wave propagation patterns across different media configurations, offering crucial applications in seismological studies and geological exploration. The code structure modularizes core components including wavefield initialization, boundary condition handling, and dispersion relationship calculations for enhanced maintainability and scalability.