In actual application, the energy utilization rate of underwater glider directly affects the total voyage range. When underwater glider is used for executing exploration mission for a fixed point, the position that the glider resurfaces should be accurate enough. In this paper, we employ a multi-objective optimization method to determine the control parameters values that can maximize the position accuracy that the glider resurfaces and the energy utilization rate simultaneously. Especially, the optimization of this paper considers the effect of uncertain input errors. The control parameters include the net buoyancy adjustment amount and the movable mass block translation amount. The input errors include the control parameters errors, the motion depth error and the current. Based on the dynamic model of an underwater glider, we propose the calculation model and evaluation flow that are used for analyzing the glider position accuracy and energy utilization rate, considering the effect of uncertain input errors. Besides, a combinatorial experimental design method is proposed to calculate the performance evaluation parameters under different control parameters values. Then the radial basis function neural network is employed to establish the surrogate models of performance evaluation parameters to participate in the optimization calculation, which can improve the optimization efficiency. After optimization calculation based on the non-dominated sorting genetic algorithm II, we obtain a Pareto optimal set consisting of 257 sets of non-dominated solutions. Finally, the selection rule of optimal control parameters values is given, and the optimization results are validated under 3 sets of solutions. This research may be valuable for the improvement of the glider work quality.