1. A computer-implemented method for fault removal of one or more faults in a subsurface geological model in order to populate the model with desired information, comprising: (a) ordering the one or more faults in reverse chronological order, and selecting the first fault; (b) removing the selected fault by iteratively solving, using a computer, an optimal control problem wherein Laplace's equation is solved for an optimal set of rigid body transformations and boundary displacement vectors on the fault's surface; (c) repeating (b) as necessary to remove any remaining faults, one at a time, in reverse chronological order, resulting in a transformation of the geological model from faulted space to continuous space; (d) populating said geological model with faults removed with selected grids or surfaces or physical property values; (e) computing a mapping of the populated geological model back to the faulted space.