A reactive power optimization method for integrated transmission and distribution networks related to a field of operation and control technology of an electric power system is provided. The reactive power optimization method includes: establishing a reactive power optimization model for a transmission and distribution network consisting of a transmission network and a plurality of distribution networks, in which the reactive power optimization model includes an objective function and a plurality of constraints; performing a second order cone relaxation on a non-convex constraint of a plurality of distribution network constraints of the plurality of constraints; and solving the reactive power optimization model by using a generalized Benders decomposition method so as to control each generator in the transmission network and each generator in the plurality of distribution networks.