Research on a Carbon Emission Calculation Model and Method for an Underground Fully Mechanized Mining Process

: With the ratiﬁcation of the Paris Agreement at the Paris Climate Conference, reducing carbon emissions has become a global interest. Coal is one of the main industries causing carbon emissions; thus, quantifying carbon emissions from coal mining is an important step in reducing these emissions. Firstly, based on the life cycle idea, in this paper, we deﬁne the Carbon Emission Boundary of the fully mechanized coal mining method. Secondly, the carbon emission accounting model (B-R model) of fully mechanized coal mining is established, which includes the total amount of carbon emissions and the carbon emissions of each mining link during the mining process. The Fifth-II mining area of the Jinda Coal Mine in Tengzhou City is taken as an example. We collect the relevant data on carbon emissions in the mining process of the Jinda Coal Mine, and the B-R model is used to obtain the carbon emissions in the mining process of this mining area. Finally, the feasibility of the B-R model is further veriﬁed according to the international authoritative carbon emission IPCC calculation method and the China Coal Production Enterprises Greenhouse Gas Emissions Accounting Methodology and Reporting Guide . The results show that the B-R model in this paper is feasible and that the greatest amount of carbon emissions arises from the coal breaking link and coal transportation, which provides a basis for other coal mines to calculate carbon emissions. The B-R model lays a foundation for coal mines to formulate a carbon emission reduction system.


Introduction
Due to the development of the economy, technology, and population growth, human consumption of coal is increasing day by day, which is causing damage to the carbon balance of the environment [1][2][3].In recent years, the Chinese Communist Party and the government have attached great importance to environmental issues, and in September 2020, they proposed to adopt stronger policies and measures to achieve the "double carbon target" [4][5][6].The double carbon target, which takes into account carbon peaking and carbon neutrality, requires that China strives to achieve carbon peaking by 2030 and carbon neutrality by 2060.According to statistics, coal accounted for 56.7% of China's energy mix in 2020, and coal mining, as the upstream link in coal consumption, plays an important role in the overall carbon emissions generated by coal consumption [7].With the continuous development of technology, the underground coal mining process has gradually been mechanized and developed into fully mechanized coal mining technology, which is now widely used [8][9][10].At present, there is very little research on carbon emission accounting for fully mechanized coal mines, and there is no accounting model to quantify the carbon emissions of each process in a fully mechanized working face.Firstly, this paper uses life cycle thinking to analyze the carbon emission sources of each process step in the mining process.Secondly, this paper establishes a carbon emission accounting model for each process step and finally uses real cases to verify the feasibility of the model and calculate the carbon emissions of each step in the example.This paper is of great significance for the realization of low carbon fully mechanized coal mining.

Fully Mechanized Carbon Emission Boundary Establishment
This paper uses the life cycle idea to determine the carbon emission boundary of the fully mechanized process under different geological conditions.The carbon emission boundary cannot be determined without the identification of carbon emission sources [11,12].The purpose of this study was to calculate the amount of CO 2 released by each source of carbon emissions throughout the underground coal mining process, so the carbon emissions were calculated using functional units of CO 2 which has the advantage of allowing all parts of the mining system to be analyzed and compared on the same basis.A uniform unit of tons (t) is used in this paper.

Carbon Emission Source Determination
The main components of the CO 2 emission boundary of underground coal mining include: energy use, escape of greenhouse gases, and spontaneous combustion of coal gangue and raw coal [13].The CO 2 emission boundary also contains some auxiliary elements: the drainage system, ventilation system, energy consumption for staff daily work, office area consumption, transportation system, etc. (as shown in Table 1).

Carbon Dioxide Emission Boundary Establishment
According to the information in Table 1, the carbon dioxide emission boundary of coal mining based on life cycle thinking is shown in Figure 1, combined with the specific coal mining processes and the preparation stage, operation stage, and closing stage of the coal mining process.

