Next Article in Journal
Study towards Integrating Secondary Forests into Future Land Use Development in Singapore for Biodiversity Objectives
Next Article in Special Issue
Open-Pit Mine Truck Dispatching System Based on Dynamic Ore Blending Decisions
Previous Article in Journal
Regional Differences in Tourism Eco-Efficiency in the Beijing–Tianjin–Hebei Region: Based on Data from 13 Cities
Previous Article in Special Issue
Slope Stability Prediction Method Based on Intelligent Optimization and Machine Learning Algorithms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Production Strategy Optimization of Integrated Exploitation for Multiple Deposits Considering Carbon Quota

1
School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China
2
Key Laboratory of Ministry of Education for Efficient Mining and Safety of Metal Mines, University of Science and Technology Beijing, Beijing 100083, China
3
Anhui Mine IOT and Security Monitoring Technology Key Laboratory, Hefei 230601, China
4
Shandong Gold Group Co., Ltd, Jinan 250101, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(4), 2917; https://doi.org/10.3390/su15042917
Submission received: 11 January 2023 / Revised: 28 January 2023 / Accepted: 1 February 2023 / Published: 6 February 2023
(This article belongs to the Special Issue Advances in Intelligent and Sustainable Mining)

Abstract

:
Nowadays, the mining industry actively advocates and practices the concept of green and integrated exploitation to realize the sustainable development of resources with low-carbon emissions. The certain carbon quota for mining companies limits the production capacity and resource utilization efficiency. The integrated exploitation of multiple deposits could coordinate resource allocation and operation facilities, which would reduce capital expenditure and operating costs for the mining company from a systematic perspective. In this condition, some deposits located nearby could be treated as one entity to make plans and optimize. An optimization framework is proposed based on the analysis of the characteristics and advantages of integrated exploitation. A new mathematical programming model is presented to optimize production capacity and extracted ore grade for each deposit considering constraints of maximum and minimum mining capacity, extracted ore grade and concentrated ore grade requirement, and metal output target, which has a significant influence on the economic benefit and resource utilization rate for a mining company. The model is verified using the data collected from three deposits of a gold mining company in China to demonstrate its ability to optimize the allocation of production capacity and improve the technical and economic effect of mining under the limitation of carbon quota. The sensitivity analysis of some key parameters is carried out to generate a series of integrated exploitation schemes under different production and operation conditions, which is useful for the mining company to make decisions in different situations.

1. Introduction

High carbon emissions represent a serious problem for the mining industry worldwide [1,2]. A large number of carbon emissions produced by the exhaustive exploitation of mineral resources bring a burden on the environment. Especially in the context of global warming, mining companies have to consider energy conservation and emission reduction in the mining process with the impact of carbon quota and ever-increasing carbon tax [3,4,5].
The large-scale exploitation of mineral resources leads to serious carbon emissions, land destruction, and environmental damage in the process of mining, which is not conducive to the sustainable exploitation of mineral resources [6,7]. The independent exploitation of each deposit is prone to the problems of rapid resource consumption, the high energy consumption of equipment, and low production efficiency. The mineral resources exploitation method is changing from the independent exploitation of a single deposit to the integrated exploitation of multiple deposits simultaneously under the guidance of efficient resource exploitation and low-carbon economy strategy.
Integrated exploitation methods could accelerate collaborative production among multiple deposits and realize the large-scale extraction and processing of mineral resources [8]. The complex interaction influence among multiple deposits should be considered comprehensively compared with the conventional exploitation of a single deposit. The mining method, production organization, operation cost, and extracted ore grade of each orebody differ due to the different characteristics of each deposit. All deposits supply mineral resources together, and they compete for investment and equipment as well. Therefore, the optimal utilization of multiple deposits could be in the way of production capacity allocation, extracted ore grade control, and mining sequence arrangement to maximize the overall economic benefit [9,10].
Each mining company has a certain carbon quota. The advantage of a high-benefit orebody could not be exploited if the carbon quota was allocated to each deposit equally. Therefore, the scientific allocation of production capacity for each deposit under the constraint of a carbon quota to obtain an optimal economic decision scheme in the integrated exploitation of mineral resources is very important for mining companies [11,12].
A theoretical system of integrated exploitation for mineral resources was proposed with optimization methods and models, which facilitates the optimization of a comprehensive mining system and has practical application in multiple deposits [13,14]. Several methods and models have been presented to optimize the mine planning and production process in the context of carbon emission [15,16]. Ivanova et al. analyzed the carbon footprint generated by the increase in carbon emissions in mining regions due to coal mining and production. It was of great significance to reduce the negative influence of mining on the environment [17]. Valderrama et al. presented a model that integrated ore grade and carbon cost in the whole mining process for the first time to reduce greenhouse gas emissions and operating costs. Furthermore, the proposed model was verified in an iron mine, which suggested that the ore grade had a direct effect on costs and carbon emissions [18]. Azhar et al. proposed a new mine planning approach to realize maximum profits when obeying the given carbon emission target. The efficient approach considered resource constraints in mining and promoted carbon-aware mining practices [19].
Low-carbon technologies effectively obtained practical applications in the integrated exploitation of mineral resources. Boubault and Maïzi pointed out that the widespread application of low-carbon technologies would lead to dynamic adjustments in mineral product demand and supply against a low-carbon-economy background [20]. Large amounts of tailings and waste generated from the integrated exploitation of deposits are still valuable with the development of new technology. The expansion of low-carbon technologies and production increase the recycling of resources. The efficient utilization of low-grade ores and the optimization of low-carbon mining have become the main goals of sustainable mining for deposits [21].
Mathematical programming has been used widely to calculate the optimum value of the objective function on the basis that the decision variables are under certain constraints [22]. It has the characteristics of multi-disciplinary integration, rigorous theory, and extensive application [23]. Currently, mathematical programming has obtained some practice in resource allocation and production planning. Several resource allocation models are applied to optimize mine production and improve economic efficiency. A fuzzy linear programming optimization model of multiple deposit mining systems was presented to seek the optimal plan of production, which was applied in the Bauxite Basin Niksic with five deposits [24]. Rahmanpour and Osanloo formulated a fuzzy linear programming model to optimize short-term mine production planning. The optimal scheme was provided to determine the amount of material that needs to be mined from each site [25]. Gligorić et al. implemented a mathematical model in a small hypothetical lead–zinc mineral deposit, which demonstrated that it could define the sequence of mining of different areas by quantifying the metal price and operating costs of mining [26]. Ahmadi and Shahabi developed a grade optimization model to optimize mine production planning. The output results of the model figured out the optimum grade was beneficial to maximizing profits, which gave a feasible vision to determine the cutoff grade [27].
Additionally, some scholars established multi-objective programming to handle mineral resource allocation and management problems. The improvement of comprehensive optimization for the mine production process was advocated to promote production efficiency and achieve economic and environmental objectives [28,29,30]. However, most optimization models do not consider the constraint of carbon quota and production capacity for the integrated exploitation of mineral resources. Profits are closely related to the resource allocation of all deposits under the limitation of carbon emission. Therefore, the overall economic benefit would not be optimal if the production coordination and competition among multiple deposits are ignored.
In the context of carbon quota, the production strategy optimization of integrated exploitation for multiple deposits is meaningful for a mining company. The integrated exploitation theory and mathematical planning method are introduced to propose an optimization framework and modeling assumptions. Then, this paper presents a new mathematical model for optimizing the allocation of production capacity and technical and economic indexes for multiple deposits. Constraints of maximum and minimum production capacity, extracted ore grade and concentrated ore grade requirement, and metal output target are considered. The model is verified in three deposits of a gold mining company in China. Eventually, sensitivity analysis is conducted to support dynamic decision making for multiple deposits. The detailed research content is shown in Figure 1.

