Next Article in Journal
Logging Identification and Distribution of Bauxite in the Southwest Ordos Basin
Previous Article in Journal
Provenance, Depositional Environment, and Paleoclimatic Conditions of a Near-Source Fan Delta: A Case Study of the Permian Jiamuhe Formation in the Shawan Sag, Junggar Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Objective Approach for Optimizing the Layout of Additional Boreholes in Mineral Exploration

by
Amin Hossein-Morshedy
1,
Farzaneh Khorram
2,3,4,* and
Xavier Emery
3,4
1
Department of Mining and Metallurgical Engineering, Yazd University, Yazd 8915818411, Iran
2
Facultad de Ingeniería, Universidad Andres Bello, Santiago 8370134, Chile
3
Department of Mining Engineering, Universidad de Chile, Santiago 8370448, Chile
4
Advanced Mining Technology Center, Universidad de Chile, Santiago 8370448, Chile
*
Author to whom correspondence should be addressed.
Submission received: 14 August 2023 / Revised: 19 September 2023 / Accepted: 23 September 2023 / Published: 25 September 2023

Abstract

:
Accurate subsurface exploration requires an optimal network of boreholes. This paper proposes a multi-objective approach to optimize the layout of additional exploratory boreholes. In order to illustrate this approach, geochemical analyses of core samples at the eastern part of the Kahang copper deposit, Central Iran, were used. A measure of the grade uncertainty (kriging standard deviation) and a confidence measure on the ore/waste classification were first calculated by implementing ordinary and indicator kriging. An ore value function was then determined to measure the total value of each block by considering the grades of all the effective variables and their ore membership degree derived from a fuzzy treatment of the grades. Finally, a misclassification cost is defined for each block based on the expected economic effects of ore loss and waste dilution. As a result, an index for the selection of additional boreholes was introduced in order to maximize the kriging standard deviation, the ore misclassification cost, and the ore value and to minimize the confidence measure. Applied to the Kahang copper deposit, this index allowed the prioritization of areas for infill sampling, leading to the recommendation for eight vertical and two directional additional boreholes.

1. Introduction

The spatial sampling design describes the sample number and its spatial pattern in a selected area [1]. In the additional sampling stage, we focus on the questions of where to locate additional samples and how many samples are needed to end this stage when spatial data are collected. This situation occurs in various disciplines, including mineral and oil exploration, environmental studies, oceanography, agriculture, meteorology, hydrology, and forestry [2,3]. Optimal additional samples must include valuable information and significantly improve the spatial distribution modeling of the regionalized variables under study [4]. In case of sampling shortage, there is a lack of sufficient knowledge about the target; in the reverse situation, sampling redundancy leads to increased exploration costs [3,4,5]. Borehole drilling is a common tool for the purpose of core sampling for subsurface investigation, in particular, in mineral exploration [5,6,7]. The amount of boreholes and the accuracy of subsequent mineral resource estimations increase during mineral exploration steps.
One of the first attempts to optimize spatial sampling designs is associated with a probabilistic geometry approach [8,9]. The probability of intersection between the target and the exploration network was calculated as a function of the target geometry and its relative orientation concerning the directional and dimensional properties of the exploration network. The target assumed ellipses because the shape of the surface projection of many natural resource targets can be approximated by an ellipse [9,10,11,12]. The optimum spacing of exploratory boreholes was evaluated by maximizing the expected gross drilling return (GDR) [13].
A second solution for optimal sampling design is geostatistical error management. Originally, the foundation of this approach is the reduction of the estimation error measured by a kriging variance or a kriging standard deviation, which can be either a global error [14,15] or an average or a maximum local error [16,17,18]. The samples are either ranked one at a time or as a set of samples are selected (e.g., via the definition of a regular sampling design and a sampling mesh) based on minimizing the kriging variance or kriging standard deviation and the number of samples. These methods have been extended by the use of simulation and optimization algorithms [19,20,21,22,23,24,25].
In the past decades, many other geostatistical objective functions have been proposed as an alternative to the kriging variance for optimal sampling design, such as the weighted kriging variance [26,27,28], interpolation variance [29,30,31], combined variance [31,32,33], conditional variance [34], information entropy [31,35], value of information [36], efficacy of information [37], GET (grade-estimation error-thickness) function [10,38], cross-validation error [39], interquartile range [40], probability interval widths [34,41], probability of classification error [42], probability of threshold exceedance [43], expected ore value [44], expected cost of classification errors [45], or increase of indicated and measured mineral resource categories [46]. These objective functions account for local ore grade variability, expected grade, expected productivity, expected profit, and/or reduction of uncertainty, unlike the kriging variance that only depends on the spatial correlation (variogram) and geometric configuration of the data. However, most of them focus on a single objective for optimizing a spatial sampling design in circumstances where this design should account for a multiplicity of factors.
Geologists, engineers, and investors are not only interested in reducing the uncertainty on the unknown ore grade values but also in how much the additional samples increase the confidence in the ore-waste classification, given some cut-off grades on the main products and by-products, how much they reduce the cost of ore loss (ore material wrongly sent to dump) or of waste dilution (waste material wrongly sent to the processing plant), or how much they increase the potential economic value of the ore material. This is the main motivation of this paper, which proposes a simple and meaningful multi-objective approach for infill-sampling design.
The paper is outlined as follows: Section 2 presents the eastern Kahang porphyry deposit, which is selected for the study, and the available data. The proposed methodology, which mixes statistical, geostatistical, and fuzzy logic considerations, is then described in Section 3 (as summarized in Figure 1), and the results are illustrated and discussed in Section 4. Conclusions follow in Section 5.

2. Case Study: Eastern Kahang Deposit

Deposit