Carbon Dioxide Emission Boundary Establishment
According to the information in Table 1, the carbon dioxide emission boundary of coal mining based on life cycle thinking is shown in Figure 1, combined with the specific coal mining processes and the preparation stage, operation stage, and closing stage of the coal mining process.Based on the CO boundary, this paper will focus on the energy use aspects of the system, which are as follows: (1) Preparation stage.The preparation stage includes the up-and-down hill mining area type, panel area type, strip type, and other preparation methods, in which the main materials used are concrete, steel, etc., and also related pioneering machines (road header, excavator, etc.).Such carbon emissions can be calculated by carbon emission factors.
(2) Operation stage.The operation phase includes the entire coal mining process: coal breaking, coal loading, coal transportation, roof support, and treatment of the coal mine gob area.The entire process generates large amounts of carbon dioxide.CO is produced by the use of machines, electricity, escape, and spontaneous combustion under different geological conditions.The machines often used in the process are coal mining machines, scraper conveyors, and hydraulic pillars, which use electricity and produce some carbon emissions.Lubricants are also applied to the machines, which also produce small amounts of CO emissions [14][15][16].(3) Final stage.A lot of equipment and machinery are dismantled or removed during this phase and the net is laid in the mine.All of these construction phases generate some greenhouse gas emissions.Based on the CO 2 boundary, this paper will focus on the energy use aspects of the system, which are as follows: (1) Preparation stage.The preparation stage includes the up-and-down hill mining area type, panel area type, strip type, and other preparation methods, in which the main materials used are concrete, steel, etc., and also related pioneering machines (road header, excavator, etc.).Such carbon emissions can be calculated by carbon emission factors.(2) Operation stage.The operation phase includes the entire coal mining process: coal breaking, coal loading, coal transportation, roof support, and treatment of the coal mine gob area.The entire process generates large amounts of carbon dioxide.CO 2 is produced by the use of machines, electricity, escape, and spontaneous combustion under different geological conditions.The machines often used in the process are coal mining machines, scraper conveyors, and hydraulic pillars, which use electricity and produce some carbon emissions.Lubricants are also applied to the machines, which also produce small amounts of CO 2 emissions [14][15][16].(3) Final stage.A lot of equipment and machinery are dismantled or removed during this phase and the net is laid in the mine.All of these construction phases generate some greenhouse gas emissions.

Accounting for Carbon Emissions from Underground Coal Mining
Through the above-mentioned analysis of CO 2 emission boundary establishment, combined with the actual situation of underground coal mining, this paper simplifies the establishment of the model.
(1) The preparation stage mainly involves coal mine tunnel layout and mine construction.
There are many factors to consider in coal mine tunnel layout: roadway design, policy, funding, geology, and labor.The greenhouse gas emissions from these factors are negligibly small.The main sources of carbon emissions considered for mine construction are machine use and fugitive emissions.There are overlaps between mine construction and coal mining machine use and the carbon emissions from fugitive emissions are low.Therefore, they are not included in the calculation of carbon emissions in this paper.(2) Emissions from the finishing phase of netting, evacuation, or removal of the equipment are small and difficult to quantify, therefore, they are not included in the calculation.

Accounting for Carbon Emissions from Fully Mechanized Processes
Based on the above analysis, the B-R model is established in this paper.The B-R model proposed in this paper is an accounting model for the carbon emissions generated in the fully mechanized processes.
The fully mechanized mining process consists of four specific processes: falling and mining coal, transporting coal, roof support, and treatment of the coal mine gob area [17].The geological conditions of a coal mine are an important factor affecting all aspects of underground coal mining.Different geological conditions can lead to different mining methods [13] and thus have an impact on carbon emissions, therefore, this paper needed to analyze the four processes based on geological conditions before building the model in order to ensure the correct source of carbon emissions.The idea was to account for the total carbon emissions and then find the carbon equivalent emissions per ton of coal mined.The total carbon emission accounting model (B-R model) for the entire fully mechanized mining process was based on: where E refers to total carbon emissions from underground fully mechanized mining processes (t/t); E f refers to carbon emissions from falling and mining coal (t/t); E t refers to carbon emissions from transporting coal (t/t); E s refers to carbon emissions from roof support (t/t); and E g refers to carbon emissions from the treatment of the coal mine gob area (t/t).

Carbon Emissions from Falling and Mining Coal
As shown in Table 2, this paper analyzed the sources of carbon emissions from falling and mining coal.In general, the shearer completes the cutting and loading of coal, and on the working face, the coal miner cuts the top coal in the front roller and the bottom coal in the rear roller [17,18].The advantages of this method are safe operation, less coal dust, and a good coal loading effect, therefore, it is widely used.Spray pumping stations are also used in this process.The spraying pumping stations provide a power source for the coal breaking and loading process where the dust needs to be sprayed.Electricity is used extensively in the processes of falling and mining coal.Falling coal causes the seam to fracture, which in turn leads to the release of stored gas from the surrounding rock seam; this is known as spillage.The fugitive gases are essentially methane, therefore, the calculation is based on methane.Lubricants are used as auxiliary supplies in the operation of machines and are omitted in this article due to their small quantity.In this paper, a classification of the calculated electricity carbon emissions of the B-R model is presented.There are two sources of electricity for coal mines: self-generation and purchased electricity, and different sources of electricity have different impacts on carbon emissions.