2. Integrated Exploitation of Mineral Resources Considering Carbon Quota

Reserve information from a gold mining company in China is collected, and the production process of each deposit is determined. Furthermore, the advantage of the integrated exploitation of resources is analyzed in detail to maximize the competitive advantage of multiple deposits. Simultaneously, a decision-making framework is constructed to guide the optimization of the enterprise’s production planning.

2.1. Background

A gold mining company in China has abundant geological resources, a high mechanization degree, advanced mining technology, and a management level, which owns three deposits: I, II, and III. All deposits have been designed and mined for several years. The carbon emission has been increasing sharply as the mining scale of each deposit has continued to expand. Today, the advantages of rich geological resources are no longer prominent as the shallow mineral resources have been rapidly consumed. The gold mining company has formulated a strategy-oriented integrated exploitation plan to relieve the pressure of resource sustainable development and reduce carbon emissions.
Table 1 shows the reserves of three deposits at the end of 2021. The reserve of deposit I is the highest. Deposit III has the lowest reserves. The average grade of deposit II is the highest, and deposit III has the lowest grade.
From the perspective of the overall production process, three deposits have a complete production process including mining, mineral processing, and smelting. Each deposit has its mining processes and ore processing independently, while the three deposits share one smelting plant to improve the overall recovery efficiency. The extracted ore is processed in the respective mill and then transported to the shared smelting plant, as shown in Figure 2, which implies that the deposits differ in terms of resource reserves, production capacity, gold throughput, mining cost, and processing cost.

2.2. Advantage Analysis of Integrated Exploitation of Resources

The integrated exploitation of resources is typical of large-scale, integrated planning and manufacturing management, the coexistence of internal synergy and competition, and decision-making consistency [31,32,33]. The integrated exploitation of multiple deposits is regarded as carrying out production decision making from a global perspective instead of only considering the economic benefit of individual deposits, which is consistent with sustainable mineral resource management and development in the context of a low-carbon economy [34]. Figure 3 demonstrates the mine planning optimization strategy for integrated exploitation in multiple deposits. The strategy is conducive to the coordinated management of all deposits, the realization of the optimum goal for integrated exploitation, and the advantage of large-scale performance.
There are obvious synergetic and competitive relationships between different deposits in the mode of integrated exploitation. The carbon quota for all deposits is fixed, which results in the overall production capacity being limited [35,36]. The optimal allocation of production capacity for each deposit needs to take the technical and economic conditions into account. In this condition, the overall economic benefit for the mining company could be improved under the carbon emission constraint. The advantage of the integrated exploitation of resources is analyzed from three aspects shown as follows [37,38,39,40]:
(1)
Resource advantage
Resource-based enterprises should have sufficient resource conditions in order to survive and develop better in fierce market competition. The advantage for enterprises with richer resources is more obvious when the demand for resources is greater in the market. Mining companies could make full use of the resource advantage to avoid the waste of resources by way of the integrated exploitation of resources. The ore from high-grade deposits could be mixed with the ore from lower-grade deposits in smelt processing.
(2)
Cost advantage
The integrated exploitation of mineral resources takes overall economic efficiency as the target. The mining company should attempt to optimize production to realize the optimal economic benefit of multiple deposits, which is the core competitiveness of integrated exploitation.
(3)
Scale performance advantage
Considering the coexistence of production coordination and competition, the advantage of production and the technical and economic indexes of each deposit should be fully utilized to improve the comprehensive utilization efficiency of resources and the overall economic benefit of multiple deposits. In addition, the competitive advantage of economies with large-scale mining should be brought into play.

