Based on the parameters of the HL-2A experiment, the effect of energetic particles (EPs) on non-resonant high-order harmonics energetic particle modes (EPMs) with qmin>1 is investigated in the present work. Hybrid kinetic-magnetohydrodynamic nonlinear code M3D-K is performed to simulate the linear properties and the nonlinear evolution of the non-resonant EPM during neutral beam injection (NBI). To deeply understand the physical mechanism of interaction resonant between energetic-ions and non-resonant EPM, this work compares the effects of passing energetic particles and trapped energetic particles on the non-resonant EPM instabilities. It is numerically identified that EPs’ effects on high n harmonics (m/n = 2/2, 3/3, 4/4) instability are more obvious than the m/n = 1/1 mode. Furthermore, the effects of energetic particles injection energy, the minimum safety factor qmin , toroidal rotation and beam ion distribution on the features of high n harmonics are also investigated specifically. Toroidal rotation is found to suppress high n harmonics, which is more obvious for the modes driven by trapped particles. Nonlinear simulation results show that these non-resonant high n harmonics can induce larger energetic ion transport, which may affect the plasma confinement performance.