Drawing mechanisms for top coal in longwall top coal caving (LTCC): A review of two decades of literature

This review details the state of the art in research on top coal drawing mechanisms in Longwall top coal caving (LTCC) by examining the relevant literature over the last two decades. It starts with an introduction of the brief history and basic procedures of LTCC. The framework of research on the drawing mechanism, basic concepts, and some theoretical models of LTCC are detailed in sect. research framework of top coal drawing mechanism. The authors note that the Top coal drawbody (TCD), Top coal boundary (TCB) and Top coal recovery ratio (TCRR) are key factors in the drawing mechanism. The Body–boundary–ratio (BBR) research system has been the classic framework for research over the last 20 years. The modified Bergmark–Roos model, which considers the effects of the supporting rear canopy, flowing velocity of top coal, and its shape factor, is optimal for characterizing the TCD. A 3D model to describe the TCB that considers the thicknesses of the coal seam and roof strata is reviewed. In sect. physical testing and numerical simulation, the physical tests and numerical simulations in the literature are classified for ease of bibliographical review, and classic conclusions regarding the drawing mechanism of top coal are presented and discussed with elaborate illustrations and descriptions. The deflection of the TCD is noted, and is caused by the shape of the rear canopy. The inclined coal seam always induces a larger TCD, and a deflection in the TCD has also been observed in it. The effects of the drawing sequence and drawing interval on the TCRR are reviewed, where a long drawing interval is found to lead to significant loss of top coal. Its flowing behavior and velocity distribution are also presented. Sect. practical applications of drawing mechanisms for LTCC mines 4 summarizes over 10 cases where the TCRR of LTCC mines improved due to the guidance of the drawing mechanism. The final section provides a summary of the work here and some open questions. Prospective investigations are highlighted to give researchers guidance on promising issues in future research on LTCC.


