A fast and efficient triangulation algorithm is presented in this paper. The scattered data on any 2D shape is triangulated by this algorithm. Triangular mesh optimization is done at the same time with mesh creating. The resulted mesh is fit to Delaunay criteria. The algorithm complexity is about linear to the number of points. Application in oil and geology exploration, the algorithm can deal with large scale scattered data including complicated faults.