We propose a true amplitude solution to the seismic imaging problem. We derive a diagonal scaling approach for the normal operator approximation in the curvelet domain. This is based on the theorem which states that curvelets remain approximately invariant under the action of the normal operator. We use curvelets as essential tools for both approximation and inversion. We also exploit the theorem which states that curvelet-domain approximation should be smooth in phase space by enforcing smoothness of curvelet coefficients in angle and space domain.We analyze our method using a reverse time migration-demigration code, simulating the acoustic wave equation on different synthetic models. Our method produces a good resolution with reflecting dips and reproduces true amplitude reflectors and compensates for incomplete illumination in seismic images.

Keywords