A novel approach to identify the crystal orientation of turbine blades with anisotropy materials is proposed. Based on enhanced mode basis, with the main advantages of its efficiency,accuracy and general applicability, the blade vibration mode of each order is linearly constructed by several specified mode shapes, which are obtained from the considered turbine blade with specified crystal orientations correspondingly. Then, a surrogate model based on Kriging method is introduced for constructing the condensed perturbed matrix of stiffness in order to improve the efficiency even further. The constructed surrogate model allows to perform the modal analysis of turbine blades with arbitrary crystal orientations in higher efficiency, due to the fact that the elements of condensed perturbed matrix of stiffness are considered in construction of the surrogate model rather than concerning the perturbation of all the elements of the initial stiffness matrix for the blade. Genetic algorithm is finally employed to optimize the defined fitness functions in order to identify the crystal orientation angles of turbine blades. Several corresponding examples demonstrated the accuracy, efficiency and general applicability of the proposed method.