Numerical Modelling of Wave Velocity Structure of Subduction Zones