Thickness of coal seam
A thin seam shearer is suitable for mining coal seams with minimum seam thickness below 1.3 m.A medium-thick coal seam mining machine is suitable for mining coal seams with a thickness from 1.3 to 3.5 m, and is the mainstream shearer, with many varieties and wide uses.A thick seam shearer is suitable for mining coal seam with a highest seam thickness of more than 3.5 m; the model basically was derived from a medium-thick seam shearer and the maximum cutting height can reach 4.5~5 m at present.

Dip angle of the coal seam
According to the applicable coal seam, the inclination angle can be divided into three types of shearers: gently inclined seam, pitching seam, and steeply pitching seam.A gently inclined seam mining machine is applicable to coal seams with inclination angle below 25 • , and most of the shearers belong to this category.A pitching seam mining machine is suitable for coal seams with inclination angle of 25 to 45 • , and this kind of shearer is generally modified on the basis of a slow-slope coal seam mining machine, in the process of development of a steeply pitching seam shearer.

Spray pumping stations
Dust spraying Geological requirements are not obvious.

Spillage
Methane spill Geological requirements are not obvious.
Energies 2022, 15, 2871 5 of 15 E f of the B-R model is expressed as: where N refers to the power of the spray pumping stations (kw); M refers to the power of the shearer (kw); α refers to carbon emission factors for electricity (t/kwh); S CH 4 refers to the total amount of fugitive CH 4 over time h (t); GWP CH 4 refers to the Global Warming Potential of CH 4 (t/t); and Q m refers to the average coal mining capacity of the shearer (t/h).In B-R, the power of spray pumping stations [18] is expressed as: where R refers to the flow rate of the spray pumps (m 3 /h); H refers to the head of the spray pump (m); η refers to the efficiency of the spray pumps.In B-R, the power of spray pumping stations is expressed as: where k b refers to the shearer cutting efficiency factor; B refers to the dinting depth of the shearer (m); D refers to the cutting height of the shearer (m); and V max refers to the maximum velocity of the shearer (m/min).In B-R, the average coal mining capacity of the shearer is expressed as: where P refers to the daily production of the working face (t); L refers to the length of the working face (m); i refers to the ratio of the coal cutting speed and haulage speed of the shearer; L m refers to the center distance between two drums (m); K refers to the average daily shearer start rate; T 1 refers to the daily production time of the working face (min); C refers to the recovery rate of the working face; T d refers to the return time of shearer (min); H refers to the mining height (m); and γ refers to the coal seam capacity (t/m 3 ).The global warming potential is a parameter of the radiative impact of a certain mass of greenhouse gases with an equal mass of CO 2 over a period of time.Generally speaking, it means that the mass of greenhouse gases other than CO 2 is converted into a uniform accounting for the mass of CO 2 .The specific GWP values for each GHG are specified in the IPCC report.

Carbon Emissions from Transporting Coal
As shown in Table 3, the main equipment of the coal transportation link of fully mechanized coal mining is the scraper conveyor to transport coal, followed by the loader conveyor.The mining equipment train is also a piece of equipment for coal transportation, which is installed on the main road or secondary road of the coal mining face.Spillage will continue to occur in the coal transportation segment in general, although it is slower than emissions from the coal seam crushing stage.
E t of the B-R model is expressed as: where F refers to the power of the scraper conveyor, (kw); W refers to the power of the band conveyor (kw); G refers to the power of the hoister (kw); and Z refers to the power of the mining equipment train (kw).In B-R, the power of each transporter and the mining equipment train [16] is expressed as: where L refers to the length of working face (m); H a refers to the vertical lifting height (m); Y a refers to the adhesive gravity of the mining equipment train, (t); Y b refers to the weight of the mining equipment train (t); ξ refers to the utilization factor of the mining equipment train; and v refers to the operating speed of the mining equipment train (km/h).ω refers to the unit basic resistance of the mining equipment train (N).In general, ω = a + bv + cv 2 ; a and b are the mechanical resistance, and c is the air resistance.

