Similar to fault zone replacement, we can also undo faulting in a seismic image to align seismic reflectors across faults. Following the theory part of the unfaulting method in Wu et al. (2016), we solve a similar regularized linear equation to get the shift vector where and coordinates (1, 2 and 3 are indices representing different spatial dimensions). Shift vector means that we can unfault the image if the sample at moves to .

If is a left point of a fault, we can use local similarity to estimate its fault slip vector , which can provide us with a simple linear equation:

This is the linear equation only applying to those samples alongside faults. For other samples away from faults, we expect unfaulting shifts to vary slowly and continuously along structural directions. This constraint can be added to the inverse problem by using a gradient operator as a regularization term (Wu et al., 2016) or by using a structure-oriented smoothing operator as a shaping operator in the framework of shaping regularization (Fomel, 2007b). We choose the latter, and solve for the shift vector using the method previously utilized by Xue et al. (2016).

After we get the shift vector , we unfault the original image to a fault-free image by doing an inverse interpolation:

Then we can estimate dip from the new image and perform predictive painting to get a painting result , which can be converted to the painting result of the original image by doing a forward interpolation:

2019-05-06