2.3. Decision-Making Framework

The optimum decision-making framework for the integrated exploitation of multiple deposits is proposed based on the advantage analysis of integrated exploitation. The optimization problem of multiple deposits’ simultaneous exploitation could be transformed into a mathematical model. Resources, economy, and technology are considered as constraints in the mathematical model to obtain the optimal solution, which eventually forms an optimal scheme of production capacity allocation for each deposit. The optimum decision-making framework is given in Figure 4, which is beneficial for all deposits.

3. Optimization Model Construction

Modeling assumptions are formed in this section. Then, the optimization model is constructed to achieve optimum production capacity allocation and improve technical and economic indexes of integrated exploitation.

3.1. Assumptions

Mathematical modeling is performed to solve the problem of production capacity allocation for multiple deposits’ integrated exploitation, which considers the factors of annual sales revenue, mining cost, mineral processing, and smelting cost, cost of smelted ore transported to a smelting plant, reserve information, technical and economical indexes, and the production scale limitation. Eventually, an optimization model of production capacity allocation will be constructed using mathematical planning methods. The basic ideas and assumptions of the model construction are as follows:
(1)
The production capacity assigned to the deposit is adjustable
There are differences in resource reserves, technical and economic levels, and production capacity among all deposits. The production capacity of one deposit is often related to the mining scale, the capital expenditure, the supply of raw materials, and the production technology. The production capacity is limited against the background of low-carbon mining. The optimum production capacity allocation scheme needs to be calculated through the optimization model instead of being determined solely from a single deposit in the mining plan, which aims to maximize the overall economic benefit of integrated exploitation for multiple deposits.
(2)
The extracted ore grade of deposits is different
The ore-matching work of deposits has a significant impact on the stability of ore quality and the improvement of resource utilization efficiency. It is necessary to optimize the extracted ore grade of deposits to effectively ensure that the ore production meets the requirements.
(3)
Sustainable exploitation is required for deposits
In order to ensure the continuity of the production for all deposits and the stable supply of mineral products, the deposit should not be interrupted for one single year, nor should a single deposit be over-mined or almost stopped. Therefore, the certain carbon quota and actual resource reserves of each deposit should be considered to limit production capacity in the optimization model.
(4)
Set constraints to meet production requirements of deposits
Constraints such as maximum and minimum mining capacity, extracted ore grade range, average concentrated ore grade requirement, overall production scale, and metal output target of each deposit are set for the optimization model to meet the production requirements of deposits. The technical and economic indexes associated with the integrated exploitation of multiple deposits are considered due to the actual production situation.
(5)
Consider interaction relations between deposits
The interaction between multiple deposits is closely linked to the overall economic benefit in the context of integrated exploitation. Competition relations between deposits to compete resources needs to be considered when constructing the optimization model.

3.2. Modeling Variables and Parameters

Table 2 represents sets in the model. Table 3 and Table 4 express the variables of the mathematical model. The parameters in the model are shown in Table 5.

3.3. Objective Function

The objective function is given by Equation (1), which seeks to maximize the overall economic benefit and the optimal allocation of production capacity of multiple deposits. Simultaneously, it may reduce the total cost and relieve the pressure of resource continuity. Equation (1) is obtained by subtracting the total cost of mining, ore processing and smelting, and the transportation cost of concentrated ore from the total sales revenues of the metal product.
M a x   S = P i = 1 I X i g i r i h i i = 1 I C i X i i = 1 I Z i X i g i r i h i a i

3.4. Constraints

Equation (2) represents the total annual metal output of the integrated exploitation for multiple deposits. The total annual metal output of multiple deposits may reach a certain value in order to meet production requirements while constraining the total metal output under a certain carbon quota of deposits.
Equations (3) and (4) ensure that it is not possible to exploit more ores than the production capacity of deposits. Equation (3) illustrates that setting a minimum production capacity of multiple deposits is to avoid enterprise layoffs. The production capacity of the respective deposits is limited and influenced by the scale, capital expenditure, raw material supply, and technical level of mines. In other words, the annual operating volume of the deposits cannot exceed the maximum productive capacity constraints.
Equations (5) and (6) control the extracted ore grade of deposits so that the fluctuation range of the deposits’ extracted ore grade is in accordance with the actual resource conditions. It ensures the stability of ore quality and improves resource utilization efficiency.
The average concentrated ore grade constraint of all deposits can be evaluated and restrained using Equations (7)–(9), which could ensure that the average concentrated ore grade from ore processing plants to the shared smelter each year is within a reasonable interval.
S .   T . i = 1 I M i = i = 1 I X i g i r i h i = Y
i = 1 I X i L X i
i = 1 I X i U X i
g i L g i
g i U g i
i = 1 I X i g i r i h i = b i = 1 I X i g i r i h i a i
b L b
b U b                   i = 1 , 2 , 3 , , I

4. Modeling Application and Result Analysis

4.1. Regression Analysis