Carbon Emissions from Roof Support
As shown in Figure 2, the main equipment used to support the roof is the self-moving hydraulic bracket, which has corresponding models for simple, complex, and extremely complex geological types according to the geological conditions.As the main carbonemitting equipment in this segment, the emulsion power pack provides energy for the hydraulic support, which is mainly used for the high-pressure fluid supply of the hydraulic bracket in the underground fully mechanized coal mining working face and is also used for spray cooling and dust removal of the coal mining machine and bracket.The use of lubricants is significant in this segment, which has an impact on carbon emissions.Negligible spillage in the support links occurs.E of the B-R model is expressed as: where P is the pressure of the emulsion power pack (MPa); Q refers to the measured flow rate of liquid from the emulsion pumps (L/min); ρ refers to the total efficiency of the emulsion power pack (value range 0.7-0.8);I refers to the supply of lubricants (L); λ refers to the carbon emission factor of the lubricant use process (kg/L); and λ is carbon emission factor of the waste process treatment of the lubricant (kg/L).
The working pressure of the emulsion pump can be calculated according to the requirements of the initial bracing force, the pulling force, and the push-conveyor force.The pressure of the emulsion power pack is expressed as: P = P P P where P refers to the working pressure of the emulsion pump calculated according to the initial support force requirement; P refers to the initial support force of the hydraulic support (kN); K refers to the correction factor of the column, K 1/cosθ, where θ refers to the inclination of the column; D1 refers to the column bore of the hydraulic support (m); a refers to the number of columns of the bracket; P refers to the working pressure of the emulsion pump calculated according to the pulling force (MPa); P refers to the pulling force of the hydraulic bracket (kN); D refers to the column bore of the pulling force (m); P refers to the working pressure of the emulsion pump, calculated according to the push-conveyor force (MPa); P refers to the push-conveyor force of hydraulic bracket (kN); d refers to the piston rod diameter of the push-pull cylinder (m); and σ refers to the pressure loss factor (1.1~1.2).The measured flow rate of the emulsion power pack is expressed as: E s of the B-R model is expressed as: where P b is the pressure of the emulsion power pack (MPa); Q b refers to the measured flow rate of liquid from the emulsion pumps (L/min); ρ refers to the total efficiency of the emulsion power pack (value range 0.7-0.8);I refers to the supply of lubricants (L); λ 1 refers to the carbon emission factor of the lubricant use process (kg/L); and λ 2 is carbon emission factor of the waste process treatment of the lubricant (kg/L).
The working pressure of the emulsion pump can be calculated according to the requirements of the initial bracing force, the pulling force, and the push-conveyor force.The pressure of the emulsion power pack is expressed as: where P i refers to the working pressure of the emulsion pump calculated according to the initial support force requirement; P 1 refers to the initial support force of the hydraulic support (kN); K refers to the correction factor of the column, K ≈ 1/cos θ, where θ refers to the inclination of the column; D 1 refers to the column bore of the hydraulic support (m); a refers to the number of columns of the bracket; P f refers to the working pressure of the emulsion pump calculated according to the pulling force (MPa); P 2 refers to the pulling force of the hydraulic bracket (kN); D 2 refers to the column bore of the pulling force (m); P c refers to the working pressure of the emulsion pump, calculated according to the pushconveyor force (MPa); P 3 refers to the push-conveyor force of hydraulic bracket (kN); d refers to the piston rod diameter of the push-pull cylinder (m); and σ refers to the pressure loss factor (1.1~1.2).The measured flow rate of the emulsion power pack is expressed as: where t refers to the fluid supply time (min); d 1 refers to the cylinder diameter of the lifting jack (m); l 1 refers to the moving distance of the moving frame jack (m); l 2 refers to the short column bore of beam support (m); d 3 refers to the column bore (m); l 3 refers to the distance of column shift frame (m); q refers to the number of columns raised and lowered; We refers to the support width of one stand (m); t 1 refers to moving rack time (min); t 2 refers to the operation adjustment time (min); V y refers to the shift speed of the bracket (m/min); and V c refers to shearer working traction speed (m/min).

