Access the full text.
Sign up today, get DeepDyve free for 14 days.
Arch. Min. Sci., Vol. 60 (2015), No 3, p. 791806 Electronic version (in color) of this paper is available: http://mining.archives.pl DO0.1515/amsc-2015-0052 JACEK MUCHA*, MONIKA WASILEWSKA-BLASZCZYK* INTERPOLATION AND SAMPLING ERRORS OF THE ASH AND SULPHUR CONTENTS IN SELECTED POLISH BITUMINOUS COAL DEPOSIT (UPPER SILESIAN COAL BASIN USCB) BLDY INTERPOLACJI I OPRÓBOWANIA ZAWARTOCI POPIOLU I SIARKI W WYTYPOWANYCH POLSKICH ZLOACH WGLA KAMIENNEGO (GÓRNOLSKIE ZAGLBIE WGLOWE) The basic sources of information on the parameters characterizing the quality of coal (i.e. its ash and sulphur contents) in the deposits of The Upper Silesian Coal Basin (Poland) are drill core sampling (the first stage of exploration) and channel sampling in mine workings (the second stage of exploration). Boreholes are irregularly spaced but provide relatively uniform coverage over an entire deposit area. Channel samples are taken regularly in mine workings, but only in the developed parts of the deposit. The present study considers selected seams of two mines. The methodology used is based on detailed geostatistical analysis, point kriging procedure and P. Gy's theory of sampling. Its purpose is: · defining and comparing geostatistical models for variability of the ash and sulphur contents for data originating from boreholes and mine workings, · predicting by means of point kriging the values of the parameters and errors of interpolation using data from boreholes at grid points where underground mine workings were later channel-sampled, · assessing the accuracy of interpolation by comparison of predicted values of parameters with real values (found by channel sampling), · evaluating the variances of total secondary sampling error (error of preparation of assay samples) and analytical error introduced by assaying of sulphur and ash, · assessing the contribution of sampling and analytical errors (global estimation error) to the interpolation errors. The authors found that the interpolation errors for ash or sulphur content are very large, with mean relative values of 35%-60%, mainly caused by the considerable natural variability, a significant role of random component of variability, and heterogeneity of spatial distribution of these characteristics. The sampling and analytical errors play a negligible role. Their values are smaller than 11% of interpolation error values. Presenting estimates of the spatial distribution of ash and sulphur contents in coal seams by means of contour maps is unreasonable if they are based on drill core sampling. Keywords: geostatistics, contour maps, sampling errors, interpolation accuracy AGH UNIVERSITY OF SCIENCE AND TECHNOLOGY, FACULTY OF GEOLOGY, GEOPHYSICS AND ENVIRONMENTAL PROTECTION, DEPARTMENT OF ECONOMIC AND MINING GEOLOGY, AL. MICKIEWICZ 30, 30-059 KRAKÓW, POLAND Prawo geologiczne i górnicze a w szczególnoci Rozporzdzenie Ministra rodowiska w sprawie dokumentacji geologicznej zloa kopaliny z dnia 22 grudnia 2011 r. naklada na dokumentatora zloa obowizek ilustrowania zmiennoci jakoci kopaliny za pomoc map. W odniesieniu do zló wgla kamiennego rozmieszczenie wartoci parametrów charakteryzujcych jego jako najczciej wizualizuje si za pomoc map izoliniowych zawartoci siarki, zawartoci popiolu lub wartoci opalowej konstruowanych w oparciu o wyniki ich interpolacji w regularnej sieci interpolacyjnej znacznie zagszczonej w stosunku do sieci rozpoznania zloa. Podstawowymi ródlami informacji o wartociach tych parametrów w warunkach zló Górnolskiego Zaglbia Wglowego s wyniki opróbowania rdzeni z otworów wiertniczych (we wczesnych fazach rozpoznania zloa) i wyniki opróbowania bruzdowego pokladów w wyrobiskach górniczych (w zaawansowanych fazach rozpoznania zloa). Otwory wiertnicze s wykonywane w nieregularnej sieci ale z reguly zapewniaj jednolite pokrycie zloa. Próby bruzdowe w wyrobiskach górniczych pobierane s na ogól ze stalym rozstawem ale charakteryzuj jedynie czci zló rozcitych wyrobiskami. W artykule przestawiono wyniki bada nad wielkoci bldów interpolacji zawartoci siarki i popiolu oraz na ich tle przedstawiono oszacowania bldów opróbowania wykonane dla dwóch pokladów kopalni Murcki (330 i 334/2) oraz jednego pokladu kopalni Brzeszcze (405/1) (Fig. 1 i 2). Metodyczn podstaw bada stanowila statystyczna i geostatystyczna analiza zmiennoci wyników opróbowa, geostatystyczna procedura interpolacji tych parametrów za pomoc zwyczajnego krigingu punktowego oraz dla okrelenia bldów przygotowania próbki do analizy chemicznej elementy teorii opróbowania Gy (Mucha i Wasilewska, 2006). Cele bada byly wielokierunkowe i obejmowaly: · okrelenie i porównanie geostatystycznych modeli zmiennoci zawartoci siarki i popiolu dla danych pochodzcych z z opróbowania rdzeni wiertniczych i wyrobisk górniczych, · prognozowanie za pomoc krigingu zwyczajnego na podstawie danych otworowych wartoci parametrów w punktach zloa opróbowanych w póniejszej fazie rozpoznania w wyrobiskach górniczych, · ocen dokladnoci interpolacji przez porównanie wartoci parametrów prognozowanych na podstawie danych otworowych z danymi z wyrobisk górniczych, które stanowily zbiór referencyjny, · ocen wariancji bldów przygotowania próbek do analizy chemicznej w oparciu o formul Gy (wzór (1)) oraz bldów oznacze zawartoci siarki i popiolu na podstawie oznacze kontrolnych, · ocen udzialu bldów opróbowania i oznacze zawartoci siarki i popiolu w bldach interpolacji. Wyniki analizy statystycznej nie pozwalaj z jednym wyjtkiem (zawarto siarki w pokladzie 405/1 kopalni Brzeszcze) na odrzucenie hipotezy o identycznoci rozkladów prawdopodobiestwa parametrów dla danych z otworów wiertniczych i wyrobisk górniczych z ryzykiem bldu mniejszym od 5% (Fig. 3). Wyniki geostatystycznego modelowania zmiennoci s niejednoznaczne i pokazuj w dwóch przypadkach (zawarto siarki w pokladzie 405/1 kopalni Brzeszcze i zawartoci popiolu w pokladzie 334/2) drastyczne rónice w postaci modeli dla danych otworowych i górniczych. (Fig. 4). W pozostalych przypadkach modele dla obu rodzajów danych wykazuj zadowalajc zgodno. Interpolacja wartoci parametrów na podstawie danych otworowych za pomoc procedury krigingu punktowego charakteryzuje si nisk dokladnoci ze rednimi wzgldnymi (absolutnymi) bldami rzdu 35%-60% (Tab. 1). Przyczyn niskiej dokladnoci interpolacji naley upatrywa przede wszystkim w do duej zmiennoci zawartoci siarki i popiolu (ze wspólczynnikami zmiennoci od 38% do 77%) (Fig. 3) oraz w znaczcym udzialem losowego skladnika ich zmiennoci reprezentowanego na modelach semiwariogramów przez silnie zaznaczon warto parametru modeli C0 (nugget effect) (Fig. 4). Znajduje to potwierdzenie w mapach izoliniowych sporzdzonych metod krigingu dla zawartoci siarki (Fig. 5) i zawartoci popiolu w wglu (Fig. 6) dla pokladu 330 w kopalni Murcki gdzie rezultaty interpolacji dokonanej metod krigingu punktowego s najdokladniejsze. Przebiegi izolinii dla danych otworowych i z wyrobisk górniczych wykazuj zauwaalne rónice, podobnie jak i prognozowane bldy interpolacji. Wysokie, rzeczywiste rednie bldy interpolacji zawartoci siarki (35%) i zawartoci popiolu (47%) dowodz, e moliwoci wiarygodnej predykcji wartoci tych parametrów w punktach zloa na podstawie danych otworowych s dalece niewystarczajce dla planowania eksploatacji uredniajcej. Bldy opróbowania okrelone w oparciu o teoretyczn formul Gy (wzór 1, Fig. 7) oraz bldy oznacze zawartoci siarki i popiolu ocenione na podstawie oznacze kontrolnych maja marginalne znaczenie. Wariancja tych bldów stanowi zaledwie od 3% do 11% redniego kwadratowego bldu interpolacji (Fig. 8). Tak wic nie maj one zauwaalnego wplywu na dokladno interpolacji. Wysokie wartoci bldów interpolacji zawartoci siarki i popiolu czyni bezzasadnym ilustrowanie rozmieszczenia wartoci tych parametrów w pokladach za pomoc map izoliniowych skonstruowanych na podstawie rzadko rozmieszczonych danych otworowych. Slowa kluczowe: geostatystyka, mapy izoliniowe, bldy opróbowania, dokladno interpolacji 1. Introduction The most common way of presenting the spatial distribution of various quality characteristics for coal seams in the Upper Silesian Coal Basin (USCB) in Poland includes contour maps. The maps are based upon interpolation of the parameters determined in the nodes of a square grid. In initial stages of exploration, contour maps of deposits are prepared on the basis of drilling results. They allow prediction of such coal quality characteristics as sulphur and ash contents in the parts of deposits being developed. The reliability of these maps is controlled by the accuracy (or, conversely, by the errors) of interpolation, and depends mainly upon the natural variability of the deposit and the density of the borehole grid, but is also affected by errors of sampling and assaying. Usually, no numerical estimate of the reliability of interpolation is presented. 2. Aim and scope of investigations Investigations were carried out to establish the reliability of contour maps showing the amounts of sulphur and ash in selected coal seams of the USCB. They particularly deal with: · assessment of accuracy of the interpolation made on the basis of drill core sampling, · verification of the variances of interpolation errors predicted using the kriging method by comparing them with the variances of real interpolation errors, · quantitative, theoretical evaluation of the variance of errors made during the comminution and preparation for assaying of samples taken from coal seams (total secondary sampling error) and analytical (assaying) error, · assessment of the contribution of sampling and analytical errors to interpolation errors. 3. Materials The investigations were based on the results of assaying for total sulphur (Str ) and ash (Ar ) in samples collected from drill cores and underground workings in coal seams No's 330 and 334/2 in the Murcki Mine as well as assaying for total sulphur (Str ) in samples collected from underground development and exploitation workings in coal seam No. 405/1 in the Brzeszcze Mine (Fig. 1). It must be emphasized that pyrite constitutes main source of sulphur in coal seams. These seams and mines were selected because there is sufficient amount of data for reliable statistical and geostatistical operations to be made. In Polish mines of bituminous coal, systematic seam sampling for assessing coal quality parameters is a rare practice and the sampling grid is usually irregular. Fig. 2 presents the locations of samples collected from drill cores and mine workings, whose sulphur and ash assays were used by the authors. Fig. 1. Localization of the Murcki and Brzeszcze mines (in red) within the Upper Silesian Coal Basin (according to PIG-PIB Warszawa), http://geoportal.pgi.gov.pl/css/zrozum_z/images/jkc/Rys_JKCw1.jpg) 4. Methods The sets of the numerical data obtained were divided into two groups. The first group is data that served as a basis for interpolation (results of sampling of drill cores in the Murcki Mine and of development workings in the Brzeszcze Mine), while the other group is data used in assessment of the accuracy of interpolation (results of sampling of exploitation workings in both mines). The two groups of data were statistically analysed, including preparation of histograms and calculation of basic variability measures, i.e. arithmetic mean, variance and coefficient of variation (Fig. 3). The authors have compared empirical distributions using the KolmogorovSmirnov test and compared medians of distributions by applying the Mann-Whitney test (Fig. 3). The structure of variability of the sulphur and ash contents was studied using Matheron geostatistics (1963) by means of sample semivariograms that express the dependence between Fig. 2. Locations of the boreholes and underground sampling points. 1 sampling points in development workings (Brzeszcze Mine), 2 sampling points in exploitation workings (Murcki and Brzeszcze mines), 3 boreholes sampled, 4 area of mining exploration the mean squared differences of the parameter values and the average distance between sampling points. The omnidirectional semivariograms calculated for each of the seams, separately for the data obtained from drill core sampling and those from underground sampling, along with theoretical fitting models are presented in Fig. 4. The reliability of the models was checked by cross-validation (Armstrong, 1998). The values of the deposit parameters analysed were interpolated using the geostatistical point kriging method (Isaaks, Srivastava 1989). Its advantage stems from the fact that it estimates the value at a given point as a weighted average and provides an estimate of the magnitude of the interpolation error, considering in both cases the model of the semivariogram assumed, the positions of the sampling points for which the value of the data is known, and the distances of those points from the point at which the value is estimated. The sulphur and ash contents were interpolated for each of the seam sampling points in underground exploitation workings using data from the eight closest boreholes (Murcki Mine) or the eight nearest points of sampling in development workings (Brzeszcze Mine). The differences between the values interpolated and the values later measured in the workings were used as a measure of interpolation accuracy. Assuming a high accuracy of sampling and assaying, these differences may be regarded as real interpolation errors. The average sizes of the real interpolation errors were compared with estimates of the reliability predicted using point kriging. The accuracy of interpolations is displayed on contour maps of the data analysed. These contour maps also show predicted and real errors of interpolation (Fig. 5, Fig. 6). The magnitude of the variance of total secondary sampling error was estimated using the theoretical formula of Gy (Gy 1983; Petersen, Minkkinen, Esbensen 2005), considering the flow sheet for analytical sample preparation (Fig. 7). The preparation of the assay portion involves two stages of sample size reduction with simultaneous reduction of the mass from 25 kg to 1 g (ash assaying) or 0.25 g (sulphur assaying). The relative variance of the fundamental error (FE) that appears during sample mass reduction was evaluated with a simplified formula (Mucha, Szuwarzyski, 2004). 2 FE where: 1 1 3 m f g li di Qi Q (1) N -- number of stages of size reduction (N = 2), Qi -- mass of the sample at the "i-th" stage of size reduction [g] (Q0 = 25,000g, Q1 = 150g, Q2 = 1 g for ash or 0.25g for sulphur), m -- mineral constitution factor [g/cm3], for sulphur: m = 241.5 g/cm3 (mean pyrite content = 2% i.e. mean sulphur content 1%, specific densities of pyrite = 5.0 Mg/m3 and coal = 1.4 Mg/m3) and for ash: m = 15.6 g/cm3 (mean ash content = 15%, specific density of ash = 3.0 Mg/m3), f -- grain shape factor (f = 0.5), g -- granulometric distribution factor (g = 0.5), di -- diameter of the largest grains of the material at the "i-th" stage of size reduction [cm] (d1 = 0.6 cm, d2 = 0.02 cm), li -- factor of liberation of the mineral of interest; li ( dl 1.5 ) (François-Bongarçon, di dl for dl di, Gy 2000), here it was used more conservative formula li di where dl liberation size for pyrite = 0.002 cm hence for sulphur l1 = 0.058 and l2 = 0.316; the liberation size for ash is not known, therefore the safe value l = 1 was assumed. The maximum relative variance of the total sampling errors [TSE], comprising selection and preparation errors, was obtained as a doubled variance of fundamental error [FE] what provides a safe estimate (Whateley & Scott, 2006). The variance of the chemical assaying error [AE] was assumed to be the same as for control assays for sulphur and ash carried out on two coal samples (internal standards) in 29 chemical laboratories of different coal mines within the USCB. Variance of the global estimation error including combined total sampling and assaying errors was evaluated from the relation: 2 [GE ] 2 [TSE ] 2 [AE ] . The variance of the total secondary sampling error and assaying error (global estimation error, according to Gy 1983, Whateley & Scott, 2006), relative to the mean squared error of interpolation are shown in Fig. 8. 5. Results 5.1. Statistical comparison of the data sets from drilling and mining explorations The coefficients of variation of sulphur content range from 40% to 61%. Those of ash content cover the wider range of 38% to 77%. Both ranges indicate generally high variability of the parameters analysed. The distributions of sulphur and ash contents show strong positive asymmetry (Fig. 3). Standardized coefficients of asymmetry and of kurtosis for these distribu- Fig. 3. Histograms of ash and sulphur contents for the data sets used in interpolation (upper section) and in estimation of interpolation accuracy (lower section). N number of data, z arithmetic mean, S2 variance, v coefficient of variation, K-S KolmogorovSmirnov test to compare the distribution of the two samples, M-W Mann-Whitney test to compare the medians of the two samples, [+] there is not a statistically significant difference between the two distributions or the medians at the 95% confidence level, [-] there is a statistically significant difference between the two distributions or the medians at the 95% confidence level, I development workings, II exploitation workings tions, except for the ash content in the seam No. 330 of the Murcki Mine, are higher than 2. These values indicate significant departures of the empirical distributions from normality. This means that the conditions for applying Snedecor's F-test for comparison of variances and Student's t-test for comparison of mean values have not been fulfilled for these data sets. Therefore, the estimation of statistical similarity of both types of data sets have been made by comparing the empirical distributions of the parameters analysed (Kolmogorov-Smirnov test) and their medians (Mann-Whitney test) (Fig. 3). These tests do not provide a basis for rejecting (at the 95% confidence level) the hypothesis that the distributions or medians are identical for the sets of sulphur and ash contents in drill cores and underground workings. Statistically significant differences in distributions of the parameters analysed were only noted for the sulphur content of samples taken from underground development and exploitation workings in the Brzeszcze Mine. 5.2. Geostatistical models of parameter variability Omnidirectional (isotropic) sample semivariograms were approximated using spherical models and, in single cases, the Gauss model and the linear model (Journel & Huijbregts, 1978) (Fig. 4). Cross-validation indicates that the models of semivariograms were properly selected (mean standardized differences between predicted and real values of parameters range from -0.08 to 0.06, while their standard deviations range from 0.98 to 1.05). The fitted semivariograms show a substantial variation in the random component (termed "nugget effect" by Royle 1991) from 21% to 75% for the sulphur content and from 12% to 82% for the ash content. Geostatistical models fitted using the data from exploration drillings (wider range of calculations) and from underground workings (shorter range of calculations) are consistent (i.e. they indicate similar semivariogram at intermediate ranges) for the sulphur content in seams 330 and 334/2 in the Murcki Mine and for the ash content in seam 330 of the same mine. In the remaining two cases they reveal a different structure of variability of both sulphur and ash for the two ranges of observations, possibly caused by differences between types and areas of exploration (Fig. 4). 5.3. Accuracy of interpolation of the sulphur and ash contents Assessment of interpolation precision was based on the data obtained from sampling of underground workings (reference set). It reveals generally poor reliability of interpolation, expressed as a poor ability to predict the results of future channel sampling. The accuracy of interpolation has been quantitatively characterized by some measures shown in Table 1. On the basis of those measures, the following conclusions have been drawn. · The differences between the predictions of geostatistical models fitted using data from drilling exploration and future measurements in underground workings are not consistent with the standard errors predicted by those geostatistical models. The values of cross validation statistics ( /s) deviate significantly from the ideal value ( = 0/s = 1). This was particularly marked for ash content, but the deviation was sufficiently large that predictions of sulphur content must be regarded as untrustworthy. · Interpolated values of sulphur and ash are often quite biased, with systematic errors ( R), reaching 18% (see Table 1) in the cases analysed here. Mean relative errors, i.e. absolute ( AR) and squared ( IR), of the parameter estimation reach high values from 35% to 60% for sampling in underground workings. Fig. 4. Semivariograms and geostatistical models of the parameters analysed for the data from drilling and mining exploration (notation of model formulae according to Deutsch and Journel 1992). 1 semivariogram for the data from drilling exploration, 2 semivariogram for the data from mining exploration (exploitation workings), 3 semivariogram for the data from mining exploration (development workings), 4 semivariogram model for the data from mining exploration, 5 semivariogram model for the data from drilling exploration, C0 variance of the local variability (nugget variance), C variance of the spatial variability, a the range of the semivariogram, NE nugget effect · Estimates of the variance of interpolation errors predicted by point kriging procedure ( KR) are inaccurate by amounts ranging from 34 to 124%. They can drastically differ from the average variance of real interpolation errors. · Coefficients of linear correlation r between the data and the predictions made using point kriging range from -0.13 to 0.41, which is regarded as surprisingly small. · Assuming the normal distribution of errors, the absolute differences between the data and the values predicted using the kriging method should be smaller than the standard error predicted by kriging 68% of the time, and smaller than twice the kriging standard error 95% of the time. However, the numbers in Table 1 (N) reveal, except for sulphur in the seam No. 330 of the Murcki Mine, that the differences calculated are not in accordance with this expectation. TABLE 1 Summary of the quality of interpolation for total sulphur content (Str ) and ash (Ar ) contents by point kriging Mine seam No. parameter 1 Crossvalidation /s (eq. 2-4) 2 R [%] (eq. 5) 3 AR [%] [Lq-Uq] (eq. 6) 4 IR [%] (eq. 7) 5 KR [%] [Lq-Uq] 6 N(%) < 1K N(%) < 2K 8 r(z zK*) 7 Brzeszcze, 405/1 Str 0.25/0.71 58.5 [19.3-83.6] 84.6 [49.6-110.0] -0.13 86.7/100 Murcki, 330 Str Murcki, 334/2 Str Murcki, 330 Ar Murcki, 334/2 Ar -0.05/0.94 -0.29/0.96 0.50/1.27 0.02/0.63 -2.4 -16.9 14.0 1.6 35.3 [11.8-45.0] 49.0 [17.1-55.7] 47.1 [12.5-47.4] 59.7 [20.5-82.6] 37.6 58.9 38.1 58.4 47.1 [32.4-54.2] 77.7 [48.1-91.1] 33.9 [21.7-38.5] 124.2 [70.0-166.9] 0.32 -0.01 0.41 0.12 69.0/96.4 74.1/98.1 50.6/87.1 88.0/100 /s measures of the cross-validation procedure: mean of the standardized estimation errors/standard deviation of the standardized estimation errors, R mean relative interpolation error, AR mean relative absolute interpola tion error, IR mean squared relative interpolation error, KR mean relative kriging standard error (predicted * interpolation error), r(z zK* ) coefficient of linear correlation between the predicted (zK ) and the observed (z) values of the parameter at sampling points, N(%) < 1K, N(%) < 2K percentage of interpolation errors lower than: kriging standard error / twice kriging standard error The statistical parameters (Tab. 1) were calculated using following equations: i 1 N N ziK zi iK (2) 1 N (3) s 1 N (i )2 (4) R AR * (ziK zi ) z 1 N 1 N (5) N * ziK zi zi ziK zi (6) IR (7) where: * ziK -- the predicted with point kriging values of the parameter at sampling point "i", zi -- the observed values of the parameter at sampling point "i". The investigations of interpolation accuracy described above have also been visualized using contour maps of the sulphur content (Fig. 5) and the ash content (Fig. 6) for seam No. 330 of the Murcki Mine, for which the results of interpolation are the most accurate. * Fig. 5. Kriging contour maps: A sulphur content (zK ), B predicted relative interpolation errors ( KR), ); D histogram of relative interpolation errors ( ); C real relative absolute interpolation errors ( AR AR seam No. 330, Murcki Mine. 1 boreholes, 2 underground sampling points, 3 contours of sulphur contents and relative kriging errors based on results of channel sampling in mine workings, 4 contours of sulphur contents and relative kriging errors based on drilling data, 5 contours of absolute relative interpolation errors, N number of data, mean value of a parameter, Lq lower quartile, Uq upper quartile z * Fig. 6. Kriging contour maps: A ash content (zK ), B predicted relative interpolation errors ( KR), C real relative absolute interpolation errors ( AR); D histogram of relative absolute interpolation errors ( AR); seam No. 330, Murcki Mine. 1 boreholes, 2 underground sampling points, 3 contours of ash contents and relative kriging errors based onr results of channel sampling in mine workings, 4 contours of ash contents and relative kriging errors based on drilling data, 5 contours of absolute relative interpolation errors, N number of data, mean value of a parameter, Lq lower quartile, Uq upper quartile z The contour maps based on point kriging show spatial distributions rather different for the data obtained from drilling and the data from underground workings (Fig. 5A and 6A). Comparison of kriging using the data from boreholes and the data from underground workings (Fig. 5B and 6B) shows substantial discrepancy between the locations of contour lines. Accuracy of interpolation has been presented as a contour map of relative absolute interpolation errors within the area of mining exploration (Fig. 5C and 6C) and histograms of these errors (Fig. 5D and 6D). Large interpolation errors for the sulphur content and extremely large errors for the ash content suggest that point prediction of these quantities is insufficiently accurate to be useful in mine planning. 5.4. Total sampling and assaying errors The relative total sampling errors [TSE] evaluated from P.Gys formula are identical for both characteristics of quality of bituminous coal and amount to 10.6% (Fig. 7). The relative standard errors of assaying [AE] estimated from control assays on two coal samples used as internal standards are distinctly lower, amounting to 1.9% for the ash content and 6.1% for the sulphur content. Fig. 7. Flow sheet of preparation of analytical samples. FE, TSE fundamental and total sampling errors, AE analytical error, GE global estimation error, Qi mass of the sample at the i-th stage of comminution, di maximum diameter of grains at the i-th stage of comminution The global estimation errors [GE] can be evaluated as 10.8% for ash content and 12.2% for sulphur content. The assessed accuracy of sampling and analytical procedures can be recognized as acceptable in mining practice. 5.5. Impact of total sampling and analytical errors on the accuracy of interpolation The total variance of sampling and analytical errors can be added to the variance of parameters in the small scale of observation (microvariability) to determine the value of the nugget variance. In equations of semivariogram models, the nugget variance is the free term C0 (Fig. 4). For the data analysed, the variances of sampling and analytical errors contribute in 3-35% to the nugget variances established for the data from drilling exploration. Therefore, they should not measurably affect interpolation errors. This conclusion is corroborated by the diagrams in Fig. 8, showing mean contributions of the combined variance of total sampling and analytical errors (global estimation error) to the mean squared interpolation errors. Their proportionate contribution is 3 to 11%. Hence total sampling and analytical errors can be treated as having only a marginal impact on prediction errors. Fig. 8. Contribution (%) of variance of global estimation error to the mean squared error of interpolation for sulphur (Str ) and ash (Ar ) contents. ( IR)2 mean squared error of interpolation, 2[GE] variance of global estimation error 6. Summary and conclusions The reliability of the predicted spatial distributions of the sulphur and ash contents represented as the contour maps for the selected seams of the Upper Silesian Coal Basin (Poland) is very low. Prediction of these quantities using point kriging is generally of low accuracy, particularly for the ash content. This prediction is burdened with a large mean relative absolute interpolation error, ranging from 35 to 60%. Quite often, the interpolated values are biased. Additionally, the predicted variance of interpolation errors can substantially differ from the observed variance of interpolation errors. It must be stressed that the authors have obtained similar results using nonlinear geostatistical methods such as the indicator and probability kriging. The variance of sample preparation errors (i.e. errors of comminution and preparation for assaying of samples collected in underground workings or from drill cores) and the variance of errors of assaying are not major problems for the accuracy of interpolation. The variances of these types of errors do not exceed 11% of the variance of interpolation errors. This means that improving the accuracy of preparation of analytical samples and of assaying would not significantly improve the quality of interpolation. The low interpolation accuracy can be explained, first of all, by enough high variability of the sulphur and ash contents in coal seams and considerable contribution of the random component of variability. An additional factor lowering the accuracy of interpolation may include heterogeneous spatial distribution of the sulphur and ash contents in the seams, statistically resulting in significant differences of empirical distributions or medians between various parts of the seams. Considering the high variability of sulphur and ash in other seams and mines of the USCB, it seems very likely that the authors conclusion of low interpolation accuracy can be projected onto the whole area of the USCB. Summarising, the contour maps of the spatial distribution of sulphur and ash contents in coal seams predicted in the early stages of exploration on the basis of drilling data have no predictive value at all. These maps do not provide the sound and reliable prediction of the sulphur and ash values in the parts of a deposit that will be mined in the future. The research project was financed from the AGH University of Science and Technology grant No. 22.214.171.1240.
Archives of Mining Sciences – de Gruyter
Published: Sep 1, 2015
Access the full text.
Sign up today, get DeepDyve free for 14 days.