This paper proposes a novel multi-searcher optimization (MSO) algorithm for the optimal energy dispatch (OED) of combined heat and power-thermal-wind-photovoltaic systems. The available power of wind turbine (WT) units and photovoltaic (PV) units is approximated with the probability density functions of wind speed and solar irradiance, respectively. The chaos theory is used to implement a wide global search, which can effectively avoid a low-quality local optimum for OED. Besides, a double-layer searcher is designed to guarantee fast convergence to a high-quality optimal solution. Finally, three benchmark functions and an energy system with 27 units are used for testing the performance of the MSO compared with nine other frequently used heuristic algorithms. The simulation results demonstrate that the proposed technique not only can solve the highly nonlinear, non-smooth, and non-convex OED problem of an energy system, but can also achieve a superior performance for the convergence speed and the optimum quality.