Jet pump efficiency heavily relies on the geometrical parameters of the pump design and parameter global optimization in the full variable space is still a big challenge. This paper proposed a global optimization method for annular jet pump design combining computational fluid dynamics (CFD) simulation, the Kriging approximate model and experimental data. The suction angle, the flow ratio, the diffusion angle, and the area ratio are selected as the design variables for optimization. The optimal space filling design (OSF) method is used to generate sampling points from the design space of the four design variables. The optimization method solves the constrained optimization problem with a given head ratio by building the functional relationship established by the Kriging model between efficiency and design parameters, which makes the method more applicable. The design result shows that the annular jet pump efficiency is predicted well by the Kriging model; m is a key variable affecting the annular jet pump efficiency. As the area ratio m decreases, the mixing effect at the suction chamber outlet can be improved, but the frictional resistance increases.