A global variance reduction(GVR)method based on the SPN method is proposed.First,the global multi-group cross-sections are obtained by Monte Carlo(MC)global homogenization.Then,the SP3 equation is solved to obtain the global flux distribution.Finally,the global weight windows are approximated by the global flux distribution,and the GVR simulation is performed.This GVR method is implemented as an automatic process in the RMC code.The SP3-coupled GVR method was tested on a modified version of C5G7 benchmark with a thick-ened water shield.The results show that the SP3-coupled GVR method can improve the efficiency of MC criticality calculation.