The annual production capacity range of the three deposits is 2.3 × 106 t to 2.9 × 106 t, 2.1 × 106 t to 2.5 × 106 t, and 1.6 × 106 t to 2.2 × 106 t, respectively. In the process of the integrated exploitation of mineral resources, the change in the extracted ore grade will affect the change in other technical and economic indexes of each deposit, which may have a significant impact on the benefit and the resource utilization rate as well. Therefore, the establishment of a regression analysis of indexes based on historical data is considered. The regression equations of the three deposits’ indexes are derived by gathering actual production data of three deposits from 2020 to 2021, excluding anomalies.
(1)
Recovery rate of ore processing and extracted ore grade
The least square method is used for polynomial fitting. Equations (10)–(12) demonstrate the correlations between the recovery rate of ore processing and extracted ore grade from the regression results for deposits I, II, and III, which are shown in Figure 5. g i represents the extracted ore grade of each deposit, and the unit is g/t. r i represents the recovery rate of ore processing, and the unit is %.
r 1 = 1.529 g 1 2 + 8.386 g 1 + 84.093
r 2 = 0.160 g 2 2 0.121 g 2 + 93.922
r 3 = 0.369 g 3 2 1.056 g 3 + 95.780
(2)
Concentrated ore grade and extracted ore grade
The least square method is used for polynomial fitting. Equations (13)–(15) demonstrate the correlations between concentrated ore grade and extracted ore grade from the regression results for deposits I, II, and III, which are shown in Figure 6. g i represents the extracted ore grade of each deposit, and the unit is g/t. a i represents the concentrated ore grade of the deposit, and the unit is g/t.
a 1 = 22.465 g 1 2 73.646 g 1 + 99.848
a 2 = 20.270 g 2 2 63.127 g 2 + 102.484
a 3 = 24.233 g 3 2 + 111.998 g 3 61.893
(3)
Operation cost and extracted ore grade
The least square method is used for polynomial fitting. Equations (16)–(18) demonstrate the correlations between operation cost and extracted ore grade from the regression results for deposits I, II, and III, which are shown in Figure 7. The currency utilized is Chinese RMB (¥). g i represents the extracted ore grade of each deposit, and the unit is g/t. C i represents the operating cost of each deposit, and the unit is RMB/t.
C 1 = 28.647 g 1 2 209.491 g 1 + 638.504
C 2 = 106.573 g 2 2 559.193 g 2 + 1050.944
C 3 = 98.927 g 3 2 485.634 g 3 + 915.038
(4)
Recovery rate of smelting and concentrated ore grade
Exploring change rules of the smelting recovery rate and concentrated ore grade of deposits is conducive to defining the limited smelting recovery rate and feeding back the reasonable concentrated ore grade of the ore processing plant. It is of great significance for determining the reasonable concentrated ore grade of deposits, which may improve the smelting recovery rate and reduce the operation cost.
The least square method is used for polynomial fitting. Equation (19) demonstrates the correlations between the recovery rate of smelting and concentrated ore grade from the regression results for deposits, which are shown in Figure 8. b represents the concentrated ore grade of deposits, and the unit is g/t. h represents the recovery rate of smelting, and the unit is %.
h = 0.00214 b 2 + 0.310 b + 86.392

4.2. Parameters

To implement the optimization model, some parameters such as gold price, total gold throughput, transportation cost of per ton, limitation of extracted ore grade, and concentrated ore grade are collected as shown in Table 6.

4.3. Result Analysis

Lingo is used as a tool to calculate the mathematical optimization model with corresponding parameters. The global options function is applied to obtain the optimal solution of the model. The optimization model is a nonlinear programming mathematical model, and the global optimal solution is finally reached. The optimal results are obtained and shown in Table 7, Table 8 and Table 9.
Table 10 and Table 11 summarize the comparative analysis of production capacity and technical and economic indexes of each gold deposit before and after the optimization. The comparison of production capacity and the economic benefit of integrated exploitation under certain carbon quota is shown in Figure 9 and Figure 10. The total gold throughput in the optimization model is set to 13 t per year, at which time, the annual gold production capacity allocation of the three deposits shown is relatively balanced and stable. The comparison analysis displays that the optimized scheme achieves the goal of maximizing the overall economic benefit.
Generally, the optimal decision-making scheme for each deposit under carbon quota is calculated based on the optimization model, which could improve the comprehensive utilization efficiency of mineral resources while reducing the costs of mining, ore processing, and smelting under the constraint of the total gold throughput.

4.4. Sensitivity Analysis

The optimization model obtains the decision-making scheme of production capacity allocation for multiple deposits from the perspective of best condition in technology and economy. Changes in each index have different effects on the results obtained from the optimization model. Therefore, sensitivity analysis can be carried out to guide dynamic decision making by changing each parameter.
Figure 11 shows the results of the gold price sensitivity analysis. Price changes from −30% to 30% are performed to evaluate the overall economic benefit of the three deposits. It can be seen that the gold price and the overall economic benefit are positively correlated, while the change in the gold price has low sensitivity to the model decision variables. The break-even point occurs when the gold price falls to 157 RMB/g to 158 RMB/g range. Due to the gold market price changing from 150 RMB/g to 430 RMB/g in the last ten years, a positive benefit is always obtained.
Figure 12 illustrates the results of gold throughput sensitivity analysis. It can be observed that the variation in gold throughput has the highest sensitivity to the overall economic benefit and average concentrated ore grade, followed by the sensitivity to the extracted ore grade, and the sensitivity to the recovery rate of ore processing and concentrated ore grade is lower. When the total gold throughput of three deposits is constrained to around 14t, the production capacity allocation of three deposits is relatively balanced, and the average concentrated ore grade could reach the maximum simultaneously.
The results of the average concentrated ore grade sensitivity analysis are presented in Figure 13. It can be seen that the variation in the average concentrated ore grade has the strongest sensitivity to the overall economic benefit. The sensitivity of the extracted ore grade for deposits I and II is higher than that of deposit III. As the average grade of the concentrated ore increases, the production capacity allocation of deposits I and III moves closer to deposit II, which implies that the allocation results tend to be balanced and stable.
Decision-making schemes of integrated exploitation for multiple deposits under different production conditions can be formed based on the above sensitivity analysis. In addition, suggestions for dynamic decision making regarding mines can be obtained. Scientific decision making regarding multiple deposits’ integrated exploitation should be made according to key factors such as actual production conditions, overall production strategic objectives, and technical and economic level of deposits. Changes to resource conditions and the market should be considered to determine the optimal allocation of production capacity under the carbon quota. Increasing the overall gold output of multiple deposits appropriately is expected to improve the comprehensive utilization efficiency of resources and maximize the economic benefit under the condition that the gold price keeps rising and the average concentrated ore grade of multiple deposits is guaranteed to meet the requirements.

