Multivariate geostatistical estimation of coal quality in a multilayered deposit