Carbon Emissions from the Treatment of the Coal Mine Gob Area
As shown in Figure 3, there are four main methods to deal with the coal mine gob area: the total fallen method, support method, slow sinking method, and filling method [19].The carbon emissions from the support method and the slow sinking method are too small and difficult to quantify, therefore, they are not included in the calculation of carbon emissions from the process of treatment of the coal mine gob area.Explosives are used in the total collapse method to perform blasting, which generates some carbon emissions.Equipment is used to fill the extraction area in the filling method.Treatment of the coal mine gob includes the extraction of gas from the extraction zone, which causes some gas spillage, which also contributes to the carbon emissions.The process of the treatment of the coal mine gob will result in some spontaneous combustion.
Energies 2022, 15, x FOR PEER REVIEW 9 of 16 where t refers to the fluid supply time (min); d refers to the cylinder diameter of the lifting jack (m); l refers to the moving distance of the moving frame jack (m); l refers to the short column bore of beam support (m); d refers to the column bore (m); l refers to the distance of column shift frame (m); q refers to the number of columns raised and lowered; We refers to the support width of one stand (m); t refers to moving rack time (min); t refers to the operation adjustment time (min); V refers to the shift speed of the bracket (m/min); and V refers to shearer working traction speed (m/min).

Carbon Emissions from the Treatment of the Coal Mine Gob Area
As shown in Figure 3, there are four main methods to deal with the coal mine gob area: the total fallen method, support method, slow sinking method, and filling method [19].The carbon emissions from the support method and the slow sinking method are too small and difficult to quantify, therefore, they are not included in the calculation of carbon emissions from the process of treatment of the coal mine gob area.Explosives are used in the total collapse method to perform blasting, which generates some carbon emissions.Equipment is used to fill the extraction area in the filling method.Treatment of the coal mine gob includes the extraction of gas from the extraction zone, which causes some gas spillage, which also contributes to the carbon emissions.The process of the treatment of the coal mine gob will result in some spontaneous combustion.(1) Treatment of the coal mine gob area in the total fallen method.
E of the B-R model in the total fallen method is expressed as: (1) Treatment of the coal mine gob area in the total fallen method.
E g of the B-R model in the total fallen method is expressed as: where θ i refers to the carbon emission factors for explosives i (t/t); U i refers to the amount of explosives i (t); NCV j refers to the net calorific value of j spontaneously combusted raw materials j (typically gangue or raw coal) (TJ/Gg); CR j refers to the carbon content of spontaneous combustion raw materials j (kg/TJ); CL j refers to the quantity of spontaneous combustion raw materials j (t); β refers to the carbon oxidation factor (usually the default is 1); S CH 4 refers to the total amount of gas emissions (t); GWP CH 4 refers to the gas GWP (CO 2 t/t); F q refers to mine winch traction (kN); V q refers to the winch speed for towing heavy vehicles (m/s); and γ refers to the transmission efficiency of the winch (0.8).
A carbon balance method is used to calculate the carbon emission factors for explosives.The specific steps are as follows [20][21][22][23]: 1 Calculating the chemical equation for the explosive.The chemical formula of explosives is mainly composed of four elements: carbon, hydrogen, oxygen, and nitrogen.As coal mines mostly use a mixture of explosives, calculations are based on a mixture of explosives.
a. Let the chemical equation of mixed explosives be C c H h O o N n ; b.Assuming that the mixed explosive is composed of n components, the chemical equation of each component is , and the proportion of component i is ϑ i (%). 2 Through the chemical equation of the above explosive, the carbon emission factor was calculated.The formula is expressed as: where c i refers to the number of moles of carbon elements in the i-th group (1 mol); h i refers to the number of moles of hydrogen elements in the i-th group (1 mol); o i refers to the number of moles of oxygen elements in the i-th group (1 mol); and n i refers to the number of moles of nitrogen elements in the i-th group (1 mol).
(2) Treatment of the coal mine gob area in the filling method.
Filling methods to treat goaf can be divided into two categories: (similar) paste filling technology and (super) highwater filling technology.Both technical methods will produce fugitive carbon emissions during the process. 1 (Similar) Paste filling technology Paste filling is used to process several materials into a mixture using a conveying system to enter the gob to complete the filling, as shown in Figure 4.
Energies 2022, 15, x FOR PEER REVIEW 10 of 16 where θ refers to the carbon emission factors for explosives i (t/t); U refers to the amount of explosives i (t); NCV refers to the net calorific value of j spontaneously combusted raw materials j (typically gangue or raw coal) (TJ/Gg); CR refers to the carbon content of spontaneous combustion raw materials j (kg/TJ); CL refers to the quantity of spontaneous combustion raw materials j (t); β refers to the carbon oxidation factor (usually the default is 1); S refers to the total amount of gas emissions (t); GWP refers to the gas GWP (CO t/t); F refers to mine winch traction (kN); V refers to the winch speed for towing heavy vehicles (m/s); and γ refers to the transmission efficiency of the winch (0.8).
A carbon balance method is used to calculate the carbon emission factors for explosives.The specific steps are as follows [20][21][22][23]: ① Calculating the chemical equation for the explosive.The chemical formula of explosives is mainly composed of four elements: carbon, hydrogen, oxygen, and nitrogen.As coal mines mostly use a mixture of explosives, calculations are based on a mixture of explosives.
a. Let the chemical equation of mixed explosives be C H O N ; b.Assuming that the mixed explosive is composed of n components, the chemical equation of each component is C H O N (i ∈ n,), and the proportion of component i is ϑ (%).
② Through the chemical equation of the above explosive, the carbon emission factor was calculated.The formula is expressed as: where c refers to the number of moles of carbon elements in the i-th group (1 mol); h refers to the number of moles of hydrogen elements in the i-th group (1 mol); o refers to the number of moles of oxygen elements in the i-th group (1 mol); and n refers to the number of moles of nitrogen elements in the i-th group (1 mol).
(2) Treatment of the coal mine gob area in the filling method.
Filling methods to treat goaf can be divided into two categories: (similar) paste filling technology and (super) highwater filling technology.Both technical methods will produce fugitive carbon emissions during the process.
① (Similar) Paste filling technology Paste filling is used to process several materials into a mixture using a conveying system to enter the gob to complete the filling, as shown in Figure 4.  E g of the B-R model in paste filling technology is expressed as: Energies 2022, 15, 2871 10 of 15 where ρ refers to the energy factor of the crusher; Q p refers to crusher productivity (t/h); r p refers to the crushed product particle size coefficient; δ a refers to the belt conveying efficiency; δ b refers to the crusher efficiency; N P refers to the criterion number of the blender power; ρ d refers to the system density; V j refers to stirring speed (r/s); d j refers to the paddle diameter (m); k 1 refers to the pump station motor reserve factor; q b refers to the flow rate conveyed by the pumping station (L/s); p k refers to the total head of filling paste conveying (m); ε v refers to pump station volumetric efficiency; and ε j refers to the total mechanical efficiency. 2 (Super) Highwater filling technology The flow of high-water filling is shown in Figure 5.
where ρ refers to the energy factor of the crusher; Q refers to crusher productivity (t/h); r refers to the crushed product particle size coefficient; δ refers to the belt conveying efficiency; δ refers to the crusher efficiency; N refers to the criterion number of the blender power; ρ refers to the system density; V refers to stirring speed (r/s); d refers to the paddle diameter (m); k refers to the pump station motor reserve factor; q refers to the flow rate conveyed by the pumping station (L/s); p refers to the total head of filling paste conveying (m); ε refers to pump station volumetric efficiency; and ε refers to the total mechanical efficiency.② (Super) Highwater filling technology The flow of high-water filling is shown in Figure 5. E of the B-R model in highwater filling technology is expressed as: where Q refers to the spiral machine conveying capacity (t/h); ω refers to the spiral machine resistance coefficient; L refers to conveying length (m); H refers to vertical height in inclined arrangement (m); D refers to the spiral diameter (m).