5. Conclusions

The strategy of mineral resources’ integrated exploitation is considered to reduce carbon tax generated by carbon emissions. The allocation of production capacity for multiple deposits could be optimized to obtain more benefits for a mining company. The optimization decision-making framework and modeling assumptions are proposed based on the analysis of the characteristics and advantages of integrated exploitation.
An optimization model for the integrated exploitation of multiple deposits is constructed and applied to three deposits under a gold mining company. Due to integrated exploitation under the carbon emission constraints, the decision-making scheme generated using the mathematical programming method outperforms that obtained considering a single deposit independently. The optimal solution of production capacity allocation and extracted ore grade for each deposit is obtained according to the optimal scheme, which provides a new method for research on the optimization of mineral resources’ integrated exploitation.
The implementation of the model in multiple gold deposits could improve the overall economic benefit and reduce the costs of mining, ore processing, and smelting. Profits of whole gold deposits increase from 2.454 × 109 RMB to 3.414 × 109 RMB. Although the overall gold throughput of multiple deposits decreases under a carbon quota, the optimal scheme implemented by the optimization model improves the technical and economic indexes of integrated exploitation, which highlights the ability to improve the comprehensive efficiency of resource utilization and profits.
The results demonstrate that the mathematical model in this paper is reasonable and feasible. The results obtained enable enterprises to allocate production capacity more purposefully and make optimal production planning. Low-carbon mining is taken into account innovatively to make mining companies more sustainable. A sensitivity analysis forms different optimal plans for the integrated exploitation of mineral resources from different aspects, which might guide the dynamic decision making of deposits.
The optimization method and model have practical applicability to solve the production capacity allocation problem of multiple deposits during the integrated exploitation of resources. The optimization model could be applied in resource exploitation in multiple similar deposits. However, parameters should be updated based on the specific situation, and regression equations have to be fitted according to the real data from the mining company. The proposed optimization model does not consider the cyclic utilization of resources. At the same time, the complex correlations between indexes of multiple deposits are simplified. Further research could be conducted in this direction.

Author Contributions

