Vibroacoustography by means of ultrasound is an imaging method that was developed to achieve higher resolutions while avoiding the drawbacks of scattering and stronger attenuation. High frequency waves that show a strongly preferred direction of propagation are sent into the medium. Therefore, we make use of a paraxial approach to arrive at a system of PDEs that involve space dependent parameters. In this talk, we will deal with the modeling and inverse problem for vibroacoustography.