Next Article in Journal
Paleoenvironment of Mesoproterozoic Gaoyuzhuang and Wumishan Formations, North China: New Insights from Geochemistry and Carbon and Oxygen Isotopes of Dolostones
Next Article in Special Issue
The Roof Safety under Large Mining Height Working Face: A Numerical and Theoretical Study
Previous Article in Journal
Low-Temperature Hydrothermal Synthesis of Chromian Spinel from Fe-Cr Hydroxides Using a Flow-Through Reactor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geostatistical Modeling of Overburden Lithofacies to Optimize Continuous Mining in the Ptolemais Lignite Mines, Greece

by
Konstantinos Modis
1,*,
Daphne Sideri
2,
Christos Roumpos
3,*,
Hélène Binet
4,
Francis Pavloudakis
5 and
Nikolaos Paraskevis
3
1
School of Mining and Metallurgical Engineering, National Technical University of Athens, 15780 Zοgrafou, Greece
2
Engineering Department, Eydap SA, 15772 Athens, Greece
3
Department of Mining Engineering and Closure Planning, Public Power Corporation of Greece, 10432 Athens, Greece
4
Business Development Department, Geovariances SA, 77210 Avon, France
5
Mineral Resources Engineering Department, School of Engineering, University of Western Macedonia, 50100 Kozani, Greece
*
Authors to whom correspondence should be addressed.
Minerals 2022, 12(9), 1109; https://doi.org/10.3390/min12091109
Submission received: 7 August 2022 / Revised: 27 August 2022 / Accepted: 28 August 2022 / Published: 30 August 2022

Abstract

:
Lignite production in Greece is implemented mainly by the Public Power Corporation (PPC), with the higher production being in the Lignite Center of Western Macedonia. A continuous surface mining method is used in order to satisfy the high production needs combined with the necessity for selective mining; however, the occasional appearance of hard rock formations in the South Field mine overburdens was critical for the adoption of a discontinuous auxiliary method of rock mass removal, at these places, by explosives and large shovels. Furthermore, to minimize the delay of changing the machinery arrangements when a hard rock formation is met, an a priori knowledge of the spatial distribution of these rock masses would be catalytic. In this work, a plurigaussian simulation model of the overburden geological formations is developed in the South Field mine. This model could be used as a guide to schedule and optimize the overburden removal process. Validation of the model was affected in two ways: by direct comparing estimated to real cross-sections as observed on mine slopes or by correlating PPC’s recorded volumetric results to the average simulated hard rock percentages.

1. Introduction

