<p style='text-indent:20px;'>We analyze mathematically the acoustic imaging modality using bubbles as contrast agents. These bubbles are modeled by mass densities and bulk moduli enjoying contrasting scales. These contrasting scales allow them to resonate at certain incident frequencies. We consider two types of such contrasts. In the first one, the bubbles are light with small bulk modulus, as compared to the ones of the background, so that they generate the Minnaert resonance (corresponding to a local surface wave). In the second one, the bubbles have moderate mass density but still with small bulk modulus so that they generate a sequence of resonances (corresponding to local body waves). <p style='text-indent:20px;'>We propose to use as measurements the far-fields collected before and after injecting a bubble, set at a given location point in the target domain, generated at a band of incident frequencies and at a fixed <i>single backscattering direction</i>. Then, we scan the target domain with such bubbles and collect the corresponding far-fields. The goal is to reconstruct both the, variable, mass density and bulk modulus of the background in the target region. <p style='text-indent:20px;'>1.We show that, for each fixed used bubble, the contrasted far-fields reach their maximum value at, incident, frequencies close to the Minnaert resonance (or the body-wave resonances depending on the types of bubbles we use). Hence, we can reconstruct this resonance from our data. The explicit dependence of these resonances in terms of the background mass density of the background allows us to recover it, i.e. the mass density, in a straightforward way. <p style='text-indent:20px;'>2.In addition, this measured contrasted far-fields allow us to recover the total field at the location points of the bubbles (i.e. the total field in the absence of the bubbles). A numerical differentiation argument, for instance, allows us to recover the bulk modulus of the targeted region as well.