Photovoltaic(PV)power is a typical intermittent energy source

and shows appreciably correlative between adjacent PV power stations. For the probabilistic load flow in the system with PV power stations

the correlation of multiple PV powers is important for analyzing the influence of PV power on the power system. This paper proposes a probabilistic load flow method based on the correlation of multiple PV power stations. The output correlation of multiple PV power stations is modeled by the Pair Copula function

which considers the various dependency structures of pairs of PV power stations. Simultaneously

the existing probabilistic load flow method is improved

and the digit interlacing technology is integrated into quasi-Monte Carlo simulation method

which effectively improves the convergence rate. Taking the history data of PV power stations in Texas

USA for example

the effectiveness of the proposed method is verified by the simulation in the IEEE 30-bus and IEEE 118-bus systems.