1. A computational method for restoring a seismic fault for data representing a seismic formation, comprising the steps of: receiving original seismic data including a formation volume containing a fault surface; identifying a fault surface in the seismic data, wherein the fault surface intersects and thereby divides a formation volume into at least two sections; identifying regions in the formation volume intersecting the fault surface; decomposing the formation volume into finite elements, wherein the finite elements comprise a simple geometric shape; applying an integral equation to the finite elements, wherein elastic mechanical properties are assumed, wherein the integral equation is derived from the virtual work principle; and evaluating the integral equation, wherein the integrals are reduced to a plurality of linear equations parameterized with displacement and traction values at each vertex of the finite elements.