Case Introduction
The Fifth-II mining area of the Jinda Coal Mine is located in Tengzhou City, Shandong Province, China (as shown in Table 4).The average east-west width of the Fifth-II mining area is about 0.66 km, the average north-south length is about 2.9 km, and the total area is about 1.94 km 2 .E g of the B-R model in highwater filling technology is expressed as: where Q g refers to the spiral machine conveying capacity (t/h); ω g refers to the spiral machine resistance coefficient; L g refers to conveying length (m); H g refers to vertical height in inclined arrangement (m); D g refers to the spiral diameter (m).

Case Introduction
The Fifth-II mining area of the Jinda Coal Mine is located in Tengzhou City, Shandong Province, China (as shown in Table 4).The average east-west width of the Fifth-II mining area is about 0.66 km, the average north-south length is about 2.9 km, and the total area is about 1.94 km 2 .The Fifth-II mining area is a wide and gentle monoclinic structure, and the dip angle of the coal bed is 2-10 • , with an average of 6 • .The faults with a drop of less than 5 m in this mining area are relatively developed, and they are all normal faults.According to the above geological conditions, the Fifth-II mining area adopts the fully mechanized coal mining method (toward longwall mining).In the calculation process, the selected time was 24 h working time, and four digits after the decimal point were reserved.In the Fifth-II mining area, the roof management of the working face adopts the full caving method and adopts hydraulic support to control the roof.

