Geological formations represent complex multiphase systems comprising elastic solid matrices and fluid-saturated pore networks. While Biot’s classical poroelastic theory has been widely adopted for wave propagation modeling, accumulated physical inconsistencies in describing dynamic filtration processes have motivated our numerical implementation of a Symmetric Hyperbolic Thermodynamically Compatible (HTC) model. This approach maintains rigorous physical validity across the complete spectrum of phase compositions in heterogeneous media. Our research focuses on optimizing parallel computation strategies for large-scale 3D wave field simulations in realistic poroelastic environments. We present a comparative analysis of two fundamental parallelization paradigms - distributed memory (MPI) and GPU-accelerated (CUDA) approaches - evaluating their computational efficiency through scalability tests. The results of numerical experiments are presented and discussed.