Surface approximation with faults: application to geophysical surfaces