Abstract: In view of the challenges brought by the uncertainty of renewable energy output and load to the economic stability of power system, a robust optimization strategy for the interactive game between multiple virtual power plants (VPP) and distribution networks is proposed. First, in the context of electricity markets and reactive auxiliary service markets, considering the market clearing strategy based on price quota curves and the point-to-point energy trading strategy between VPPs, the operation models of distribution networks and VPPs, including reactive power compensation mechanisms, are established respectively. Subsequently, the characteristic scenario probability uncertainty set for PV output and load is constructed using Kullback-Leibler (KL) divergence, and optimization models for distribution network and VPP operation are established respectively. The column and constraint generation (C&CG) algorithm is employed to address the issue. On this basis, the interactive game framework of multi-VPPs and the distribution network is established, with the electricity price and the on-grid power determined by VPP as interactive variables. Finally, the effectiveness and economic benefits of the proposed strategy are verified through a simulation example on an improved urban power distribution system.