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
(1, 2 and 3 are indices representing different spatial dimensions).
means that we can unfault the image if the sample at moves to
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: