A multi-objective design exploration for a three-element airfoil consisting of a slat, a main wing, and a flap was carried out. The lift curve improvement is important to design high-lift system, thus design has to be performed with considered multi-angle. The objective functions considered here are to maximize the lift coefficient at landing and near-stall conditions simultaneously. Genetic algorithm is used as an optimizer. Although it has the advantage of global exploration, its computational cost is expensive. To reduce the computational cost, the Kriging surrogate model, which was constructed based on several sample designs, is introduced. The solution space was explored based on the maximization of expected improvement (EI) value corresponding to objective functions on the Kriging surrogate models. The improvement of Kriging surrogate model and the exploration of the optimum can be advanced at the same time by maximizing EI value. In this study, a total of 90 sample points are evaluated using the Reynolds-averaged Navier-Stokes simulation for the construction of the Kriging surrogate model. Through the present exploration process, several designs were obtained with better performance than the baseline setting in each objective function. To obtain the information of the design space, functional analysis of variance, which is one of the data mining techniques showing the effect of each design variable on the objectives, is applied. Main effects of the design variables are calculated to recognize which design variable has the effect on the objective functions. This result suggests that the gap and the deflection of the flap have a remarkable effect on each objective function and the gap of the slat has an effect at near-stall condition.