A good resolution of the subsurface image helps to reduce uncertainties, which will guide a good interpretation of the reservoir rocks, spotting possibly an oil & gas accumulation. In this context Least-Squares Reverse Time Migration ( LSRTM ) plays an important role in attenuating migration artifacts, helping to focus the reflectors and improve spatial resolution. In this work it has been used the LSRTM acoustic case using a state-of-art finite-difference (FD) code-generation, Devito, which allows to solve the wave equation, using symbolic expressions. The LSRTM using Devito has been performed in two 2D models, which show its particularities and difficulties. In one of the examples, a synthetic velocity model from the B úzios field was used to validate the LSRTM. Also, some piece of codes has been provided to show the way Devito handles with the wave equation and how easily it is coded.