With the increase in the proportion of multiple renewable energy sources, power electronics equipment and new loads, power systems are gradually evolving towards the integration of multi-energy, multi-network and multi-subject affected by more stochastic excitation with greater intensity. There is a problem of establishing an effective stochastic dynamic model and algorithm under different stochastic excitation intensities. A Milstein-Euler predictor-corrector method for a nonlinear and linearized stochastic dynamic model of a power system is constructed to numerically discretize the models. The optimal threshold model of stochastic excitation intensity for linearizing the nonlinear stochastic dynamic model is proposed to obtain the corresponding linearization threshold condition. The simulation results of one-machine infinite-bus (OMIB) systems show the correctness and rationality of the predictor-corrector method and the linearization threshold condition for the power system stochastic dynamic model. This study provides a reference for stochastic modelling and efficient simulation of power systems with multiple stochastic excitations and has important application value for stability judgment and security evaluation.