Software, X.Q.; formal analysis, B.B.; resources, C.F. and Y.Z.; writing—original draft, Y.G.; Writing—review and editing, J.H.; supervision, G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Key R&D Program of China (no. 2022YFC2903905), the National Natural Science Foundation of China (no. 52074022), and Anhui Mine IoT and Security Monitoring Technology Key Laboratory (Grant No. 2109Y-09-03).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The lingo 18.0 optimization software was used in this work to solve the nonlinear programming model.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Franzen, A.; Mader, S. Consumption-based versus production-based accounting of CO2 emissions: Is there evidence for carbon leakage? Environ. Sci. Policy 2018, 84, 34–40. [Google Scholar] [CrossRef]
  2. Azadi, M.; Northey, S.A.; Ali, S.H.; Edraki, M. Transparency on greenhouse gas emissions from mining to enable climate change mitigation. Nat. Geosci. 2020, 13, 100–104. [Google Scholar] [CrossRef]
  3. Cox, B.; Innis, S.; Kunz, N.C.; Steen, J. The mining industry as a net beneficiary of a global tax on carbon emissions. Commun. Earth Environ. 2022, 3, 1–8. [Google Scholar] [CrossRef]
  4. Dincer, H.; Yüksel, S.; Karakuş, H.; Kalkavan, H. The importance of carbon emission of manufacturing companies on sustainable economic development. In Globalization, Income Distribution and Sustainable Development; Emerald Publishing Limited: Bingley, UK, 2022; pp. 81–93. [Google Scholar]
  5. Sam, U.; Allan, T.; Steffen, H. Gold mining greenhouse gas emissions, abatement measures, and the impact of a carbon price. J. Clean. Prod. 2022, 340, 130851. [Google Scholar]
  6. Lagos, G.; Peters, D.; Videla, A.; Jara, J.J. The effect of mine aging on the evolution of environmental footprint indicators in the Chilean copper mining industry 2001–2015. J. Clean. Prod. 2018, 174, 389–400. [Google Scholar] [CrossRef]
  7. Zafeiriou, E.; Spinthiropoulos, K.; Tsanaktsidis, C.; Garefalakis, S.; Panitsidis, K.; Garefalakis, A.; Arabatzis, G. Energy and mineral resources exploitation in the deliginitization era: The case of greek peripheries. Energies 2022, 15, 4732. [Google Scholar] [CrossRef]
  8. Daju, R.D. Comprehensive mineral-exploration/exploitation of ore deposits, with some new techniques, value-addition and creation of wealth from waste. J. Appl. Geochem. 2019, 21, 411–427. [Google Scholar]
  9. Tkačíková, L. More economical exploitation of mineral resources. Int. Multidiscip. Sci. GeoConf. SGEM 2018, 18, 593–599. [Google Scholar]
  10. Moore, K.R.; Whyte, N.; Roberts, D.; Allwood, J.; Leal-Ayala, D.R.; Bertrand, G.; Bloodworth, A.J. The re-direction of small deposit mining: Technological solutions for raw materials supply security in a whole systems context. Resour. Conserv. Recycl. X 2020, 7, 100040. [Google Scholar] [CrossRef]
  11. Marklund, P.O.; Samakovlis, E. What is driving the EU burden-sharing agreement: Efficiency or equity? J. Environ. Manag. 2007, 85, 317–329. [Google Scholar] [CrossRef]
  12. Nasih, M.; Harymawan, I.; Paramitasari, Y.I.; Handayani, A. Carbon emissions, firm size, and corporate governance structure: Evidence from the mining and agricultural industries in Indonesia. Sustainability 2019, 11, 2483. [Google Scholar] [CrossRef]
  13. Tyuleneva, T. Problems and prospects of regional mining industry digitalization. EDP Sci. 2020, 174, 04019. [Google Scholar] [CrossRef]
  14. Zhernov, E.; Nekhoda, E.; Rusak, I. Digitalization in a resource-extracting region: From cluster to business ecosystem. EDP Sci. 2021, 315, 04009. [Google Scholar] [CrossRef]
  15. Sandeman, T.; Fricke, C.; Bodon, P.; Stanford, C. Integrating optimization and simulation-A comparison of two case studies in mine planning. In Proceedings of the 2010 Winter Simulation Conference IEEE, Baltimore, MD, USA, 5–8 December 2010; pp. 1898–1910. [Google Scholar]
  16. Musingwini, C. Optimization in underground mine planning-developments and opportunities. J. S. Afr. Inst. Min. Metall. 2016, 116, 809–820. [Google Scholar] [CrossRef]
  17. Ivanova, S.; Vesnina, A.; Fotina, N.; Prosekov, A. An overview of carbon footprint of coal mining to curtail greenhouse gas emissions. Sustainability 2022, 14, 15135. [Google Scholar] [CrossRef]
  18. Valderrama, C.V.; Santibanez-González, E.; Pimentel, B.; Candia-Véjar, A.; Canales-Bustos, L. Designing an environmental supply chain network in the mining industry to reduce carbon emissions. J. Clean. Prod. 2020, 254, 119688. [Google Scholar] [CrossRef]
  19. Azhar, N.A.B.; Gunawan, A.; Cheng, S.-F.; Leonardi, E. A carbon-aware planning framework for production scheduling in mining. In Lecture Notes in Computer Science, Proceedings of the International Conference on Computational Logistics, Barcelona, Spain, 21–23 September 2022; Springer: Cham, Switzerland, 2022; pp. 441–456. [Google Scholar]
  20. Boubault, A.; Maïzi, N. Devising mineral resource supply pathways to a low-carbon electricity generation by 2100. Resources 2019, 8, 33. [Google Scholar] [CrossRef]
  21. Markevych, K.; Maistro, S.; Koval, V.; Paliukh, V. Mining sustainability and circular economy in the context of economic security in Ukraine. Min. Miner. Depos. 2022, 16, 101–113. [Google Scholar] [CrossRef]
  22. Khazaei, S.; Pourrahimian, Y. Mathematical programming application in sublevel caving production scheduling. Mining 2021, 1, 180–191. [Google Scholar] [CrossRef]
  23. Newman, A.M.; Rubio, E.; Caro, R.; Weintraub, A.; Eurek, K. A review of operations research in mine planning. Interfaces 2010, 40, 222–245. [Google Scholar] [CrossRef]
  24. Vujic, S.; Benovic, T.; Miljanovic, I.; Hudej, M.; Milutinovic, A.; Pavlovic, P. Fuzzy linear model for production optimization of mining systems with multiple entities. Int. J. Miner. Metall. Mater. 2011, 18, 633–637. [Google Scholar] [CrossRef]
  25. Rahmanpour, M.; Osanloo, M. Application of fuzzy linear programming for short-term planning and quality control in mine complexes. J. S. Afr. Inst. Min. Metall. 2017, 117, 684–694. [Google Scholar] [CrossRef]
  26. Gligorić, M.V.; Gligorić, Z.M.; Beljić, Č.R.; Lutovac, S.M.; Damnjanović, V.M. Long-term room and pillar mine production planning based on fuzzy 0-1 linear programming and multicriteria clustering algorithm with uncertainty. Math. Probl. Eng. 2019, 2019, 3078234. [Google Scholar] [CrossRef]
  27. Ahmadi, M.R.; Shahabi, R.S. Cutoff grade optimization in open pit mines using genetic algorithm. Resour. Policy 2018, 55, 184–191. [Google Scholar] [CrossRef]
  28. Montiel, L.; Dimitrakopoulos, R. Optimizing mining complexes with multiple processing and transportation alternatives: An uncertainty-based approach. Eur. J. Oper. Res. 2015, 247, 166–178. [Google Scholar] [CrossRef]
  29. Canales-Bustos, L.; Santibañez-González, E.; Candia-Véjar, A. A multi-objective optimization model for the design of an effective decarbonized supply chain in mining. Int. J. Prod. Econ. 2017, 193, 449–464. [Google Scholar] [CrossRef]
  30. Gupta, P.; Mehlawat, M.K.; Aggarwal, U.; Charles, V. An integrated AHP-DEA multi-objective optimization model for sustainable transportation in mining industry. Resour. Policy 2018, 74, 101180. [Google Scholar] [CrossRef]
  31. Benndorf, J.; Yueksel, C.; Shishvan, M.; Rosenberg, H.; Thielemann, T.; Mittmann, R.; Lohsträter, O.; Lindig, M.; Minnecker, C.; Donner, R.; et al. RTRO–coal: Real-time resource-reconciliation and optimization for exploitation of coal deposits. Minerals 2015, 5, 546–569. [Google Scholar] [CrossRef]
  32. Helen, L.; Wilson, H.N.; Macdonald, E.K.; Clarke, B. Resource integration, value creation and value destruction in collective consumption contexts. J. Bus. Res. 2019, 103, 173–185. [Google Scholar]
  33. Lomakina, N. Complex exploitation of the Far East mineral resources: Effects new model of the macroregional development. EDP Sci. 2020, 192, 03006. [Google Scholar] [CrossRef]
  34. Smol, M.; Marcinek, P.; Duda, J.; Szołdrowska, D. Importance of sustainable mineral resource management in implementing the circular economy (CE) model and the European green deal strategy. Resources 2020, 9, 55. [Google Scholar] [CrossRef]
  35. Sauer, P.C.; Seuring, S. Sustainable supply chain management for minerals. J. Clean. Prod. 2017, 151, 235–249. [Google Scholar] [CrossRef]
  36. Tabelin, C.B.; Dallas, J.; Casanova, S.; Pelech, T.; Bournival, G.; Saydam, S.; Canbulat, I. Towards a low-carbon society: A review of lithium resource availability, challenges and innovations in mining, extraction and recycling, and future perspectives. Miner. Eng. 2021, 163, 106743. [Google Scholar] [CrossRef]
  37. Luger, J.; Raisch, S.; Schimmer, M. Dynamic balancing of exploration and exploitation: The contingent benefits of ambidexterity. Organ. Sci. 2018, 29, 449–470. [Google Scholar] [CrossRef] [Green Version]
  38. Poufinas, T.; Galanos, G.; Papadimitriou, P. The Competitiveness of Small and Medium Enterprises in Adverse Economic Environments. Theor. Econ. Lett. 2018, 8, 2788–2802. [Google Scholar] [CrossRef]
  39. Ruokonen, E.; Temmes, A. The approaches of strategic environmental management used by mining companies in Finland. J. Clean. Prod. 2019, 210, 466–476. [Google Scholar] [CrossRef]
  40. Pochtovyuk, A.; Latyshev, K.; Zaika, K. Resource potential as a key factor of Ukrainian mining enterprises competitive position. EDP Sci. 2021, 91, 01017. [Google Scholar] [CrossRef]
