Aiming at the development of the application of massive and heterogeneous LEO constellations, the influence of typical perturbation forces on the orbit prediction of LEO satellites at different altitudes is analyzed. Firstly, the main perturbation force model of LEO satellites is introduced. Then experiments for analyzing the influence of orbit prediction on LEO satellites at different altitudes are designed. Finally, the influence of typical perturbation forces on the orbit prediction under different prediction durations by numerical analysis is assessed systematically, and the orbit prediction of LEO satellites with different non-spherical gravitational field orders, ocean tide orders and atmospheric density models were compared and calculated. The research results show that different orders of non-spherical gravitational fields should be considered when the altitude of satellites is different, the reasonable gravitational field orders for LEO satellites at 500km, 1000km, 1500km and 2000km under 1-hour forecast time are 140x140, 80x80, 60x60 and 40x40 orders, the reasonable orders of ocean tide perturbation for 1500km and 2000km satellites are 20x20 orders. Under the 24-hour forecast time, atmospheric resistance has a great influence on the precision of LEO satellites at orbit altitude below 1,000km. The forecast deviation of the LEO satellites orbit at 500km reaches 8751.576m. The difference of orbit precision between different atmospheric density models is in the meter level. The orbit deviation by the Standard Atmosphere 1976 and NRLMSISE 2000 atmospheric models respectively reaches 25.494m. Relevant conclusions provide support for the optimization and selection of dynamical models in high precision orbit determination and fast orbit prediction of LEO satellites.
