Prediction of the Adaptability of Using Continuous Extraction and Continuous Backfill Mining Method to Sequestrate CO2-A Case Study

The consumption of coal resources has caused an increase in CO2 emissions. A scientific concept that can realize CO2 sequestration, the harmless treatment of solid wastes, and coal extraction under buildings, railways, and water bodies (BRW) is proposed. First, a novel CO2 mineralized filling body (CMFB) is developed by employing CO2 gas, fly ash, silicate additives, and cement. It is then injected into the mined-out mining roadways (MRs) of the continuous extracting and continuous backfill (CECB) mining method to ameliorate the overburden migration and thus extract the coal body under the BRW. The AHP-fuzzy comprehensive evaluation method was employed to construct a prediction model for the suitability of this concept. Subsequently, the evaluation model is generalized and applied to the Yu-Shen mining area. Each indicator affecting adaptability is plotted on a thematic map, and the corresponding membership degree is determined. The aptness for 400 boreholes distributed in the entire area was determined and a zoning map which divides the whole area into good, moderate, slightly poor, and extremely poor suitability was drawn. This paper puts forward a mathematical model for predicting the suitability of using CECB and CMFB to sequestrate CO2. Research results can provide references for determining the site of CO2 sequestration under the premise of maximizing the economic and ecological benefits, which is conducive to constructing ecological, green, and sustainable coal mines.


Introduction
Coal consumption accounts for 70% of primary energy consumption, making significant contributions to China's economic development [1][2][3]. However, large-scale and high-intensity coal extraction has resulted in overburden breakage, water-conducting fracture development, and surface subsidence. If water diversion channels between the underground aquifer and the mined-out area occur, large amounts of water will flow into the gob, leading to the water inrush or even the sand inrush [4][5][6][7][8]. In addition, the original balance state of recharge, runoff, and discharge of the aquifer will be broken, and a drastic decline in water level and a large range of water resource loss will be inevitable [9,10]. Furthermore, severe mining disturbance causes ground subsidence, contributing to the damage of the surface buildings and structures. The relocation of villages is therefore unavoidable [11][12][13][14].
On the other hand, the utilization of coal resources has produced large amounts of CO 2 gas [15,16], indicating the importance of CO 2 sequestration. In-depth research on predicting the adaptability of CO 2 sequestration has been carried out at home and abroad.