Calculation Results of Carbon Emissions in the Fifth-II Mining Area Based on the B-R Model
The coal mining capacity of the Fifth-II mining area is 59.47 t/h, and the electricity consumption depends on the purchase of electricity from the grid.Tables 5 and 6 show the carbon emissions calculated according to the actual situation of the Fifth-II mining areas using the B-R model established above.

Result and Discussion
The total carbon emission in the five mining areas is 41,670.953tons, and the carbon emissions per ton of coal mining amount to 29.196 tons.As Figure 6 shows, in the whole process of coal mining, carbon emissions are mainly concentrated in the coal drop and coal transportation links, accounting for 29% and 38% of the total carbon emissions, respectively, the other links account for 33%.Relevant administration departments can start from these two links to reduce greenhouse gas emissions.The most effective way is to control the use of electricity by making secondary use of coal gas or methane gas generated in the coal mines and using it for power generation through technological innovation.The main source of carbon emissions in the Fifth-II mining areas is from equipment use according to the calculation results.Shearers, emulsification pumps, and scraper conveyors are the main high carbon emission equipment.If carbon emissions are reduced from the perspective of the equipment, these three kinds of equipment can be used to formulate relevant indicators to regulate the equipment.The main source of carbon emissions in the Fifth-II mining areas is from equipment use according to the calculation results.Shearers, emulsification pumps, and scraper conveyors are the main high carbon emission equipment.If carbon emissions are reduced from the perspective of the equipment, these three kinds of equipment can be used to formulate relevant indicators to regulate the equipment.

Model Feasibility Verification
In this paper, two common carbon emission accounting methods are used to validate the B-R model based on relevant data from the Fifth-II mining areas: the IPCC method and the China Coal Production Enterprises Greenhouse Gas Emissions Accounting Methodology and Reporting Guide.
The IPCC method is an internationally authoritative carbon accounting model.The basic idea of the IPCC method is shown in Figure 7.According to the relevant content of the "IPCC Energy Volume", the main sources of carbon emissions from underground coal mining are: energy consumption, dynamite, spillover, and spontaneous combustion.The IPCC model is as follows: where E ipcc refers to CO 2 emissions calculated by the IPCC model (t); co t 2j refers to total CO 2 emissions of fuel j in time t; 44  12 refers to the molecular weight of carbon dioxide; EF t j refers to the energy consumption of fuel j (t); OF j refers to the CO 2 emission factor; and O j refers to the carbon oxidation rate of j.

E
= E fuel + E flare + E fugitive • GWP fugitive + E electricity + E heat (24) where E refers to the total GHG emissions of an enterprise (t); E fuel refers to the carbon emissions from fuel combustion (t); E flare refers to carbon emissions from flare combustion (t); E fugitive refers to carbon emissions from GHG fugitive (t); GWP fugitive refers to the Global Warming Potential of GHG; E electricity refers to carbon emissions from net purchased electricity (t); and E heat refers to carbon emissions from net purchased heat (t).E GHG = E fuel + E flare + E fugitive •GWP fugitive + E electricity + E heat (24) where E GHG refers to the total GHG emissions of an enterprise (t); E fuel refers to the carbon emissions from fuel combustion (t); E flare refers to carbon emissions from flare combustion (t); E fugitive refers to carbon emissions from GHG fugitive (t); GWP fugitive refers to the Global Warming Potential of GHG; E electricity refers to carbon emissions from net purchased electricity (t); and E heat refers to carbon emissions from net purchased heat (t).
where E refers to the total GHG emissions of an enterprise (t); E fuel refers to the carbon emissions from fuel combustion (t); E flare refers to carbon emissions from flare combustion (t); E fugitive refers to carbon emissions from GHG fugitive (t); GWP fugitive refers to the Global Warming Potential of GHG; E electricity refers to carbon emissions from net purchased electricity (t); and E heat refers to carbon emissions from net purchased heat (t).In this paper, the carbon emissions of the Fifth-II mining areas were calculated according to the IPCC method and Method G, and the specific results are shown in Table 7.In this paper, the carbon emissions of the Fifth-II mining areas were calculated according to the IPCC method and Method G, and the specific results are shown in Table 7.The IPCC method calculates that the total carbon emission of the Fifth-II mining areas is about 41,675 tons.The usual international error range is ±5%.A difference of 5 tons is within a reasonable range from the 41,670 tons calculated by the B-R model established in this paper.Method G calculated that the total carbon emissions of the Fifth-II mining areas are about 40,560 tons.The difference of 1110 tons is within a reasonable range from the other method's finding of 41,670 tons.It can be concluded that the B-R model established in this paper is feasible.

