For the purpose of improving voltage support capacity of receiving-end power grid and reducing the risk of high voltage direct current (HVDC) continuous commutation failure, this paper presents an optimization method for the operation of multi-dynamic reactive power devices. Firstly, considering the key characteristics of voltage recovery speed and extinction angle of HVDC after faults, the evaluation index of voltage resilience for receiving-end power grid is introduced to measure the voltage recovery ability. Then, taking the maximum voltage resilience as the optimization objective, a mathematical optimization model is established. In this model, the reactive power output of multi-dynamic reactive power devices is treated as the optimization variable, while the steady and transient constraints of power grid are included at the same time. In the end, particle swarm optimization (PSO) method is used for calculating in the proposed model. To testify the correctness of the model, Suzhou southern power grid is utilized to simulated. The result shows that it effectively improves the recovery speed of grid voltage after short-circuit fault and reduces the number of HVDC continuous commutation failure.