Figure 1. Research content.
Figure 1. Research content.
Sustainability 15 02917 g001
Figure 2. Production process of these three deposits.
Figure 2. Production process of these three deposits.
Sustainability 15 02917 g002
Figure 3. Optimization strategy of integrated exploitation for multiple deposits.
Figure 3. Optimization strategy of integrated exploitation for multiple deposits.
Sustainability 15 02917 g003
Figure 4. Decision-making framework of integrated exploitation for multiple deposits.
Figure 4. Decision-making framework of integrated exploitation for multiple deposits.
Sustainability 15 02917 g004
Figure 5. Regression results of the recovery rate of ore processing and extracted ore grade. (a) Deposit I; (b) deposit II; (c) deposit III.
Figure 5. Regression results of the recovery rate of ore processing and extracted ore grade. (a) Deposit I; (b) deposit II; (c) deposit III.
Sustainability 15 02917 g005
Figure 6. Regression results of concentrated ore grade and extracted ore grade. (a) Deposit I; (b) deposit II; (c) deposit III.
Figure 6. Regression results of concentrated ore grade and extracted ore grade. (a) Deposit I; (b) deposit II; (c) deposit III.
Sustainability 15 02917 g006
Figure 7. Regression results of operation cost and extracted ore grade. (a) Deposit I; (b) deposit II; (c) deposit III.
Figure 7. Regression results of operation cost and extracted ore grade. (a) Deposit I; (b) deposit II; (c) deposit III.
Sustainability 15 02917 g007
Figure 8. Regression results of recovery rate of smelting and concentrated ore grade.
Figure 8. Regression results of recovery rate of smelting and concentrated ore grade.
Sustainability 15 02917 g008
Figure 9. Comparison analysis of production capacity allocation before and after optimization.
Figure 9. Comparison analysis of production capacity allocation before and after optimization.
Sustainability 15 02917 g009
Figure 10. Comparison analysis of economic benefit before and after optimization.
Figure 10. Comparison analysis of economic benefit before and after optimization.
Sustainability 15 02917 g010
Figure 11. Sensitivity analysis results with different gold prices. (a) Overall economic benefit; (b) extracted ore grade; (c) production capacity allocation.
Figure 11. Sensitivity analysis results with different gold prices. (a) Overall economic benefit; (b) extracted ore grade; (c) production capacity allocation.
Sustainability 15 02917 g011
Figure 12. Sensitivity analysis results of gold throughput. (a) Overall economic benefit and average concentrated ore grade; (b) extracted ore grade of each deposit with different gold throughput; (c) production capacity allocation of each deposit with different gold throughput; (d) recovery rate of ore processing of each deposit with different gold throughput; (e) concentrated ore grade of each deposit with different gold throughput.
Figure 12. Sensitivity analysis results of gold throughput. (a) Overall economic benefit and average concentrated ore grade; (b) extracted ore grade of each deposit with different gold throughput; (c) production capacity allocation of each deposit with different gold throughput; (d) recovery rate of ore processing of each deposit with different gold throughput; (e) concentrated ore grade of each deposit with different gold throughput.
Sustainability 15 02917 g012
Figure 13. Sensitivity analysis results of average concentrated ore grade. (a) Overall economic benefit with different average concentrated ore grade; (b) extracted ore grade of each deposit with different average concentrated ore grade; (c) production capacity allocation of each deposit with different average concentrated ore grade.
Figure 13. Sensitivity analysis results of average concentrated ore grade. (a) Overall economic benefit with different average concentrated ore grade; (b) extracted ore grade of each deposit with different average concentrated ore grade; (c) production capacity allocation of each deposit with different average concentrated ore grade.
Sustainability 15 02917 g013
Table 1. Statistics of reserve information from three deposits (up to 2021).
Table 1. Statistics of reserve information from three deposits (up to 2021).
DepositOre Tonnage (104 t)Grade (g/t)Metal Tonnage (t)
I43022.84122.237
II9683.5434.286
III4232.8111.873
Table 2. Set list.
Table 2. Set list.
SetDescription
i Set of multiple deposits
Table 3. Decision variables.
Table 3. Decision variables.
VariablesDescription
X i Annual ore processing capacity of deposit i
g i Extracted ore grade of deposit i
Table 4. Other variables.
Table 4. Other variables.
VariablesDescription
r i Recovery rate of ore processing for deposit i
h i Recovery rate of ore smelting for deposit i
a i Concentrated ore grade of deposit i
C i Total costs of mining, ore processing, and smelting per ton of deposit i
M i Metal output of deposit i
S Overall economic benefit
Table 5. Parameters.
Table 5. Parameters.
ParametersDescription
I Total number of deposits
P Price of mineral products
Y Total annual metal output of multiple deposits
b Average concentrated ore grade from ore processing plants to the shared smelting plant
Z i Transportation cost per ton of concentrated ore for deposit i
L X i Minimum production capacity of deposit i
U X i Maximum production capacity of deposit i
L g i Lower limit of extracted ore grade of deposit i
U g i Upper limit of extracted ore grade of deposit i
L   b Lower limit of average concentrated ore grade
U   b Upper limit of average concentrated ore grade
Table 6. Parameter information collection.
Table 6. Parameter information collection.
DescriptionUnitValueExplanation
Gold priceRMB/g420Statistics from iiMedia Think Tank show that the higher yield of gold price mainly fluctuates around 420 RMB/g in recent years.
Total gold throughputt13The total annual gold throughput of three deposits is 13t under certain carbon emission quota of deposits.
Transportation cost of ore per tonRMB/t14.6According to the distance between the respective ore processing plant and the shared smelter, the transportation cost of ore per ton can be estimated.
4.2
6.3
Range of extracted ore gradeg/t1.59~2.11According to the previously extracted ore grade data and the regression analysis of the above indexes, the extracted ore grade intervals of three deposits are estimated.
1.50~2.28
1.58~2.36
Range of concentrated ore gradeg/t51~58Ensure the overall concentrated ore quality of three deposits.
Table 7. Statistics of extracted ore grade and production capacity.
Table 7. Statistics of extracted ore grade and production capacity.
DepositExtracted Ore Grade (g/t)Ore Processing Capacity (104 t)Gold Throughput (kg)
I2.112304474
II2.282104391
III2.361894135
Table 8. Statistics of technical indexes.
Table 8. Statistics of technical indexes.
DepositRecovery Rate of Ore Processing (%)Concentrated Ore Grade (g/t)Average Concentrated Ore Grade (g/t)Recovery Rate of Smelting (%)
I94.9844.4756.3897.07
II94.4863.93
III95.3467.45
Table 9. Statistics of economic indexes.
Table 9. Statistics of economic indexes.
DepositCost Per Ton (RMB/t)Economic Benefit (108 RMB)General Economic Benefit (108 RMB)
I32411.324634.1395
II33011.5092
III32011.3058
Table 10. Comparative analysis of overall gold deposits before and after optimization.
Table 10. Comparative analysis of overall gold deposits before and after optimization.
SortTotal Ore Tonnage (104 t)Total Gold Throughput (t)Total Cost (108 RMB)General Economic Benefit (108 RMB)
Before11041941.078424.5465
After6291320.430034.1395
Table 11. Comparative analysis of technical and economic indexes of each gold deposit before and after optimization.
Table 11. Comparative analysis of technical and economic indexes of each gold deposit before and after optimization.
SortExtracted Ore Grade (g/t)Recovery Rate of Ore Processing (%)Concentrated Ore Grade (g/t)Recovery Rate of Smelting (%)
IIIIIIIIIIIIIIIIIIIIIIII
Before1.941.732.0994.6394.2295.1338.7853.2165.6095.6997.2297.67
After2.112.282.3694.9894.4895.3444.4769.9367.4597.0797.0797.07
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

