Abstract: The paper introduces the elastic wave equation of 2D three components first-order staggered-grid stress and velocity, simulates wave field of dilatational source, vertical concentrated force source, shear source in fracture media, analyzes wave propagation characteristics with different azimuth angle. The result shows that fast and slow shear wave will be observed when seismic wave propagates in fracture media. This work will provide a strong basis for further analysis of shear wave splitting characteristics and actual seismic data.