In this paper, we propose a true amplitude solution to the seismic imaging problem. We derive a diagonal scaling approach in curvelet domainfor the normal operator approximation, which is based on the fact that curvelets remain approxiately invariant under the action of operator. We use curvelets as essential tools for both approximation and inversion. Our method produces a 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 BP 2004 synthetic dataset.

Keywords