Greece was considered a large producer of soft brown coal, producing about 45–70 Mt of lignite annually, mainly at the Lignite Center of Western Macedonia; however, as a result of the European CO2 emissions trading scheme and the increase in lignite price, the annual lignite production during the period 2016–2018 ranged to 35 Mt, while in the period 2019–2021, it was further reduced to almost half (17 Mt), due to the phase-out of lignite-fired power plants. The Lignite Center of Western Macedonia currently operates with two mines, namely Mavropigi and South Field, in the Ptolemais area. The remaining coal reserves in the two mines have been estimated at approximately 300 Mt.
The increasing requirements for sustainable power generation make it necessary to optimize the corresponding production process from a technical, economic, and environmental point of view. In this context, the appropriate use of lignite deposits can contribute to sustainability.
Greek lignite deposits have a multiple-layered geological structure consisting of several coal seams, which are separated mainly by calcareous and argillaceous waste beds. The necessity of selective mining due to the above stratigraphic deposit structure, combined with the requirements for high production rates, was the reason for the application of the continuous surface mining method for more than 60 years [1]; however, surface exploitation of multiple-layered lignite deposits often faces the problem of the presence of hard formations, with a particularly unfavorable effect on the productive and economical sizes. In these cases, the removal of hard formations is carried out discontinuously by trucks after the application of drilling and blasting operations to loosen the material [2]. In addition, similar semi-hard formations within the overburden strata that cannot be excavated by bucket wheel excavators (BWEs) are loaded by shovels without blasting; therefore, identifying the spatial distribution of hard formations within the mine area is crucial for the appropriate long and short-term mine planning and design activities. A recent research in the area was focused on the application of geophysical methods to detect the hard formations on the stope [3,4], but this approach can only be useful in very short-term planning.
The underlying idea of this work is to create a geological model of the burden formations overlying the lignite beds, so as to be able to optimally program the different operations for overburden removal according to the rock type and, thus, save significant time and cost. In order to model the patterns of such complex geological structures, conceptual and deterministic models are not generally as flexible and realistic to represent the internal geometry [5] adequately. For this reason, stochastic models such as geostatistical estimation and simulation are most commonly used [6,7,8,9], particularly for the characterization of lithofacies distribution.
Geostatistical estimation is used to interpolate the value of a given attribute at a given location by minimizing the error and bias of the estimate. The use of estimation methods for obtaining facies reconstruction has not been as widespread in the literature as the use of simulation methods. Johnson et al. and Ritzi et al. [10,11] applied a categorical estimation method (Indicator Kriging: IK) to obtain hydrofacies reconstructions. Falivene et al. [12] used IK to reconstruct facies distribution in a fine-grain alluvial fan. Moreover, Falivene et al. [5] compared visually and statistically several facies reconstruction methods (Truncated Inverse Distance Weighting: TIDW, Truncated Kriging: TK, Indicator Inverse Distance Weighting: IIDW, and Indicator Kriging, among others) applied to a heterogeneous coal seam. Other studies have focused on the estimation of continuous parameters that can directly be related to facies (mud fraction in Flach et al., [13]; grain-size compositions in Koike et al. [14]; or results of geotechnical cone penetration tests in Lafuerza et al. [15]. This approach is conceptually similar to continuous methods for facies reconstruction (TK and TIDW).
The scarcity of facies estimation in the literature (compared to simulation facies models) is due to the smoothing effect of estimation methods [16,17,18,19], which results in more homogeneous distributions compared to reality, thus limiting their predictive use. Moreover, as the density of the dataset for facies heterogeneity decreases, the smoothing effect of facies reconstructions increases.
Geostatistical simulations, on the other hand, aim to generate equiprobable images of the real situation, obeying the variability of the dataset and the probability distribution function [20,21]. Simulated images exhibit the complexity and the dispersion of the facies as equally probable scenarios of the internal geological structure. Among the simulation techniques for several facies, Plurigaussian Simulation (PGS) allows the incorporation of geological concepts into stochastic simulations [22]. This unique feature is difficult to establish in other similar simulation techniques. The geological rule is not only derived from a statistical analysis of the borehole data, but it also allows the inclusion of other forms of geological and empirical knowledge. We further refer to the comparison between PGS used in this paper vs. other similar simulation methods in Section 2.3. In this line, the study of Liu et al. faces the problem of the thick loose coal seam by applying numerical simulations to predict the geomechanical behavior of the overburden materials [23]. Oggeri et al. [24] applied a rational approach to select the most suitable methods for excavation, transportation, and disposal of thick formations in open pit mining operations and showed some crucial parameters that should be taken into account for their thorough management. In the Big Bell Gold Mine in the Murchison mineral field (Perth, Western Australia), the PGS techniques were used to simulate the lithotypes of the area [25]. The PGS and IK methods were applied in a mining exploration area in southern Ecuador, where only surficial lithological and geochemical data samples were available. It has arrived that PGS reaches a higher accuracy compared with the IK [26]. In addition, PGS was applied in a Chilean porphyry copper deposit to simulate five lithological facies spatially, and it is considered a simple and mathematically consistent method, as it gives information regarding the topology of lithology, geological knowledge, as well as the spatial continuity of the lithology [27]. Stochastic modeling based on PGS to obtain numerical models of the ore body spatial variability was applied in the Sungun porphyry copper deposit in Iran [28].
The objective of this study is to identify the spatial distribution of the hard formations within the mine area. In this framework, the PGS methodology is used to simulate the geology of the overburden formations in the Ptolemais lignite mines, aiming to determine the average volume of hard formation lithofacies in broader mining segments. Based on these simulation results, the cost of hiring discontinuous mining equipment in these areas can be calculated and optimized. From those mentioned above, it is clear that the majority of the scientific publications focus on the use of geostatistics for the reserve’s estimation and the ore deposits’ geological investigation. In contrast, to our knowledge, no prior studies have focused on hard rock overburden management in mining areas.

2. Materials and Methods

2.1. Geology of the Study Area

The South Field Mine (SFM) forms the southern boundary of the total Ptolemais mining area (Figure 1). It extends approximately 8 km towards the northwest and the southeast. Its SW–NE extension ranges between approximately 5 km in the northwest and 4 km in the southeast.
Approximately 30% of the overburden formations consist of hard and semi-hard material that cannot be excavated by BWE continuous mining system [29]. The hard formations in lignite mines consist of fine and coarse clastic sediments such as clays, marls, gravel, and conglomerates with embedded hard layers of sandstones, cemented conglomerates, and mudstones [30,31]. The average specific weight of the overburden is 19.62 kN/m3 (2 ton/m3), and the bulking factor ranges from 1.4 to 1.5. The average thickness of overburden material is 90 m.
The series of the overburden rocks of the SFM can be distinguished according to the rock coloration and way of formation into the following groups:
Red to brown clastic sediments: This group consists mainly of clays, conglomerates, and limestone gravels. The thickness of the whole formation is about 25 m, while the average thickness of hard material is approximately 10 m.
Gray to yellow clastic sediments: This group consists of clay, sandstones, sand, conglomerates, and siltstones. The average thickness of the whole formation is 25 m, while the average thickness of the hard material is 11 m.
Green-gray clay, sand, and silt sediments: The thickness of this formation varies from 25 m to 50 m. No hard material has been found in this formation.
The extent and the distribution of hard and semi-hard formations in the overburden strata are presented in the geological section of Figure 2.
Water outlets from the coarsely clastic sediments have been observed in the opencast mine. Furthermore, these layers, which principally consist of conglomerates, sandstone, and calcareous marl, locally incorporate consolidated intercalations, seriously impeding mining operations [30].

2.2. Dataset and Discretization

PPC created an extended borehole network throughout the SFM by drilling 777 exploratory boreholes (Figure 3).
The borehole dataset contains information on 51,783 core samples. Sample lengths range from 0.01 to 213 m, with an average value of 2.5 m. and a standard deviation of 6.8 m. The large number of samples in each borehole and the diversity of sample lengths are due to the strong multi-seam character of the deposit and the variation of the quality characteristics between the seams [1]. For simplicity reasons, the geological descriptions of the samples were grouped into five lithotypes: hard rock formations (HD), clay (AL), marl (MG), sand with gravels (SG), and lignite (CO) (Figure 4). The (HD) formations consist of conglomerates, claystone, sandstone, and marly limestone; the (AL) of clay, soil, silt, and losses; the (MG) of marls and chalk; the (SG) of sand and gravels; the (CO) of coal (lignite). The aforementioned geological formations are interpreted in the photo of SFM presented in Figure 5. The spatial distribution of HD facies causes significant problems and delays since bucket wheel excavators cannot dig the entire face of the bench (see in the photo of Figure 5 the traces of bucket teeth stop in the lower marly limestone seam).
The overburden area can then be defined between the earth’s surface and the geological rooftop of the lignite orebody, as seen in Figure 6. This domain is divided into 58,609 horizontal parallelepiped blocks sizing 100 m × 100 m × 10 m each. For functionality reasons, the drills’ data were discretized in the 3D overburden grid, and a lithofacies value was assigned to each block crossed by a borehole. This operation was accomplished in two steps. First, the boreholes were cut into regular samples of 10 m, a size that corresponds to the depth of the basic grid cell. Then, after compositing, these regularized samples were migrated to the closest grid block so a single lithotype value would be assigned to the blocks crossed by a well. These blocks will, this way, become the conditioning data in the lithofacies simulation process.

2.3. Methods

Geological modeling is different from the usual geostatistical applications in mining, in the sense that it confronts the need to simulate categorical instead of numerical variables, such as the geological facies. The usual alternatives to treat this problem, when in covariance-based geostatistics, are to use the traditional Sequential Indicator Simulation (SIS) or the PGS, the second being most appealing and mathematically rigorous [32,33,34]. On the other hand, Multiple Point Geostatistics (MPG) is more appropriate to capture very complex and channelized facies [22] and is preferentially applied in under- or over-informed models and usually in groundwater or petroleum exploration [35].
The idea behind SIS is to simulate each category of the study variable in every location with an indicator, i.e., 0 or 1, thus having as many indicators as the categories. PGS, instead, can use only one artificial (hidden) continuous variable that is conveniently truncated to as many intervals as the real categories are.
More rigorously, in PGS, the idea is to represent and model the lithofacies using one or more hidden continuous Gaussian Random Functions (GRFs) at every point of the study area and then to convert them back to facies using a coding function:
Consider a vector of standardized (GRF) Z ( x ) = ( Z 1 ( x ) , , Z q ( x ) ) where x R 3 . Let ( D 1 , , D k ) be a measurable partition of R q into k disjoint subdomains. A categorical random field F with k categories F i ,   i = 1 , , k is then obtained by letting
x R 3 ,   F ( x ) = i ,   if   and   only   if   Z ( x ) D i
while the indicator random field for each facies F i is defined as:
x R 3 , I F i ( x ) = { 1   if   I ( x ) = i 0   otherwise
The proper selection of the partition ( D 1 , , D k ) is very important in the arrangement of categories, as we will show below, because it defines the permissible contacts and characterizes the spatial distribution of facies [34].
Then, one can approximate the proportion of a particular facies F i at point x by the probability of having these facies F i at that point x :
P F i ( x ) = P { ( Z 1 ( x ) , , Z q ( x ) ) D i } = D i g Σ ( Z 1 ( x ) , , Z q ( x ) )   dZ 1 dZ q
where g Σ ( Z 1 ( x ) , , Z q ( x ) ) is the q-variate Gaussian density of Z ( x ) and Σ its covariance matrix.
In order to compute P F i it is enough to know the D i ’s of the chosen partition, but the inverse is, in general, not true because there can be many solutions to the problem. This means that one has a degree of freedom here to decide the spatial arrangement of D i ’s. and it is here where the empirical geological knowledge can conveniently enter the scene under the “geological rule”. To ease the solution, the partitioning of the Gaussian space is affected into cuboids defined by their projections on the Gaussian axes. The endpoints of the projections on each axis are then considered thresholds. In practical applications, the proportions of facies are known experimentally from core data by calculating the Vertical Proportion Curves (VPCs), so the k × 2 q thresholds should be determined together with the definition of the geological rule. This can be achieved by numerically inverting Equation (3).
After the definition of the geological rule that describes the sets associated with each category, the second integrant to be determined for the complete specification of the plurigaussian model is the multivariate variogram functions of the underlying GRFs. These variograms represent the structural characteristics of the physical phenomenon, such as spatial continuity and anisotropies; however, the inference of this variogram model is one of the major difficulties associated with the use of the plurigaussian model [8]. This is because, while the available information restricts the experimental indicator variograms and cross variograms of the facies, the required variogram models are ones of the underlying GRFs, which are hidden by the truncation procedure; therefore, these models, which represent the spatial continuity of facies, were usually determined by a trial and error procedure, using, for example, the program VMODEL [36]. In this work, we use an algorithm based on the Pairwise Likelihood (PL) maximization principle to directly estimate the variograms of the underlying GRFs [37].
After the definition of the model variogram of Z ( x ) , the Gibbs sampler is used to generate Gaussian values at the sample points that have the appropriate covariance and belong to the right intervals. Then, a conditional simulation is used to produce several realizations conditional to the Gaussian values at the above points. At the end of the process, the Gaussian values of the realizations are converted to facies using the truncation rule.

3. Results and Discussion

3.1. Structural Analysis and Spatial Modeling

First, care should be taken so that the simulation results are consistent with the deposition process, meaning that calculations have to be performed inside layers conforming to the deposition beds in order to attend the stratifying order. An approach to secure these geological restrictions is to work in a flattened domain while the flattening process confronts the sedimentogenesis model. In the present work, the Isatis.neo™ algorithm was used to unfold the working grid parallel to the top surface (Figure 7).
Next, the proportions of each lithotype have to be estimated in each block of the flattened grid. These proportions are important ingredients of the modeling process because, as mentioned in the previous paragraph, they approximate the probability of having certain phases at a specified block, and so they constitute the soft constraints of the plurigaussian model. The proportions are estimated from the borehole data in three dimensions (3D) using an interpolation approach based on Stochastic Partial Differential Equations (Figure 7). Finally, the lithological rule, which embodies the empirical geological knowledge concerning lithofacies contacts, complementary to the proportions, allows the calculation of the phase’s probabilities, as shown in Figure 8.
After the estimation of the proportions and the lithological rule, the last ingredient for the complete determination of the plurigaussian model is the variograms of the two underlying GRFs. In the present study, these variograms are calculated directly from Gaussian values assigned to the lithotype samples, as previously mentioned [37]. The variogram of the first GRF is a spherical scheme with a horizontal range of 300 m, and an anisotropic vertical range of 30m. The variogram of the second GRF is an exponential model with sill = 1, a horizontal range of 400 m, and an anisotropic vertical range of 200 m (Figure 9).
The main part of the plurigaussian simulation is performed in four steps: First, a Gibbs sampler was used to conditionally simulate the two GRFs on the data points by assessing the right Gaussian values to the experimental lithofacies. Second, the two GRFs were unconditionally simulated at each block of the flattened working grid. Then, the realizations of the second simulation were conditioned on the results of the first step. Finally, the Gaussian outputs of the last simulations were inverted to lithofacies by using the thresholds calculated by the combination of the lithotypes proportions and the lithological rule.
The last thing to perform is to copy the simulation results from the unfolded grid to the original domain by using the inverse geometrical transformation. An overview of the resulting plurigaussian model realization is shown in Figure 10.

3.2. Validation of the Results

Unlike estimation, where the best map is sought under a certain optimality criterion (for kriging, it is the minimization of the mean square error), in geostatistical simulations, the main concern is to provide maps that reproduce the spatial variability of the study variable. As a result, different realizations are produced, and there is not a single best realization, since they are all considered equally likely to occur. This distribution of outcomes is interpreted as the uncertainty of the variable and thus allows us to explore the risk associated with the produced model of reality. As previously stated, in our case, this risk is connected to the average cost of discontinuously removing the hard rock volumes of overburden formations.
Validation of simulation results cannot rely on a point-wise comparison but rather on examining the average match of estimated images to reality. In our case, we examine the similarity of our model to reality in two available experimental cross sections, AA’ and BB’, the traces of which are presented in Figure 1 and Figure 2.
Figure 11and Figure 12 show that the shapes and relative volumes of the overburden formations are, in general, similar to the realistic geological cross-sections.
A more straightforward way, though, to validate the model results is to compare the volume of the assessed hard rock formations to the actual one. In this framework, PPC has available accounting volumetric results in quite a large area on the northern part of the orebody, as shown in Figure 13 and Figure 14. In order to calculate these actual results in the ArcGIS environment, the top and bottom surface elevation models were elaborated with value ranges from 670 m. a.s.l. to 740 m. a.s.l. The top surface is considered the natural topographic relief, whereas the bottom surface is considered the hard formations excavated surface up to a specific terrace. By using the “surface difference” tool from the ArcGIS 3D Analyst Toolbox, the volume calculation was affected inside the defined volume calculation polygon. The total area amounts to 8,500,000 m3, while the company’s reported hard rock volume in this area amounts to 1,700,000 m3.
In order to check these results with the model output, the top and bottom surfaces were entered into the grid to define the calculation domain. Finally, a number of 84 grid blocks were found inside this part of the model. In a total of 100 simulations, an average percentage of 18 blocks was characterized as HD, which means that a total of 1,800,000 m3 are expected to be met here. This is quite in accordance with the company’s reported results, which, as previously mentioned, amount to 1,700,000 m3 in the examined area.
Considering the improvement of the results obtained by the suggested method, it could be noted that this approach can offer a more accurate prediction regarding the spatial distribution of various rock lithofacies in broader mining sectors of coal or other surface mines. In addition, complicated structures of the overburden material can be successfully modeled and simulated more efficiently than indicator kriging or other similar categorical simulation methods stated in the introduction. Furthermore, the analysis can be extended to other types of rocks, such as water-permeable rocks.

4. Conclusions and Future Research

Plurigaussian simulations can be successfully applied to the overburden lithofacies of the South Field Mine by using the plethora of accessible core data with geological descriptions.
The available field data, such as cross-sections drawn on mine slopes and actual production data from specific areas of non-continuous mining equipment operations, can be used to calibrate the model parameters, such as the unfolding of the working domain and the regularization of core samples.
Validation of model results showed that the estimated thickness and quantities of the hard rock formations in broader mining sectors agree with the available PPC accounting volumetric data.
Thus, model predictions can be successfully used for the optimization and scheduling of the mining process by minimization of the delays from changing the machinery arrangements when hard rock formations are encountered. This approach can be incorporated into the mine planning and design procedures, emphasizing the mining cost minimization considering a better mine production scheduling. In this context, the improved spatial interpolation of hard rock formations could support better planning and scheduling of non-continuous mining equipment in specific mining sectors and also the strategic planning of continuous surface mining equipment (e.g., bucket wheel excavators). Significantly, if the developed methodology could be incorporated from the first stages of the strategic mining plan, it could help to achieve the mining sustainability goals (e.g., environmental by minimizing the impacts of non-continuous mining, or economic ones by minimizing the mine production costs, also considering the delays from changing the mining machinery arrangements).
Possible limitations of the applied method include the need for accurate interpretation of borehole samples and a suitable definition of the lithological rule based on geological experience.
Future research could focus on examining multiple-point statistics, an alternative method for lithofacies modeling and prediction of hard rock spatial distribution.

Author Contributions

Conceptualization, K.M., D.S. and C.R.; methodology, K.M., D.S. and H.B.; validation, D.S., C.R., K.M., F.P. and H.B.; formal analysis, K.M., C.R., F.P., N.P. and D.S.; data curation, D.S., C.R., N.P, F.P. and K.M.; writing—original draft preparation, K.M., C.R. and N.P.; writing—review and editing, D.S., C.R., F.P., N.P. and K.M.; visualization, K.M., C.R., F.P., N.P. and D.S.; supervision, K.M., C.R. and D.S.; project administration, C.R. and K.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Agni Patra and Aikaterini Servou, geologists, for their assistance in the pre-processing of the provided data and their contribution to the interpretation of the geological formations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sideri, D.; Roumpos, C.; Pavloudakis, F.; Paraskevis, N.; Modis, K. Multivariate Geostatistical Modeling of Lower Calorific Value in Multi-Seam Coal Deposits. Appl. Sci. 2020, 10, 6208. [Google Scholar] [CrossRef]
  2. Agioutantis, Z.G.; Panagiotou, G.N.; Kavouridis, C. Evaluation and Analysis of Blasting Procedures for Removing Hard Formations at the South Field Lignite Mine. In Proceedings of the 17th International Mining Congress and Exhibition of Turkey, Ankara, Turkey, 19–22 June 2001. [Google Scholar]
  3. Vafidis, A.; Kritikakis, G.; Economou, N.; Galetakis, M.; Soultana, A.; Michalakopoulos, T.; Apostolopoulos, G.; Roumpos, C.; Pavloudakis, F. Assessment of Geophysical Methods for the Detection of Hard Rock Formations in Lignite Mines Employing Bucket Wheel Excavators. The Case of South Field Lignite Mine, Macedonia, Greece. Górnictwo Odkrywkowe 2018, 4, 23–30. [Google Scholar]
  4. Roumpos, C.; Sideri, D.; Pavloudakis, F.; Kolovos, N.; Michalakopoulos, T.; Apostolopoulos, G.; Galetakis, M.; Vafidis, A. Possibilities for Improving Work Efficiency of Continuous Surface Mining Systems Operating in Rocks with Excessive Digging Resistance. Górnictwo Odkrywkowe 2018, 4, 31–41. [Google Scholar]
  5. Falivene, O.; Cabrera, L.; Muñoz, J.A.; Arbués, P.; Fernández, O.; Sáez, A. Statistical Grid-Based Facies Reconstruction and Modelling for Sedimentary Bodies. Alluvial-Palustrine and Turbiditic Examples. Geol. Acta Int. Earth Sci. J. 2007, 5, 199–230. [Google Scholar]
  6. Bierkens, M.F.P.; Burrough, P.A. The Indicator Approach to Categorical Soil Data: I. Theory. J. Soil Sci. 1993, 44, 361–368. [Google Scholar] [CrossRef]
  7. Deutsch, C.V. A Sequential Indicator Simulation Program for Categorical Variables with Point and Block Data: BlockSIS. Comput. Geosci. 2006, 32, 1669–1681. [Google Scholar] [CrossRef]
  8. Mariethoz, G.; Renard, P.; Cornaton, F.; Jaquet, O. Truncated Plurigaussian Simulations to Characterize Aquifer Heterogeneity. Ground Water 2009, 47, 13–24. [Google Scholar] [CrossRef]
  9. Stafleu, J.; Maljers, D.; Gunnink, J.L.; Menkovic, A.; Busschers, F.S. 3D Modelling of the Shallow Subsurface of Zeeland, the Netherlands. Neth. J. Geosci. 2011, 90, 293–310. [Google Scholar] [CrossRef]
  10. Johnson, N.M.; Dreiss, S.J. Hydrostratigraphic Interpretation Using Indicator Geostatistics. Water Resour. Res. 1989, 25, 2501–2510. [Google Scholar] [CrossRef]
  11. Ritzi, R.W.; Dominic, D.F.; Brown, N.R.; Kausch, K.W.; McAlenney, P.J.; Basial, M.J. Hydrofacies Distribution and Correlation in the Miami Valley Aquifer System. Water Resour. Res. 1995, 31, 3271–3281. [Google Scholar] [CrossRef]
  12. Falivene, O.; Cabrera, L.; Sáez, A. Large to Intermediate-Scale Aquifer Heterogeneity in Fine-Grain Dominated Alluvial Fans (Cenozoic As Pontes Basin, Northwestern Spain): Insight Based on Three-Dimensional Geostatistical Reconstruction. Hydrogeol. J. 2007, 15, 861–876. [Google Scholar] [CrossRef]
  13. Flach, G.; Hamm, L.L.; Harris, M.K.; Thayer, P.A.; Haselow, J.S.; Smits, A.D. Characterizing Hydrogeologic Heterogeneity Using Lithologic Data; Savannah River Site (SRS): Aiken, SC, USA, 1997; p. 565245.
  14. Koike, K.; Shiraishi, Y.; Verdeja, E.; Fujimura, K. Three-Dimensional Interpolation and Lithofacies Analysis of Granular Composition Data for Earthquake-Engineering Characterization of Shallow Soil. Math. Geol. 1998, 30, 733–759. [Google Scholar] [CrossRef]
  15. Lafuerza, S.; Canals, M.; Casamor, J.L.; Devincenzi, J.M. Characterization of Deltaic Sediment Bodies Based on in Situ CPT/CPTU Profiles: A Case Study on the Llobregat Delta Plain, Barcelona, Spain. Mar. Geol. 2005, 222–223, 497–510. [Google Scholar] [CrossRef]
  16. Hohn, M.E. An Introduction to Applied Geostatistics. Comput. Geosci. 1991, 17, 471–473. [Google Scholar] [CrossRef]
  17. Olea, R.A.; Pawlowsky, V. Compensating for Estimation Smoothing in Kriging. Math. Geol. 1996, 28, 407–417. [Google Scholar] [CrossRef]
  18. Journel, A.G.; Kyriakidis, P.C.; Mao, S. Correcting the Smoothing Effect of Estimators: A Spectral Postprocessor. Math. Geol. 2000, 32, 787–813. [Google Scholar] [CrossRef]
  19. Yamamoto, J.K. Correcting the Smoothing Effect of Ordinary Kriging Estimates. Math. Geol. 2005, 37, 69–94. [Google Scholar] [CrossRef]
  20. Deutch, D.V.; Journel, A.G. GSLIB: Geostatistical Software Library and User’s Guide, 2nd ed.; Oxford University Press: New York, NY, USA, 1997. [Google Scholar]
  21. Stein, M.L.; Chiles, J.-P.; Delfiner, P. Geostatistics: Modeling Spatial Uncertainty. J. Am. Stat. Assoc. 2000, 95, 335. [Google Scholar] [CrossRef]
  22. Hashemi, S.; Javaherian, A.; Ataee-pour, M.; Khoshdel, H. Two-Point versus Multiple-Point Geostatistics: The Ability of Geostatistical Methods to Capture Complex Geobodies and Their Facies Associations—an Application to a Channelized Carbonate Reservoir, Southwest Iran. J. Geophys. Eng. 2014, 11, 065002. [Google Scholar] [CrossRef]
  23. Liu, W.; Pang, L.; Xu, B.; Sun, X. Study on Overburden Failure Characteristics in Deep Thick Loose Seam and Thick Coal Seam Mining. Geomat. Nat. Hazards Risk 2020, 11, 632–653. [Google Scholar] [CrossRef]
  24. Oggeri, C.; Fenoglio, T.M.; Godio, A.; Vinai, R. Overburden Management in Open Pits: Options and Limits in Large Limestone Quarries. Int. J. Min. Sci. Technol. 2019, 29, 217–228. [Google Scholar] [CrossRef]
  25. Dunn, R. Plurigaussian Simulation of Rocktypes Using Data from a Gold Mine in Western Australia. Master’s Thesis, Edith Cowan University, Joondalup, WA, Australia, 2011. [Google Scholar]
  26. Guartán, J.A.; Emery, X. Predictive Lithological Mapping Based on Geostatistical Joint Modeling of Lithology and Geochemical Element Concentrations. J. Geochem. Explor. 2021, 227, 106810. [Google Scholar] [CrossRef]
  27. Emery, X. Simulation of Geological Domains Using the Plurigaussian Model: New Developments and Computer Programs. Comput. Geosci. 2007, 33, 1189–1201. [Google Scholar] [CrossRef]
  28. Talebi, H.; Asghari, O.; Emery, X. Application of Plurigaussian Simulation to Delineate the Layout of Alteration Domains in Sungun Copper Deposit. Open Geosci. 2013, 5, 514–522. [Google Scholar] [CrossRef]
  29. Kavouridis, K.; Galetakis, M.; Roumpos, C.; Pavloudakis, F. Methods and Technological Improvements for the Efficient Removal of the Overburden Hard Rock Formations at South Field Lignite Mine, Ptolemais, Greece. Ann. Univ. Petroşani Min. Eng. 2008, 9, 132–148. [Google Scholar]
  30. Anastopoulos, J.C.; Koukouzas, C.N. Economic geology of the southern part of Ptolemais lignite basin I.G.S.R. Geol. Geophys. Res. 1972, XVI, 1–189. [Google Scholar]
  31. Kavouridis, K.; Agioutantis, Z. The Impact of a Large-Scale Failure of the External Waste Dump to Theoperations at the South Field Mine. In Proceedings of the 2nd International Conference on Advances in Mineral Resources Management and Environmental Geotechnology, Hania, Greece, 25–27 September 2006; pp. 183–190. [Google Scholar]
  32. Emery, X. Properties and Limitations of Sequential Indicator Simulation. Stoch. Environ. Res. Risk Assess. 2004, 18, 414–424. [Google Scholar] [CrossRef]
  33. Modis, K.; Sideri, D. Geostatistical Simulation of Hydrofacies Heterogeneity of the West Thessaly Aquifer Systems in Greece. Nat. Resour. Res. 2013, 22, 123–138. [Google Scholar] [CrossRef]
  34. Armstrong, M.; Galli, A.; Beucher, H.; Loc’h, G.; Renard, D.; Doligez, B.; Eschard, R.; Geffroy, F. Plurigaussian Simulations in Geosciences; Springer: Berlin/Heidelberg, Germany, 2011; ISBN 978-3-642-19606-5. [Google Scholar]
  35. Mariethoz, G. When Should We Use Multiple-Point Geostatistics? In Handbook of Mathematical Geosciences; Daya Sagar, B.S., Cheng, Q., Agterberg, F., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 645–653. ISBN 978-3-319-78998-9. [Google Scholar]
  36. Emery, X.; Silva, D.A. Conditional Co-Simulation of Continuous and Categorical Variables for Geostatistical Applications. Comput. Geosci. 2009, 35, 1234–1246. [Google Scholar] [CrossRef]
  37. Desassis, N.; Renard, D.; Beucher, H.; Petiteau, S.; Freulon, X. A Pairwise Likelihood Approach for the Empirical Estimation of the Underlyingvariograms in the Plurigaussian Models. arXiv 2015, arXiv:1510.02668. [Google Scholar] [CrossRef]
Figure 1. South Field Mine and relative land use. Location of cross-sections (A–A’ and B–B’).
Figure 1. South Field Mine and relative land use. Location of cross-sections (A–A’ and B–B’).
Minerals 12 01109 g001
Figure 2. Geological setting of the South Field Mine, included the location of cross-sections (A–A’ and B–B’).
Figure 2. Geological setting of the South Field Mine, included the location of cross-sections (A–A’ and B–B’).
Minerals 12 01109 g002
Figure 3. South Field Mine, drills, and location of cross sections and photo.
Figure 3. South Field Mine, drills, and location of cross sections and photo.
Minerals 12 01109 g003
Figure 4. Boreholes with regularized data and the lignite rooftop surface.
Figure 4. Boreholes with regularized data and the lignite rooftop surface.
Minerals 12 01109 g004
Figure 5. Interpretation of geological formations (the location of the photo is shown in Figure 3).
Figure 5. Interpretation of geological formations (the location of the photo is shown in Figure 3).
Minerals 12 01109 g005
Figure 6. Discretization of the study domain.
Figure 6. Discretization of the study domain.
Minerals 12 01109 g006
Figure 7. Unfolded domain with HD proportions (%).
Figure 7. Unfolded domain with HD proportions (%).
Minerals 12 01109 g007
Figure 8. Lithological rule.
Figure 8. Lithological rule.
Minerals 12 01109 g008
Figure 9. Experimental and fitted anisotropic variogram models for the two GRFs: (a) a spherical scheme and (b) an exponential model.
Figure 9. Experimental and fitted anisotropic variogram models for the two GRFs: (a) a spherical scheme and (b) an exponential model.
Minerals 12 01109 g009
Figure 10. Overview of the PGS model.
Figure 10. Overview of the PGS model.
Minerals 12 01109 g010
Figure 11. Comparison between (a) a realization of the modeland (b) reality in Section B–B’. The bold frame shows the model area that corresponds to the real image. In addition, note that the model does not include the lignite beds shown in black color and their interburden and underburden waste layers on Section B–B’.
Figure 11. Comparison between (a) a realization of the modeland (b) reality in Section B–B’. The bold frame shows the model area that corresponds to the real image. In addition, note that the model does not include the lignite beds shown in black color and their interburden and underburden waste layers on Section B–B’.
Minerals 12 01109 g011
Figure 12. Comparison between (a) a realization of the model and (b) reality in Section A–A’. The bold frame shows the model area that corresponds to the real image. In addition, note that the model does not include the lignite beds shown in black color and their interburden and underburden waste layers on Section A–A’.
Figure 12. Comparison between (a) a realization of the model and (b) reality in Section A–A’. The bold frame shows the model area that corresponds to the real image. In addition, note that the model does not include the lignite beds shown in black color and their interburden and underburden waste layers on Section A–A’.
Minerals 12 01109 g012
Figure 13. Volume calculation polygon and location of cross-section B–B’.
Figure 13. Volume calculation polygon and location of cross-section B–B’.
Minerals 12 01109 g013
Figure 14. Overburden model blocks inside the test area. Only the bottom surface is shown.
Figure 14. Overburden model blocks inside the test area. Only the bottom surface is shown.
Minerals 12 01109 g014
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Modis, K.; Sideri, D.; Roumpos, C.; Binet, H.; Pavloudakis, F.; Paraskevis, N. Geostatistical Modeling of Overburden Lithofacies to Optimize Continuous Mining in the Ptolemais Lignite Mines, Greece. Minerals 2022, 12, 1109. https://doi.org/10.3390/min12091109

AMA Style

Modis K, Sideri D, Roumpos C, Binet H, Pavloudakis F, Paraskevis N. Geostatistical Modeling of Overburden Lithofacies to Optimize Continuous Mining in the Ptolemais Lignite Mines, Greece. Minerals. 2022; 12(9):1109. https://doi.org/10.3390/min12091109

Chicago/Turabian Style

Modis, Konstantinos, Daphne Sideri, Christos Roumpos, Hélène Binet, Francis Pavloudakis, and Nikolaos Paraskevis. 2022. "Geostatistical Modeling of Overburden Lithofacies to Optimize Continuous Mining in the Ptolemais Lignite Mines, Greece" Minerals 12, no. 9: 1109. https://doi.org/10.3390/min12091109

APA Style

Modis, K., Sideri, D., Roumpos, C., Binet, H., Pavloudakis, F., & Paraskevis, N. (2022). Geostatistical Modeling of Overburden Lithofacies to Optimize Continuous Mining in the Ptolemais Lignite Mines, Greece. Minerals, 12(9), 1109. https://doi.org/10.3390/min12091109

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop