Figure 4 shows the example for unfaulting method. The image (Figure 4a) contains 6 faults, and the fault slip of the second and last faults is relatively large. We first solve a regularized inverse problem based on equation 4 using shaping regularization to obtain the shift vector (Figures 4c and 4d), and then unfault the original image to an image shown in Figure 4e. After unfaulting, the seismic events have been aligned across faults. Then we estimate dip from the new image and carry out predictive painting to get the result shown in Figure 4f. To get the predictive painting result of the original image, we interpolate the result in Figure 4f back to the original coordinates, and get the result in Figure 4g, where the sharp changes caused by fault displacements can be clearly observed. Figure 4h shows the several picked horizons overlaid on the original image. The consistency between the yellow curves and the true horizons verifies the effectiveness of the unfaulting method.