In this paper we explore a unified cosmological framework combining variable generalized Chaplygin gas (VGCG) and gravitationally–induced adiabatic matter creation cosmology within the context of $f(Q)$ gravity. Theoretical and statistical analyses are performed with two functional forms, $f(Q)=n_{1} Q^{k}$ and $f(Q)=n_{0}+n_{1}Q^{k}$ , where $n_{0}$ , $n_{1}$ and $k$ are dimensionless free parameters. Using the latest observational data sets we perform a robust statistical analysis to constrain the model’s parameters. A statistical comparison, including goodness-of-fit with standard $\Lambda $ CDM model is performed. The results show that the Bayesian inference favors the complex models while the model selection criterion, like Akaike information criteria (AIC) and Bayesian information criteria (BIC) favor the simpler $\Lambda $ CDM model. The results also show a smooth transition from deceleration to acceleration around the redshift $z \approx 0.6$ with effective equation of state parameter remains in the quintessence regime. The present value of Hubble parameter align closely with Planck 2018 measurements.