The Kahang porphyry copper deposit is located 73 km NE of Isfahan in Central Iran. This porphyry system is situated on the Urumieh-Dokhtar magmatic arc, which extends from NW to SE and contains the main copper-bearing regions in Iran, such as Sarcheshmeh, Sungun, Meiduk, Darehzar, Sarkuh, and Daralu [47,48,49]. The Kahang ore consists of three main parts: eastern, central, and western; the eastern body has been drilled for a detailed exploration. It mainly consists of andesite, andesite porphyry, dacite porphyry, quartz monzonite, diorite, andesitic breccia, and quartz veins with iron oxide fillings. Andesitic lavas and Eocene volcanic breccias with propylitic alteration as the oldest igneous unit have a high potential for copper–molybdenum mineralization. Dacitic units with porphyry texture are accompanied by quartz-sericite, argillic, and quartz–goethite alterations intruded into andesite rocks. Semi-deep sections of andesite porphyry are intruded by quartz monzonite and diorite stocks with phyllic and propylitic alterations leading to mineralization. The alteration zoning from center to margin includes silicification, phyllic, argillic, and propylitic, respectively. In this porphyry system, the ore minerals comprise chalcocite, covelite, chalcopyrite, pyrite, malachite, chalcantite, magnetite, hematite, jarosite, and goethite; these Fe-hydroxides are concerned with the central part of alteration. The studied area is affected by two fault systems with NE–SW and NW–SE orientations [50,51,52]. In Figure 2, the geological map with the location of the eastern part of the Kahang deposit is depicted.
In the eastern region of Kahang, 32 vertical and directional boreholes were drilled in the various mineral exploration phases, totaling about 15,550 m of drilling. The depth of these boreholes varies from 100 to 700 m, and their dip fluctuates within a range of 60 to 90 degrees. Cores samples were assayed for Cu, Mo, Ag, and Au and were composited to a length of 2 m, which corresponds to the mode of the assayed core lengths. Most often (>95% of the time), the original core length was in the range of 1 to 3 m; core samples larger than 3 m had negligible grades for all the geochemical elements, so their decomposition into 2 m composites was not an issue.
The Kahang deposit has a significant value with copper and molybdenum, which constitute the main products of interest. Their three-dimensional distribution in the composite samples is illustrated in Figure 3. Gold and silver anomalies were also detected in limited areas that are reasonable prospects for eventual economic extraction and may, therefore, be considered by-products. The statistical parameters of Cu, Mo, Ag, and Au are shown in Table 1. The mean values are 1800, 34, 1.17, and 0.01 ppm, respectively.
Geologists and engineers of Kahang’s project defined thresholds (cut-offs) for copper (1500 ppm), molybdenum (400 ppm), silver (20 ppm), and gold (0.5 ppm), which have been applied to codify the geochemical grades into indicators. These thresholds are considered in the definition of the recoverable mineral resources subject to an eventual economic extraction by open pit mining.

3. Methodology

3.1. Directional Variogram Analysis

The composite samples are used for variogram analysis. Three-dimensional experimental variograms are first computed to analyze the spatial distributions of Cu, Mo, Ag, and Au for 24 directions with varying azimuths and inclinations. In the first instance, theoretical variogram models are fitted to the directional experimental variograms in the identified main direction of continuity, which depends on the variable under consideration, with a range of approximately 70 m. The models consider a nugget effect and an exponential structure for Cu, Mo, and Ag or a Gaussian structure for Au. As shown in Figure 4, the relative nugget effect for Cu and Mo is approximately 15 % , and for Ag and Au, it increases to 30 % . The fitting of variogram models in the directions of lower continuity will be explained in the next subsection.
In the same way, indicator variograms are calculated based on the indicator data of Cu, Mo, Ag, and Au in the same 24 directions, and theoretical models are fitted in the direction of maximum continuity, involving a nugget effect and an exponential structure for Cu or a spherical structure for Mo, Ag, and Au. The correlation ranges of these models vary from 55 to 100 m. The relative nugget effect of the Ag model is larger than its continuous structure, yet for Cu, Mo, and Au, the spatial structure is considerable (Figure 5).

3.2. Spatial Anisotropy Modeling

Anisotropy can affect the geostatistical estimation procedure and should, therefore, be modeled as accurately as possible [38,53]. In the Kahang deposit, structural features and geological formations with various rock units cause anisotropy of grade distributions. The characteristics of the direction and dimensions of geostatistical anisotropy ellipsoids are the basis for interpreting spatial changes and identifying geological controls. In this research, as explained in the previous subsection, variograms were computed for 24 directions. The anisotropy parameters, which consist of two anisotropy factors (ratio of major axis range to semi-major axis range and of the major axis range to minor axis range) and three rotation angles (azimuth, plunge, and dip) were determined for the four variables based on the experimental directional variograms of the original grades and indicators (Table 2).

3.3. Ordinary Kriging

Ordinary kriging was applied to estimate the three-dimensional spatial distribution of the ore grades (Cu, Mo, Ag, and Au), based on the variogram analysis results. The horizontal dimensions of the grid were determined by considering the anisotropy modeling and the borehole layout. The size of the elementary blocks was set to 25 m × 25 m × 12.5 m. The block height was associated with mine excavation planning and with the vertical distribution of the ore and waste thicknesses [54,55].
The kriging variance can be used as a criterion for choosing additional samples [17,18,25,56,57,58]. Additional samples are commonly designed away from observation locations or where the kriging variance is maximized [22,55]. Here, the ordinary block kriging standard deviation ( σ o b k ) was calculated and normalized ( σ ^ o b k ) within a range of [0; 1] for each ore variable:
σ ^ o b k = σ o b k min ( σ o b k ) max ( σ o b k ) min ( σ o b k ) ,
where min ( σ o b k ) and max ( σ o b k ) stand for the minimal and maximal kriging standard deviations calculated over all the blocks, respectively. Then, an average of σ ^ o b k for the ore variables (denoted as KSD, the acronym for kriging standard deviation, hereinafter) is calculated in each block (Equation (2)):
K S D = ( σ ^ o b k ) C u + ( σ ^ o b k ) M o + ( σ ^ o b k ) A g + ( σ ^ o b k ) A u 4 .
Drilling in areas where KSD reaches high values will be selected as one of the objectives to plan the configuration of additional boreholes.

3.4. Indicator Kriging

The output of indicator kriging represents an estimate of the probability for the indicator variable to be equal to 1 at the designated point. In other words, it approximates the probability that the point has a grade greater than the threshold. One application of indicator kriging is ore-waste classification [59]. The confidence of the classification result can be measured by the relative difference between the ore and waste probabilities (denoted as P o and P w ):
Confidence   measure ( CM ) = P o P w P o if P o > P w P w P o P w if P w > P o 0 if P o = P w .
Such a confidence measure ranges between 0 and 1. It is 0 when the probabilities of ore and waste are equal, and the sample or the block cannot be specifically assigned to ore or waste. In contrast, when the two probabilities are completely discriminated against each other, i.e., when P o 1 and P w 0 or when P o 0 and P w 1 , the confidence measure is equal to 1.
The confidence measure provides an approach to making decisions about the uncertainty of assigning blocks to the ore or waste classes. In a detailed exploration stage, additional boreholes can be realized in areas where the confidence measure reaches low values. In a multi-element deposit such as eastern Kahang, this measure is computed for each ore element and then averaged in each block.

3.5. Introducing Fuzzy Thresholds and the Ore Membership Degree

