1. Sinopec Geophysical Research Institute, Nanjing 211103, China;
2. College of Geo-exploration Science and Technology, Jilin University, Changchun 130026, China

Aiming at solving the problem of low precision and low computation efficiency in conventional simulation algorithm of the complex structure seismic exploration, the authors propose a high precision spectral element simulation based on element by element algorithm. Based on strong and weak forms of wave equation, the authors derivate the spectral element solution process of second-order seismic wave equation, followed by calculating the separate small-scale element structure using element by element algorithm without global matrix, which improves the computational efficiency and reduces storage space consumption. Numerical simulation of a complicated orebody model proves that this method gets high precision, small numerical dispersion and high computation efficiency.