The damage of wind turbines suffered from lightning strikes has been a key issue for the safe and reliable operation of wind farms. Accurate determination of the annual lightning flash number to a wind turbine is essential for designing proper lightning protection measures. The interaction of downward and multiple upward leaders (MULs) is studied in this paper, which also considers the stochastic nature and branched behaviour of the lightning attachment phenomenon. Firstly, with an improved stochastic lightning model, the relationship among the striking distance, the height of wind turbines and the return stroke current is established. Moreover, a modified method for predicting the annual lightning flash number strikes to a wind turbine is proposed. The simulation results show that the striking distance and the collection area not only depend on the return stroke current, but also on the height and blade angle. Besides, the comparations between the calculation results and field statistics indicate that the conventional electric geometry model is not satisfied with the need of lightning protection for wind turbines. Note that there is only a difference of 4% between the modified method and the observation value. These quantitative researches provide the guidance of calculation results for the optimization of the blade lightning protection system (LPS) design.