Abstract: To address the randomness of distributed photovoltaic (PV) output and load during the power supply restoration process in distribution networks, a robust island division method considering the flexible operation mode of distributed PV is proposed. Firstly, network reconstruction technology is incorporated into the constraints of island division in the distribution network, and an island division model that integrates distribution network reconstruction is established. In light of this, a new island division model for distribution networks is established based on the flexible operation mode of distributed PV, considering the black start mode, P-Q operation mode, and P-V operation mode of distributed PV generators. Subsequently, with the randomness of distributed PV output and load taken into account, the probability distribution of its characteristic scenarios is constrained using the 1-norm and ∞-norm. Additionally, a robust partitioning optimization model for distribution network islands is developed, considering the flexible operation mode of distributed PV. Finally, through the comparative analysis of numerical examples, it is confirmed that the proposed method can effectively enhance the load recovery capability and operational stability of the distribution network.