It is well known that the global strength of multi-layered composite structures strongly depends on the quality of the adhesion between its constituent elements. Imperfections along interfaces of adhesion can strongly compromise structure s performance. The characterization of such defects is a very difficult task. The main goal of this study is the development of an analytic-numerical method to simulate the acoustic field resulting from the interaction between ultrasonic waves and imperfect interfaces, helping in selection of parameters for ultra-sonic inspecting methods. The Quasistatic-approximation (QSA), introduced by Thompson in 1982, is applied together with the perturbation method to allow modelling of interfacial localized aws. A solution algorithm for the problem is developed with the aid of the invariant embedding method. It is applicable to solve wave propagation problems in arbitrarily anisotropic layered plates and it is stable for high frequencies. Three simulations of multi-layered plates immersed in acoustic uid are presented as illustration of the application of the developed method. Results of those simulations indicate the frequencies and angles of incidence where the scattering effects, which allow the characterization of localized defects, are more significant.