CO 2 Sequestration Using CECB Mining Method
Prior to CECB mining, the mining panel is divided into many MRs, and then the MRs are allocated into different mining phases (generally [3][4][5]. The three-dimensional sketch map of the MRs arrangement of the CECB method is shown in Figure 1. The schematic map of the processes of this method is shown in Figure 2. With due consideration of the ventilation and excavation as well as the roof subsidence caused by CECB mining, the length of MR is generally less than 150 m and its width usually ranges from 4.0 to 7.0 m. To facilitate the turning of the continuous shearer and the hydraulic drilling rig, the angle α between the MR and the headgate ranges from 40 to 60 • . Due to the restriction on the height of the mining equipment, such as continuous shearers, etc., the height of MR is generally lower than 5.5 m. This method avoids the problem of uncoordinated operations between mining and filling during the longwall filling mining method, which is not conducive to conducting the filling process of CMFB. To effectively control the movement of the overlying strata and avoid surface subsidence, each MR is sealed and backfilled with CMFB immediately after the coal body is extracted. The MRs in the second phase will be extracted and backfilled after all the MRs in the first phase are mined and backfilled. Finally, the coal body in the entire mining panel is completely substituted with the CMFB, and the CO 2 gas can be sequestered in the form of solid CMFB. The replacement of carbon and coal can be realized by using the CECB mining method and CMFB.

Preparation of CMFB Specimens
In this paper, the proportions of low-calcium fly ash and the curing time are selected as the variables while the silicate additives, as well as the flow rate and concentration of CO2, gas are invariant. In this context, the control variable method is employed to design the testing schemes. The strength characteristics of CMFB samples are subject to the watersolid ratio. The ratio of the water weight to the total CMFB of 30% is invariable. In other words, the water-solid ratio is set to a constant of 3:7. The schemes consist of four groups, i.e., FA 55, FA 65, FA 75, and FA 85. The symbols FA 55, FA 65, FA 75, and FA 85 refer to the ratio of fly ash to the total amount of fly ash and cement is 55%, 65%, 75%, and 85%, respectively. For mechanic testing of UCS, the setting times are set to 3 d, 7 d, 14 d, 28 d, and 56 d, respectively. Three specimens are prepared for each curing time and each fly ash content. The number of the specimens for compressive testing is therefore at least 4 × 5 × 3 = 60. The testing scheme is designed and listed in Table 1.

Filling bodies in Phase Ⅱ
Filling bodies in Phase Ⅲ Filling bodies in Phase Ⅳ The sealing wall of MRs

Preparation of CMFB Specimens
In this paper, the proportions of low-calcium fly ash and the curing time are selected as the variables while the silicate additives, as well as the flow rate and concentration of CO2, gas are invariant. In this context, the control variable method is employed to design the testing schemes. The strength characteristics of CMFB samples are subject to the watersolid ratio. The ratio of the water weight to the total CMFB of 30% is invariable. In other words, the water-solid ratio is set to a constant of 3:7. The schemes consist of four groups, i.e., FA 55, FA 65, FA 75, and FA 85. The symbols FA 55, FA 65, FA 75, and FA 85 refer to the ratio of fly ash to the total amount of fly ash and cement is 55%, 65%, 75%, and 85%, respectively. For mechanic testing of UCS, the setting times are set to 3 d, 7 d, 14 d, 28 d, and 56 d, respectively. Three specimens are prepared for each curing time and each fly ash content. The number of the specimens for compressive testing is therefore at least 4 × 5 × 3 = 60. The testing scheme is designed and listed in Table 1.

Filling bodies in Phase Ⅱ
Filling bodies in Phase Ⅲ Filling bodies in Phase Ⅳ The sealing wall of MRs

Preparation of CMFB Specimens
In this paper, the proportions of low-calcium fly ash and the curing time are selected as the variables while the silicate additives, as well as the flow rate and concentration of CO 2 , gas are invariant. In this context, the control variable method is employed to design the testing schemes. The strength characteristics of CMFB samples are subject to the water-solid ratio. The ratio of the water weight to the total CMFB of 30% is invariable. In other words, the water-solid ratio is set to a constant of 3:7. The schemes consist of four groups, i.e., FA 55, FA 65, FA 75, and FA 85. The symbols FA 55, FA 65, FA 75, and FA 85 refer to the ratio of fly ash to the total amount of fly ash and cement is 55%, 65%, 75%, and 85%, respectively. For mechanic testing of UCS, the setting times are set to 3 d, 7 d, 14 d, 28 d, and 56 d, respectively. Three specimens are prepared for each curing time and each fly ash content. The number of the specimens for compressive testing is therefore at least 4 × 5 × 3 = 60. The testing scheme is designed and listed in Table 1. Both the development of CMFB specimens and the mechanical test are carried out under the regulation of the national standard (GBT17671-1999). The detailed process of CMFB sample preparation and mechanical experiment is shown in Figure 3.  Both the development of CMFB specimens and the mechanical test are carried out under the regulation of the national standard (GBT17671-1999). The detailed process of CMFB sample preparation and mechanical experiment is shown in Figure 3. First, weigh the fly ash and cement and mix them evenly. Meanwhile, weigh and put the silicate additive into a container and keep stirring until all the particles are dissolved. The silicate additive was dissolved in tap water by stirring at 600 rpm for 15 min at room temperature, as shown in Figure 3a-c. Subsequently, the dissolved solution is First, weigh the fly ash and cement and mix them evenly. Meanwhile, weigh and put the silicate additive into a container and keep stirring until all the particles are dissolved. The silicate additive was dissolved in tap water by stirring at 600 rpm for 15 min at room temperature, as shown in Figure 3a-c. Subsequently, the dissolved solution is immediately poured into the homogeneous mixture of fly ash and cement, and continued stirring until a uniform paste is obtained. The paste was then poured into the reactor, and the motor started and stirred at 500 rpm while CO 2 gas was introduced at a constant rate of 1 L/min under atmospheric conditions. The entire stirring time lasted 20 min, and then the fresh slurry of CMFB is achieved, as shown in Figure 3d-f. Then, the filling slurry was decanted into standard moulds with a size of ϕ×h of 50 mm × 100 mm for UCS experiments. Twelve to fifteen hours later, the samples are taken out and put into the standard curing box. Its humidity is adjusted to approximately 98% and the temperature remains at 20 • C. Samples of CMFB are cured for 3, 7, 14, 28, and 56 days, respectively. Afterwards, the specimens are taken out of the curing box in batches. The dimensions are calibrated to make the specimens standard ones for the UCS test. After that, the UCS experiments were implemented, as shown in Figure 3g-o.
The diagrammatic drawing of the chemical reaction of CO 2 gas and fresh slurry is shown in Figure 4. immediately poured into the homogeneous mixture of fly ash and cement, and continued stirring until a uniform paste is obtained. The paste was then poured into the reactor, and the motor started and stirred at 500 rpm while CO2 gas was introduced at a constant rate of 1 L/min under atmospheric conditions. The entire stirring time lasted 20 min, and then the fresh slurry of CMFB is achieved, as shown in Figure 3d-f. Then, the filling slurry was decanted into standard moulds with a size of φ×h of 50 mm × 100 mm for UCS experiments. Twelve to fifteen hours later, the samples are taken out and put into the standard curing box. Its humidity is adjusted to approximately 98% and the temperature remains at 20 °C. Samples of CMFB are cured for 3, 7, 14, 28, and 56 days, respectively. Afterwards, the specimens are taken out of the curing box in batches. The dimensions are calibrated to make the specimens standard ones for the UCS test. After that, the UCS experiments were implemented, as shown in Figure 3g-o. The diagrammatic drawing of the chemical reaction of CO2 gas and fresh slurry is shown in Figure 4.

Testing Results of UCS and Elastic Modulus of CMFB
The changing law of UCS of CMFB corresponding to various proportions of fly ash and curing times are shown in Figure 5a,b, respectively. Scholars working on the development of paste filling bodies generally believe that increasing the proportion of fly ash will result in low backfill strength. The low hydration reaction and the fewer hydration products caused by the low fly ash content may account Curing time (d) Fly ash content：55% Fly ash content：65% Fly ash content：75% Fly ash content：85%

Testing Results of UCS and Elastic Modulus of CMFB
The changing law of UCS of CMFB corresponding to various proportions of fly ash and curing times are shown in Figure 5a,b, respectively. immediately poured into the homogeneous mixture of fly ash and cement, and continued stirring until a uniform paste is obtained. The paste was then poured into the reactor, and the motor started and stirred at 500 rpm while CO2 gas was introduced at a constant rate of 1 L/min under atmospheric conditions. The entire stirring time lasted 20 min, and then the fresh slurry of CMFB is achieved, as shown in Figure 3d-f. Then, the filling slurry was decanted into standard moulds with a size of φ×h of 50 mm × 100 mm for UCS experiments. Twelve to fifteen hours later, the samples are taken out and put into the standard curing box. Its humidity is adjusted to approximately 98% and the temperature remains at 20 °C. Samples of CMFB are cured for 3, 7, 14, 28, and 56 days, respectively. Afterwards, the specimens are taken out of the curing box in batches. The dimensions are calibrated to make the specimens standard ones for the UCS test. After that, the UCS experiments were implemented, as shown in Figure 3g-o. The diagrammatic drawing of the chemical reaction of CO2 gas and fresh slurry is shown in Figure 4.

Testing Results of UCS and Elastic Modulus of CMFB
The changing law of UCS of CMFB corresponding to various proportions of fly ash and curing times are shown in Figure 5a,b, respectively. Scholars working on the development of paste filling bodies generally believe that increasing the proportion of fly ash will result in low backfill strength. The low hydration reaction and the fewer hydration products caused by the low fly ash content may account Curing time (d) Fly ash content：55% Fly ash content：65% Fly ash content：75% Fly ash content：85% reaction and the fewer hydration products caused by the low fly ash content may account for this phenomenon [27]. In this paper, it is interesting that the UCS rises first and then drops with the increasing fly ash content, reaching the peak at 75% fly ash content. This trend is almost the opposite to that of cemented backfill without CO 2 . Note that, in the early stage of the test, we have developed a sample of CMFB with 95% fly ash. However, the UCS of 95% CMFB is too small to bear the preloading, and the creep phenomenon occurs during the UCS test.
The UCS of CMFB with 75% fly ash is 1.982, 2.356, 3.832, 4.162, and 7.499 MPa, respectively, when the curing time is 3, 7, 14, 28, and 56 d, respectively. The UCS of CMFB with 85% fly ash content is lower than that of CMFB whose fly ash content is 75% and 65% when the setting time is invariable. Therefore, the influence of the addition of CO 2 gas and silicate additives on the strength of the CMFB is significant.
According to the results for UCS of CMFB from the laboratory, the elastic modulus is determined by calculating the slope of the linear elastic stage of the stress-strain curve: where E f is the elastic modulus of CMFB; σ a and σ b refers to the stress of the starting and ending point of the linear elastic stage, respectively; ε a and ε b represent the strain of the origin and the end of the linear elastic stage, respectively. The results of the Young's modulus calculation of CMFB with a curing time of 56 d are shown in Figure 6. for this phenomenon [27]. In this paper, it is interesting that the UCS rises first and then drops with the increasing fly ash content, reaching the peak at 75% fly ash content. This trend is almost the opposite to that of cemented backfill without CO2. Note that, in the early stage of the test, we have developed a sample of CMFB with 95% fly ash. However, the UCS of 95% CMFB is too small to bear the preloading, and the creep phenomenon occurs during the UCS test.
The UCS of CMFB with 75% fly ash is 1.982, 2.356, 3.832, 4.162, and 7.499 MPa, respectively, when the curing time is 3, 7, 14, 28, and 56 d, respectively. The UCS of CMFB with 85% fly ash content is lower than that of CMFB whose fly ash content is 75% and 65% when the setting time is invariable. Therefore, the influence of the addition of CO2 gas and silicate additives on the strength of the CMFB is significant.
According to the results for UCS of CMFB from the laboratory, the elastic modulus is determined by calculating the slope of the linear elastic stage of the stress-strain curve: where Ef is the elastic modulus of CMFB; a σ and b σ refers to the stress of the starting and ending point of the linear elastic stage, respectively; a ε and b ε represent the strain of the origin and the end of the linear elastic stage, respectively. The results of the Young's modulus calculation of CMFB with a curing time of 56 d are shown in Figure 6.

Determination of Indicators and Membership Functions
The indicators affecting the suitability of using CECB to sequestrate CO2 are complex and fuzzy, reaching good agreement with the characteristics and connotations of fuzzy mathematics. The combination of AHP with the fuzzy comprehensive evaluation method is capable of quantifying the factors and making them more concrete and accurate.
The AHP-fuzzy comprehensive evaluation method was employed to construct a triple-leveled structural model and obtain the distribution of the weights. The mathematical model includes six factors as secondary indicators, i.e., roof and floor system, BRW system, geological structures system, mining block system, coal seam system, and filling body system. Moreover, nineteen influencing factors, including the strength of the immediate roof, the distance from the overlying aquifer to the coal seam, the fold complexity

Determination of Indicators and Membership Functions
The indicators affecting the suitability of using CECB to sequestrate CO 2 are complex and fuzzy, reaching good agreement with the characteristics and connotations of fuzzy mathematics. The combination of AHP with the fuzzy comprehensive evaluation method is capable of quantifying the factors and making them more concrete and accurate.
The AHP-fuzzy comprehensive evaluation method was employed to construct a tripleleveled structural model and obtain the distribution of the weights. The mathematical model includes six factors as secondary indicators, i.e., roof and floor system, BRW system, geological structures system, mining block system, coal seam system, and filling body system. Moreover, nineteen influencing factors, including the strength of the immediate roof, the distance from the overlying aquifer to the coal seam, the fold complexity coefficient, the coal reserves of the mining panel, the seam thickness, and the filling body strength, etc., were selected as tertiary indicators, as shown in Figure 7. Furthermore, the membership function of each indicator was analyzed and obtained based on the field practice of the CECB mining method, the UCS test of CMFB, and the summary of the current literature. coefficient, the coal reserves of the mining panel, the seam thickness, and the filling body strength, etc., were selected as tertiary indicators, as shown in Figure 7. Furthermore, the membership function of each indicator was analyzed and obtained based on the field practice of the CECB mining method, the UCS test of CMFB, and the summary of the current literature. Figure 7. The hierarchy structure for evaluating the suitability of using CECB to sequestrate CO2.

The Strength of the Basic Roof
During the process of using CECB to sequestrate CO2, the strength of the basic roof should be considered since the soft roof may cause large amounts of roof subsidence despite the temporary support of the MR roof being adopted. Meanwhile, each mining equipment has a requirement for the height of MR. Therefore, the height of MR may be lower than that of the continuous miner and make it difficult to pass through, which will lower the adaptability of the concept. Therefore, the method of CO2 sequestration employing CECB can be carried out better in the coal mines with a hard roof. The firmness coefficient of the basic roof was selected as the evaluation standard.
The membership function of the firmness coefficient of the basic roof is as follows: Figure 7. The hierarchy structure for evaluating the suitability of using CECB to sequestrate CO 2 .

The Strength of the Basic Roof
During the process of using CECB to sequestrate CO 2 , the strength of the basic roof should be considered since the soft roof may cause large amounts of roof subsidence despite the temporary support of the MR roof being adopted. Meanwhile, each mining equipment has a requirement for the height of MR. Therefore, the height of MR may be lower than that of the continuous miner and make it difficult to pass through, which will lower the adaptability of the concept. Therefore, the method of CO 2 sequestration employing CECB can be carried out better in the coal mines with a hard roof. The firmness coefficient of the basic roof was selected as the evaluation standard.
The membership function of the firmness coefficient of the basic roof is as follows:

The Thickness of the Basic Roof
The basic roof controls the overburden migration and has an indispensable influence on the realization of using CECB to sequestrate CO 2 . The ratio of the thickness of the basic roof to that of the coal seam was selected as the evaluation standard.
where n is the number of boreholes in the mining block; H ri is the thickness of the basic roof of the i-th borehole, m; H ci is the thickness of the coal seam of the i-th borehole, m.

The Strength of the Immediate Floor
The strength of the immediate floor has a great influence on equipment transportation and worker passage. The argillaceous soft floor suffers from argillization, disintegration, expansion, and fragmentation soaked by water. The CMFB may mix with the disintegrated floor, contributing to a decline in strength. The transportation of mining equipment may be seriously affected, and the production efficiency will be reduced. Consequently, countermeasures such as immediate floor grouting reinforcement should be taken to make it feasible for CO 2 sequestration in the CECB mining face. The firmness coefficient of the immediate floor was selected as the index.

The Required Protection Grade of Surface Buildings
Mining-caused surface subsidence usually exerts a significant influence on BRW. There are five surface deformation indexes, i.e., vertical and horizontal displacement, inclination, curvature, and horizontal deformation. Different types of buildings have various abilities to anti-surface deformation, and when it exceeds the critical value, the buildings will be damaged. The corresponding damage levels and standards of the buildings regulated by state are listed in Table 2. Table 2. State classification criteria of the buildings damage grade [12,49].

Grade
Horizontal Deformation (mm/m) Curvature (mm/m 2 ) Tilt (mm/m) With due consideration of the economic and ecological benefits, the CO 2 sequestration utilizing the CECB method is more suitable for extracting coal resources under buildings whose protection class is higher than grade II. For the coal seams trapped under the buildings, which need grade III protection class, partial mining, i.e., strip mining, room and pillar mining, and partial backfill mining methods, are the most suitable ones that can be employed to extract coal resources under BRW due to their low production cost. Additionally, for buildings that need a grade IV protection class, the coordinated mining and slice mining methods can meet the requirements. The CO 2 sequestration using CECB can also be realized when the protective grade is III or IV, while the costs may be higher.

The Distance from Overlying Aquifer to Coal Seam
The height of the water-conducting fractured zone of the CECB method depends on the mining height, the CMFB strength, the overlying strata lithology, etc. During the process of generalizing and applying the CECB, it was found that after all MRs were extracted and backfilled, the height of the water-flowing fractured zone is approximately 5 times the mining height. Moreover, to prevent the overlying water bodies from penetrating into the mined-out area and thus resulting in a water inrush, the thickness of the protective zone between the top boundary of the fractured zone and the overlying aquifer was analyzed. The strength of various types of rock was employed to quantize the water-resisting indicator of the protection zone. The impermeability decreases with the increasing rock strength of the protection zone. Generally, a clay stratum with a thickness of more than 12 m can effectively preserve the water of the overburdened aquifer from percolating into the goaf [12]. Hence, if the water-resistant indicator of clay is set to 1, the values of the waterresisting index of other strata with various lithologies range from 0 to 1. Then the formula for the protective zone thickness can be expressed as follows: where H p is the thickness of the protective zone (m); H c is the thickness of the clay stratum (m); k is the water-resistant indicator of the protective zone. If the distance is greater than the sum of the height of the water-flowing fractured zone and the thickness of the protective zone, the water-resisting layers can totally block the hydraulic interaction between the aquifer and the mined-out area. If the distance is only greater than the height of the water-conducting fracture zone, then the groundwater only permeates slowly in the overlying layers between the aquifer and the mined-out areas, and no water inrush accident will occur. However, when it is less than the height of the water-conducting fracture zone, the water-flowing fractures will reach the aquifer, contributing to a large-scale loss in groundwater resources and a decline in the water table. In this context, the concept is totally incapable of achieving its purpose of CO 2 sequestration and water preservation under BRW.
where d denotes the distance between the overlying aquifer and the coal seam, m; H f denotes the height of the water conductive fractured zone, m; H p denotes the thickness of the protective zone.

Water Yield Property of the Overlying Aquifer
The water yield of the aquifer is closely related to the thickness and the permeability coefficient of the aquifer. The larger the thickness and permeability, the higher its water storage and recharge capacity, and thus the higher the water-richness. With due consideration of the mining costs, the high water-richness means the high adaptability of CO 2 sequestration using the CECB method since it is capable of controlling overburden migration and maintaining the stability and integrity of the aquifuge beneath the aquifer using the developed CMFB. According to the classification standard of aquifer water yield, the water abundance can be divided into three levels according to the unit water inflow of the borehole, as shown in Table 3. Table 3. Classification standard of water abundance of overlying aquifer [12].

Grade of Water Yield Property Uniform Drawdown Unit Flow (L·s
The faults destroy the continuity and integrity of the overburden. The overlying strata near the medium and small faults are generally broken, which is prone to causing support break-off during the process of CO 2 sequestration by the CECB. The large faults will disturb the original mining arrangement of MRs and thus lower the mining efficiency. The coal seam in the vicinity of the large faults is usually left unmined as the permanent protective coal pillar, leading to a low recovery rate of coal resources and a low storage rate of CO 2 .
The fault strength was selected as the evaluation indicator affecting the appropriateness of the concept. The formula of fault strength is as follows: where l i refers to the length of the i-th fault (m); h i refers to the fault throw of the i-th fault (m); A refers to the area of the block (m 2 ); k refers to the number of faults in the block. The membership function of fault strength is:

The Characteristic of Fold
There is high residual stress and elastic energy in the folded rock mass. When sequestrating CO 2 in the MRs of the CECB, the elastic energy accumulated in the coal and rock mass is further released, which may lead to serious deformation of the MRs and even destructive rock burst of the surrounding rock. The development and complexity of the fold have an influence on the mining layout and mining sequence of the method. The fold complexity coefficient P was selected as the indicator.
where ∆h is the contour height difference between the lowest and highest coal seam in the CECB mining panel, m; w is the change value of the trend of the contour line expressed by radian, rad; l is the horizontal contour distance between the lowest and highest coal seam in the CECB mining panel, m; s is the acreage of the CECB mining block, km 2 .

The Strength Index of the Collapsed Pole
The collapse column destroys the continuity and integrity of the coal seam, which is to the disadvantage of the layout and extraction of MRs. The strength index of the collapse column k c was chosen as the evaluation standard.
where S i is the acreage of the affected area of the i-th collapse column in the CECB mining block, km 2 ; S is the acreage of the CECB mining block, km 2 ; η i is the ratio of the destruction thickness of the i-th collapse column to the mining thickness. The deployment and layout of the CECB method is flexible, which is suitable for extracting irregular coal seams and sequestrating CO 2 . However, for the extremely irregular mining block, with the increase in irregularity, it is difficult to optimize the layout of the MRs and the mining sequence, and thus mining efficiency and CO 2 sequestration rate will be decreased. Taking the irregularity P ir of the CECB mining block as the evaluation index: where P ir is the degree of irregularity of the CECB mining block; m is the number of irregular zones in the CECB mining block, km 2 ; A i is the number of edges of the i-th irregular zone.

The Coal Reserves of the CECB Mining Block
It is necessary for the concept of using CECB to sequestrate CO 2 to construct filling stations, collect CO 2 gas and set up the filling pipelines to transport CMFB. The costs for recovering coal and sequestrating CO 2 are higher than those of other general backfill mining methods. Subsequently, the concept is better employed to extract large amounts of coal resources under BRW, or large pieces of boundary coal pillar, and it is inappropriate to extract a small amount of barrier pillar between panels. According to the former field practice of the CECB mining method, the maximum daily output can reach 3000 t, and the annual output is more than 1000,000 t.
where C r is coal reserves of the CECB mining block, ten thousand tons. The CECB mining method is appropriate for shallowly buried coal seams since the deeper depth means greater in-site stress of the mining stope. Under the action of the high static pressure, the MRs will suffer from large deformations and the roof support and equipment passage will be difficult, thereby the suitability grade will decrease. The membership function of the buried depth of the coal seam b is as follows:

The Thickness of Coal Seam
The thickness of the coal seam, varying from 2 to 3.5 m, is the most suitable one for CO 2 sequestration using the CECB mining method. When the thickness is less than 1.5 m, the equipment transportation and operations are limited, which greatly affects the working efficiency and thus lowers its suitability. Furthermore, coal seams, with a thickness exceeding 4.5 m are inappropriate for the method since they exceed the mining height of the most commonly used continuous miners. If the CECB mining method is applied to an extremely thick coal seam, part of the top coal must be discarded, contributing to the low filling volume and CO 2 sequestration rate.

The Variation Coefficient of the Thickness of Coal Seam
The variation coefficient is a quantitative index to appraise the deviation degree of the thickness of a coal seam from the average thickness. The higher the variation coefficient, the more unfavorable to the application of CO 2 sequestration using the CECB method. The calculation formula for the variation coefficient is as follows: where γ c refers to the variation coefficient of the thickness of the coal seam; S refers to the standard deviation of the variation in coal thickness.
where n is the number of boreholes; x i is the measured coal thickness of the i-th borehole, m; x is the average thickness of the coal seam in the CECB mining panel, m. Mining panels, with a variation coefficient of less than 30%, are completely suitable for the application of the concept since the narrow change of coal seam hosting. For mining panels whose coefficients range between 30% and 50%, most of the coal seams in the mining block can be mined and the rest can be extracted by taking countermeasures. When the coefficient varies from 50% to 70%, extracting coal and storing CO 2 using the CECB method is neither economically desirable nor technically feasible. If the coefficients are greater than 70%, the coal seams are usually chicken coop shaped, lens shaped or lentil shaped, and are discontinuous and scattered. The coal resources are totally unrecoverable using the CECB method and the adaptability grade is therefore equal to 0.

The Dip Angle of the Coal Seam
Generally, the continuous miner and crawler-type support in CECB mining are selfmoving equipment, which is appropriate for the coal seam with a small dip angle. If the dip exceeds 10 • , it is difficult for the equipment to move automatically. Additionally, the filling slurry of CMFB tends to flow down to the bottom of the MRs, which goes against the CMFB connecting with the roof. Therefore, the adaptability degree equals 0. The CECB mining face is suitable to be arranged in the flat coal seam whose dip angle is less than 8 • , especially in the ones with a dip ranging from 1 • to 3 • . Coal seams, with a dip angle of less than 1 • , are not conducive to the self-flow of the fresh slurry of CMFB and thus have a negative influence on the filling process. The average dip angle in the mining block is taken as the evaluation index.
where n is the number of boreholes; α i is the dip angel of the i-th borehole, m.

The Coefficient of Dirt Band
The MRs in the CECB mining face are generally suitable to be arranged in the coal seam with no or less dirt band. The pure coal seam can contribute to the CECB method, giving full play to its advantages of fast extraction and fast CO 2 sequestration. On the contrary, the higher the proportion of dirt band, the greater the resistance of continuous miners to coal cutting, contributing to the concept's low mining efficiency and adaptability. The coefficient of dirt band is chosen as the indicator: where n is the number of boreholes; X i is the thickness of the coal seam of the i-th borehole, m; G i refers to the thickness of the dirt band of the i-th borehole, m.
4.5.6. The Strength of Coal Seam The influence of the strength of the coal seam on the appropriateness is manifested primarily in two aspects. On the one hand, a high-strength coal body will affect the cutting speed and mining efficiency of the continuous miner, leading to a low rate of CO 2 sequestration. In contrast, the low strength may lead to the instability and rib spalling of the isolated coal pillar since the CECB method adopts skip mining. Moreover, continuous shear, shuttle cars, crawler supports, and other equipment have relatively strict requirements on the dimension of MR. The soft coal seam is prone to causing great deformation of MR and its size cannot be guaranteed. The comprehensive strength of coal seam f c is used as an indicator.
where k g denotes the coefficient of the dirt band; R c denotes the UCS of the coal body, MPa; R e denotes the strength of the dirt band, MPa.
4.6. CMFB System 4.6.1. The UCS of the CMFB The CMFB substitutes the coal body to support the roof and restrains overburden movement and breakage. The UCS of CMFB plays a decisive role in mitigating the deformation of the overburden. It determines whether the CO 2 sequestration using CECB can successfully be implemented under BRW. At present, the filling materials are primary paste backfill, high water backfill, and gangue backfill. The development and application of CMFB have not been reported up to now. The UCS of CMFB is manipulated by material ratios, water-solid ratio, the setting time, etc. As Figure 5 shows, the UCS of CMFB with a curing time of 56 d remains stable and is therefore utilized as the evaluation index.

The Compression Rate of the CMFB
Diverse compression ratios have various effects on restricting the overburden's movement and deformation. The smaller the compression ratio, the better the control effect of overburden migration. The compression rate of the paste filling is roughly less than 4%. The compressibility of ultra-high water filling material is generally less than 0.3% due to 97% water content. In contrast, the compression rate of the gangue filling body, ranging from 10% to 40%, is much larger than that of the paste filling body and ultra-high water filling body owing to the large void between the gangue particles. It plays a minor role in reducing surface subsidence.
The Young's modulus of the CMFB is a significant index reflecting its supporting capacity for the overlying strata. The CMFB with a large elastic modulus is capable of effectively alleviating the migration and fracture development of the overburden and thus ameliorating the mining-induced impact on the underground aquifer and surface buildings. On the contrary, the CMFB with a small elastic modulus has a higher compression rate under the same load than that of the larger one, leading to the decreasing suitability of using CECB to sequestrate CO 2 . The compression rate can be calculated by: The membership function of the compression rate of CMFB is:

Mathematical Modeling
The AHP-fuzzy comprehensive evaluation method was employed to evaluate the suitability of using CECB to sequestrate CO 2 . Based on the analysis of the factors influencing the applicability, six sub-factors and nineteen triple-leveled factors were identified and determined. In the AHP model, A represents the overall goal, B denotes the evaluation standard, and C is the sub-standard. The judgment matrix was constructed using this hierarchy to realize the numerical expression and quantification of the appropriateness. The influencing factors were put into a discourse domain U, U = {u 1 , u 2 , u 3 , u 4 , u 5 , u 6 }, Let V denote the discourse domain of applicability. The grade of adaptability was represented by the comprehensive evaluation value, Φ, which is the membership degree of U in the fuzzy subset V. The fuzzy subset V was defined by Equation (34) [26]: The comprehensive evaluation value Φ can be derived via the following equation: where u i (u i ) is the membership degree of the i-th influencing factor and w i is the weight of the i-th indicator. The fuzzy subset V of the adaptability of CO 2 sequestration in the CECB mining face can be divided into four grades, as listed in Table 4. The applicability is therefore quantified and embodied by calculating the value Φ.
The method is not totally appropriate for the colliery due to the complex engineering and hydrogeological conditions. Several countermeasures are necessary, and fewer profits can be made.
Due to the extremely adverse conditions, the CECB is totally unsuitable for CO 2 sequestration no matter what countermeasures will be taken. Although the colliery can extract a part of the coal resources and store small amounts of CO 2 gas using the CECB mining method, the filling scale cannot meet the requirements. The investment is far greater than the economic and environmental benefits, and thus the profits cannot meet the disbursements.

Weight Distribution
The 1-9 and its reciprocal scaling method proposed by Thomas. L. Saaty was employed to determine the value of the judgment matrix elements and thus reflect the relative importance of each element. a ij denotes the relative importance of the i-th indicator to the jth indicator. Experts and scholars engaged in CECB mining, cemented backfill development, coal extraction under BRW, and CO 2 sequestration were invited to assign the weights of various factors in the AHP model to evaluate the applicability of using CECB to sequestrate CO 2 by developing CMFB. The following judgment matrixes were determined according to the relative weight of the experts: The roof and floor system whose matrix is W B 1 ∼C was taken as an example. Its maximum eigenvalue is 3.0652, and the corresponding eigenvector W = [0.2461, 0.4923, 0.2616]. Then Formulae (36) and (37) were used for consistency test.
where C.I. denotes the consistency indicator; λ max denotes the largest eigenvalue; n denotes the number of the factors in the matrix W B 1 ∼C .
C.R. = (C.I.)/(R.I.), (37) where C.R. is the consistency ratio and R.I. is the average consistency index. If C.R. ≥ 0.1, the relative weights in the matrix should be reassigned. Otherwise, the weight distribution is reasonable, and the calculation results can be taken as the final results. The C.R. of the matrix is 0.0627, which is much less than 0.1, suggesting the relative weights given by the experts are acceptable.
The weight of the roof and floor system is 0.1037. Therefore, the weight of the three third-level factors in the roof and floor system is 0.0255, 0.0511, and 0.0271, respectively.
The consistency tests of the other six matrices were conducted in the same way and the results are listed in Table 5. After finishing the consistency test, the weight distribution of all the influencing factors were obtained and are listed in Table 6.

Generalization and Application of the Prediction Model
The Yu-Shen coal area, covering a total area of 5160 km 2 , is situated in the north of Yulin City, Shaanxi Province, China. The total in situ coal reserve in this area is greater than 30 billion tons. In the early phase of coal mining, owing to the high-intensity and large-scale mining activities, the water level of the Salawusu Formation aquifer dropped drastically over a large area, and the number of springs and rivers declined sharply. This phenomenon seriously jeopardizes domestic and ecological water use for local residents and the environment. Additionally, the surface buildings, structures, and ecosystem were damaged to a large extent due to mining-caused surface subsidence, forcing the relocation of the villages in the vicinity of the many collieries in the coal area. Backfill mining is thus necessary for the coal area. On the other hand, large amounts of coal-based solid wastes such as fly ash and gangue are produced and accumulated in the process of coal extraction and coal dressing. Meanwhile, the coal burning and consumption in the thermal power plant in the Yu-Shen mining area generates a quantity of CO 2 gas to generate electricity. Therefore, the coal area is taken as an instance to predict the suitability of CO 2 sequestration using CECB and CMFB. It is noted that the prospective investigation is only a feasibility analysis for the field application of the concept in the future. The concept has not been employed in any colliery in the Yu-Shen mining area.
Prior to employing CECB to sequestrate CO 2 by means of developing CMFB in the mining area, the adaptability of this method for various collieries (400 boreholes) is predicted by using the proposed mathematical evaluation model. Subsequently, the thematic map of the contour map of Φ is plotted, to provide a decision basis for the planners and engineers. The primary calculating process of the comprehensive evaluation value Φ of the applicability of using CECB to sequestrate CO 2 is shown in Figure 8. First, plotting the thematic map of a single indicator affects the suitability of the concept. Then, determining the membership degree of each factor based on comprehensively analyzing the thematic map and summarizing the relevant literature. Third, based on the weight distribution of each influencing factor listed in Table 6 as well as the membership degree of each factor, calculate the value of Φ by using Equation (35). According to the borehole data, the basic roof is dominantly off-white and light gray middle-fine grained feldspar sandstone, debris-feldspar, and calcareous sandstone, followed by gray black siltstone. The UCS of the basic roof is more than 80 MPa. The ratio of the thickness of the basic roof to the first-mined coal seam varies greatly, with an average Figure 8. Process for evaluating and predicting the adaptability of using CECB to sequestrate CO 2 .

Roof and Floor Conditions in the Yu-Shen Mining Area
According to the borehole data, the basic roof is dominantly off-white and light gray middle-fine grained feldspar sandstone, debris-feldspar, and calcareous sandstone, followed by gray black siltstone. The UCS of the basic roof is more than 80 MPa. The ratio of the thickness of the basic roof to the first-mined coal seam varies greatly, with an average of eight. The immediate floor is primarily composed of sandy mudstone, mudstone, and carbonaceous mudstone, with lenticular marl distributed sporadically. The UCS of the floor ranges from 24 to 66 MPa, with an average of 45 MPa.

Thematic Maps of BWR System
Due to the fragile ecological environment of the whole mining area, the protection level of surface buildings is regarded as grade I. Moreover, there are four types of overlying aquifers in the Yu-Shen mining area, namely, the confined porous bedrock aquifer, the burnt rock phreatic aquifer, and the Salawusu phreatic aquifer, and the unconsolidated porous phreatic aquifer, with increasing elevation. The Salawusu Formation aquifer, featuring a large thickness, shallowly buried depth, and wide distribution, is the main aquifer in the coal area which should be taken into account when predicting the adaptability. Based on the collected borehole data, the contour maps of the distance from the overlying aquifer to the first-mined coal seam and the contour map of the height of the water conductive fractured zone, were drawn by employing the Kriging interpolation method, as shown in Figure 9a,b respectively. The contour map of the aquifuge (clay and red soil), whose water-resisting index equals 1 is shown in Figure 9c. The bedrock of the entire Yu-Shen mining area is composed of sandstone, mudstone, and sandy mudstone, and the waterresisting index of the bedrock is 0.4. Using Equation (3), the thicknesses of the protection zone of 400 boreholes in the Yu-Shen mining area were obtained, and the contour map was plotted by using Surfer software, as shown in Figure 9d. Furthermore, the areas of high water abundance of the Sarawusu Formation aquifer cover a total area of 504 km 2 , which represents 10 percent of the entire Yu-Shen mining area. The total acreage of the moderate and low water-rich regions is 1911 km 2 and 1919 km 2 , respectively, accounting for approximately 37% of the mining area. The zoning map of the water-richness of the overlying aquifer is shown in Figure 10.
contour map was plotted by using Surfer software, as shown in Figure 9d. Furthermore, the areas of high water abundance of the Sarawusu Formation aquifer cover a total area of 504 km 2 , which represents 10 percent of the entire Yu-Shen mining area. The total acreage of the moderate and low water-rich regions is 1911 km 2 and 1919 km 2 , respectively, accounting for approximately 37% of the mining area. The zoning map of the water-richness of the overlying aquifer is shown in Figure 10.

Thematic Map of Geological Structures System
The geological exploitation results show that there are seven faults in the entire Yu-Shen coal area and the geological structure is simple. The fault distribution and the fault strength are shown in Figure 11. Generally, we tend to agree that there are no obvious folds or collapse column structures in the entire coal area. The shape of the CECB mining panel is thought of as quadrangular, so the irregular rate is 1. According to the preliminary colliery planning before production, the annual output of a working face is usually more than 400,000 tons.

Thematic Maps of Coal Seam System
The buried depth of the first-mined coal seam varies from 0 to 650 m, with a thickness ranging from 0 to 12 m and the dip angle varying from 1° to 3°. The contour maps of the buried depth and the thickness of the first-mined coal seam were plotted, respectively, as shown in Figure 12a,b. The coal seam thickness shows narrow variation in the CECB  Figure 10. Zoning map of the water-richness of the overlying Salawusu Formation aquifer.

Thematic Map of Geological Structures System
The geological exploitation results show that there are seven faults in the entire Yu-Shen coal area and the geological structure is simple. The fault distribution and the fault strength are shown in Figure 11. Generally, we tend to agree that there are no obvious folds or collapse column structures in the entire coal area.

Thematic Map of Geological Structures System
The geological exploitation results show that there are seven faults in the entire Yu-Shen coal area and the geological structure is simple. The fault distribution and the fault strength are shown in Figure 11. Generally, we tend to agree that there are no obvious folds or collapse column structures in the entire coal area. The shape of the CECB mining panel is thought of as quadrangular, so the irregular rate is 1. According to the preliminary colliery planning before production, the annual output of a working face is usually more than 400,000 tons.

Thematic Maps of Coal Seam System
The buried depth of the first-mined coal seam varies from 0 to 650 m, with a thickness ranging from 0 to 12 m and the dip angle varying from 1° to 3°. The contour maps of the buried depth and the thickness of the first-mined coal seam were plotted, respectively, as shown in Figure 12a,b. The coal seam thickness shows narrow variation in the CECB  The shape of the CECB mining panel is thought of as quadrangular, so the irregular rate is 1. According to the preliminary colliery planning before production, the annual output of a working face is usually more than 400,000 tons.

Thematic Maps of Coal Seam System
The buried depth of the first-mined coal seam varies from 0 to 650 m, with a thickness ranging from 0 to 12 m and the dip angle varying from 1 • to 3 • . The contour maps of the buried depth and the thickness of the first-mined coal seam were plotted, respectively, as shown in Figure 12a,b. The coal seam thickness shows narrow variation in the CECB mining panel, and the variation coefficient is, therefore, less than 30%. The first-mined coal seam is generally free of gangue. Only a small part of the coal seam contains 1 or 2 layers of gangue, indicating that the coefficient of the dirt band is less than 10%. The first-mined coal seam belongs to long-flame coal or non-caking coal, and the Proctor coefficient is from 2 to 2.5.
Minerals 2022, 12, x FOR PEER REVIEW 24 of 30 mining panel, and the variation coefficient is, therefore, less than 30%. The first-mined coal seam is generally free of gangue. Only a small part of the coal seam contains 1 or 2 layers of gangue, indicating that the coefficient of the dirt band is less than 10%. The first-mined coal seam belongs to long-flame coal or non-caking coal, and the Proctor coefficient is from 2 to 2.5.
(a) (b) Taking CMFB with a curing time of 56 d and a fly ash content of 75% as an example, its UCS is 7.499 MPa and its elastic modulus is 1.5 GPa. According to the coal seam burial depth of different boreholes, the compressibility of CMFB corresponding to 400 boreholes is calculated by using Equation (31) and the contour map is drawn by the Kriging interpolation method, as shown in Figure 13.

Thematic Map of CMFB System
Taking CMFB with a curing time of 56 d and a fly ash content of 75% as an example, its UCS is 7.499 MPa and its elastic modulus is 1.5 GPa. According to the coal seam burial depth of different boreholes, the compressibility of CMFB corresponding to 400 boreholes is calculated by using Equation (31) and the contour map is drawn by the Kriging interpolation method, as shown in Figure 13. mining panel, and the variation coefficient is, therefore, less than 30%. The first-mined coal seam is generally free of gangue. Only a small part of the coal seam contains 1 or 2 layers of gangue, indicating that the coefficient of the dirt band is less than 10%. The first-mined coal seam belongs to long-flame coal or non-caking coal, and the Proctor coefficient is from 2 to 2.5.
(a) (b) Taking CMFB with a curing time of 56 d and a fly ash content of 75% as an example, its UCS is 7.499 MPa and its elastic modulus is 1.5 GPa. According to the coal seam burial depth of different boreholes, the compressibility of CMFB corresponding to 400 boreholes is calculated by using Equation (31) and the contour map is drawn by the Kriging interpolation method, as shown in Figure 13.

Prediction Results of the Suitability of Using CECB to Sequestrate CO 2
Based on the calculation results of the comprehensive evaluation value Φ of 400 boreholes, the contour map of the appropriateness degree Φ in the Yu-Shen mining area was plotted using the Kriging method, as shown in Figure 14.

Prediction Results of the Suitability of Using CECB to Sequestrate CO2
Based on the calculation results of the comprehensive evaluation value Φ of 400 boreholes, the contour map of the appropriateness degree Φ in the Yu-Shen mining area was plotted using the Kriging method, as shown in Figure 14. The grades of the adaptability of various boreholes in the northwest area are generally higher than those of the southeast overall. There are three grade I regions, all of which are geographically fragmented in the northwest, southwest, and eastern areas of the Yu-Shen coal area. Good coal seam occurrence and its roof and floor conditions, large distance from the aquifer and coal seam, and high water-richness are typical features of these regions. These three areas can employ the CECB mining method to sequestrate CO2 directly without adopting other countermeasures, and building protection, water body preservation, and conservation can be realized. The appropriateness of the western border and the middle and northeast parts of the mining area belong to grade II. The collieries in these regions can sequestrate CO2 utilizing the CECB method by taking minor supplementary measurements, and coordinated development among village protection, water body preservation, CO2 sequestration, and treatment of coal-based solid waste can be obtained.
The regions, with the degree of appropriateness ranging between grade III and grade IV, are geologically dispersed in the vicinity of the area where the geological structures exist while there is no Salawusu Formation aquifer occurrence. For collieries whose suitability belongs to grade III and grade IV, it is not suitable for the CECB method to extract coal resources and sequestrate CO2 by means of CMFB under BRW since the economic and ecological benefits are too minor or even lower than the investment. The evaluation results of our research are consistent with other evaluation methods, indicating the prediction model is scientific, reasonable, and reliable.

Discussion
(1) In this paper, a novel method for CO2 sequestration was proposed by developing CMFB and injecting it into the mined-out MRs of CECB mining under BRW. The   The grades of the adaptability of various boreholes in the northwest area are generally higher than those of the southeast overall. There are three grade I regions, all of which are geographically fragmented in the northwest, southwest, and eastern areas of the Yu-Shen coal area. Good coal seam occurrence and its roof and floor conditions, large distance from the aquifer and coal seam, and high water-richness are typical features of these regions. These three areas can employ the CECB mining method to sequestrate CO 2 directly without adopting other countermeasures, and building protection, water body preservation, and conservation can be realized. The appropriateness of the western border and the middle and northeast parts of the mining area belong to grade II. The collieries in these regions can sequestrate CO 2 utilizing the CECB method by taking minor supplementary measurements, and coordinated development among village protection, water body preservation, CO 2 sequestration, and treatment of coal-based solid waste can be obtained.
The regions, with the degree of appropriateness ranging between grade III and grade IV, are geologically dispersed in the vicinity of the area where the geological structures exist while there is no Salawusu Formation aquifer occurrence. For collieries whose suitability belongs to grade III and grade IV, it is not suitable for the CECB method to extract coal resources and sequestrate CO 2 by means of CMFB under BRW since the economic and ecological benefits are too minor or even lower than the investment. The evaluation results of our research are consistent with other evaluation methods, indicating the prediction model is scientific, reasonable, and reliable.

Discussion
(1) In this paper, a novel method for CO 2 sequestration was proposed by developing CMFB and injecting it into the mined-out MRs of CECB mining under BRW. The adaptability of this concept is investigated, to generalize and apply this method in the future. For common CO 2 sequestration methods, the supercritical CO 2 is injected into the subsurface structures such as abandoned oil reservoirs, deep saline aquifers, un-minable coal seams and so forth [55][56][57]. However, CO 2 may leak and escape into the air again while being confronted with large disturbances of the overlying strata. Therefore, developing CMFB and introducing it into the mined-out area by the CECB mining method is the most effective approach to sequestrate CO 2 . As far as economic benefits are concerned, the investment in our method is far less than the conventional methods for CO 2 sequestration. What is known to us all is that special facilities such as pressurizing equipment which pressurizes CO 2 gas to turn it into supercritical CO 2 are necessary for conventional CO 2 storage. The extra apparatus indicate high cost and the economic benefits are therefore limited. In contrast, the CECB is usually employed to extract the coal body under BRW [58,59]. The filling body here is necessary and the CMFB is developed only by adhering CO 2 gas to the fresh slurry of the filling body. The CMFB share the filling transportation pipeline and other facilities with traditional cement filling materials. The extra investment for installing the tank to store the CO 2 gas is too low to ignore. As illustrated in the paper, the process of developing CMFB is similar to the common filling body development. The reaction rate of CO 2 gas and fly ash, and coal-based solid wastes, at ambient temperature and pressure, is rapid by adding catalyzers such as silicate additives, which makes it feasible for the large-scale industrialization of CO 2 sequestration. The time cost of CO 2 sequestration using CMFB and CECB is of great minority. To sum up, compared with conventional CO 2 storage, the method proposed in our paper has advantages in both economic benefits and time cost over other traditional methods for CO 2 sequestration.
(2) The CO 2 mineralization process with solid wastes should be carried out continuously in a space where fresh slurry cannot seep out. Meanwhile, the MR of the CECB mining method can tackle this problem easily since it features excellent leakage-isolated properties. In addition, the mining disturbance at the storage site of CMFB should be mitigated as much as possible since the excessive disturbance may cause the instability of CMFB. The span of MR is narrower than that of longwall mining, contributing to the restricted mined-out space and thus maintaining the stability of the mining stope. It is noted that this paper is only a preliminary study of using CECB to sequestrate CO 2 by means of developing CMFB and injecting it into the mined-out area. In the process of field practice of this method, there are many problems waiting to be resolved, such as large-scale capture of CO 2 gas, CMFB transportation, and the weakening mechanism of CMFB soaked with water.
(3) To avoid the mutual restrictions between extracting and filling, these two operations are separated into two different positions in the CECB mining method. This is essentially different from the operations of traditional longwall backfill mining. The CMFB can maintain stability before reaching its designed strength owing to the limited exposure of each MR, which leads to restricted roof subsidence. It is necessary to investigate the suitability of using other backfill mining methods to sequestrate CO 2 by developing CMFB. For instance, the mining processes of longwall backfill mining are essentially different from those of CECB mining. The extracting and backfill operations are not totally separated from each other, giving rise to mutual influence, and thus it is unsuitable to be employed to sequestrate CO 2 . In addition, the CMFB in CECB mining stays in three-dimensional stress, while the CMFB in longwall backfill mining is in uniaxial stress. This is apt to result in CMFB instability since its three-dimensional strength is far greater than that of the uniaxial one. Furthermore, if the slurry of CMFB is injected into the flexible bags in longwall backfill mining, they may suffer from instability owing to weathering caused by torn flexible materials.
(4) The mathematical model for predicting the adaptability of using CECB to sequestrate CO 2 can provide references to the generalization and application of this method. In future research, the model will be employed in other mining areas that are still in preparation or exploration and have not yet been put into production. For instance, it can be utilized to evaluate the adaptability of the Ningdong, Shendong, and Huanglong mining areas, since they share similar engineering and hydrogeological conditions with the Yu-Shen mining area. Whereas it is inappropriate for the model to be applied in the Xinjiang mining area, owing to its extremely thick coal seam (>100 m) and totally different geological and engineering conditions. To accurately predict the suitability and hence promote its application in a larger area, the reidentification of influencing factors and the redistribution of weights in the prediction model should be carried out by the scholars engaged in coal extraction in the Xinjiang coal area.
(5) To analyze a complex problem, the AHP-fuzzy comprehensive evaluation method has been widely used. Generally, this research usually takes the average membership degree of one influencing factor in various regions to represent the whole evaluation object, indicating the inaccuracy and irrationality of the prediction results. The innovation point of our research is that we evaluated the aptness of using CECB to sequestrate CO 2 in the entire mining area and plotted the zoning map of the adaptability grade based on the evaluation results of 400 drillholes in various locations in the Yu-Shen mining area. The research findings can offer a prospective reference for the generalization and application of employing CECB to realize CO 2 sequestration in the future.

Conclusions
The conclusions drawn from the research are as follows: (1) A novel concept of CO 2 sequestration by developing CMFB together with injecting it into the mined-out area of CECB mining is proposed. Currently, the filling materials are primarily cemented backfill, solid backfill (crushed gangue), and high-water content backfill. However, none of these aforementioned filling materials takes CO 2 sequestration into consideration. In this paper, industrialized CO 2 sequestration is incorporated into the framework of green mining and the "trinity" of low-carbon, sustainable, and ecologically protective mining is proposed. This is, large-scale and permanent CO 2 sequestration, harmless treatment of coal-based solid wastes, and coal extraction under BRW, which is innovative compared with the previous research; (2) The CMFB was developed by making CO 2 react with coal-based solid wastes at ambient temperature and pressure. Fly ash was employed as aggregate, while CO 2 gas, silicate additives, and cement were taken as accessories. The UCS of CMFB with various curing times and different fly ash contents was tested indoors. The results show that, as the fly ash content increases, the UCS ascends first and then decreases, reaching the peak when the fly ash content is 75%. Furthermore, the Young's modulus of the CMFBs of various setting times and material ratios is determined, so as to determine the membership degrees of the two indicators in the CMFB system; (3) The AHP-fuzzy comprehensive evaluation method was employed to construct a triple-leveled hierarchy model for predicting the adaptability of using CECB to sequester CO 2 . Six sub-indicators and nineteen tertiary indicators were chosen as the factors, and the weight distribution and membership functions of the influencing factors were analyzed and determined. The results show that, among the secondary influence indicators, the BWR system, with a weight of 0.3262, is the most important one, followed by the geological structures system and the coal seam system, whose weights are both 0.2073. The distance from the overlying aquifer to the coal seam, with a weight of 0.1631, is the critical tertiary indicator; (4) The prediction model was generalized and employed to evaluate its applicability in the Yu-Shen mining area. The thematic map of a single factor was drawn to determine the membership degree of each indicator. Then, the adaptability of 400 boreholes in the entire coal area was calculated, and the Kriging interpolation method was utilized to plot the zoning map of suitability in the whole mining area. The prediction results can provide a theoretical basis and practical reference for the field application of developing CMFB and using CECB to sequestrate CO 2 gas in the future.

Conflicts of Interest:
The authors declare no conflict of interest. the ratio of the fly ash to the total of the fly ash and the cement is 55% FA65 the ratio of the fly ash to the total of the fly ash and the cement is 65% FA75 the ratio of the fly ash to the total of the fly ash and the cement is 75% FA85 the ratio of the fly ash to the total of the fly ash and the cement is 85%