1. A method for enhancing three dimensional fault attributes within a set of three dimensional seismic data, the method comprising: converting the set of three dimensional seismic data from a spatial domain to a spatial-frequency domain using a Fourier transform, the set of three dimensional seismic data including an array of points identifying the three dimensional fault attributes; convolving the set of three dimensional seismic data with a plurality of three dimensional orientation filters in the spatial-frequency domain, the plurality of three dimensional filters forming an orientation filter array; converting the convolution result back to the three dimensional spatial domain by inverse Fourier transform; calculating an energy for each orientation filter using the converted convolution result; identifying an orientation filter of the orientation filter array having a response with a maximum energy for each point in the seismic data; associating an orientation of the orientation filter having the maximum energy with the corresponding energy and the corresponding point for each point in the seismic data; and enhancing a continuity of the three-dimensional fault attributes in the seismic data using the orientations associated with points in the seismic data.