Variational Bayesian inference for spatial negative binomial count data models with unobserved heterogeneity