In this paper, we have proposed a true amplitude solution to the seismic imaging problem. We derive a diagonal scaling approach for the normal operator approximation in curvelet domain, which is based on the fact that curvelets remain approximately invariant under the action of operator. We use curvelets as essential tools for both approximation and inversion.
Our method produces a very good resolution of the conflicting dips reflectors, reproduces true reflector amplitudes and compensates for incomplete illumination in seismic images. Our method is tested with a reverse-time ’wave equation’ migration code simulating the acoustic wave equation on different synthetic models.