Data-driven modeling of heterogeneous viscoelastic biofilms