A novel probabilistic small signal stability analysis (PSSSA) method considering load correlation is proposed in this paper. The superiority Latin hypercube sampling (LHS) technique combined with Monte Carlo simulation (MCS) is utilized to investigate the probabilistic small signal stability of power system in presence of load correlation. LHS helps to reduce the sampling size, meanwhile guarantees the accuracy and robustness of the solutions. The correlation coefficient matrix is adopted to represent the correlations between loads. Simulation results of the two-area, four-machine system prove that the proposed method is an efficient and robust sampling method. Simulation results of the 16-machine, 68-bus test system indicate that load correlation has a significant impact on the probabilistic analysis result of the critical oscillation mode under a certain degree of load uncertainty.