A sharp threshold categorizes data into two classes. To reflect uncertainty, one can define a fuzzy threshold instead of a sharp one. The triangular fuzzy number (TFN) is a proper choice for designing a fuzzy threshold due to its simplicity and its membership being one of the most commonly applied in different scientific disciplines [44,56,60,61,62,63]. It defines a ‘fuzzy membership function’ μ f , e for each ore element e (e = Cu, Mo, Ag, or Au):
μ f , e ( g e ) = 0 if g e < a 1 g e a 1 a 2 a 1 if a 1 g e < a 2 a 3 g e a 3 a 2 if a 2 g e < a 3 0 if g e a 3 ,
where g e stands for the estimated grade of element e, and a 1 , a 2 , and a 3 are thresholds (that vary from one element to another) such that a 1 a 2 a 3 .
In the usual approach to defining this fuzzy membership function, the most promising value, a 2 , is equal to the sharp threshold, and the limiting thresholds ( a 1 and a 3 ) are determined based on equal data frequency on both sides of a 2 . In this work, a different approach is proposed, in which the most promising value, the lower, and the upper thresholds are calculated by comparing the classification of actual and estimated values obtained by leave-one-out cross-validation.
Considering the indicator threshold when plotting the scatter of actual values against their estimates creates four combinations of sub-areas as follows: true positive (TP), false negative (FN), true negative (TN), and false positive (FP) [64]. True positive refers to ore data that were correctly labeled by a classifier, and a true negative is the set of waste data that was correctly labeled by the classifier. In the false positive category, the output of estimation is ore, while the data are, in fact, waste. False negative refers to ore data that were misclassified as waste [59,65]. A useful ranking tool for analyzing the performance of a classifier is the confusion matrix, where each row corresponds to the data in an estimated class, while each column corresponds to the data in an actual class [66,67]. The correct classification rate (CCR) can be calculated by summing the correct decisions in the diagonal of the confusion matrix and dividing this sum by the total number of tests [68]:
CCR = Count ( TP + TN ) Count ( TP + TN + FP + FN ) .
For defining the most promising value of the fuzzy membership function, different modes of classification are considered based on a constant sharp threshold for the actual value and dynamic thresholds for the estimated value chosen in a tolerance interval around the sharp threshold. The CCR is calculated for each dynamic threshold to find the maximum CCR and corresponding threshold for allocating a promising value ( a 2 ) to the fuzzy membership function. The lower and upper thresholds ( a 1 and a 3 ) are computed by averaging the estimated values in the FN and FP categories, respectively (Figure 6).
Finally, the ‘ore membership function’ ( μ o , e ) of element e (e = Cu, Mo, Ag, or Au) is introduced based on the thresholds of the fuzzy membership function μ f , e so defined:
μ o , e ( g e ) = 0 if g e < a 1 g e a 1 2 ( a 2 a 1 ) if a 1 g e < a 2 g e + a 3 2 a 2 2 ( a 3 a 2 ) if a 2 g e < a 3 1 if g e a 3 ,
with g e standing for the estimated grade of element e. The ore membership function equals 0 for an estimated grade less than or equal to a 1 . It reaches 0.5 for an estimated grade equal to a 2 and increases to 1 for an estimated grade greater than or equal to a 3 (Figure 6).
As a result of this procedure, the promising values of Cu, Mo, Ag, and Au are found to be 1600, 390, 20, and 0.5 ppm, respectively. In the fuzzy membership function of Cu, the right side is larger than the other side: the average estimated grade of the false positives differs more from the promising value than the average estimated grade of the false negatives. The reverse situation occurs for Ag and Au, for which the left side of the fuzzy membership function is larger than its right side, whereas the fuzzy membership function of Mo is described by a symmetric TFN. Accordingly, the lower threshold a 1 , most promising value a 2 , and upper threshold a 3 of Cu, Mo, Ag, and Au (expressed in ppm) are determined as: ( 1124 , 1600 , 2342 ) , ( 256 , 390 , 538 ) , ( 9 , 20 , 20 ) , and ( 0.03 , 0.5 , 0.5 ) , respectively (Figure 7).

3.6. Ore Value Modeling

The ‘ore value’ function is applied to transform the estimated grade into a relative value of ore variables such as Cu, Mo, Ag, and Au based on the different thresholds. This function determines the total value of each block by considering all the effective variables. It is calculated by the sum of a weight factor ( w e ) times the ore membership degree ( μ o , e ( g e ) ) times the estimated grade in ppm ( g e ) for all the ore elements (e = Cu, Mo, Ag, or Au) in each block. The weight factor for each variable depends on the expected recovery and selling price of the element [44]:
Ore   value ( OV ) = e = 1 4 w e × μ o , e ( g e ) × g e ,
with w C u = 0.155 , w M o = 0.493 , w A g = 18.8 , and w A u = 1030 .

3.7. Misclassification Cost

A misclassification cost can be introduced when there are different types of classification errors and the costs they bring are different [69]. In the assessment of mineral resources and mineral reserves, errors have an impact on the metal recovery and production tonnage: if ore blocks are inadvertently classified as waste, the ore is lost; conversely, if waste blocks are included in ore, the total tonnage increases but the average grade decreases [57,70]. The economic impact of ore loss is the income lost when the ore is sent to the waste dumps and never recovered. This is calculated using the ore value (Equation (7)). The dilution economic effect is imposed on the extra cost of mining and processing waste that is treated as ore [71].
According to the previous statements, the ‘misclassification cost’ (MC) can be defined as follows:
Misclassification   cost   ( MC ) = e = 1 4 w ( 1 μ o , e ( g e ) ) × I ( g e , t e ) + e = 1 4 w e × μ o , e ( g e ) × g e × ( 1 I ( g e , t e ) ) ,
where w is a weight factor that depends on the mining and processing cost per block (in the present case, w = 180 ), t e is the cut-off for element e ( t C u = 1500 ppm, t M o = 400 ppm, t A g = 20 ppm, and t A u = 0.5 ppm), g e is the estimated grade of element e, and I is an indicator function that determines whether a block is classified as ore or as waste: I ( g e , t e ) = 1 if g e t e , 0 otherwise. One goal of detailed exploration is to reduce such a misclassification cost.

3.8. Our Proposal: Multi-Objective Function Index

To design the additional boreholes, an index is introduced to detect locations with a high potential [44,72]. This index has four sub-objectives to identify new boreholes that can improve the accuracy of classification and estimation of mineral resources and reserves:
  • maximize the kriging standard deviation (KSD), reflecting a high uncertainty in the true unknown grades;
  • minimize confidence measure (CM), reflecting a high uncertainty in the ore-waste classification;
  • maximize ore value (OV), reflecting a high potential economic value;
  • maximize misclassification cost (MC), reflecting a high impact of ore loss and waste dilution.
The 3D block models of these four sub-objectives are illustrated in Figure 8a–d, while Figure 8e depicts a 2D map of the average ore value along the vertical column of blocks.
It is common to standardize each term in the same range for a simultaneous application of these objectives in the index. In the present study, to introduce and maximize the index in terms of a fraction, its numerator includes parameters to be maximized (kriging standard deviation, ore value, and misclassification cost), while the confidence measure is situated in the denominator:
Index = KSD × OV × MC CM .

