Abstract: This paper proposes a robust optimization method for active power and reactive power collaboration in active distribution networks considering wind power uncertainty and energy storage flexibility. Firstly, the active power and reactive power cooperative optimization model of active distribution network considering the adjustable power characteristics of energy storage devices in four quadrants is established. Secondly, the linearized representation of AC power flow equation and energy storage devices constraints are given. Then, by using the box set to describe the uncertainty of wind power output, the active power and reactive power cooperative robust optimization model of active distribution network is established. Through the dual transformation, the main problem and subproblem of two-stage robust optimization in the form of mixed integer linear programming are obtained, and the column constraint generation algorithm is used to solve the interactive and iterative problem. Finally, the IEEE 33 node system is taken as an example for simulation verification. The simulation results show that the proposed method can overcome the influence of high non-convexity and nonlinearity of the model on the convergence accuracy and solving efficiency, and make full use of the reactive power regulation ability of the energy storage devices to realize the robust optimal operation of active power and reactive power coordination in the active distribution network.