Using population of models to investigate and quantify gas production in a spatially heterogeneous coal seam gas field