4. Results and Discussion

4.1. Application of the Index to Define Additional Boreholes

Figure 9a displays the 3D spatial distribution of the index within the block model. The index values are averaged along a vertical column of blocks assigned to the block centers’ longitude and latitude and plotted on a 2D map. Then, eight vertical boreholes are located based on large amounts of the index values in the 2D map (Figure 9c). Two potential directional boreholes are also candidates based on the fact that they intersect areas exhibiting large index values that are not located along the same vertical line. The locations and directions of these boreholes are shown in Figure 9b.
The index average and its variations are surveyed along the imaginary axis of the boreholes to determine the length and priority of each additional borehole. A change in the index average tending to zero or becoming negative, indicating that the index value remains constant or trends downward, is the main criterion to calculate the length of the borehole. The priority of boreholes is defined based on comparing the index average value in the calculated length. In Figure 10, the index average and change of index average are shown with a solid smooth line and dash line, respectively, and the crossed dash line (black) represents the borehole length.
Finally, the dimensional and directional properties, such as the header coordinates (x, y, z), length, azimuth, and dip of ten proposed additional boreholes with their priorities, are presented in Table 3.

4.2. Discussion

One focus for discussion relates to the simultaneous selection of additional boreholes, which may introduce some interference effects due to possible information redundancy between these boreholes: drilling ‘twinned’ holes (i.e., two or more holes next to each other) may decrease the target index just a little more than drilling a single hole. Therefore, an extra criterion to select the additional borehole should be a minimal spacing between them. In the present case, the boreholes with priorities 1 and 2 are distant 42.6 m from each other, while the borehole with priority 10 is distant 45.9 m from borehole 4 near the collar (elevation 2275 m) and 28 m from borehole 3 near elevation 2157 m. The remaining additional boreholes are distant more than 55 m from each other, a distance for which the grade variograms (Figure 4) and the indicator variograms (Figure 5) reach more than 90 % of their sills (or 80 % of the sill in the case of the Mo indicator); this means that, except for the twinned holes 1-2, 3-10, and 4-10, there is practically no interference between the selected boreholes. Concerning these twinned holes, their locations may be adjusted, e.g., by moving borehole 2 slightly (say, 10 m) to the south-east in order to decrease its redundancy with borehole 1 while still drilling an area with a large index value (Figure 9c). Moving borehole 10 may not be desirable, as the interest for this borehole mainly lies in its deepest part (Figure 9b), for which the distance with boreholes 3 and 4 is greater than 55 m; hence, the information redundancy is minimal.
As a second focus for discussion, it is seen that the proposed index (Equation (9)) gives the same importance to the four parameters KSD, CM, OV, and MC: an increase of, say, 10 % of KSD has the same effect on the index as an increase of 10 % of OV or MC, or a decrease of 10 % of CM. One could generalize this index to weight these parameters differently, at the cost of greater complexity, in the following fashion:
Index = KSD α × OV β × MC γ CM δ ,
with α , β , γ , δ > 0 . Equivalently, ln ( Index ) = α ln ( KSD ) + β ln ( OV ) + γ ln ( MC ) δ ln ( CM ) . This way, one can assign more importance to one parameter to the detriment of another, depending on the decision-makers’ experience and risk preferences.
The indices (9) and (10) may also be extended to incorporate more sub-objectives in order to account for auxiliary geoscientific information of the subsurface, if available. In particular, lithological, structural, alteration and topographical maps, hyperspectral imaging, as well as ground or airborne magnetic, electrical, gravimetric, and seismic surveys, are often used for mineral potential mapping and are therefore good candidates for the definition of a ‘geological potential’ parameter relevant to the definition of exploratory boreholes (in addition to KSD, CM, OV, and MC). The reader is referred to [73,74,75,76,77,78] for applications of geophysical surveys to the exploration of volcanogenic massive sulfide, iron oxide copper gold, and porphyry ore deposits and to the sampling design problem.
However, it should be noted that the geological and geophysical layers are subject to some degree of subjectivity and to possible cognitive biases, as they result from human decisions (e.g., the choice of a conceptual model of the mineral deposit) or from indirect, incomplete and imperfect information (e.g., geophysical inversion, or interpolated information due to incomplete spatial coverage), which raises the question of the confidence to place in these layers [79,80]. This is the reason why they are used in the prospective stage to propose the primary boreholes in and around the detected anomalies, and their weight is often reduced against existing drilling data when designing the layout of infill boreholes for detailed exploration. According to the level of exploration studies of the eastern Kahang deposit, the geological and geophysical layers cannot play a significant role, and the existing drilling network is more prominent in the design of complementary boreholes.

5. Conclusions

A new approach has been presented to determine the layout of additional boreholes in mineral exploration projects. The proposed method was applied to the eastern part of the Kahang copper deposit in Central Iran. The locations and lengths of eight vertical and two directional complementary drillings were suggested for further exploration of this area. These complementary drillings were selected using the following steps:
  • The average kriging standard deviation quantifies the uncertainty in the true grades and was obtained based on ordinary kriging results.
  • The confidence measure expresses the certainty of the decision to assign blocks to the ore or waste classes. This measure was calculated by a discrimination rate between ore and waste probabilities, which were obtained based on indicator kriging results.
  • Fuzzy thresholds were recommended instead of a sharp threshold. The actual and kriging estimated grades were classified to calculate a confusion matrix and correct classification rate (CCR), which were applied to define fuzzy thresholds for Cu, Mo, Ag, and Au. The center of the triangular fuzzy number (TFN) was determined by maximizing CCR, and the lower and upper limits of the TFN were determined by the average of estimated grades in the false positive (FP) and false negative (FN) categories, respectively. The TFNs of Cu, Mo, Ag, and Au were obtained as ( 1124 , 1600 , 2342 ) ; ( 256 , 390 , 538 ) ; ( 9 , 20 , 20 ) ; and ( 0.03 , 0.5 , 0.5 ) . The ore membership function was then modeled by a quasi-trapezoidal membership function based on the TFN center and limits.
  • The ore value function represents the total value of each block by calculating the sum of a weight factor times the ore membership degree times the estimated grade for all the ore elements. It measures the economic potential of the block.
  • The misclassification cost includes the ore loss and dilution economic impact computed by the income lost when sending ore to the dump and the extra cost of mining and processing waste material.
  • An index for the selection of additional boreholes was finally defined by multiplying the kriging standard deviation, ore value, and misclassification cost and dividing by the confidence measure. By surveying the index average and its change along the axis of the additional boreholes, the properties of eight vertical and two directional additional boreholes were proposed together with their priorities.

Author Contributions

