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).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 single backscattering direction. 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.