Discussion - Short Scale Spatial Variability Of Sulfur In A Coal Seam

Discussion by A. Unal Barbaro et al. (1990) implemented a tedious study aimed at delineating the short-scale spatial variability of sulfur in a coal seam. It is not possible, however, to extend their conclusions to any other coal seam or use their in any other fashion, because the background geology was not presented. Nor were the conclusions accompanied by geological interpretations in the paper. In addition, an unfortunate printing error occurred where the total sulfur (determined by High Temperature Combustion) variogram in Fig.3 was duplicated in Fig. 2 instead of the seam height variogram. A more serious error, however, has been committed in the definition of the variable, seam height. What is defined as seam height by Barbaro et al. (1990) is, in essence, mining height, and is not a regionalized variable that should be studied by geostatistical methods directly. Despite the fact that the variable, seam height, is not defined in the paper explicitly, the following quotation discloses the fallacy: "All roof rock that exceeded 1.8 m (6 ft) from the floor was not taken because the longwall was operated to not mine more than 1.8 in (6 ft) unless the coal seam height exceeded 1.8 m (6 ft)." From this definition, the seam height, and therefore the sampling height, is equal to the actual coal thickness, if the coal thickness is greater than 1.8 m (6 ft). Otherwise, it is equal to the sum of the coal thickness plus the thickness of the roof rock that complements the thickness to 1.8 m (6 ft). In the latter case, the seam height is constant and equal to 1.8 m (6 ft). It is possible to conduct a variogram study on a pool of samples that are realizations of two different variables. But the conclusions derived would not belong to any one of the two variables uniquely and, therefore, do not possess any significance. Geostatistical analysis is irrelevant for the sum of multiple regionalized variables formed by arbitrary selections. In a two-seam setting, for example, the mining height, as defined by the thickness of one seam at one location and the thickness of both seams at another location (due to quality and/or minimum thickness considerations perhaps), should not be used in the calculation of one common variogram. The two seams should be modeled separately. They can then be combined according to the specific purposes of the study. On the other hand, if a constant is added to a regionalized variable (to incorporate dilution perhaps), the variogram of the new variable will not change. Barbaro et al. (1990), surprisingly, does not give any geological interpretation of their results despite the fact that most of them can be explained by the origins of sulfur in a coal seam. The presentation of the results of a geostatistical study with no reference to the geology of the deposit is uninformative. It may also be misleading for the potential users of geostatistics, It is not unusual to find nuggets of pyrite in coal seams. In such cases, pyritic sulfur will probably display a spatial structure for only a very small distance that will appear in the experimental variogram as no spatial correlation. This very well known phenomenon is called the pure nugget effect in geostatistics (Journel and Huijbregts, 1978) and perhaps can explain the lack of correlation found for pyritic sulfur content. The lack of correlation found for the total sulfur content may also be explained in the same way because the total sulfur content is dominated by the pyritic sulfur content in this case study. One should notice, however, that the situation may completely be reversed after cleaning the coal. Not all of the inorganic sulfur should be expected to be in the form of pyrite nuggets in a coal seam. It may also be disseminated in the coal seam and it is expected that it follows a certain spatial structure. However, an existing spatial structure may be masked by including a part of the roof rock, rich in sulfides, into the seam thickness in an arbitrary fashion because areas having a sandstone roof sometimes are known to show a higher sulfur content due to the downward percolation of solutions rich in iron sulfides (Clark, 1979). Plants use sulfur in their growth processes. Much of this sulfur is bound organically during peat accumulation and coal formation (Cecil et al., 1978.) This suggests a spatial structure of some sort for the organic sulfur. However, it is not possible to test this hypothesis because the results obtained by the authors for the organic sulfur content are not given in this paper. For this reason, the conclusion that simple average of the nearby samples would provide the best unbiased estimates is questionable for organic sulfur and is not based on any substantial supporting evidence. It is suspected that no spatial structure was detected and this was due to high sampling and laboratory analysis errors. Before concluding that sulfur variability in the seam at the location of study was random, a more detailed study for the disseminated non-pyritic sulfur should have been conducted, not for the sake of scientific curiosity only, but also due to its utmost importance with regard to coal cleaning and emission control. Pyritic sulfur can be cleaned to a considerable extent, whereas organic sulfur can not, making the emission control strategies highly dependent on the spatial distribution of the organic sulfur (Knudsen, 1981). In the light of these facts, one wonders why Barbaro et al. (1990) did not present the results of their study for the sulfate and organic sulfur content. References Barbaro, R.W.. et al., 1990, "Short-Scale Spatial Variability of Sulfur in a Coal Seam,' Mining Engineering, Vol. 42, No. 11, pp. 1267-1269. Cecil, C.B., et al., 1978, "Geology of Ccontaminants in coal," report prepared for Environmental Protection Agency, North Carolina, 123 pp. Clark, W.J.. 1979, "An interfluve model of the upper Freeport coal Bed in part of western Pennsylvania," unpublished MS thesis, University of South Carolina, 57 pp. Journel, A.G., and [Huijbregls], Ch. J., 1978, Mining Geostatistics, Academic Press, London, 600 pp. Knudsen, H.P., 1981. "Development of a Conditional simulation model of a coal deposit," unpublished PhD dissertation, The University of Arizona, 109 pp.