The residual stress introduced by laser shock peening (LSP) is one of the most important factors in improving metallic fatigue life. The shock wave pressure has considerable influence on residual stress distribution, which is affected by the distribution of laser energy. In this work, a titanium alloy is treated by LSP with flat-top and Gaussian laser beams, and the effects of spatial energy distribution on residual stress are investigated. Firstly, a 3D finite element model (FEM) is developed to predict residual stress with different spatial energy distribution, and the predicted residual stress is validated by experimental data. Secondly, three kinds of pulse energies, 3 J, 4 J and 5 J, are chosen to study the difference of residual stress introduced by flat-top and Gaussian laser beams. Lastly, the effect mechanism of spatial energy distribution on residual stress is revealed.