Conclusions
In this paper, the carbon emission boundary of fully mechanized mining was determined based on the idea of the life cycle, and a carbon emission calculation model based on each link of fully mechanized mining is proposed.For the fully mechanized mining process in the area, the following conclusions were drawn: (1) Among the four links of fully mechanized mining, the two links of coal transportation and coal drop contribute to a relatively high proportion of carbon emissions.The carbon emission sources were analyzed, and the main carbon emission source was equipment.(2) The IPCC method and The Methodology and Reporting Guide for Greenhouse Gas Emissions of Chinese Coal Producers were introduced to verify the B-R model, and we found that the error of the calculation results was small and the B-R model is feasible.(3) In order to better reduce carbon emissions, coal mines should control the use of electricity and improve the ability to transform and upgrade equipment.The secondary utilization of coal gas or methane gas generated inside the coal mine for power generation through technological innovation would be useful.

Further Work
(1) The carbon emission boundary in this paper assumes that, in the process of coal mining, there is no evaluation and calculation of related auxiliary systems.In the

Figure 1 .
Figure 1.The carbon dioxide emission boundary of coal mining.

Figure 1 .
Figure 1.The carbon dioxide emission boundary of coal mining.

16 Figure 2 .
Figure 2. Analysis of the sources of carbon emissions from roof support.

Figure 2 .
Figure 2. Analysis of the sources of carbon emissions from roof support.

Figure 3 .
Figure 3. Analysis of the sources of carbon emissions from the coal mine gob area.

Figure 3 .
Figure 3. Analysis of the sources of carbon emissions from the coal mine gob area.

Figure 4 .
Figure 4. Process of paste filling technology E of the B-R model in paste filling technology is expressed as:

Figure 4 .
Figure 4. Process of paste filling technology.

Figure 7 .
Figure 7. Flow chart of the IPCC method.The specific ideas of the "Greenhouse Gas Emissions Accounting Methodology and Reporting Guide for China's Coal Producers" (referred to as Method G) are shown in Figure 8.The Methodology and Reporting Guide for Greenhouse Gas Emissions of Chinese Coal Producers was commissioned by the China Development and Reform Commission and prepared by the National Center for Strategic Research and International Cooperation to address climate change.The compilation team drew on the research results and practical experience of GHG accounting reports of relevant enterprises at home and abroad and completed the compilation after field research and in-depth study.The Methodology and Reporting Guide on GHG Emissions of Chinese Coal Producers considers that there are five main sources of carbon emissions in the coal mining process: fuel combustion, flare combustion, GHG fugitive, net purchased electricity, and net purchased heat.

Figure 7 .
Figure 7. Flow chart of the IPCC method.The specific ideas of the "Greenhouse Gas Emissions Accounting Methodology and Reporting Guide for China's Coal Producers" (referred to as Method G) are shown in Figure 8.The Methodology and Reporting Guide for Greenhouse Gas Emissions of Chinese Coal Producers was commissioned by the China Development and Reform Commission and prepared by the National Center for Strategic Research and International Cooperation to address climate change.The compilation team drew on the research results and practical experience of GHG accounting reports of relevant enterprises at home and abroad and completed the compilation after field research and in-depth study.The Methodology and Reporting Guide on GHG Emissions of Chinese Coal Producers considers that there are five main sources of carbon emissions in the coal mining process: fuel combustion, flare combustion, GHG fugitive, net purchased electricity, and net purchased heat.

Figure 8 .
Figure 8. Flow chart of Method G.

Figure 8 .
Figure 8. Flow chart of Method G.

Table 1 .
Carbon emission source analysis of coal mines.

Table 2 .
Analysis of sources of carbon emissions from falling and mining coal.

Table 3 .
Analysis of the sources of carbon emissions from transporting coal.

Table 4 .
The coordinates of the inflection point in the fifth-II mining area.

Table 4 .
The coordinates of the inflection point in the fifth-II mining area.

Table 5 .
Summary of carbon emissions from major carbon emission sources in the mining areas.

Table 6 .
Carbon emissions in all stages of fully mechanized mining.

Table 7 .
Results of IPCC, Method G, and the B-R model in the Fifth-II mining area.