1. A method for identifying faults in seismic data, comprising the steps of: obtaining said seismic data; calculating a coherence cube of said seismic data; performing the following steps, using at least one processing device, for a plurality of two-dimensional seismic sections of the coherence cube: (i) applying a threshold to the coherence cube to obtain a binary image representation comprising one or more continuities; (ii) identifying edges of continuity areas in the binary image representation to identify changes in said one or more continuities as fault point candidates; (iii) identifying fault points in the obtained seismic data based on a fault confidence value above a specific threshold, indicating a likelihood that a given point is part of a fault; (iv) creating one or more fault segments from the identified fault points; and (v) joining one or more fault segments into one or more fault lines using one or more of geological and geometrical constraints; and generating, using the at least one processing device, one or more three-dimensional fault surfaces from said one or more fault lines in said plurality of two-dimensional seismic sections.