With distributed generation (DG) being continuously connected into distribution networks, the stochastic and fluctuating nature of its power generation brings ever more problems than before, such as increasing operating costs and frequent voltage violations. However, existing robust scheduling methods of flexible resources tend to make rather conservative decisions, resulting in high operation costs. In view of this, a two-stage robust optimal scheduling method for flexible distribution networks is proposed in this paper, based on the pairwise convex hull (PWCH) uncertainty set. A two-stage robust scheduling model is first formulated considering coordination among on-load tap changers, energy storage systems and flexible distribution switches. In the first stage, the temporal correlated OLTCs and energy storage systems are globally scheduled using day-ahead forecasted DG outputs. In the second stage, FDSs are scheduled in real time in each time period based on the first-stage decisions and accurate short-term forecasted DG outputs. The spatial correlation and uncertainties of the outputs of multiple DGs are modeled based on the PWCH, such that the decision conservativeness can be reduced by cutting regions in the box with low probability of occurrence. The improved column-and-constraint generation algorithm is then used to solve the robust optimization model. Through alternating iterations of auxiliary variables and dual variables, the nonconvex bilinear terms induced by the PWCH are eliminated, and the subproblem is significantly accelerated. Test results on the 33-bus distribution system and a realistic 104-bus distribution system validate that the proposed PWCH-based method can obtain much less conservative scheduling schemes than using the box uncertainty set.