Rayleigh Surface Wave Forward Modeling in Bi-directional Media

This program implements forward modeling of Rayleigh surface waves in bi-directional media. The core functionality simulates seismic wave propagation through subsurface geological formations. The implementation utilizes Perfectly Matched Layer (PML) boundary conditions to effectively mimic real-world wave transmission scenarios, significantly enhancing simulation accuracy by minimizing artificial boundary reflections. The algorithm employs finite-difference time-domain (FDTD) methods to solve elastic wave equations, with optimized matrix operations ensuring high computational efficiency. Key functions include wavefield initialization, velocity-stress relation calculations, and PML absorption coefficient implementations. The program can simulate complex wave propagation patterns within reasonable timeframes, making it suitable for both research and practical applications. If you are interested in seismic wave propagation simulation, we encourage you to download and test this program. We welcome valuable feedback and suggestions to further improve its functionality and performance optimization. The code structure allows for easy modification of medium parameters and boundary conditions to accommodate various geological scenarios.