A method and a system for simulating magnetic resonance echo-planar imaging artifacts. Firstly, for K-space artifacts, K-space data are restored through normal magnetic resonance images, and the K-space data are modified pertinently, and then images with artifacts are reconstructed; for susceptibility artifacts, a susceptibility model is constructed through normal magnetic resonance images, and the magnetic field distribution is reconstructed, and then the images with distortion artifacts are reconstructed. According to the present disclosure, a large number of artifact data sets with different artifact types and artifact degrees can be quickly created through a small number of normal images, thus laying a foundation for the research of identifying artifacts, eliminating or weakening artifacts. A simulation algorithm is designed according to the principle of generation of EPI sequence artifacts, and the obtained images such as stripe artifacts, Moer artifacts, Nyquist artifacts, susceptibility artifacts and the like have good scientificity, accuracy and interpretability.