1. INTRODUCTION
⌅In 2022, global coal consumption hit an unprecedented level of 8,415 Mt, marking a 4% rise. This increase was predominantly driven by countries like China, where coal continues to account for over 60% of electricity production. Despite efforts to diversify energy sources, coal demand in 2023 was expected to see a modest uptick of 1.4%, setting another record at approximately 8,536 Mt (1Coal2023. Analysis and forecast to 2026. Retrieved fromhttps://iea.blob.core.windows.net/assets/a72a7ffa-c5f2-4ed8-a2bf-eb035931d95c/Coal_2023.pdf, 2bp Energy Outlook. 2023 edition. Retrieved fromhttps://www.bp.com/content/dam/bp/business-sites/en/global/corporate/pdfs/energy-economics/energy-outlook/bp-energy-outlook-2023.pdf). However, the downside of coal production is the substantial waste it generates, ranging from 10-15% by weight (3LiJ, WangJ. 2019. Comprehensive utilization and environmental risks of coal gangue: A review. J. Clean. Prod.239:117946. 10.1016/j.jclepro.2019.117946). This could result in more than 150 gigatons of additional waste globally, causing adverse ecological impacts on water, air, and soil. In Morocco, the closure of the Jerada coal mine in 2001 has left around 40 million tons of coal waste, leading to issues like acid mine drainage affecting groundwater. The interaction with rainwater triggers acid mine drainage (AMD), a phenomenon resulting from the oxidation of pyrite (Fe2S) present in the coal waste. This AMD negatively affects groundwater, emphasizing the urgent need for effective environmental management and remediation strategies to mitigate the ecological consequences of coal mining activities. In general, the recycling of mining waste aligns with the principles of a circular economy, emphasizing sustainable resource utilization (4Tayebi-KhoramiM, EdrakiM, CorderG, GolevA. 2019. Re-thinking mining waste through an integrative approach led by circular economy aspirations. J. Minerals, 9(5):286. 10.3390/min9050286.). Specifically, the recycling of coal gangue aims to exploit its mineral content, but also the coal remains it contains. China has successfully utilized coal gangue as a low-quality energy source for power generation (5YuJL, MengFR, LiXC, TahmasebiA. 2012. Power generation from coal gangue in China: Current status and development. J. Open J. Adv. Mater.550:443-446. 10.4028/www.scientific.net/AMR.550-553.443). In general, the major components of coal gangue are silica, alumina, iron oxide, other oxides and carbon remains, hence a possibility of recycling in civil engineering, or energy production. Numerous scientific studies have explored its applications in cement (6Shuang-xiZ. 2009. Study on the reaction degree of calcined coal gangue powder in blended cement by selective solution method. J. Procedia Earth and Planetary Science, 1(1):634-639. 10.1016/j.proeps.2009.09.100), building materials (7AddouR, HannawiK, AgbodjanWP, ZenasniM. 2015. Caractérisation des déchets stériles de charbon de la mine de Jerada (Est du Maroc), en vue de l’élaboration d’un éco-matériau en génie civil. Mécanique et Electrique pour l’Energie (CMEEE 2015).), bricks (8XuH, SongW, CaoW, ShaoG, LuH, YangD, ChenD, ZhangR. 2017. Utilization of coal gangue for the production of brick. J. Material cycles and waste management. 19:1270-1278. 10.1007/s10163-016-0521-0) and ceramics (9YangM, GuoZ, DengY, XingX, QiuK, LongJ, LiJ. 2012. Preparation of CaO–Al2O3–SiO2 glass ceramics from coal gangue. Int. J. Miner. Process, 102:112-115. 10.1016/j.minpro.2011.11.004).
Moreover, cement, a vital strategic material, holds a significant global production scale. In 2020, global cement consumption reached 4.17 billion tons, making concrete the second most consumed product worldwide after water, as reported by Cembureau and the United Nations (10Cembureau and the United Nations, Cementing the European Green Deal. Reaching climate neutrality along the cement and concrete value chain by 2050. Retrieved fromhttps://cembureau.eu/media/kuxd32gi/cembureau-2050-roadmap_final-version_web.pdf). Understanding various stages in cement production is crucial. A study by Zhenning et al (11ZhaoH, ZhangN, WangHJ. 2014. Power consumption prediction modeling of cement manufacturing based on the improved multiple non-linear regression algorithm. Appl. Mech. Mater.687:5185-5189. https://www.scientific.net/AMM.687-691.5185) outlines the electricity consumption distribution in a typical cement plant: extraction and crushing (5%), grinding raw materials into powder (24%), homogenization (6%), clinker burning and cooling (22%), clinker grinding (38%), and packaging and loading (5%). The cement industry is highly energy-intensive, requiring 3.2-6.3 GJ per ton of clinker (12RahmanA, RasulMG, KhanMMK, SharmaS. 2015. Recent development on the uses of alternative fuels in cement manufacturing process. J. Fuel.145:84-99. 10.1016/j.fuel.2014.12.029). By 2050, energy costs for cement manufacturing are projected to represent 20-40% of total production expenses. Currently, electricity accounts for 12% of this energy consumption, but it’s expected to increase to over 20% by that year (13KermeliK, EdelenboschOY, Crijns-GrausW, van RuijvenBJ, MimaS, van VuurenDP, WorrellE. 2019. The scope for better industry representation in long-term energy models: Modeling the cement industry. J. Appl. Energy.240:964-985. 10.1016/j.apenergy.2019.01.252-14International Energy Agency. 2023. World energy outlook. Retrieved fromhttps://iea.blob.core.windows.net/assets/86ede39e-4436-42d7-ba2a-edf61467e070/WorldEnergyOutlook2023.pdf). The enormous quantities of cement produced offer the possibility of recycling industrial and construction waste. This includes the potential for recycling and recovering coal gangue, especially considering its composition containing key cement oxides.
In a previous study by Belkheiri et al. (15BelkheiriaD, DiouriA, TaibiM, SassiO, ArideJ. 2015. Recycling of Moroccan coal gangue in the elaboration of a Portland clinker. J. Mater. Environ. Sci.6(6):1570-1577.), Moroccan coal gangue was used to produce a Portland clinker. Their study involved clinkerization of a mixture of limestone-gangue at 18.5% (wt), resulting in a satisfactory Portland clinker denoted LG. The correct composition of the gangue in the limestone-gangue mixture (LG) was determined through computer simulation, aligning 11 parameters of LG with a reference raw cement plant (RR). These parameters included four oxides (CaO, SiO2, Al2O3, Fe2O3), three moduli (LSF, SR, AF), and the four phases of the clinker. This current work is a continuation of the previous study, with a focus on the energy gain achieved through the combustion of the remaining coal within the cement raw limestone-gangue mixture. DSC characterization of the LG mixture revealed a combustion phenomenon occurring around 300°C-600°C. This exothermic effect has the potential to contribute to reducing energy consumption during the heat treatments of the clinker process. The ultimate goal is to achieve total recycling and recover coal waste while producing a Portland limestone-gangue clinker. The mineral components become part of the clinker, and the combustion of the organic part during the cement process could lead to a reduction in energy costs, contributing to environmental, social and economic sustainability.
2. MATERIALS AND METHODS
⌅2.1. Used materials
⌅The coal gangue sample was collected from various points in the backfill of Jerada's old mine in Morocco, while the limestone and raw meal originated from a cement plant. The materials and their abbreviations used are presented in Table 1.
For thermal studies, a pure coal sample is considered as a reference, denoted as PC. All materials were finely ground to a diameter of approximately 160 µm for consistency in the study.
2.2. Composition of materials
⌅The elemental and oxide compositions of coal gangue (CG) and pure coal (PC) used in this study are derived from characterizations reported in the previous studies (15-18). The detailed compositions are presented in Table 2 and Table 3, respectively. Additionally, Table 4 provides the percentages of the studied cementitious materials, including the limestone reference (LR), cementitious mix Limestone-Gangue (LGx), and the cement raw meal reference (RR).
The mineralogical composition of the materials was determined through Rietveld quantitative X-ray diffraction (XRD) analysis. The crystalline phases identified in each material are as follows:
CG is formed by quartz 60.29%, illite ( K+,H3O+) Al2Si3AlO10(OH)2 :19.2%, muscovite Si3Al2.88Fe0.12KO12H2 :10.12%, clinochlore Si3Al2Mg5O18H8 : 6.8 % and gypsum 3.5%. LR is formed by calcite 80.26 %, dolomite : 12.2% and quartz 6.76%. RR is formed by calcite 63.40 %, dolomite 9.63 %, illite K0.7 Al2.1 (Si, Al)4 O10 (OH)2 : 13.13%, quartz 11.42% and gypsum 0.68%.
2.3. Methods
⌅We will adopt the following methodology
- Thermal characterizations of Coal Gangue (CG) and Limestone-Gangue (LG) will qualitatively illustrate the energy contribution of CG during the heating of the cement mixture LG. The Differential Scanning Calorimetry analysis (DSC) and thermos-gravimetric analysis (TG) were performed using DSC/TG type SDTQ 600 with a heating rate of 10°C/min up to 1000°C, and using 100ml/min of dry air free of CO2 as carrier gas.
- During the progressive heating process leading to the production of clinker, coal dispersed within the limestone-gangue mixture undergoes combustion. The combustion behavior of coal in this specific context is distinct from that in ambient air. The combustion observed during the progressive heating in the clinker production process will be assimilated to the combustion recorded in the Differential Scanning Calorimetry (DSC) analysis. The heat value derived from the combustion observed in the DSC analysis will be utilized to quantify and understand the energy contribution of coal during the clinker production.
- Based on the composition analyses, we will calculate the masses mi of the different substances i contained in the 1kg of each clinker from the two cementitious mixtures LG and RR, these clinkers are noted LG* and RR* respectively.
- We calculate the difference in thermochemical energy required to produce 1 kg of each of the LG* and RR* clinkers; minimizing the errors due to the approximations used.
- Energies for heating compounds: : formula (F1) where mi is the mass of compound i and Cp,i is the specific heat capacity at constant pressure.
- Energies related to chemical reactions denoted: where vk < 0 if Ak is a reactant and vk > 0 if Ak is a product; at the temperature T0, and pressure P°=1bar, the standard enthalpy of reaction is: : formula (F2), where is the standard enthalpy of formation for the compound Ak, at T0. At another temperature T, we calculate.
The standard enthalpy of reaction by using the following formula:
The thermodynamic data, including specific heat capacities (Cp) and enthalpies of formation, will be presented in Table 6.
- Computational calculations are conducted using Python software.
| Comp | a-2 | a-1 | a-0.5 | a0 | a1 | a2 | a3 | a4 | ∆m | M | Hmol | References |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| calcite(T<500) | -1,11E+09 | 2,08E+07 | 0,00E+00 | -1,59E+05 | 6,33E+02 | -1,38E+00 | 1,56E-03 | -7,25E-07 | 2,74E-02 | 1,00E-01 | -1,21E+03 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| calcite | -2,15E+07 | 0,00E+00 | 0,00E+00 | 9,95E+02 | 2,71E-01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 2,74E-02 | 1,00E-01 | -1,21E+03 | (26BlancP, VieillardP, GailhanouH, GaboreauS, GaucherE, FialipsCI, MadeB, GiffautE. 2015. A generalized model for predicting the thermodynamic properties of clay minerals. Am. J. Sci.315(8):734-780. 10.2475/08.2015.02) |
| illite | -1,81E+07 | 0,00E+00 | 0,00E+00 | 8,63E+02 | 5,97E-01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | -1,64E-01 | 3,89E-01 | -5,97E+03 | (26BlancP, VieillardP, GailhanouH, GaboreauS, GaucherE, FialipsCI, MadeB, GiffautE. 2015. A generalized model for predicting the thermodynamic properties of clay minerals. Am. J. Sci.315(8):734-780. 10.2475/08.2015.02) |
| muscovite | -1,76E+07 | 0,00E+00 | 0,00E+00 | 8,42E+02 | 5,82E-01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 2,22E-02 | 4,02E-01 | -5,97E+03 | (26BlancP, VieillardP, GailhanouH, GaboreauS, GaucherE, FialipsCI, MadeB, GiffautE. 2015. A generalized model for predicting the thermodynamic properties of clay minerals. Am. J. Sci.315(8):734-780. 10.2475/08.2015.02) |
| clinochlore | -3,20E+07 | 0,00E+00 | 0,00E+00 | 1,14E+03 | 5,09E-01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 1,49E-02 | 5,56E-01 | -8,91E+03 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| quartz-a | -8,00E+07 | 1,00E+06 | 0,00E+00 | -4,36E+03 | 1,03E+01 | -1,20E-03 | -1,50E-05 | 1,15E-08 | 3,94E-02 | 6,01E-02 | -9,11E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| quartzb | 3,21E+06 | 0,00E+00 | 0,00E+00 | 9,74E+02 | 1,72E-01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 3,94E-02 | 6,01E-02 | -9,11E+02 | (27GottschalkM. 1997. Internally Consistent Thermodynamic Data for rock-forming minerals in the system Si02-Ti02-Al203-Fe203-CaO-MgO-FeO-K20-Na20-H20-C02. Eur. J. Mineral.9:175-223cp dolomite.) |
| dolomite | 1,50E+05 | 0,00E+00 | -3,55E+04 | 2,97E+03 | -9,09E-01 | 4,18E-04 | 0,00E+00 | 0,00E+00 | 4,04E-03 | 1,84E-01 | -2,33E+03 | (28LeśniakB, ŁukaszS, JakubinaG. 2013. Institute for chemical processing of coal, zabrze, poland. The determination of the specific heat capacity of coal based on literature data. Chemik. 67(6):560-571.,29EisermannW, JohnsonP, Conge, WL. 1980. Estimating thermodynamic properties of coal, char, tar and ash. Fuel Process. Technol.3(1):39-53. 10.1016/0378-3820(80)90022-3) |
| coal | 0,00E+00 | 0,00E+00 | 0,00E+00 | -8,91E+01 | 3,67E+00 | -1,63E-03 | 0,00E+00 | 0,00E+00 | 3,16E-02 | 1,80E-02 | 0,00E+00 | (30LothenbachB, MatscheiT, MöschnerG, GlasserFP. 2008. Thermodynamic modelling of the effect of temperature on the hydration and porosity of Portland cement. Cem. Concr. Res.38(1):1-18. 10.1016/j.cemconres.2007.08.017,31ThoenenT, KulikD. 2003. Nagra/PSI chemical thermodynamic database 01/01 for the GEM-Selektor (V. 2-PSI) geochemical modeling code. PSI, Villingen. 10.1524/ract.2002.90.9-11_2002.805.,32HummelW, BernerU, CurtiE, PearsonFJ, ThoenenT. 2002. Chemical thermodynamic data base (01/01). Nagra/PSI Chemical Thermodynamic Data Base, 1(01):565. DOI: 10.1524/ract.2002.90.9-11_2002.805) |
| gypsum | 0,00E+00 | 0,00E+00 | 0,00E+00 | 5,29E+02 | 1,85E+00 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 1,06E-03 | 1,72E-01 | -2,02E+03 | (30LothenbachB, MatscheiT, MöschnerG, GlasserFP. 2008. Thermodynamic modelling of the effect of temperature on the hydration and porosity of Portland cement. Cem. Concr. Res.38(1):1-18. 10.1016/j.cemconres.2007.08.017,31ThoenenT, KulikD. 2003. Nagra/PSI chemical thermodynamic database 01/01 for the GEM-Selektor (V. 2-PSI) geochemical modeling code. PSI, Villingen. 10.1524/ract.2002.90.9-11_2002.805.,32HummelW, BernerU, CurtiE, PearsonFJ, ThoenenT. 2002. Chemical thermodynamic data base (01/01). Nagra/PSI Chemical Thermodynamic Data Base, 1(01):565. DOI: 10.1524/ract.2002.90.9-11_2002.805) |
| plaster | 0,00E+00 | 0,00E+00 | 0,00E+00 | 3,25E+03 | 1,13E+01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 9,42E-04 | 1,43E-01 | -1,58E+03 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| anhydrite | -1,82E+07 | 0,00E+00 | 0,00E+00 | 8,30E+02 | 3,58E-01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 8,37E-04 | 1,36E-01 | -1,43E+03 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| T | -1,78E+05 | 2,83E+04 | 0,00E+00 | 6,21E+01 | 2,00E+00 | -2,77E-03 | 1,89E-06 | -5,08E-10 | 0,00E+00 | 7,99E-02 | -3,05E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| P | -1,76E+06 | -2,69E+03 | 0,00E+00 | 4,09E+02 | 1,69E+00 | -1,82E-03 | 9,63E-07 | -2,03E-10 | 0,00E+00 | 1,42E-01 | -1,12E+03 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| C | -2,17E+07 | 0,00E+00 | 0,00E+00 | 1,07E+03 | -2,91E-01 | 1,92E-04 | -3,08E-08 | 0,00E+00 | 2,51E-03 | 5,60E-02 | -6,35E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| A | -4,92E+07 | 0,00E+00 | 0,00E+00 | 1,20E+03 | 6,74E-02 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 1,79E-02 | 1,02E-01 | -1,68E+03 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| F | -9,33E+06 | 0,00E+00 | 0,00E+00 | 6,14E+02 | 4,88E-01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | -7,42E-03 | 1,60E-01 | -8,24E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| F>600K | -9,33E+06 | 0,00E+00 | 0,00E+00 | 6,14E+02 | 4,88E-01 | 0,00E+00 | 0,00E+00 | 0,00E+00 | -7,42E-03 | 1,60E-01 | -8,24E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| M | -2,43E+07 | 0,00E+00 | 0,00E+00 | 1,15E+03 | 2,26E-01 | -1,01E-04 | 2,60E-08 | 0,00E+00 | -8,78E-03 | 4,03E-02 | -6,02E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| K | 0,00E+00 | 0,00E+00 | 0,00E+00 | 1,06E+03 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 3,23E-04 | 9,42E-02 | -3,62E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| N | 3,42E+07 | -3,40E+05 | 0,00E+00 | 2,23E+03 | -1,65E+00 | 2,07E-03 | -1,08E-06 | 2,34E-10 | 0,00E+00 | 6,20E-02 | -4,15E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| CO2 | 2,22E+07 | -3,38E+05 | 0,00E+00 | 1,57E+03 | -1,74E-02 | 9,19E-07 | -3,57E-10 | 1,20E-13 | 7,92E-02 | 4,40E-02 | -3,94E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| H2O | 4,78E+08 | -1,11E+06 | 0,00E+00 | 2,15E+03 | 1,06E+00 | -3,16E-04 | 4,35E-08 | -2,23E-12 | 1,21E-02 | 1,80E-02 | -2,42E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| SO2 | -6,89E+06 | 1,18E+05 | 0,00E+00 | -3,06E+02 | 2,86E+00 | -3,26E-03 | 1,88E-06 | -4,37E-10 | 1,68E-03 | 6,40E-02 | -2,97E+02 | (33HaneinT, GlasserFP, BannermanMN.2020. Thermodynamic data for cement clinkering. J. Cement and Concrete Research. 132:106043. 10.1016/j.cemconres.2020.106043) |
| C3S | -2,86E+05 | 0,00E+00 | -1,21E+04 | 1,46E+03 | -1,02E-02 | 0,00E+00 | 0,00E+00 | 0,00E+00 | -9,26E-04 | 2,28E-01 | -2,93E+03 | (33HaneinT, GlasserFP, BannermanMN.2020. Thermodynamic data for cement clinkering. J. Cement and Concrete Research. 132:106043. 10.1016/j.cemconres.2020.106043) |
| C2Sg | -1,37E+07 | 0,00E+00 | 9,63E+03 | 0,00E+00 | 1,24E+00 | -4,74E-04 | 0,00E+00 | 0,00E+00 | -1,19E-03 | 1,72E-01 | -2,32E+03 | (33HaneinT, GlasserFP, BannermanMN.2020. Thermodynamic data for cement clinkering. J. Cement and Concrete Research. 132:106043. 10.1016/j.cemconres.2020.106043) |
| C2Sap | 0,00E+00 | 0,00E+00 | 0,00E+00 | 9,40E+02 | 0,00E+00 | 1,10E-04 | 0,00E+00 | 0,00E+00 | -1,19E-03 | 1,72E-01 | -2,31E+03 | (33HaneinT, GlasserFP, BannermanMN.2020. Thermodynamic data for cement clinkering. J. Cement and Concrete Research. 132:106043. 10.1016/j.cemconres.2020.106043, 34HaasJrJL, RobinsonJrGR, HemingwayBS. 1981. Thermodynamic tabulations for selected phases in the system CaO-Al2O3-SiO2-H2 at 101.325 kPa (1 atm) between 273.15 and 1800 K. J. Phys. Chem. Ref. Data.10(3):575-670. 10.1063/1.555645) |
| C3A | 1,86E+07 | 0,00E+00 | 0,00E+00 | 9,65E+02 | 1,77E-02 | 0,00E+00 | 0,00E+00 | 0,00E+00 | 3,74E-03 | 2,70E-01 | -3,59E+03 | (33HaneinT, GlasserFP, BannermanMN.2020. Thermodynamic data for cement clinkering. J. Cement and Concrete Research. 132:106043. 10.1016/j.cemconres.2020.106043, 35BonnicksonKR. 1955. High temperature heat contents of aluminates of calcium and magnesium. J. Phys. Chem.59(3):220-221. 10.1021/j150525a006.) |
| C4AF | 1,12E+07 | 0,00E+00 | 0,00E+00 | 7,70E+02 | 3,74E-02 | 0,00E+00 | 0,00E+00 | 0,00E+00 | -1,11E-03 | 4,86E-01 | -5,08E+03 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| FeO | -1,38E+07 | -4,20E+04 | 0,00E+00 | 6,39E+02 | -1,15E-01 | 5,07E-05 | -7,86E-09 | 4,21E-13 | 0,00E+00 | 7,18E-02 | -2,72E+02 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| O2 | -8,90E+06 | 1,26E+05 | 0,00E+00 | 2,91E+02 | 1,12E+00 | -1,78E-04 | -5,26E-07 | 2,70E-10 | 5,53E-02 | 3,20E-02 | 0,00E+00 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| N2 | 6,56E+06 | -1,13E+05 | 0,00E+00 | 1,81E+03 | -2,53E+00 | 4,11E-03 | -2,86E-06 | 7,48E-10 | 1,94E-01 | 2,80E-02 | 0,00E+00 | (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) |
| H2Ol | -2690906 | 0 | 485,00939 | 1,1057984 | -0,000848 | 0 | 0 | 0 | 0,00247 | 0,018 | -241,826 | (19TahaY, BenzaazouaM, HakkouR, MansoriM. 2017. Coal mine wastes recycling for coal recovery and eco-friendly bricks production. J. Miner. Eng.107:123-138. 10.1016/j.mineng.2016.09.001) |
3. RESULTS AND DISCUSSION
⌅3.1. About the combustion of coal gangue and coal
⌅Large heaps of waste coal are found at various points in the coal mine, they were generated by the industrial exploitation of the mine during 1930-2000. Figure 1 shows the largest heap, with around 20Mt of coal mining waste rocks. The distinct black color of this coal waste signifies a notable presence of coal content. Upon closer examination, stages of heating a Coal Gangue (CG) sample are illustrated in Figure 1 and it is observed that the dark color disappears at about 620 °C. This was also observed by TGA characterization, with more than 10% of loss of mass due to the combustion of coal residues. The thermal characterizations of CG show an exothermic contribution attributed to the combustion of the coal. According to a mine staff report, the calorific value of gangue (CVCG) is determined to be 3.767 MJ/kg. This value aligns with the literature, where reported values fall within the range of 2.307 to 8.309 MJ/kg (18TahaY, ElghaliA, DerhyM, AmraniM, HakkouR, Benzaazoua, M. 2023. Towards an integrated approach for zero coal mine waste storage: solutions based on materials circularity and sustainable resource governance. Mineral Processing and Extractive Metallurgy Review. 44(6):375-388. 10.1080/08827508.2022.2084733); and the calorific value of this anthracite coal is approximately 31.8 MJ.kg-1. From these values, we can deduce an estimation of 11.8% coal in the coal gangue. Based on the ratio of carbon atoms in the gangue and in the coal reported in Table1, we can estimate a percentage of 12.4% coal in the CG. This estimation aligns with findings from a study by Taha et al. (18TahaY, ElghaliA, DerhyM, AmraniM, HakkouR, Benzaazoua, M. 2023. Towards an integrated approach for zero coal mine waste storage: solutions based on materials circularity and sustainable resource governance. Mineral Processing and Extractive Metallurgy Review. 44(6):375-388. 10.1080/08827508.2022.2084733, 19TahaY, BenzaazouaM, HakkouR, MansoriM. 2017. Coal mine wastes recycling for coal recovery and eco-friendly bricks production. J. Miner. Eng.107:123-138. 10.1016/j.mineng.2016.09.001), where the content of coal in coal gangue is reported to be within the range of 10-15%. Moreover, the study suggests that coal extraction through flotation yields a percentage higher than 6.09%, with an extraction efficiency of 60%.
3.2. Thermal characterizations
⌅3.2.1. DSC characterization of coal gangue CG compared to the pure coal PC
⌅The DSC analyses of coal gangue CG compared to the pure coal PC are presented by in Figure 2; the heating rate was 10 °C/min under air.
The DSC analysis of coal gangue shows an endothermic region occurring between 100°C-150°C attributed to the evaporation of free water. Additionally, two dehydration processes, specifically from gypsum to plaster and then to anhydrite:
and
They are theoretically expected to occur at 106°C and 153°C, respectively: Comparing these theoretical values to those reported in literature (20Kuntze, R.A.2009. Gypsum: Connecting Science and Technology. ASTM MNL 67, ASTM International.), the temperature for the first dehydration aligns closely, with literature reporting 105°C. However, for the second dehydration, the literature reports a slightly higher temperature of 175°C.
At 580°C, the coal gangue (CG) undergoes a polymorphic transformation of quartz: SiO2α → SiO2β.
At 580°C, we have the polymorphic transformation of quartz: SiO2α → SiO2β. Then, in the range, an exothermal effect is observed, mainly attributed to the combustion of coal.
Using the DSC characterization (10°C/min, in air) for CG, we obtain the combustion enthalpy of CG is HCG = - 63.74 J.g-1; this value is used in the thermochemical balance calculations for the cementitious mixture LG. It's important to note that this value represents a minimum estimate, as the observed exothermic effect persists across multiple heating and cooling cycles, indicating incomplete combustion. This persistence is demonstrated in Figure 2(b), where the PC sample exhibits the observable exothermic effect even after five cycles of heating and cooling.
3.2.2. Thermal characterizations of materials CG and LG.
⌅Characterizations of materials LG and CG was done using thermogravimetric TGA and thermo-differential TDA analyses. The graphical representations of these analyses are illustrated in Figure 3.
In the temperature interval of 300-600°C, we observe the presence of the exothermic effect due, mainly, to the combustion of coal in the coal gangue, this effect is also noticeable in the limestone-gangue cementitious mixture LG. At about 850°C (P=1bar), we observe a main endothermic peak attributed to the calcination of the limestone, in the LG mixture:
At this temperature the standard free enthalpy of this equilibrium verifies the condition ∆ca G° (Tca) =0: ∆ca G° (T) = 178.2-0.1601.T (kJ.mol-1), Tca=1113K(840°C).
Similarly, the temperature at which dolomite decomposes is determined, considering a hypothetical single-stage decomposition:
The standard free enthalpy of this equilibrium is:
This reaction occurs at Tdo~=831K (558°C). Adopting a two-step decomposition, as reported in some literatures (21GunasekaranS, AnbalaganG. 2007. Thermal decomposition of natural dolomite. J. B. Mater. Sci.30:339-344. 10.1007/s12034-007-0056-z); the process involves the initial formation of MgO, CO2, and calcite at the temperature of 750°C, followed by the subsequent decarbonation of calcite around 900°C.
A DSC characterization has been done under argon atmosphere for the mixtures LG and RR; therefore, the thermal effects excluding combustions can be seen; thus, we will consider that the combustion phenomenon, in the LG mixture, is predominant in front of other reactions taking place in the interval 300-600°C. The comparative results are reported in Figure 4.
The thermal behavior of the two mixtures, LG and RR, exhibits notable similarities. Two endothermic phenomena are observed:
- Water Losses: Involving both free water and the dehydration of gypsum.
- Calcination of Calcite: The endothermic process where calcium carbonate (calcite) is converted into calcium oxide and carbon dioxide
It's important to note that limestone typically contains dolomite (calcium magnesium carbonate), which undergoes calcination around 500°C. Additionally, around 550°C, the decomposition of clays initiates through dehydrations, dehydroxylations, and decompositions (22GualtieriAF, FerrariS. 2006. Kinetics of illite dehydroxylation. J. Phys.Chem.Miner, 33:490-501. 10.1007/s00269-006-0092-z). These reactions continue until approximately 900°C.
The main thermochemical events involve the heating of all substances, and the main reactions are: decomposition of calcite and dolomite, combustion of coal residues, decomposition of clays and formation of the four phases of clinker. In addition to these primary events, there are, also, minor thermal occurrences like the alpha-beta transition of silica and dehydration of gypsum, and polymorphism in phases of clinker. In the literature, these main physicochemical phenomena take place in temperature intervals (23KaanteeU, ZevenhovenR, BackmanR, MikkoH. 2003. Modelling a cement manufacturing process to study possible impacts of alternative fuels. Ciments Betons Platres Chaux. 46-52.). Due to the complexity of these physicochemical phenomena, detailed studies are challenging, and approximations become necessary. To facilitate analysis, relevant temperatures Tr=1,2… are introduced, corresponding to the key reactions. For instance, the combustion phenomena of coal begins at 727K as showed in thermal characterization and in literature (24MohalikNK, MandalS, RaySK, KhanAM, MishraD, PandeyJK. 2022. TGA/DSC study to characterise and classify coal seams conforming to susceptibility towards spontaneous combustion. Int. J. Min. Sci. Technol.32(1):75-88.10.1016/j.ijmst.2021.12.002); the αβ transformation of quartz occurs at 848K. T0=298K is the initial temperature and T7=1723K is the final temperature of clinkerization.
We study mass and energy balances during the production of LG* and RR* clinkers, from raw materials. These balances are examined for different stages as presented in the Table 5:
Stage I: the first stage of decomposition.
Stage II: the second stage of clinkerization.
Overall balance:
In the energy balance analysis, the focus is on the algebraic determination of energies required to produce 1 kg of clinkers denoted LG* and RR*. We are interested in the difference of the energies , in this way we will reduce the impact of approximations done. Moreover, since the limestone LR is the same for both raw feeds, certain terms in the balance are inherently similar.
We present the balance involving the main thermochemical transformations. These include the heating energies for the main compounds in the temperature intervals considered (Formula F1), as well as the enthalpies of the chemical reactions listed in Table 5 (Formulas F2 and F3). The heats of change of physical state or of polymorphism are included with those of the reactions presented in Table 5.
Where , and , and are respectively the differences in energies of: reactions and heating of solids s, liquids l and gases g in stages I and II.
The terms of heating are of the form: in general, for a given substance i, the massic heat capacity is:
The data on the specific heat capacity Cp,i shown in Table 5 is given in literature, mainly from the NASA/TP-2002-211556 report (25McBrideBJ. 2002. NASA Glenn coefficients for calculating thermodynamic properties of individual species. National Aeronautics and Space Administration, John H. Glenn Research Center at Lewis Field.) and other sources (26-36). The specific heat capacity of coal is reported based on the work of Lesniak et al and Eisermann et al (28LeśniakB, ŁukaszS, JakubinaG. 2013. Institute for chemical processing of coal, zabrze, poland. The determination of the specific heat capacity of coal based on literature data. Chemik. 67(6):560-571., 29EisermannW, JohnsonP, Conge, WL. 1980. Estimating thermodynamic properties of coal, char, tar and ash. Fuel Process. Technol.3(1):39-53. 10.1016/0378-3820(80)90022-3), it is formed by two contributions of fixed carbon FC (fraction mC, CpC) and volatile matter (fraction 1- mC, CpV). In our case, for anthracite, the value is FC=93 and VM =5 (37ChellaiH, EssamoudR. RjimatiE. 2011. Le bassin houiller de Jerada (Chaîne des Horsts, Maroc oriental) / The Jerada Coal Basin (Horst Chain, Eastern Morocco). 1(556-564):331-335.).
Using formula (F1), we calculate the heating of all solids, liquids and gases; the gases CO2, H2O, SO2, O2, and N2 are involved in the reactions of calcinations of calcite and dolomite, dehydrations of gypsum, decomposition of clays and combustion of coal.
3.3. Mass and energy balances for the production of 1 kg LG* and RR* clinkers
⌅XRF and XRD data are cross-referenced to give the mass fractions of the compounds in the three materials CG, LR and RR. Amorphous compounds, not visible, are considered, and total oxide content from XRD is adjusted to match XRF values.
3.3.1. Massic compositions of coal gangue CG, limestone LR, and cement mixtures LG and RR
⌅The coal gangue consists of various compounds, including quartz, illite, muscovite, clinochlore, gypsum, ferric oxide (F), and sulfur dioxide (SO2). These compounds are derived, in part, from pyrite, which makes up approximately 2-5% of the coal (38DarmaneY, AlaouiA, KitaneS, BennajahM, DaramyA, CherkaouiM. 2009. Recycling the slagheap of an old coal mine (Morocco). Sep. Purif. Technol.68(1):125-128. 10.1016/j.seppur.2009.04.026). In the analysis of coal gangue (CG), adjustments are made to align the total amount of oxides ‘CSAMKS’ given by XRD (94.51%), to that given by XRF (79.24%). The experimental percentage of coal in the gangue is estimated using the loss on ignition L.O.I, which is formed by dehydration of the three water-bound in clays and gypsum (2.39 %). This estimation suggests a maximum coal content approaching 11.85%, consistent with previous estimates. For different samples of gangue, the losses on ignition varied between 10.6 and 14.22, corresponding to a coal content between 8.8 and 11.8%. For the mass balance, the study uses the conservative assumption of 11.85 % coal in GC.
The gases released during combustion include water vapor, carbon dioxide, and sulfur dioxide. In this study, we use the experimental amounts of CO2 and SO2 estimated by Taha et al: 0.243 kg CO2 and 0.00628 kg SO2 per kg of gangue and we obtain the amount of water vapor released from an approximate ‘formula’ of the anthracite. Coal can be modeled by the formula ‘CxHyOzStNu’. Jerada's anthracite analysis indicates volatile matter (VM) at about 5%, fixed carbon (FC) at approximately 93%, and ash (ASH) at 5.5
A model proposed by Shen et al. (39ShenJ, ZhuS, LiuX, ZhangH, Tan, J. 2010. The prediction of elemental composition of biomass based on proximate analysis. Energy Convers. Manag.51(5):983-987. 10.1016/j.enconman.2009.11.039) and reported in Lawal et al’s study (40LawalAI, AladejareAE, OnifadeM, BadaS, IdrisMA. 2021. Predictions of elemental composition of coal and biomass from their proximate analyses using ANFIS, ANN and MLR. Int. J. Coal Sci. Technol.8:124-140.) provides an approximate formula (‘C1H 1.158O 0.418 NuSt’). This model predicts that the combustion of 1 ton of CG produces 62.1 kg H2O and 263 kg CO2. The optimal limestone-gangue mixture retained has a composition of 17.1% coal gangue. The massic composition of coal gangue is given in Table 7.
We complete the determination of the mass fractions of the main compounds forming the LR and RR materials. The total amount of carbon dioxide from calcite and dolomite, and from the XRD analysis we obtain the ratios of calcite to dolomite and to quartz: 6.58 and 11.9 respectively. For the LR material, the total 96.089% in CMC given by XRD is adjusted to 74.865% given by XRF; and for the RR material the total of L.O.I is composed of carbon dioxide and water of clay, maintaining the same ratio, 6.58, of CC to CMC. All the results are reported in Table 7 and used to determine the energies involved in the comparative energy balance of the two cement mixtures LG, and RR.
3.3.2. Energy balance
⌅The masses of the gases I: CO2, H2O, SO2, O2 and N2 concerned by this phase I are respectively: 66.14, 16.81, 1.682, 55.04, 192.7 gram for 1 kg of clinker. For the energy balance, the solids in this phase are denotes solids I.
During the clinkering phase, the oxide compositions (XRF) of the two mixtures are used to calculate the saturation factor lime moduli LSF=100C/ (2.8S+1.65A+0.35F), the ratio SR=S/ (A+F) as well as the four phases of the clinker and possibly, the percentage of free lime Cf. In order to calculate compositions of the four phases, we use the Taylor’s matrix (41ShimSH, LeeTH, YangSJ, NoorNBM, KimJHJ. 2021. Calculation of cement composition using a new model compared to the bogue model. Materials. 14(16):4663. 10.3390/ma14164663) leading to the system of equations:
After calculating the total mass of the four phases (solids II1), we assumed that the difference with the mass of the C-S-A-F oxides was that of the free lime, in order to adjust the total mass. The theoritical percentage of free lime can be approximated using Fundal’s empirical equation :% CaO(free)=0.31 (LSF-100)+2.18 (SR-1.8)+(0.73 Q45 + 0.33 C125 + 0.34 Aq). The total amount of C2S formed at 1173K will be distributed between the final masses of C2S and that of C3S which will form at 1573 K; these two masses are obtained by the Taylor model above.
In Table 8, we give the massic composition of these four phases, anhydrite, minor oxides and the gases. All the results are used to determine the energies involved in the comparative energy balance of the two cement mixtures LG and RR.
Finally, in Table 9, we present the energy balance of the reactions and the heating, then we deduce the overall balance; these results are generated by the program Python.
4. Conclusions
⌅This study highlights the energy contribution of coal combustion as a prime material for clinker production. Regarding energy considerations, the following conclusions can be drawn:
- 1. The combustion of the coal contained in the gangue is more difficult than that in the open air, because this combustion takes place within the LG cement mixture; for this reason, the value of the enthalpy of combustion, during the DSC, can be used to estimate the energy input, but the heating rate is very high in the cement industry where the calcination temperature rises to 900°C in a few minutes. The efficiency of this energy input is improved by favorable kinetics resulting from the small particle size and consequently larger surface area.
- 2. Coal combustion generates heat within the limestone-gangue mixture itself during the first phase of clinker’s elaboration effectively improving calcination; this internal heat is added to that received by conduction and convection. On the other hand, using the approximation of an adiabatic transformation, it is possible to present the energy gain as a lowering of the working temperature.
- 3. Incorporating coal waste in cement production offers economic and ecological benefits, potentially making the industry more sustainable. As this practice gains traction, it's crucial to address associated social challenges. Effective collaboration among industry stakeholders, policymakers, and communities is key for ensuring both sustainability and social responsibility in cement production.