1. A method to create a fault surface from a three-dimensional seismic data volume, comprising: constructing an initial fault surface in three-dimensions, said surface containing at least two fault sticks, the fault sticks being from the same fault but from different slices of the seismic data volume, each fault stick being defined by at least two fault nodes; and reconstructing the initial fault surface to fit the three-dimensional seismic data, said reconstruction using an iterative evolution of a deformable surface model of the fault surface, said evolution being based on smoothness of the fault surface and a fault-indicating parameter of each location on the fault surface.