The approaches presently disclosed provide for fault-interpretation in a seismic volume with computer assistance, allowing automatic or semi-automatic determination of a fault surface and associated displacement across the fault. The present fault interpretation approach uses pattern matching algorithms and does not require prior interpretation of the stratigraphic horizons. In certain implementations the fault interpretation approach estimates the 3D fault surface as part of a joint fault surface location and displacement optimization process.