Abstract: With the progress of the transformation of rural clean heating and the construction of urban distributed photovoltaics, a large number of rural microgrids have appeared in the distribution networks. However, the strong randomness and volatility of distributed photovoltaics also bring many challenges to the operation of the power systems and make the power grids face such problems as severe power flow fluctuations and voltage violations. Firstly, a decomposition and coordination model is constructed for distribution networks and microgrids considering heat pumps, photovoltaics and energy storage systems , in which the heat pumps are described by the equivalent state space model. Subsequently, an improved Benders decomposition algorithm is proposed to solve the active and reactive power coordination problem of distribution networks and microgrids. Finally, numerical tests are carried out on the D141-M4 system, which has verified the validity of the proposed method in reducing the peak load and alleviating voltage violations on the premise of ensuring the users' comfort.