Introduction
Longwall top coal caving (LTCC) is a special type of mining method commonly used to extract coal seams thicker than 4.5 m ). This method originated in France and Eastern Europe in the 1950s for extracting coal pillars from marginal areas of coalfields . LTCC is mainly used in China these days, and has been applied in a few cases in Australia, Vietnam, Russia, Turkey, and Slovenia (Le et al. 2017(Le et al. , 2019Song et al. 2020;Jangara and Ozturk 2021). According to the China National Bureau of Statistics, the country produced 3.75 billion tons of raw coal in 2019, and mined the largest amount of it using LTCC in the world.
With 30 years of development, LTCC has become a representative mining technology in China (Wang 2014(Wang , 2018. Figure 1a illustrates the basic operations of LTCC. The coal seam in LTCC consists of two parts: 1) bottom coal seam, and 2) top coal. The bottom part is undercut by a shearer located on an armored face conveyor. This is a typical layout in conventional longwall mining. Top coal is fractured by the effect of mining-induced stress (Alehossein and Poulsen 2010;Sajjad et al. 2019;Hao et al. 2020), repeated movement of the shield beams (Xie and Zhao 2009a;Wang et al. 2015c), or is artificially broken through hydraulic fracturing (Huang et al. 2015(Huang et al. , 2018 or pre-blasting (Lai et al. 2014). The Drawing opening (DO) is embedded into the rear of the hydraulic support in LTCC. Once the DO is opened, the fractured top coal flows to the Armored face conveyor (AFC) through the DO. This process is often called top coal drawing. From the viewpoint of mechanics, the top coal drawing process is dominated by gravity and frictional forces. This indicates that research on the top coal drawing mechanism should be conducted within the scope of gravitational flow. In metal mines, block caving is often used. It is an underground mining method that uses gravity to exploit massive, steeply dipping ore bodies located at a considerable depth. Many studies on top coal drawing mechanisms often use gravitational flow in block caving for reference (Wang et al. 2015d;Jin et al. 2016Jin et al. , 2017Song et al. 2018;Song and Konietzky 2019). This study reviews and summarizes of top coal drawing mechanisms in LTCC excluding those based on gravitational flow in block caving. The reader interested in research on gravitational flow can refer to Melo et al. (2007Melo et al. ( , 2008Melo et al. ( , 2009, Castro et al. (2016), and Sánchez et al. (2019). The primary objective of this work is to present a comprehensive account of the state of the art in research on top coal drawing mechanisms, and herein lies its originality. The authors use the literature over the last 20 years to provide a global perspective of published articles and highlight potential areas for future research on LTCC. The remainder of this paper is organized as follows: Sect. 2 summarizes the framework of research on top coal drawing mechanisms. Some theoretical models, basic concepts, and theories are summarized and discussed. Section 3 provides a comprehensive summary of the relevant physical testing and numerical simulations. The corresponding experimental and numerical results are discussed in detail rather than merely listed. Section 4 presents more than 10 practical cases of the application of the drawing mechanism. Section 5 contains a discussion of outstanding questions and gray areas as an incentive for future work on LTCC. The current limitations and prevalent perspectives are also noted.
2 Research framework of top coal drawing mechanism 2.1 Basic conceptions

Top coal drawbody (TCD)
As shown in Fig. 1, the Top coal drawbody (TCD) refers to the drawn out top coal blocks before the DO closes (termination of the drawing procedure)-see the red zones in Fig. 1. According to the practical principle applied to LTCC mines, the DO should be immediately closed when the gangue (waste rocks above the top coal) flows out. For current, the close of DO is mainly dependent on the experiences of miner, however, the intelligent control of DO has been applied in a few cases in China. The cameras set at the end of the rear canopy will automatically recognize the gangue by using the picture analysis algorithm, so the close of the drawing opening can be realized automatically, this technology has now applied in some LTCC mines and will be soon used by more LTCC mines in near future. Dominated by gravity, the shape of the TCD quite is similar to the isolated extraction zone (IEZ) in block caving. However, the hydraulic drawing support (DS) used in LTCC leads to the formation of a distinct boundary in comparison with the block caving where no DS is used, and the rear canopy of the DS influences the shape of the TCB Wang et al. 2016d). Furthermore, the DS moves forward after each drawing procedure and the TCB evolves with the advance of the DS. In summary, top coal drawing is a continuous and interactive process ).

Top coal boundary (TCB)
The top coal boundary (TCB) is the continuous interface between the top coal and the gangue . Because top coal is drawn under zones of the TCB, the process is determined and controlled by the shape of the TCB. The blue line in Fig. 1 indicates the outline of the TCB in a 2D version. A 3D TCB and TCD are plotted in the right part of Fig. 1. It also shows the interrelations between the TCB and the TCD. The TCB directly impacts and restricts the volume of the TCD, which impacts the shape of TCB in the next drawing procedure in turn.

Top coal recovery ratio (TCRR)
The top coal recovery ratio (TCRR) refers to the percentage of recovered (drawn out) top coal divided by the total top coal in designated zones. The TCRR is an effective index to characterize the loss of top coal. A higher TCRR can yield significant economy benefit and little resource loss.

Top coal stress distribution (TCSD)
The stress distribution of top coal refers to the complex contact forces between the fractured blocks of top coal. Interlocking between the blocks always occurs when the localized contact force is large. Arching then occurs, and drawing is terminated due to the effects of jamming and arching.

Top coal flowing behavior (TCFB) and drawing efficiency (TCDE)
Top coal flowing behavior (TCFB) represents features of the flow of the top coal and gangue. It includes the flow trajectory and velocity of top coal. Top coal drawing efficiency (TCDE) is an index to evaluate the extent of fluency of the top coal blocks. The more fluent is the drawing, the higher is the efficiency in LTCC.

Key issues in top coal drawing mechanism
Based on the relevant literature over the last two decades (2000-2020), Fig. 2 shows the key issues, questions, and expected goals in research on top coal drawing in LTCC. The research framework that has largely been used is also shown in Fig. 2. The contributions made by different researchers are summarized in Table 1 based on literature in the last two decades . It presents studies according to issue. The main contributions or conclusions are listed in the rightmost column.

Theoretical models
This section reviews the theoretical models proposed to investigate the top coal drawing mechanism, and these include models of the TCD and TCB.
2.3.1 Models of TCD 2.3.1.1 Bergmark-Roos (B-R) model The Bergmark-Roos (B-R) model is commonly used to describe gravitational flow (Rustan 2000). This model was first used in block caving mining to theoretically characterize the shape of the IEZ (Melo et al. 2007;Tao et al. 2010;Song et al. 2018). Two basic assumptions in the B-R model are as follows: 1) The blocks move in straight lines from their resting point to the DO. 2) The blocks are dominated by gravitational and frictional forces (Kuchta 2002). Figure 3a shows the coordinate system of the original B-R model, and its equations are as follows (Melo et al. 2007;Wang et al. 2016b): Figure 3a, in which h G is the maximum allowed angle of motion, indicates that the particle can move only when Drawing mechanism for top coal in longwall top coal caving (LTCC): A review… 1173 angle is smaller than h G . r max is the distance between the farthest resting point and the DO, t is drawing time, and h is the angle of motion. r D is the distance between the origin of the polar coordinates and the aperture. According to Eq. (1), the resting points with different moving angles can be calculated, following which the boundary of the TCD can be obtained by linking discrete boundary points. However, the original B-R model does not consider the effect of the rear canopy of the DS on the TCB. The modified B-R is shown in Eq.
(2), obtained after considering the effect of the shape of the rear canopy (Wang et al. 2016b, d;Zhang et al. 2019). As shown in Fig. 3b, the maximum allowed angle close to the rear canopy (near side of rear canopy) is set to 90°-b (b is the angle of the rear canopy with respect to the horizontal direction), and is different from the original B-R model, while the far side of rear canopy still uses h G as the maximum allowed angle: Even though it considers the effect of the rear canopy, Eq. (2) still has two major disadvantages  (1) The TCB exhibits a discontinuity at the point connecting the far side and the near side.
(2) The model significantly underestimates the drawing time.
The first shortcoming occurs because of non-uniform maximum allowable angles in the far and near sides. The underestimation of the drawing time is induced by ignoring the interlocking between the top coal blocks. To overcome these flaws, the model is further amended by introducing a velocity correction factor, k, see Eq. (3). In Eq. (3), k F and k N are, respectively, the velocity correction factors for the far and near sides of the rear canopy, t F and t N are the theoretical flowing times at the far and near side, respectively, and t is the actual flowing time of the top coal blocks that can be measured through the physical tests: Equation (3) can address the discontinuities of the TCD and amend the unrealistic drawing time encountered in Eq. (2). However, it still has one significant defect: Once the thickness of a coal seam is fixed, the shape of the TCD is also fixed without considering other factors, such as the shape and diameters of the top coal blocks. The model

TCRR
The manner of evolution of TCRR with respect to drawing height was revealed  TCSD A field investigation-based system for top-coal cavability assessment was proposed Wang et al. (2018b) TCSD The mechanism of force chain evolution in fractured top coal was revealed 2.3.1.2 Kinematic model The kinematic model is a classic means of characterizing the flow of granular particles, and was proposed by Nedderman and Tüzün (1979). The basic assumptions for the kinematic model are as follows: (1) The bottom of the DO is flat and the particles are loosely packed. In the 3D kinematic model, the particle velocity is shown as in Eq. (5): (2) The continuous flow of granular particles is regarded as a fluid state, and thus the movements of particles meet the continuity equation in Eq. (6): In Eqs. (5) and (6), V x , V y , and V z are particle velocities along the x, y, and z directions, respectively, and D p is the diffusion coefficient, which is proportional to particle size (Litwiniszyn 1963). The expression for the TCD is then obtained as in Eq. (7): where, Q t is the total volume of the drawn out particles within a period t, and z 0 = z 0 (r 0 ,t) are the particle coordinates when t = 0. Considering that D p is related to the particle size in Eq. (7), the kinematic model can be used to investigate the effect of size of top coal on the shape of the TCD. When r 0 = 0, the maximum height of the TCD (z 0max ) according to Eq. (8) is As is known from Eq. (8), when z 0max is fixed, Q t increases linearly with D p , which also increases linearly with the size of the top coal particles. Figure 4 shows the theoretical shapes of the TCD exposed to a fixed thickness of the coal seam and different sizes of top coal.   (2008) The abbreviations in Table 1 are provided in the Nomenclature 2.3.1.3 Stochastic medium model Litwiniszyn (1956) first proposed stochastic medium theory to investigate the flow behavior of particles using probabilistic statistics. The model has been continually improved by scholars and introduced to block caving mining (Liu and Zhang 1995;Qiao et al. 2003). In recent years, it has also been applied to theoretical research on LTCC, and can adequately describe the shape of the TCD (Yu et al. 2017;Zhu et al. 2018). The assumptions are as follows:

90-β
(1) The blocks of top coal and gangue are completely granular blocks.
(2) Top coal and gangue can flow continuously, and are treated as a random medium.
(3) Their movement is independent in different directions.
Equations (9) and (10) are the probability density equations of top coal along the directions of advance and layout of the DS, respectively.
A is average cutting coefficient of the rear canopy, and can be given by the following equations The parameter A in Eq. (11) is defined as the average cutting coefficient of the rear canopy, which can be calculated from the following method: it is assumed that the left particles moving boundary coincides with the extension line of the rear canopy, and A' in Eq. (11) is the cutting coefficient of the rear canopy. From the perspective of statistics, there is at least one particle in any layer to move down during top coal drawing. To reduce the calculation complexity and difficulty in solving the integral equation, A' is averaged along the vertical direction and A is obtained. where h is the dip angle of the rear canopy, H is the height of the seam, a and b are the coefficients along the direction of advance of the DS, and a 1 and b 1 are coefficients along direction of layout of the DS. After combining Eqs. (9) and (11), the theoretical equation of the TCD along the direction of advance of the DS can be obtained in Eq. (12), where z H is the height of the drawnout top coal: For LTCC and blocking caving, the movements of top coal and gangue are both dominated by gravity, thus some basic models for blocking caving can be also applied to the LTCC. However, attention should be paid to due to the specific features of LTCC: (1) Drawing supports are arranged in LTCC panels, with an advancing process (2) Density of the drawn top coal is smaller than that of gangue in LTCC panels, while density of the ore block is much larger.
In summary, the B-R model, kinematic model and stochastic medium model are three commonly used models to characterize TCD. The B-R model mainly considers the force of the loose top coal particles, and then obtains the drawing body equation by calculating the moving path of the particles. For kinematic model, the state of top coal is regarded as the uninterruptible fluid, and the drawing body equation is obtained through the continuity equation. For stochastic medium model, the drawing body equation is calculated by considering the flowing probability of the particle. The three methods have distinctive advantages and disadvantages. In terms of the drawing body shape, it is wide for the top part and narrow for the bottom part in B-R model. However, it is opposite in kinematic model. In terms of application range, B-R model is more suitable to establish the mechanical model.

Models of TCB
2.3.2.1 2D model Based on stochastic medium theory, a 2D theoretical equation for the TCB can be obtained by analyzing the flow trajectory of the particles (Yu et al. 2017;Zhu et al. 2018). The initial coordinates of a given particle are (x 0 , y 0 ), and the volume of the drawn out top coal is S f when the particle moves to a new location (x, y). The flow trajectory of top coal and S f can be obtained by Eqs. (13) and (14), respectively, where a and b are the coefficients along the direction of advance of the DS as documented in Eq. (11): Combining Eqs. (13) and (14), the dynamic equation of the TCB can be obtained as shown in Eq. (15): According to Eq. (15), the shape of the TCB can be obtained when subjected to different values of S f in the drawing process. When the first gangue flows out through the DO, the TCB is funnel shaped, and is also called the initial TCB. The theoretical equation of the initial TCB is  Fig. 6, curves of various colors illustrate the theoretical shapes of the TCB subjected to different values of S f . When the first drawing process is terminated, S f is 37.5 m 2 and the shape of the TCB is indicated by the black curve in Fig. 6. 2.3.2.2 3D model Based on granular mechanics, the 3D equation of the TCB is deduced by analyzing the stress acting on the particles of top coal. As shown in Fig. 7, the force analysis includes the vertical force (F N ), lateral force (KcHDS), frictional force (f), and the weight of the overlying gangue (cHDS). c is the unit weight, DS is the stressed area, H is the thickness of the overlying gangue pillar, and K is the lateral pressure coefficient.
When the drawing process is finished and the DO is closed, the particle reaches a balanced state as shown in Eq. (17): According to the Cartesian coordinate system in Fig. 7, Eq. (17) is transformed into Eq. (18).
where, m is the transfer coefficient of lateral force and l is the frictional coefficient between loose particles of top coal. Then, the equation of the TCB, Eq. (19), can be deduced through the integral of Eq. (18), where M r and M c are the thicknesses of the immediate roof and the coal seam, respectively: The 3D theoretical equation of the TCB can be derived by rotating Eq. (19): Figure 8a shows the theoretical 3D shape of the initial TCB, which fits the shape of the TCB in physical experiments by comparing the size and boundary of the DO. In addition, as shown in Fig. 8b, the theoretical equation of the TCB is also deduced when an excessive amount of gangue flows out, and the theoretical shape is different from the initial shape of the TCB such that it includes an elliptical top opening and rectangular bottom opening.
3 Physical testing and numerical simulation 3.1 Brief summary of physical testing on drawing mechanism Table 2 summarizes the physical tests documented in the literature in last two decades, and classifies them based on the geometries of the relevant models (dimension, scaling, etc.) and the research items (such as the TCRR, TCD, and TCB were considered). The physical tests were often conducted to simulate multi-round top coal drawing procedures. According to Table 2, a scaled ratio (usually varying between 1:20 and 1:40) is used to reduce the intensity of labor during the tests. Small rock fragments are often used to simulate the fractured top coal and gangue in lab testing (Huang et al. 2007;Wang et al. 2015a;Zhang et al. 2015aZhang et al. , 2018a, including limestone (Wang et al. 2019a), granite, and marble (Wang et al. 2015d, a). Their frictional coefficients are close to those of actual top coal and gangue in mines. Figure 9 shows typical devices used for top coal drawing tests (Wang et al. 2019b, a). A 3D testing device is shown in Fig. 9a, in which the DS can move forward as controlled by a screw to simulate the advancing support in mines. The corresponding results can be obtained and analyzed along the directions of advance and layout of the DS. A quasi-3D device is shown in Fig. 9b; the DS in this device cannot move forward due to restrictions along the direction of advance of the DS. Analysis is thus possible only along the direction of layout of the support.
''Marked particles'' are often used in physical tests to characterize the shape of the TCD (Wang et al. 2015d(Wang et al. , 2019aZhang et al. 2018a). Marked particles refer to numbered markers (the same materials as the top coal and gangue) that correspond to unique coordinate-related information to characterize location (Jin et al. 2017). Figure 10 illustrates some typical layouts of marked particles in physical tests (Wang et al. 2015d(Wang et al. , 2019aZhang et al. 2018a). The spacing between marked particles directly impacts the accuracy of description of the TCD. The dense layout of the marked particles usually leads to a smooth boundary of the TCD. The TCD and TCRR can be calculated according to the location information and mass of the drawn out marked particles Zhang et al. 2018a;Wang et al. 2019bWang et al. , 2020b. Marked particles are also commonly used in block caving testing (Sharrock 2008;Jin et al. 2017 Table 4 lists the numerical studies on top coal drawing mechanisms classified according to the software used, dimensions considered, and research items. Compared with the physical tests, the numerical methods exhibit specific advantages as well as defects. A comparison between the physical tests and numerical methods is presented in Table 3. Figure 11 illustrates a numerical model of an LTCC panel. The rear canopy is modeled by considering the actual shape of the DS, and the DS can continuously move forward to simulate the advancement of the DS in the mine. Top coal can flow out of the DO when it is opened, and the thicknesses of the coal layers and gangue can be modeled by entering different parameters. The velocity of advancement of the DS as well as the drawing intervals can be assigned with different values according to the researcher's demands. For in-situ LTCC mines, the top coal above the rear canopy is fractured. When DO is opened, the top coal blocks can flow to the rear scraper conveyor. Therefore, the discrete methods, such as PFC and F/DEM, can be used to investigate the drawing mechanism. For the characteristics of roof movement and deformation in LTCC panels, the finite difference or finite element software, such as FLAC, is more convenient and applicable. In addition, the coupling of CDEM and bonded models makes it possible to analyze the top coal from the intact to fractured state.

Classic conclusions based on physical testing and numerical simulations
Physical tests and numerical simulations are both effective approaches to investigate the top coal drawing mechanism in LTCC. Based on them, some classic conclusions can be Drawing mechanism for top coal in longwall top coal caving (LTCC): A review… drawn. Some conclusions or mechanisms drawn from physical tests have also been validated by numerical simulations, and these are regarded as relatively reliable results. In this section, we summarize some classical conclusions that have been verified by both physical tests and numerical simulations.

Conclusions on TCD
3.3.1.1 Effect of rear canopy of DS on TCD Figure 12 shows the effect of the rear canopy of the DS on the TCD based on experimental and simulation-based insights. As plotted in Fig. 12a, the shape of TCD is illustrated based on the drawn out marked particles in physical tests. The TCD without considering the effect of the DS (DS is not modeled) is approximately a symmetric ellipse (Wang et al. 2015d). This is similar to the shape of the IEZ in block caving mining (Melo et al. 2007;Jin et al. 2017;Song et al. 2018). When the effect of the DS is considered, it is clear that the TCD exhibits a deflection along the direction of advance of the support. To quantitatively characterize the deflection, a deflection angle, h, is proposed. It refers to the angle between the direction of the long axis of the ellipse and the plumb line, as shown in Fig. 12a. It shows that the extent of the deflection (h) decreases with an increase in the height of the TCD (Wang et al. 2015d), which indicates that the effect of the DS is more pronounced when the height of the TCD is small, and the effect of the deflection is weakened when the height is much greater than the height of the DS. A similar  phenomenon is observed in Fig. 12b based on the results of a numerical simulation. With an increase in the height of the TCD, the deflection in it decreases. The simulation shows that the TCD shows a specific ''cut-variant ellipsoid'' shape when the effect of the rear canopy of the DS is considered Wang et al. 2016d, b). Figure 12b also compares the extent of deflection of the TCD between the approaches, and shows that the deflection in both physical tests and numerical simulations exhibits similar behaviors with regard to TCD height (z coordinate). It is concluded that the rear canopy of the DS causes deflections in the TCD. Figure 12c shows the effect of the angle of the rear canopy on the TCD . As the angle increases from 45°to 72°, the deflection of the TCD is weakened. It can be concluded that a smaller value of the angle of the rear canopy corresponds to a more pronounced deflection of the TCD. Figure 13 illustrates the shapes of the TCD subjected to different dip angles of the coal seams (0-50°) based on physical testing and numerical simulations (Zhang et al. 2018b, a;Wang et al. 2019b). Figure 13a shows the effect of the dip angle on the TCD in physical tests, and the results of the corresponding numerical simulation are shown in Fig. 13b. It is clear that the TCD exhibits a similar pattern in both approaches. With an increase in the dip angle, the TCD exhibits a more pronounced deflection toward the upper side of the face of the panel. In addition, the quality (volume) of the TCD increases with the rise in seam inclination. Figure 14 summarizes the effect of the dip angle on TCD quality (volume) based on the results documented by Yang et al. (2016), Zhang et al. (2018b) and Wang et al. (2019b). It reveals that with the rise in the dip angle (at least within the range of 0-50°), the quality or volume of the TCD exhibits an approximately linear increase. One interpretation of the increased TCD quality when subjected to an inclined seam is that the inclination of the coal seam its thickness as well as the drawing height of the TCD.

Effect of sizes of top coal and DO on TCD
Many studies have reported that the sizes of top coal and DO have a significant effect on the shape of the TCD (Wang et al. 2016c(Wang et al. , 2019a(Wang et al. , 2020a. Table 5 shows the detailed size distribution of top coal in three schemes in a physical test. The percentage of top coal samples of small size is high in scheme 2 but low Melo 3, which indicates that the average size of top coal is the smallest in scheme 2 and the largest in scheme 3. In addition, the variation in the size of the DO changes the relative size of top coal (size of top coal divided by the size of DO), resulting in the different movements of top coal Wang et al. 2020b). Figure 15 illustrates the shape of the TCD subjected to three different sizes of top coal and DO. By analyzing the number and position of the marked particles drawn from different DS, the number of drawn out marked particles in the right side is larger than that in the left side. This is owing to smoother contact between the top coal and the DS, which causes the TCD to develop toward the support side. It can be concluded that when the size of the DO is fixed, the volume of the TCD linearly increases with the growth in the size of top coal. When the latter is fixed, the ratio of width to height of the TCD gradually increases with the size of the DO, which indicates that the TCD is fatter.

Conclusions on TCRR
The Top coal recovery ratio (TCRR) is an important economic indicator to evaluate the production efficiency of an LTCC mine (Wang et al. 2015b;Zhu et al. 2018). A low TCRR often results in the loss of a large amount of permanent coal. Many scholars have investigated methods to improve the TCRR by using lab testing and numerical simulations (Xie and Zhao 2009b;Wang et al. 2016d;Yang et al. 2016;Zhang and Wimmer 2018;Zhu et al. 2018;Zhang et al. 2019). This section focuses on the effects of the drawing parameters (such as drawing interval, drawing sequence, top coal size, dip angle of coal seam, etc.) on the TCRR, 3.3.2.1 Effect of drawing sequence on TCRR Horizontal coal seam The drawing sequence has a direct impact on the TCRR in a horizontal coal seam (Wang et al. 2016a, c). Table 6 presents the TCRR in the horizontal coal seam subjected to different drawing sequences, investigated using both physical tests and numerical simulations. Nine schemes were used, and the corresponding results of the TCRR are shown in Table 6. Schemes 1-3 considered five DS in a physical test. The results show that the sequential drawing (scheme 1) had the lowest TCRR (80.7%). Schemes 2 and 3 used the alternate drawing, and one with the largest interval (scheme 3) had the largest TCRR. Schemes 4-9 studied the effects of the DO size (number of DO used in drawing) on the TCRR by performing experiments and numerical simulations. They concluded that the TCRR increases with DO size. To sum up, alternate and multi-DO drawings are beneficial for increasing the TCRR. Inclined coal seam Figure 16 illustrates the effects of the drawing sequence on the TCRR in a steeply inclined (around 80°) coal seam (Wang et al. 2016d). Two drawing sequences were used: (1) roof-floor drawing, shown in Fig. 16a; (2) floor-roof drawing, shown in Fig. 16b. The losses in coal for the two plans are clearly distinct. The roof-floor drawing led to a clearly higher loss, which corresponded to a smaller TCRR than the floor-roof drawing. The simulation in Fig. 16c also confirms this conclusion. This results hint that for a steeply inclined seam, the floor-roof drawing is preferable for a higher TCRR. Figure 17 summarizes the effects of the drawing interval on the TCRR. The drawing procedures used on site are usually of three types:

Effects of drawing interval on TCRR
(1) Top coal drawing is performed after one-time support advance. (2) Top coal drawing is performed after two-time support advance. (3) Top coal drawing is performed after three-time support advance.   Figure 17 presents the TCRR versus different drawing internals at a lab testing scale based on the published literature. Some studies have conclude that ''drawing after two-time support advance'' leads to the largest TCRR (Huang et al. 2006(Huang et al. , 2008Wang et al. 2014). However, other studies have obtained different results (Wang 2008). As shown in Fig. 17, data from Wang (2008) shows that ''drawing after two-time support advance'' does not yield the largest TCRR in most cases, and ''drawing after onetime support advance'' usually leads to the largest TCRR.  The same study (Wang 2008) also examined the TCRR when subjected to distinct cutting-drawing ratios (C/D ratio: the ratio of bottom coal/top coal thickness), which indicates that the impact of the drawing interval on the TCRR is also influenced by other parameters such as the C/D ratio. However, it can least be concluded that ''drawing after three-time support advance'' always yields the smallest TCRR. To sum up, a larger drawing interval is not recommended as it may cause a smaller TCRR.

Conclusions on TCFB
3.3.3.1 Flow trajectory The flow trajectory can be precisely visualized only with the aid of a numerical simulation. Physical testing cannot trace the flow trajectory of a given particle during the drawing process. Figure 18 presents the trajectories of particles of top coal in a numerical simulation. Figure 18a shows three trajectories of different particles of top coal , where flow traces and signifies that particle trajectories subjected to the opening of at least two different DO. Figure 18b shows that once two or more drawing procedures were used (Ball 144,071 had two drawing procedures and Ball 132,082 had three), the trajectories deviated from a straight line. The Drawing mechanism for top coal in longwall top coal caving (LTCC): A review… 1189 trajectory exposed to multiple drawing procedures can be divided into several sub-trajectories, where each points to the opened DO at the give time. For example, the particles in Fig. 18b respectively have two and three sub-trajectories, from which it is clear that each sub-trajectory points to the opened DO at the time. This shows that even a complete trajectory exposed to multiple drawing procedures is not a linear line, but the sub-trajectory still confirms the assumption (linear trajectory pointing to the DO) of the B-R model. Figure 18c shows the trajectory of the particles by considering two drawing rounds and one support advance procedure (Song and Konietzky 2019). It is clear that the advance of the support causes the trajectory to develop in its direction. This indicates that the modeling of the support advance should be considered carefully in numerical simulations because it directly impacts the flow trajectory of top coal.

Flow velocity field
The distribution of the velocity field of top coal can be used to characterize flow fluency during the drawing process. The arching always occurs in areas where the flow velocity is low, and can cause a jam. Figure 19 plots the distribution of the flow velocity of top coal during the drawing procedure Zhang et al. 2018a;Song and Konietzky 2019). In Fig. 19a, the flow velocity is proportional to the length of the black arrows; a longer arrow indicates a higher flow velocity. The velocity drastically increases when the particles are closer to the DO. Furthermore, the velocity close to the rear canopy of the support is significantly higher than in any other area, and this can be attributed to the smooth contact between the top coal and the rear canopy of the support. This also leads to the deflection of the TCD presented in Sect. 3.3.1.1. Figure 19 plots the velocity field, which is indicated by different colors. A fluent flow zone (FFZ) is shown, indicating a higher flow velocity. The height of the FFZ is approximately equal to the height of the DS, which also verifies the effect of the rear canopy of the DS on the velocity field. The FFZ in Fig. 19b has a similar shape to the TCD reported in Sect. 3.3.1.1, which also exhibits a prominent deflection.  Practical applications of drawing mechanism to improve production 4 Practical applications of drawing mechanisms for LTCC mines The research system based on the framework documented in Fig. 2 can provide useful suggestions for optimization to improve production in LTCC mines. The results of physical testing and numerical simulations can be used to determine the appropriate drawing parameters, such as drawing sequences and drawing intervals, to obtain a higher TCRR. Table 7 presents the practical applications of drawing mechanisms to improve production in LTCC mines subjected to different geological characteristics based on more than 10 detailed cases investigated by our group over the last two decades. The guidance and improvements are listed in the right column, and verify that a fundamental understanding of the drawing mechanism of top coal is important for optimizing production in LTCC mines.

Summary
This review presented and classified published research on the top coal drawing mechanism in the last two decades. It offers a comprehensive insight into the research frameworks, testing schemes, and modeling methods used to study the drawing mechanism in LTCC. LTCC is an advanced and promising mining technology for exacting ultra-thick coal seams, and research on the drawing mechanism of top coal has important implications for the application of LTCC. This section presents concluding remarks and some prospects for future work in the area.
The B-R model, kinematic model, and stochastic medium model can theoretically characterize the shape of the TCD. The modified B-R model is the most appropriate one to describe the shape of the TCD. The trajectory of flow of top coal is visualized using 3D numerical simulations, and is an approximately linear line point to the DO. The linear trajectory verifies the assumption of the B-R model. The angle of the rear canopy, inclination of the coal seam, and sizes of top coal size and the DO all have direct impacts on the shape of the TCD. A larger dip angle and size of top coal will result in a larger and fatter TCD. Different drawing sequences and drawing intervals lead to distinct instances of the TCRR, alternate and multi-DO drawings are beneficial for enhancing it, drawing after three instances of DS advance yields the smallest TCRR, and an excessively large drawing interval leads to a greater loss of top coal. The velocity field of top coal is plotted, and top coal close to the rear canopy has a significantly higher flow velocity owing to the smooth contact between it and the rear canopy of the DS.

Perspectives on future research
Research on the drawing mechanism of top coal assumes that top coal is completely fractured at the outset. Such over-simplified assumptions neglect the progressive failure and fracture of top coal during extraction as well as continuous change in its size owing to secondary crushing induced by strata pressure. The precise measurement of onsite TCRR is still challenging for now; therefore, highly sensitive and portable measuring devices should be developed to precisely measure the TCRR on site. More theoretical models should be proposed to cater to more complex geological coal seams, such as seams with steep  The inclination of the coal seam will influence the shape of TCD The inclination of the coal seam determines the top boundary of the flowing field for fractured top coal. This is proportional to the height of TCD. Normally, with the increase of the seam inclination, the maximum flowing height will increase as a result. This will result in more drawn-out top coal blocks from the top part of the working face. However, the inclination will also cause the asymmetrical TCD. The TCD will incline to the upper side of the working face. In summary, the inclination of coal seam equivalently increases the height of TCD

Direct factor
Size of top coal The size of top coal will influence the shape of TCD The repose angles of top coal particles with different sizes are distinct. In a certain range, the repose angle decreases with increasing top coal sizes. According to the mathematic relationship between the repose angle and h G in Bergmark-Roos model, it can be obtained that h G gradually increases with the increase of top coal size, which results in a fatter and larger TCD. Therefore, the variation of top coal size changes the physical and mechanical properties of the loose top coal blocks, mainly refers to the change of repose angle and interlocking

Direct factor
Size of DO The size of DO will influence the shape of TCD and TCRR As the boundary condition of top coal drawing, the change of DO's size directly affects the TCD shape, and its essence is to change the range and volume of drawn-out top coal. In general, with rising the length of DO, h G and the size of top coal gradually increase, which makes the volume of TCD gradually rises. However, considering the influence of other factors, TCRR increases gradually with increasing size of DO in extra-thick coal seams. However, for the thinner coal seams, TCRR tends to decrease with increasing the size of DO, which is due to the top coal loss between supports resulted from the increase of h G

Direct factor
Sequence of top coal drawing The optimal drawing sequence should be determined according to the specific geological conditions The sequence of top coal drawing will influence the location as well as the size of DO, the variations have obvious impacts on the flowing fields of the fractured top coal

Indirect factor
Top coal drawing interval The optimal drawing interval should be determined according to the specific geological conditions The drawing interval will affect the shape of TCB, and then change the interaction between TCB and TCD. Specifically, for normal coal seams, increasing drawing interval will reduce the tangent range between TCB and TCD, and increase the volume of top coal loss in the goaf, which is not beneficial to improve TCRR. However, in the extra-thick coal seams, the increase of drawing interval is beneficial to the extraction of upper top coal and can improve TCRR in LTCC panels

Indirect factor
Drawing mechanism for top coal in longwall top coal caving (LTCC): A review… 1193 inclination and variations in thickness. Both physical testing and numerical simulations should be performed while considering the pressure of the strata during the drawing process. The development of modeling that couples the finite element method and the discrete element method would be an interesting approach in this vein.
Authors contribution The author contributions include: JW proposed the research framework and the structure of the manuscript, and reviewed the manuscript. ZS and WW conducted the literature review and wrote the manuscript. SY and JZ reviewed and revised the manuscript.
Funding This work is funded by the National Natural Science Foundation of China (Grant No. 51934008, 51674264 to Jiachen Wang;Grant No. 51974320  Data availability The data can be available from the authors upon request.
Code availability The code can be available from the authors upon request.

Declaration
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.