Conceptualization, A.H.-M.; methodology, A.H.-M., F.K. and X.E.; investigation, A.H.-M.; writing—original draft preparation, A.H.-M. and F.K.; writing—review and editing, X.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the National Agency for Research and Development of Chile through postgraduate study scholarship ANID-Subdirección de Capital Humano/Doctorado Nacional/2022-21220317 (F.K.) and grants ANID FONDECYT 1210050 and ANID PIA AFB220002 (F.K. and X.E.). The APC was funded by the National Agency for Research and Development of Chile through grant ANID FONDECYT 1210050 (F.K. and X.E.).

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge three anonymous reviewers for their insightful comments that improved the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Bengough, A.G.; Castrignano, A.; Pages, L.; Noordwijk, M. Sampling strategies, scaling and statistics. In Root Methods; Smit, A.L., Bengough, A.G., Engels, C., Noordwijk, M., Pellerin, S., Geijn, S.C., Eds.; Springer: Berlin, Germany, 2000; pp. 147–173. [Google Scholar]
  2. Bilisoly, R.L.; McKenna, S.A. Determining Optimal Location and Numbers of Sample Transects for Characterization of UXO Sites; SAND Report SAND2002-3962; Sandia National Laboratories: Albuquerque, NM, USA, 2003; 47p.
  3. Vanem, E.; Gramstad, O.; Bitner-Gregersen, E.M. A simulation study on the uncertainty of environmental contours due to sampling variability for different estimation methods. Appl. Ocean Res. 2019, 91, 101870. [Google Scholar] [CrossRef]
  4. Lin, Y.P.; Yeh, M.S.; Deng, D.P.; Wang, Y.C. Geostatistical approaches and optimal additional sampling schemes for spatial patterns and future sampling of bird diversity. Glob. Ecol. Biogeogr. 2008, 17, 175–188. [Google Scholar] [CrossRef]
  5. Rossi, M.E.; Deutsch, C.V. Mineral Resource Estimation; Springer: Dordrecht, The Netherlands, 2014; pp. 1–332. [Google Scholar]
  6. Knödel, K.; Lange, G.; Voigt, H.J. Environmental Geology: Handbook of Field Methods and Case Studies; Springer: Berlin, Germany, 2007; pp. 1–1357. [Google Scholar]
  7. Emery, X.; Séguret, S.A. Geostatistics for the Mining Industry—Applications to Porphyry Copper Deposits; CRC Press: Boca Raton, FL, USA, 2020; pp. 1–247. [Google Scholar]
  8. Savinskii, J.D. Probability Tables for Locating Elliptical Underground Masses with a Rectangular Grid; Springer: New York, NY, USA, 1965; 110p. [Google Scholar]
  9. Singer, D.A.; Wickman, F.E. Probability Tables for Locating Elliptical Targets with Square, Rectangular, and Hexagonal Point-Nets; The Pennsylvania State University: State College, PA, USA, 1969; 100p. [Google Scholar]
  10. Hossein Morshedy, A.; Memarian, H. A novel algorithm for designing the layout of additional boreholes. Ore Geol. Rev. 2015, 67, 34–42. [Google Scholar] [CrossRef]
  11. Colombera, L.; Mountney, N.P.; Russell, C.E.; Shiers, M.N.; McCaffrey, W.D. Geometry and compartmentalization of fluvial meander-belt reservoirs at the bar-form scale: Quantitative insight from outcrop, modern and subsurface analogues. Mar. Pet. Geol. 2017, 82, 35–55. [Google Scholar] [CrossRef]
  12. Oh, I.H.; Yang, S.J.; Heo, C.H.; Lee, C.; Jin, K.; Lee, J.H.; Jets, S. Applying three-dimensional modeling based on trench surveys to chromitite orebodies in the BophiVum area, northwestern Myanmar. Geosci. J. 2021, 25, 333–349. [Google Scholar] [CrossRef]
  13. Drew, L.J. Pattern drilling exploration: Optimum pattern types and hole spacings when searching for elliptical shaped targets. Math. Geol. 1979, 11, 223–254. [Google Scholar] [CrossRef]
  14. Bras, R.L.; Rodríguez-Iturbe, I. Network design for the estimation of areal mean of rainfall events. Water Resour. Res. 1976, 12, 1185–1195. [Google Scholar] [CrossRef]
  15. Scheck, D.E.; Chou, D.R. Optimum locations for exploratory drill holes. Int. J. Min. Eng. 1983, 1, 343–355. [Google Scholar] [CrossRef]
  16. McBratney, A.B.; Webster, R.; Burgess, T.M. The design of optimal sampling schemes for local estimation and mapping of regionalized variables—I: Theory and method. Comput. Geosci. 1981, 7, 331–334. [Google Scholar] [CrossRef]
  17. Olea, R.A. Sampling design optimization for spatial functions. Math. Geol. 1984, 16, 369–392. [Google Scholar] [CrossRef]
  18. Barnes, R.J. Sample design for geologic site characterization. In Geostatistics; Armstrong, M., Ed.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1989; pp. 809–822. [Google Scholar]
  19. Szidarovszky, F. Multiobjective observation network design for regionalized variables. Int. J. Min. Eng. 1983, 1, 331–342. [Google Scholar] [CrossRef]
  20. Gershon, M.; Allen, L.E.; Manley, G. Application of a new approach for drillholes location optimization. Int. J. Surf. Min. Reclam. Environ. 1998, 2, 27–31. [Google Scholar] [CrossRef]
  21. Van Groenigen, J.W.; Siderius, W.; Stein, A. Constrained optimisation of soil sampling for minimisation of the kriging variance. Geoderma 1999, 87, 239–259. [Google Scholar] [CrossRef]
  22. Delmelle, E.M.; Goovaerts, P. Second-phase sampling designs for non-stationary spatial variables. Geoderma 2009, 153, 205–216. [Google Scholar] [CrossRef] [PubMed]
  23. Fatehi, M.; Asadi Haroni, H.; Hossein Morshedy, A. Designing infill directional drilling in mineral exploration by using particle swarm optimization algorithm. Arab. J. Geosci. 2017, 10, 487. [Google Scholar] [CrossRef]
  24. Pactwa, K.; Wozniak, J.; Strempski, A. Sustainable mining–Challenge of Polish mines. Resour. Policy 2017, 74, 101269. [Google Scholar] [CrossRef]
  25. Kleijnen, J.P. Design and Analysis of Simulation Experiments; Springer International Publishing: Dordrecht, The Netherlands, 2018; pp. 1–216. [Google Scholar]
  26. Mohammadi, S.S.; Hezarkhani, A.; Tercan, A.E. Optimally locating additional drill holes in three dimensions using grade and simulated annealing. J. Geol. Soc. India 2012, 80, 700–706. [Google Scholar] [CrossRef]
  27. Jafrasteh, B.; Fathianpour, N. Optimal location of additional exploratory drillholes using a fuzzy-artificial bee colony algorithm. Arab. J. Geosci. 2017, 10, 213. [Google Scholar] [CrossRef]
  28. Dutaut, R.; Marcotte, D. A new semi-greedy approach to enhance drillhole planning. Nat. Resour. Res. 2020, 29, 3599–3612. [Google Scholar] [CrossRef]
  29. Kumral, M.; Ozer, U. Planning additional drilling campaign using two-space genetic algorithm: A game theoretical approach. Comput. Geosci. 2013, 52, 117–125. [Google Scholar] [CrossRef]
  30. Zagré, G.E.; Marcotte, D.; Gamache, M.; Guibault, F. New tabu algorithm for positioning mining drillholes with blocks uncertainty. Nat. Resour. Res. 2019, 28, 609–629. [Google Scholar] [CrossRef]
  31. Soltani-Mohammadi, S.; Safa, M.; Sohrabian, B. Reducing uncertainty in mineralization boundary by optimally locating additional drill holes through particle swarm optimization. Nat. Resour. Res. 2021, 30, 2067–2083. [Google Scholar] [CrossRef]
  32. Soltani-Mohammadi, S.; Safa, M.; Mokhtari, H. Comparison of particle swarm optimization and simulated annealing for locating additional boreholes considering combined variance minimization. Comput. Geosci. 2016, 95, 146–155. [Google Scholar] [CrossRef]
  33. Safa, M.; Soltani-Mohammadi, S. Distance function modeling in optimally locating additional boreholes. Spat. Stat. 2018, 23, 17–35. [Google Scholar] [CrossRef]
  34. Pinheiro, M.; Emery, X.; Rocha, A.M.; Miranda, T.; Lamas, L. Boreholes plans optimization methodology combining geostatistical simulation and simulated annealing. Tunn. Undergr. Space Technol. 2017, 70, 65–75. [Google Scholar] [CrossRef]
  35. Chen, Y.C.; Chiang, W.; Yeh, H.C. Rainfall network design using kriging and entropy. Hydrol. Process. 2008, 22, 340–346. [Google Scholar] [CrossRef]
  36. Soltani, S.; Hezarkhani, A. Determination of realistic and statistical value of the information gathered from exploratory drilling. Nat. Resour. Res. 2011, 20, 207–216. [Google Scholar] [CrossRef]
  37. Caers, J.; Scheidt, C.; Yin, Z.; Wang, L.; Mukerji, T.; House, K. Efficacy of information in mineral exploration drilling. Nat. Resour. Res. 2022, 31, 1157–1173. [Google Scholar] [CrossRef]
  38. Hassanipak, A.A.; Sharafodin, M. GET: A function for preferential site selection of additional borehole drilling. Explor. Min. Geol. 2004, 13, 139–146. [Google Scholar] [CrossRef]
  39. Juang, K.W.; Lee, D.Y.; Chen, Z.S. Geostatistical cross-validation for the design of additional sampling regimes in heavy-metal contaminated soils. J. Chin. Inst. Environ. Eng. 1999, 9, 89–96. [Google Scholar]
  40. Pilger, G.G.; Costa, J.F.C.L.; Koppe, J.C. Additional samples: Where they should be located. Nat. Resour. Res. 2001, 10, 197–207. [Google Scholar] [CrossRef]
  41. Hernández, J.; Emery, X. A geostatistical approach to optimize sampling designs for local forest inventories. Can. J. For. Res. 2009, 39, 1465–1474. [Google Scholar] [CrossRef]
  42. Juang, K.W.; Liao, W.J.; Liu, T.L.; Tsui, L.; Lee, D.Y. Additional sampling based on regulation threshold and kriging variance to reduce the probability of false delineation in a contaminated site. Sci. Total Environ. 2008, 389, 20–28. [Google Scholar] [CrossRef] [PubMed]
  43. Van Groenigen, J.W.; Pieters, G.; Stein, A. Optimizing spatial sampling for multivariate contamination in urban areas. Environmetrics 2000, 11, 227–244. [Google Scholar] [CrossRef]
  44. Hossein Morshedy, A.; Torabi, S.A.; Memarian, H. A new method for 3D designing of complementary exploration drilling layout based on ore value and objective functions. Arab. J. Geosci. 2015, 8, 8175–8195. [Google Scholar] [CrossRef]
  45. Aspie, D.; Barnes, R.J. Infill-sampling design and the cost of classification errors. Math. Geol. 1990, 22, 915–932. [Google Scholar] [CrossRef]
  46. Silva, D.S.; Jewbali, A.; Boisvert, J.B.; Deutsch, C.V. Drillhole placement subject to constraints for improved resource classification. CIM J. 2018, 9, 21–32. [Google Scholar] [CrossRef]
  47. Arfania, R.; Shahriari, S. Role of southeastern Sanandaj–Sirjan zone in the tectonic evolution of Zagros Orogenic Belt, Iran. Isl. Arc 2009, 18, 555–576. [Google Scholar] [CrossRef]
  48. Shafiei, B.; Haschke, M.; Shahabpour, J. Recycling of orogenic arc crust triggers porphyry Cu mineralization in Kerman Cenozoic arc rocks, southeastern Iran. Miner. Depos. 2009, 44, 265–283. [Google Scholar] [CrossRef]
  49. Haschke, M.; Ahmadian, J.; Murata, M.; McDonald, I. Copper mineralization prevented by arc-root delamination during Alpine-Himalayan collision in central Iran. Econ. Geol. 2010, 105, 855–865. [Google Scholar] [CrossRef]
  50. National Iranian Copper Industries Co. Summary Geological Report of Kahang Deposit; National Iranian Copper Industries Co.: Tehran, Iran, 2011; 23p. (In Persian) [Google Scholar]
  51. Afzal, P.; Alghalandis, Y.F.; Moarefvand, P.; Omran, N.R.; Haroni, H.A. Application of power-spectrum–volume fractal method for detecting hypogene, supergene enrichment, leached and barren zones in Kahang Cu porphyry deposit, Central Iran. J. Geochem. Explor. 2012, 112, 131–138. [Google Scholar] [CrossRef]
  52. Afshooni, S.Z.; Mirnejad, H.; Esmaeily, D.; Haroni, H.A. Mineral chemistry of hydrothermal biotite from the Kahang porphyry copper deposit (NE Isfahan), Central Province of Iran. Ore Geol. Rev. 2013, 54, 214–232. [Google Scholar] [CrossRef]
  53. Wackernagel, H. Multivariate Geostatistics: An Introduction with Applications; Springer: Berlin, Germany, 2003; pp. 1–387. [Google Scholar]
  54. Sun, T.; Chen, F.; Zhong, L.; Liu, W.; Wang, Y. GIS-based mineral prospectivity mapping using machine learning methods: A case study from Tongling ore district, eastern China. Ore Geol. Rev. 2019, 109, 26–49. [Google Scholar] [CrossRef]
  55. Joshi, D.; Paithankar, A.; Chatterjee, S.; Equeenuddin, S.M. Integrated parametric graph closure and branch-and-cut algorithm for open pit mine scheduling under uncertainty. Mining 2022, 2, 32–51. [Google Scholar] [CrossRef]
  56. Hossein Morshedy, A.; Torabi, S.A.; Memarian, H. A hybrid fuzzy zoning approach for 3-dimensional exploration geotechnical modeling: A case study at Semilan dam, southern Iran. Bull. Eng. Geol. Environ. 2019, 78, 691–708. [Google Scholar] [CrossRef]
  57. Sinclair, A.J.; Blackwell, G.H. Applied Mineral Inventory Estimation; Cambridge University Press: London, UK, 2002; pp. 1–381. [Google Scholar]
  58. Brus, D.J.; Heuvelink, G.B.M. Optimization of sample patterns for universal kriging of environmental variables. Geoderma 2007, 138, 86–95. [Google Scholar] [CrossRef]
  59. Shamseddin-Meigooni, M.; Lotfi, M.; Afzal, P.; Nezafati, N.; Razi, M.K. Application of multivariate geostatistical simulation and fractal analysis for detection of rare-earth element geochemical anomalies in the Esfordi phosphate mine, Central Iran. Geochem. Explor. Environ. Anal. 2021, 21, geochem2020-035. [Google Scholar] [CrossRef]
  60. Lee, K.H. First Course on Fuzzy Theory and Applications; Springer: Berlin, Germany, 2005; pp. 1–336. [Google Scholar]
  61. Taboada, J.; Ordóñez, C.; Saavedra, A.; Fiestras-Janeiro, G. Fuzzy expert system for economic zonation of an ornamental slate deposit. Eng. Geol. 2006, 84, 220–228. [Google Scholar] [CrossRef]
  62. Ross, T.J. Fuzzy Logic with Engineering Applications; John Wiley & Sons: Chichester, UK, 2010; pp. 1–585. [Google Scholar]
  63. Kambalimath, S.; Deka, P.C. A basic review of fuzzy logic applications in hydrology and water resources. Appl. Water Sci. 2020, 10, 191. [Google Scholar] [CrossRef]
  64. Han, J.; Kamber, M.; Pei, J. Data Mining: Concepts and Techniques; Morgan Kaufmann: Burlington, VT, USA, 2006; pp. 1–744. [Google Scholar]
  65. Moon, C.; Whateley, M.; Evans, A.M. Introduction to Mineral Exploration; Blackwell Publishing: Malden, MA, USA, 2006; pp. 1–481. [Google Scholar]
  66. Aytaç, T.; Barshan, B. Simultaneous extraction of geometry and surface properties of targets using simple infrared sensors. Opt. Eng. 2004, 43, 2437–2447. [Google Scholar] [CrossRef]
  67. Zapata, J.; Vilar, R.; Ruiz, R. An adaptive-network-based fuzzy inference system for classification of welding defects. NDT E Int. 2010, 43, 191–199. [Google Scholar] [CrossRef]
  68. Mitchell, H.B. Data Fusion: Concepts and Ideas; Springer: Berlin, Germany, 2012; pp. 1–331. [Google Scholar]
  69. Kaur, H.; Pannu, H.S.; Malhi, A.K. A systematic review on imbalanced data challenges in machine learning: Applications and solutions. ACM Comput. Surv. 2019, 52, 1–36. [Google Scholar] [CrossRef]
  70. Khorram, F.; Asghari, O.; Memarian, H.; Hossein Morshedy, A.; Emery, X.M. Adaptive multi-size block modeling for mineral resources and ore reserves evaluation. J. Min. Environ. 2021, 12, 339–350. [Google Scholar]
  71. Rogers, W.D.; Kanchibotla, S. Application of stochastic approach to predict blast movement. In Rock Fragmentation by Blasting; Singh, P.K., Sinha, A., Eds.; CRC Press: Boca Raton, FL, USA, 2013; pp. 257–265. [Google Scholar]
  72. Fatehi, M.; Asadi-Haroni, H.; Hossein Morshedy, A. 3D design of optimum complementary boreholes by integrated analysis of various exploratory data using a sequential-MADM approach. Nat. Resour. Res. 2020, 29, 1041–1061. [Google Scholar] [CrossRef]
  73. Abedi, M.; Norouzi, G.H. Integration of various geophysical data with geological and geochemical data to determine additional drilling for copper exploration. J. Appl. Geophys. 2012, 83, 35–45. [Google Scholar] [CrossRef]
  74. Kashani, S.B.M.; Abedi, M.; Norouzi, G.H. Fuzzy logic mineral potential mapping for copper exploration using multi-disciplinary geo-datasets, a case study in Seridune deposit, Iran. Earth Sci. Inform. 2016, 9, 167–181. [Google Scholar] [CrossRef]
  75. White, D.J.; Malinowski, M.; Devine, C.; Gilmore, K.; Schetselaar, E.; Pehrsson, S. Drill targeting with 3-D seismics for volcanogenic massive sulfide exploration in the Flin Flon mining camp. Econ. Geol. 2016, 111, 903–912. [Google Scholar] [CrossRef]
  76. Karbalaei-Ramezanali, A.A.; Mansouri, E.; Feizi, F. Integration of aeromagnetic geophysical data with other exploration data layers based on fuzzy AHP and C-A fractal model for Cu-porphyry potential mapping: A case study in the Fordo area, central Iran. Boll. Geofis. Teor. Appl. 2017, 58, 55–73. [Google Scholar]
  77. Feizi, F.; Karbalaei-Ramezanali, A.A.; Farhadi, S. Application of multivariate regression on magnetic data to determine further drilling site for iron exploration. Open Geosci. 2021, 13, 138–147. [Google Scholar] [CrossRef]
  78. Moreno, L.; Contreras, L.; Pavez, A.; Emery, X.; Cruz, A.; Cáceres, A. Stochastic delineation of potential sulfide exploration targets of the Sierra Norte project. In Proceedings of the 7th International Conference on Geology and Mine Planning, Online, 9–11 June 2021; Gecamin Ltda: Santiago, Chile, 2021; pp. 107–109. [Google Scholar]
  79. Bond, C.E.; Gibbs, A.D.; Shipton, Z.K.; Jones, S. What do you think this is? “Conceptual uncertainty” in geoscience interpretation. GSA Today 2007, 17, 4–10. [Google Scholar] [CrossRef]
  80. Bianchi, M.; Turner, A.K.; Lark, M.; Courrioux, G. Uncertainty in 3-D geological models. In Applied Multidimensional Geological Modeling: Informing Sustainable Human Interactions with the Shallow Subsurface; Turner, A.K., Kessler, H., van der Meulen, M.J., Eds.; John Wiley & Sons: Chichester, UK, 2021; pp. 357–382. [Google Scholar]
