Viscoacoustic Wave Numerical Simulation Method in an Auxiliary System for Irregular Surface Topography based on a Hybrid Spatial Partial Derivative Scheme