Current multiple suppression methods based on virtual events can effectively attenuate internal multiples in seismic data processing. However, challenges remain, including incomplete prediction of internal multiples and amplitude-phase mismatches between predicted and actual multiples, leading to residual artifacts after subtraction. To address these limitations, this study proposes a joint suppression strategy integrating prestack virtual events with high-precision Radon transform. The method first applies prestack virtual event construction for preliminary multiple suppression, followed by high-precision Radon transform to eliminate residual multiples, thereby minimizing subtraction-induced artifacts. Tests on synthetic models and field datasets demonstrate that compared to standalone virtual event methods, this integrated approach achieves more complete multiple removal, significantly enhancing signal-to-noise ratio and imaging quality. The method not only advances theoretical frameworks for multiple suppression but also provides a practical solution for processing seismic data in complex structural areas.