Constructing high approximation accuracy surrogate model with lower computational cost has great engineering significance. In this paper, using co-Kriging method, an efficient multi-fidelity surrogate model is constructed based on two independent high and low fidelity samples. Co-Kriging method can use a greater quantity of low-fidelity information to enhance the accuracy of a surrogate of the high-fidelity model by modeling the correlation between high and low fidelity model, thus computational cost of building surrogate model can be greatly reduced. A wing-body problem is taken as an example to compare characteristics of co-Kriging multi-fidelity (CKMF) model with traditional Kriging based multi-fidelity (KMF) model. A sampling convergence of the CKMF model and the KMF model is conducted, and an appropriate sampling design is selected through the sampling convergence analysis. The results indicate that CKMF model has higher approximation accuracy with the same high-fidelity samples, and converges at less high-fidelity samples. A wing-body drag reduction optimization design using genetic algorithm is implemented. Satisfying design results are obtained, which validate the feasibility of CKMF model in engineering design. Constructing high approximation accuracy surrogate model with lower computational cost has great engineering significance. In this paper, using co-Kriging method, an efficient multi-fidelity surrogate model is constructed based on two independent high and low fidelity samples. a greater quantity of low-fidelity information to enhance the accuracy of a surrogate of the high-fidelity model by modeling the correlation between high and low fidelity model, thus computational cost of building surrogate model can be greatly reduced. A wing-body problem is A sampling convergence of the CKMF model and the KMF model is conducted, and an appropriate sampling design is selected through the sampling convergence analysis. The results indicate that CKMF model has higher approximation accuracy with the same high-fidelity samples, and converges at less high-fidelity samples. A wing-body drag reduction optimization design using genetic algorithm is implemented. Satisfying design results are obtained, which validate the feasibility of CKMF model in engineering design.