Gu, Y.; Li, G.; Hou, J.; Fan, C.; Qiang, X.; Bai, B.; Zhang, Y. Production Strategy Optimization of Integrated Exploitation for Multiple Deposits Considering Carbon Quota. Sustainability 2023, 15, 2917. https://doi.org/10.3390/su15042917

AMA Style

Gu Y, Li G, Hou J, Fan C, Qiang X, Bai B, Zhang Y. Production Strategy Optimization of Integrated Exploitation for Multiple Deposits Considering Carbon Quota. Sustainability. 2023; 15(4):2917. https://doi.org/10.3390/su15042917

Chicago/Turabian Style

Gu, Yingyu, Guoqing Li, Jie Hou, Chunchao Fan, Xingbang Qiang, Bin Bai, and Yongfang Zhang. 2023. "Production Strategy Optimization of Integrated Exploitation for Multiple Deposits Considering Carbon Quota" Sustainability 15, no. 4: 2917. https://doi.org/10.3390/su15042917

APA Style

Gu, Y., Li, G., Hou, J., Fan, C., Qiang, X., Bai, B., & Zhang, Y. (2023). Production Strategy Optimization of Integrated Exploitation for Multiple Deposits Considering Carbon Quota. Sustainability, 15(4), 2917. https://doi.org/10.3390/su15042917

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