Identifying spatial heterogeneity of coal resource quality in a multilayer coal deposit by multivariate geostatistics