Figure 1. A flowchart of the proposed method for designing additional boreholes in mineral exploration.
Figure 1. A flowchart of the proposed method for designing additional boreholes in mineral exploration.
Minerals 13 01252 g001
Figure 2. (a) Map of Iran and the location of the eastern Kahang copper deposit (blue square); (b) geological map of the deposit. Adapted from [10].
Figure 2. (a) Map of Iran and the location of the eastern Kahang copper deposit (blue square); (b) geological map of the deposit. Adapted from [10].
Minerals 13 01252 g002
Figure 3. Three-dimensional distribution of ore grades (in ppm) in the composite samples: (a) Cu and (b) Mo.
Figure 3. Three-dimensional distribution of ore grades (in ppm) in the composite samples: (a) Cu and (b) Mo.
Minerals 13 01252 g003
Figure 4. Experimental variograms (dots) with fitted models (solid lines) in the main direction of continuity: (a) Cu, (b) Mo, (c) Ag, and (d) Au.
Figure 4. Experimental variograms (dots) with fitted models (solid lines) in the main direction of continuity: (a) Cu, (b) Mo, (c) Ag, and (d) Au.
Minerals 13 01252 g004
Figure 5. Experimental indicator variograms (dots) with fitted models (solid lines) in the main direction of continuity: (a) Cu, (b) Mo, (c) Ag, and (d) Au.
Figure 5. Experimental indicator variograms (dots) with fitted models (solid lines) in the main direction of continuity: (a) Cu, (b) Mo, (c) Ag, and (d) Au.
Minerals 13 01252 g005
Figure 6. Graphical description of the definition fuzzy and ore membership functions.
Figure 6. Graphical description of the definition fuzzy and ore membership functions.
Minerals 13 01252 g006
Figure 7. Fuzzy and ore membership functions of Cu, Mo, Ag, and Au.
Figure 7. Fuzzy and ore membership functions of Cu, Mo, Ag, and Au.
Minerals 13 01252 g007
Figure 8. The 3D modeling of (a) normalized kriging standard deviation (KSD), (b) confidence measure (CM), (c) ore value (OV), (d) misclassification cost (MC), (e) and 2D map of ore value based on averaging along vertical columns of blocks. In subfigures (c,d), the zero values of OV and MC are filtered to emphasize the mineable regions of the deposit and the regions likely to be affected by ore loss and waste dilution.
Figure 8. The 3D modeling of (a) normalized kriging standard deviation (KSD), (b) confidence measure (CM), (c) ore value (OV), (d) misclassification cost (MC), (e) and 2D map of ore value based on averaging along vertical columns of blocks. In subfigures (c,d), the zero values of OV and MC are filtered to emphasize the mineable regions of the deposit and the regions likely to be affected by ore loss and waste dilution.
Minerals 13 01252 g008
Figure 9. (a) The 3D spatial distribution of index (9), (b) location of directional boreholes with vertical cross-sections showing the index values along the boreholes, (c) and the location of additional boreholes with their priority according to the 2D index average.
Figure 9. (a) The 3D spatial distribution of index (9), (b) location of directional boreholes with vertical cross-sections showing the index values along the boreholes, (c) and the location of additional boreholes with their priority according to the 2D index average.
Minerals 13 01252 g009
Figure 10. Illustration of index average (solid red lines) and its change (dashed blue lines) along the borehole axes, with their priority (from 1 to 10). The dashed black lines indicate the proposed borehole lengths, from which the index change is close to zero.
Figure 10. Illustration of index average (solid red lines) and its change (dashed blue lines) along the borehole axes, with their priority (from 1 to 10). The dashed black lines indicate the proposed borehole lengths, from which the index change is close to zero.
Minerals 13 01252 g010
Table 1. Descriptive statistics of copper, molybdenum, silver, and gold grades in the composite samples.
Table 1. Descriptive statistics of copper, molybdenum, silver, and gold grades in the composite samples.
ParameterCu (ppm)Mo (ppm)Ag (ppm)Au (ppm)
Number of data6205576051694657
Minimum0.0000.0000.0000.000
Maximum49,200147931.80.789
Mean1800.933.761.1730.013
Mode0.0001.50.60.005
Quartile 13403.30.360.005
Median890.510.80.60.005
Quartile 3198826.381.30.015
St. deviation3209.588.731.7510.028
Table 2. Anisotropy parameters of Cu, Mo, Ag, and Au for original grades and indicators.
Table 2. Anisotropy parameters of Cu, Mo, Ag, and Au for original grades and indicators.
CaseVariableAnisotropy Factor (Major/Semimajor)Anisotropy Factor (Major/Minor)Azimuth ( )Plunge ( )Dip ( )
OriginalCu1.171.32266.947.440.5
dataMo1.351.78260.370.00.74
Ag1.311.93264.060.156.6
Au1.271.45263.861.129.7
IndicatorCu1.161.4231.142.662.3
dataMo1.271.6552.620.619.3
Ag1.121.2291.680.046.9
Au1.151.2758.619.328.5
Table 3. Dimensional and directional properties of proposed boreholes with higher priority.
Table 3. Dimensional and directional properties of proposed boreholes with higher priority.
PriorityCollar X (m)Collar Y (m)Collar Z (m)Azimuth ( )Dip ( )Length (m)Average Index Value
1638,484.93,644,90522690904602.146
2638,511.83,644,87222690903402.046
3638,638.73,644,7842272.40904802.011
4638,680.23,644,7452275.70905302.002
5638,610.73,644,7082272.10904801.985
6638,437.33,644,86122690904001.960
7638,604.93,644,5992283.129.575.75201.928
8638,437.33,644,56322870904001.925
9638,106.53,644,9112274.80902101.881
10638,694.33,644,7882284.4296.169.25101.818
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hossein-Morshedy, A.; Khorram, F.; Emery, X. A Multi-Objective Approach for Optimizing the Layout of Additional Boreholes in Mineral Exploration. Minerals 2023, 13, 1252. https://0-doi-org.brum.beds.ac.uk/10.3390/min13101252

AMA Style

Hossein-Morshedy A, Khorram F, Emery X. A Multi-Objective Approach for Optimizing the Layout of Additional Boreholes in Mineral Exploration. Minerals. 2023; 13(10):1252. https://0-doi-org.brum.beds.ac.uk/10.3390/min13101252

Chicago/Turabian Style

Hossein-Morshedy, Amin, Farzaneh Khorram, and Xavier Emery. 2023. "A Multi-Objective Approach for Optimizing the Layout of Additional Boreholes in Mineral Exploration" Minerals 13, no. 10: 1252. https://0-doi-org.brum.beds.ac.uk/10.3390/min13101252

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