Next Article in Journal
System-Level Fault Diagnosis for an Industrial Wafer Transfer Robot with Multi-Component Failure Modes
Next Article in Special Issue
Longwall Mining Automation—The Shearer Positioning Methods between the Longwall Automation Steering Committee and China University of Mining and Technology
Previous Article in Journal
Behavioral Intentions to Use Energy Efficiency Smart Solutions under the Impact of Social Influence: An Extended TAM Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Longwall Face Automation: Coal Seam Floor Cutting Path Planning Based on Multiple Hierarchical Clustering

1
School of Mechanical and Electronic Engineering, China University of Mining and Technology, Xuzhou 221116, China
2
Intelligent Mining Equipment Collaborative Innovation Center, China University of Mining and Technology, Xuzhou 221116, China
3
School of Mechanical Electronic and Information Engineering, China University of Mining and Technology, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(18), 10242; https://doi.org/10.3390/app131810242
Submission received: 4 August 2023 / Revised: 7 September 2023 / Accepted: 11 September 2023 / Published: 12 September 2023
(This article belongs to the Special Issue Advanced Intelligent Mining Technology)

Abstract

:
Space adaptability between mining equipment and coal-rock mass, to ensure the machines cut in a coal seam, is an importance technique in longwall mining automation. In order to guide the mining equipment cutting in the coal seam, a cutting path planning method based on multiple hierarchical clustering was proposed. Morphology similarity and the coplanarity measurement method were defined to evaluate the similarity of clusters. The coal seam floor series in the face-advancing direction were clustered according to the morphology similarity and coplanarity, respectively. Taking the morphology-based and coplanarity-based cluster centers as generating lines and stretching angle, respectively, the coal seam floor was reconstructed. The reconstructed floor can be regarded as the cutting path. The coal seam geological model of the 18,201 longwall face was analyzed with the proposed cutting path planning method. Comparing the reconstructed floor and original floor, the amounts of coal left and cut gangue were 1999 m3 and 1856 m3, respectively, for the segmental floor. For the case of whole floor, the amounts of coal left and cut gangue were 5642 m3 and 5463 m3, respectively. The coal loss rates only were 0.57% and 0.87% for the segmental and whole coal seam, respectively.

1. Introduction

Longwall mining, as the most efficient extracting method for underground coal mines, will remain the principal extraction method in the world in the foreseeable future [1,2]. Usually, the equipment used for working in a longwall mining face include a longwall shearer, an armored face conveyor (AFC), and more than 100 hydraulic supports, as shown in Figure 1. The shearer cuts the coal seam while traveling across the longwall face along rails attached to the AFC. The cut coal is transported by the AFC. The hydraulic supports provide temporary support of the roof material above the extracted coal seam. As the shearer moves across the coal seam, large hydraulic push rams attached to the supports advance the AFC progressively in the face-advancing direction, and thereby the longwall mining moves in a programmed sequence. Longwall mining is a coupling process between the mining equipment and coal-rock mass. In addition to the mechanical adaptability between the mining equipment and coal-rock mass (such as coal fracture mechanics induced by the shearer drum and pick cutter [3,4,5,6]), the space adaptability between them, which means that the mining equipment cuts in the coal seam, is a key issue that must be addressed in longwall mining automation.
Longwall mining automation, which can improve productivity and safety for mining personnel, is the focus of the world coal industry. Two key problems facing longwall mining automation are that of machine and face guidance to ensure the machines cut in the seam, and the face is kept straight [7,8]. Face alignment [9,10,11], produced from the shearer positioning with an inertial navigation system (INS), has been researched and applied widely in many coal mines [12,13]. It can keep the face straight for consistent production and minimal wear of the conveyor system. Ensuring the machines cut in the seam needs two technologies: horizon control and extraction control. Horizon control guides the machine to follow the seam undulations in the direction of the longwall face’s advance. Extraction control allows the machine to follow defined floor and roof profiles along the face, better following any given coal seam undulations for optimal extraction and production [14].
Enhanced horizon control (EHC), an outcome of the ARCAP landmark project, which is based on a cut model and some other information (manual input, coal interface input, and other seam information), can make a prediction of what the next floor and roof horizon should be [15]. The EHC is a memory-cut-based method [16], and will absolutely degenerate into a memory cut if the other information (manual input, coal interface input, and other seam information) becomes unavailable. Memory cutting is incapable of working for uneven coal seams [17]. Shield-data-based horizon control (SDHC) [18], based on the fact that the coal seam trend is given by the inclination of the support’s canopy, automatically adjusts the cutting drum height according to data from the roof support so that the longwall face follows the seam trend in the face-advancing direction. In fact, EHC and SDHC are not completely the same as horizon control, described in reference [14]. They include some functions of extraction control, described in reference [14]. The intrinsic and principal problem is to maintain the alignment of the shearer’s exploitation gradient with the coal seam gradient in the face-advancing direction, regardless of EHC, SDHC, and horizon control, described in reference [14]. The coal seam gradient in the face-advancing direction is the control target of horizon control.
The concept of the transparent longwall face [19,20] was proposed, which aims at autonomous cutting by the longwall shearer in the coal seam with the use of a high-precision transparent geological model [21,22,23]. But research about how to decide the cutting path of the longwall shearer according to the geological model is rarely reported. Based on the coal seam geological model, the authors of this paper proposed a coal seam floor gradient identification method using piecewise linear representation (PLR), as shown in Figure 2a [17]. This method synthesized the control resolution of the shearer range arm and the geometric characteristics of the coal seam, so that horizon control was achieved with the maximum exploitation rate within the shearer’s capacity. But this method can only analyze one cross section (such as 2–2′ in Figure 2a) in the face-advancing direction, so it is insufficient for coal seams with undulations in the face direction or different floor series in the face-advancing direction. As shown in Figure 2b,c, the floor series at different face positions have different morphologies, and the gradients in the face direction are different. The distribution of the whole coal seam gradient in the face-advancing direction is necessary for the horizon control of longwall mining.
According to the concept of the transparent longwall face, the objective of this work is to obtain a coal seam floor cutting path planning method based on the geological model. In this work, the distribution of the coal seam gradient in both the face-advancing direction and the face direction was obtained using a multiple hierarchical clustering algorithm, as shown in Figure 3. First, a morphology similarity and coplanarity measurement method was proposed to evaluate the similarity of clusters. The floor series (Fs) were clustered using a hierarchical clustering algorithm according to the morphology similarity. The cluster centers are the representative morphology of the zones covered by each cluster. By analyzing the morphology-based cluster centers using the PLR method in Ref. [17], the distribution of the coal seam gradient in the face-advancing direction can be obtained. The surfaces (Mk) composed of two neighbor floor series are clustered according to the coplanarity. The gradient of the coplanarity-based cluster center is the representative gradient of the zones covered by each cluster, so the distribution of the coal seam gradient in the face direction can be obtained. Meanwhile, taking the morphology-based and coplanarity-based cluster centers as the generating lines and stretching angle, respectively, the coal seam floor is reconstructed. And the reconstructed floor is the planned cutting path.
The rest of this paper is structured as follows. Section 2 introduces the hierarchical clustering algorithm. Section 3 describes the morphology and gradient similarity measurement method. Section 4 introduces the validity measure method of clustering in order to determine the optimal number of clusters. The cutting path planning method and its limitations are described in Section 5. Section 6 presents the results and discussion. Finally, Section 7 concludes the paper and discusses future avenues for research.

2. Hierarchical Clustering Algorithm

Clustering is the most well-known and frequently utilized unsupervised learning technique in data mining, pattern recognition, image segmentation, and so on [24,25,26,27]. The cluster algorithm is a set of methodologies for automatic classification of samples into a number of groups using a measure of similarity, so that the samples in one group are similar to each other and samples belonging to different groups are not similar to each other. Usually, there are seven categories of clustering algorithms, in particular (i) graph-based algorithms, (ii) hierarchical clustering algorithms, (iii) partitioning clustering algorithms, (iv) density-based clustering algorithms, (v) grid-based algorithms, (vi) combinational clustering algorithms, and (vii) model-based clustering algorithms [28].
Hierarchical clustering algorithms create a hierarchical decomposition of the given set of data objects. On the basis of how the hierarchical decomposition is formed, hierarchical methods can be classified into two categories: agglomerative and divisive hierarchical clustering algorithms. The agglomerative hierarchical clustering algorithm is known as a bottom-up approach. It starts with each object forming a separate group. It keeps on merging the objects or groups that are close to one another. It keeps on doing so until all of the groups are merged into one or until the termination condition holds [29]. The divisive hierarchical clustering algorithm is commonly known as the top-down approach. All objects initially belong to a single root group. Then the group is iteratively split up into sub-groups until each object is in one cluster or the termination condition takes holds [29,30]. In the paper, the agglomerative hierarchical clustering algorithm is selected because that it is usually more accurate than the divisive hierarchical clustering algorithm [30].
For both agglomerative and divisive hierarchical clustering algorithms, the termination condition usually is the number of groups or the similarity threshold value. In engineering applications, the number of clusters are not always known a priori. Different similarity measures lead to different types of clusters, and the effectiveness of the clustering algorithm depends intensively on the definition of the similarity. In the next two sections, the similarity measurement and the decision of the cluster number are discussed.

3. Morphology Similarity and Coplanarity Measurement

3.1. Morphology Similarity Measurement

Usually, the distance between two series, such as the Manhattan distance [31], Euclidean distance [32], and dynamic time warping (DTW) distance [33], is used to evaluate the similarity. The smaller the distance is, the greater the similarity is. Both the Manhattan distance and Euclidean distance are incapable of morphology recognition because these two distances are calculated from the absolute differences of vectors. In Figure 2b, The Z-axis coordinates of the coal seam floor series vary largely from 0 to about 70 m. Both the Manhattan and Euclidean distance will be very large even if the floor series have similar morphologies. DTW is usually used to calculate the shape similarity between two series with different lengths. DTW distance is computed by first finding the best warped path to align the two series. It is only because of the warped path that the position information of the floor series is eliminated. It is obvious that DTW distance is not an appropriate evaluation index for the morphology similarity of the floor series.
The gradient of the coal seam floor series is not only the key parameter describing the morphology, but also the target of horizon control. The PLR of the floor series including gradient information can be used to calculate the morphology similarity distance. Figure 4 shows the schematic of the morphology similarity measurement according to the PLR of the floor series. PLRi and PLRj are the PLR series of two floor series, Fsi and Fsj, at different positions on the longwall face. The subscripts i and j represent different X-axis values. The solid circles on each PLR series represent turning points, while the hollow circles represent alignment points that are projected from corresponding turning points on the other PLR series. The distance between Fsi and Fsj, D(Fsi, Fsj), is calculated according to the follow formula,
D ( Fs i , Fs j ) = D ( PLR i , PLR j ) = k = 1 n 1 { ( r a d i k r a d j k ) 2 + ( Δ z i k Δ z j k ) 2 }
where n is the length of the alignment PLR of the floor series. The smaller the distance is, the greater the morphology similarity is. If the distance equal zero, the two series are the same or similar.

3.2. Coplanarity Measurement

The coal seam floor is composed of k surfaces, which are formed by orderly connecting the k + 1 floor series, as shown in Figure 3. Clustering based on coplanarity automatically classifies the k surfaces into a number of clusters. Each cluster center is used to represent a floor zone of the coal seam that is covered by the cluster. In each zone, the gradient in the face direction of each cluster center is take as the extension direction of the cluster centers obtained with clustering based on morphology similarity.
Each surface sample, composed of two neighboring floor series, Fsi and Fsi + 1, is defined as
M i = Fs i Fs i + 1 = ( x i , y i 1 , z i 1 ) , ( x i , y i n , z i n ) , ( x i + 1 , y ( i + 1 ) 1 , z ( i + 1 ) 1 ) , ,   ( x i + 1 , y ( i + 1 ) n , z ( i + 1 ) n )
Figure 5 shows the schematic of the coplanarity measure between surfaces Mi and Mj. A(xiv, ziv) and B(x(i + 1)v, z(i + 1)v), with a same coordinate in the face-advancing direction yv, are the points belonged to Fsi and Fsi + 1, respectively. Line AB is the intersecting line of surface Mi and the X–Z plane at yv. The slope of Mi at yv is defined as θ v < i , i + 1 > . Analogously, the slope of Mj at yv is defined as θ v < j , j + 1 > . Points C’ and D’ are determined with the slope θ v < i , i + 1 > and the X-axis coordinates of points C and D. The collinearity between lines AB and CD, c o l v i , j , is calculated as follows,
c o l v i , j = |   ε j v   | + |   ε ( j + 1 ) v   | + | ε j v ε j v | ε j v = z ( i + 1 ) v z i v x ( i + 1 ) v x i v ( x j v x i v ) ε ( j + 1 ) v = z ( i + 1 ) v z i v x ( i + 1 ) v x i v ( x ( j + 1 ) v x i v )
The coplanarity measurement between surfaces Mi and Mj, C o l ( M i , M j ) , is defined as the collinearity summation of all sub-samples in the two surfaces,
C o p ( Μ i , M j ) = v = 1 n c o l v i , j
where n is the length of the floor series Fsi.

4. Determining the Optimal Number of Clusters

In engineering applications, the termination condition of clustering usually is determined by the number of groups or the similarity threshold value. However, it is very difficult to determine the threshold value of morphology similarity and coplanarity, so the termination condition of clustering has to be determined by the number of groups. There are some cluster validity measures that are concerned with determining the optimal number of clusters and checking the quality of clustering results. In this paper, the Dunn index (DI), Davies–Bouldin index (DBI), and CS were used to determine the optimal number of clusters.

4.1. Dunn Index (DI)

The DI, a well-established cluster validity measure, was proposed by Dunn in the 1970s [34,35]. It is defined as
DI = min ( C ) min j i ( d min ( C i , C j ) ) max 1 l K ( d i a m ( C l ) )
d min ( C i , C j ) = min s i C k , s j C l d i s t ( s i , s j )
d i a m ( C l ) = max s i , s j C l , i j d i s t ( s i , s j )
where C and s represents cluster and samples, respectively. In Equation (5), the denominator represents the maximum distance between two samples belonging to the same cluster. The numerator represents the minimum distance between two samples belonging to different two clusters. The larger the DI is, the better the clustering is.

4.2. Davies–Bouldin Index (DBI)

The DBI, which can measure compact and separate clusters, was proposed by Davies and Bouldin [36]. It is defined as
DBI = 1 K i = 1 K max i j ( C i ¯ + C j ¯ d i s t ( C i 0 , C j 0 ) )
C i ¯ = 1 | C i | l = 1 | C i | d i s t ( s l , C i 0 )
where K is the number of clusters, | C i | is the number of samples in Ci, and C i ¯ is the average distance between all samples and the cluster center sample C i 0 of cluster Ci. The lower value of C i ¯ indicates the higher similarity and compact for all samples of cluster Ci. The denominator in Equation (8), d i s t ( C i 0 , C j 0 ) , is the distance between the cluster centers of cluster Ci and Cj. So this measure is a function of the ratio of the sum of within-cluster scatter to between-cluster separation, and it uses both the clusters and their sample means. The smallest DBI indicates a valid optimal partition.

4.3. CS

Chou [37] proposed the CS validity measure of clustering, which can deal with clusters with different densities and/or sizes. It is defined as
CS = i = 1 K 1 | C i | s j C i max s l C i d i s t ( s l , s j ) i = 1 K min j K , j i d i s t ( C i 0 , C j 0 )
The denominator in Equation (10) computes the average distance between cluster centers. The numerator measures the average largest distance between two data points lying in the same cluster. The denominator measures the sum of the minimum distance between the cluster centers. So, this measure is a function of the ratio of the sum of within-cluster scatter to between-cluster separation. The smallest CS indicates a valid optimal partition.
Among the three validity measure parameters, the DI is an extremely large index. In order to compare and analyze the three parameters conveniently, DI is converted into an extremely small index by reciprocal transformation. Meanwhile, the inverse of the DI, DBI, and CS were normalized with Equation (11).
X = X X min X max X min
where X, Xmax, Xmin are the index parameter, and the maximum and minimum value of the index, respectively. X’ is the normalized value of the index parameter.

5. Cutting Path Planning Based on Clustering Algorithm

Figure 6 shows the flowchart of the agglomerative hierarchical clustering algorithm. The original data of the clustering algorithm include the floor series and its number m, and the predefined cluster number p = m − 1. According to the piecewise linear representation method proposed in reference [17], the floor series are transformed into the PLR series, and then the PLR series are aligned with each other with linear interpolation. Take each individual PLR series as a cluster, calculate the distances between each of the clusters. The pair of clusters with the minimum distance is merged into one cluster, and then update the cluster center and cluster number. Repeat steps 2 to 4 until the cluster number, K, is equal to the predefined cluster number, p. Finally calculate the clustering validity index and cluster results. Steps 1 to 5 are recursively repeated until the predefined cluster number, p, is equal to 1. The clustering validity and dendrogram were recorded during the recursive process. It should be noted that the implementations of the agglomerative hierarchical clustering based on morphology similarity and coplanarity are the same, but the distances between clusters must be calculated by Equations (1) and (4), respectively.
According to the recorded clustering validity and dendrogram, the morphology-based and conplanarity-based cluster centers and their cover zones were determined. The coal seam floor cutting path can be obtained through reconstructing the floor with the morphology-based and conplanarity-based cluster centers and their cover zones. The floor reconstructing process includes: (1) aligning the morphology-based and conplanarity-based clusters, so that the cover zones of the morphology-based cluster and conplanarity-based cluster are coincident; (2) taking the morphology-based cluster center as a generating line; (3) taking the gradient in the face direction obtained from conplanarity-based clustering as the stretching angle; and (4) moving the generating line at the stretching angle along the face direction in a cover zone.
The proposed cutting path planning method based on multiple clustering algorithms aims to guide the longwall shearer cutting in the coal seam by considering the geometry of the coal seam. So it is inapplicable to a low-thickness coal seam, in which undercut rock is unavoidable because the extracting thickness of longwall mining machines is larger than the thickness of the coal seam [38,39].

6. Results and Discussion

The coal seam geological model of the 18,201 longwall face in Shanxi province of China, as shown in Figure 2a, was used to verify the proposed cutting path planning method. The details of this geological model were described in Ref. [17]. In Figure 2a, the X-axis represents the longwall face direction, the Y-axis represents the face-advancing direction, and the Z-axis represents the vertically upward direction. Within the coal seam of the 18,201 face, there is a fault structure, which is located between 2200 m and 2300 m in the face-advancing direction. The clustering analysis of the coal seam floor was respectively performed under two cases: the segmental floor (Figure 7), located between 0 and 2100 m in the face-advancing direction, which does not include the fault structure, and the whole floor (Figure 2a), including the fault structure.

6.1. Clustering of Segmental Floor

6.1.1. Clustering Based on Morphology Similarity

Figure 8 shows the variation of the dendrogram and validity measure index during clustering based on morphology similarity. With the evolution of the clustering process, the neighbor floor series were agglomerated gradually into the same cluster, as shown in Figure 8a. The DBI increased with the decreasing of number of clusters. When the clusters were agglomerated to 36, the DBI reached the maximum. After that, it began to decrease rapidly. The DI decreased slowly until the number of clusters decreased to about 100. And then it increased slowly until the number of clusters decreased to about 5. It increased sharply when the clusters further agglomerated. The CS increased slowly when the number of clusters decreased until the cluster number reached about 120. After that, the CS remained constant at about 0.35 until the cluster number reached 7, and then it increased sharply.
Taking into account the dendrogram and the validity measure index, the number of clusters is determined as 7, as shown in Figure 8a. The morphology of the determined cluster centers is exhibited in Figure 9. According to the intuitive overall morphology, the clusters can be classified into two categories. One category includes the clusters No.1–4, and another includes the clusters No.5–7. This means that these clusters would be agglomerated into two clusters if the clustering process took a step further. This coincides with the dendrogram (Figure 8a). For each category, some observable morphology differences between the clusters are shown in Figure 9.
Table 1 lists the position of the cluster centers and cover zone of each cluster. Each cluster center is the representative morphology of the cluster zone, and can be regarded as the generating line (generatrix) of the cluster zone. The floor in this cover zone can be re-constructed by stretching the generating line with a predefined angle along X-axis. The predefined angle can be determined by the cluster centers based on coplanarity.

6.1.2. Clustering Based on Coplanarity

Figure 10 shows the variation of the dendrogram and validity measure index during clustering based on coplanarity. Similar to the morphology-based clustering, the neighbor floor series were agglomerated gradually into the same cluster with the evolution of the clustering process. During the agglomerative clustering process, the DBI increased to the maximum when the number of clusters decreased to 99. Then it began to decrease with the agglomerative clustering process. The DI decreased with the agglomerating of clusters until the number of clusters decreased to 5. Then it increased rapidly when the number of clusters further decreased. The CS remained almost constant until the number of clusters decreased to about 15. Then it increased sharply. Considering the dendrogram and validity measure index, the number of clusters was selected as 6, as shown in Figure 10a.
Figure 11 presents the gradient surface of cluster centers based on coplanarity for the segmental floor. The position, cover zone of each cluster center surface, and the gradient in the face direction determined by the center surfaces are listed in Table 2. The gradient in the face direction of the cluster center represents the floor gradient in the cover zone of the cluster. It can be regarded as the stretching angle of the generating line determined by clustering based on morphology similarity.

6.1.3. Floor Reconstructing

In order to reconstruct the floor, the cover zones of the morphology-based cluster and conplanarity-based cluster should be coincident. But, comparing Table 1 and Table 2, it is found that the cover zones of the two kinds of clusters are different. The alignment of the two kinds of clusters was performed, as shown in Figure 12, so that the morphology and face-direction gradient are the same in each alignment cluster. After aligning the two kinds of clusters, the floor was divided into 11 zones. Table 3 lists the cover zone, position, and stretching angle of the alignment clusters.
According to the generating line and stretching angle listed in Table 3, the reconstructed floor is shown in Figure 13a, which can be regarded as the planned cutting path. The coal recovery is influenced by the difference between the reconstructed floor and original floor. The amount of cut gangue is a total volume of the reconstructed floor below the original floor. The amount of coal left is the total volume of the reconstructed floor overtop the original floor. By comparing the reconstructed and original floor, the amount of coal left is 1999.52 m3, and the amount of cut gangue is 1856.47 m3. The coal loss rate, losscoal, is calculated according to Equation (12), as below
l o s s c o a l = V c o a l - l e f t V c o a l - s e a m × 100 %
where Vcoal-left is the amount of coal left, and Vcoal-seam is the volume of the coal seam, which is the volume between the roof and floor of the coal seam model. The volume of the segmental coal seam located between 0 and 2100 m in the face-advancing direction is 351,190.58 m3. Therefore, the coal loss rate only is 0.57% when exploiting the coal seam according to the planned cutting path. The planned floor cutting path can be used as the extraction control target of the longwall shearer, so that the longwall shearer can cut along these defined floor profiles for optimal extraction and production.
The distribution of the gradient in the face-advancing direction is shown in Figure 13b, which is obtained by analyzing the morphology-based cluster centers using the PLR method of Ref. [15]. This gradient information can be used as the horizon control target of the longwall shearer, so that the path of the longwall shearer can follow the coal seam undulations in the direction of the longwall face’s advance.

6.2. Clustering of Whole Floor

The variation of the dendrogram and validity measure index during whole floor clustering based on morphology similarity and coplanarity are presented in Figure 14 and Figure 15, respectively. According to the validity measure index, the optimized morphology-based and coplanarity-based cluster numbers are selected as 29 and 10, respectively. The alignment cluster number of the two kinds of clusters is 53. The reconstructed floor produced with the 53 alignment clusters is exhibited in Figure 16. By comparing the reconstructed floor and original floor, the amount of coal left is 5642.35 m3, and the amount of cut gangue is 5463.76 m3. The volume of the whole coal seam is 648,330.77 m3. The coal loss rate of the whole coal seam is 0.87%. Comparing Figure 8a and Figure 14a, the morphology-based clustering dendrograms of the segmental and whole floor are different. In the case of the whole floor, the agglomeration level is very low, in the zone of 0–80 m, in the face direction. The number of clusters is 22 in the zone. In the case of the segmental floor, the number of clusters is 2 in the same zone. The reason for the difference is that there is a deflective fault at the location of 2200 m in the face-advancing direction. The deflective fault fades away in the zone of 0–80 m in the face direction. The deflective and fading characteristics of the fault cause the dissimilar morphology of each floor series samples.

7. Conclusions

Space adaptability between mining equipment and the coal-rock mass, to ensure the machines cut in the seam, is an important technology in longwall automation. In this study, a cutting path planning method based on hierarchical clustering was proposed according to the concept of the transparent longwall face, so that the mining machines were guided in the coal seam.
(1)
Two distance functions including morphology similarity and the coplanarity measurement method were defined to evaluate the similarity of clusters. The coal seam floor series in the face-advancing direction were clustered according to the morphology similarity. The surfaces composed of two neighboring floor series were clustered according to the coplanarity. The planned cutting path was obtained by reconstructing the coal seam floor, taking the morphology-based and coplanarity-based cluster centers as generating lines and stretching angle, respectively. Meanwhile, the distribution of the coal seam gradient in both the face-advancing direction and face direction was withdrawn, which can be used as the horizon control target of the longwall shearer.
(2)
The coal seam geological model of the 18,201 longwall face was analysed using the proposed clustering method. Comparing the planned cutting floor and original floor, the amount of coal left and cut gangue was 1999 m3 and 1856 m3, respectively, for the segmental floor. For the case of the whole floor, the amount of coal left and cut gangue was 5642 m3 and 5463 m3, respectively.
(3)
The proposed cutting path planning method aims to guide the longwall shearer cutting in the coal seam by considering the geometry of the coal seam. So, it is inapplicable to the low-thickness coal seam, in which undercut rock is unavoidable because the extracting thickness of longwall mining machines is larger than the thickness of the coal seam. A cutting path planning method for low-thickness coal seams can be selected for future work.

Author Contributions

Conceptualization, Z.G., S.W. and S.G.; methodology, Z.G. and S.W.; software, J.W.; validation, J.W.; formal analysis, S.W. and Z.G.; investigation, Z.G. and S.W.; resources, J.W. and S.W.; data curation, J.W.; writing—original draft preparation, Z.G.; writing—review and editing, S.W. and S.G. All authors have read and agreed to the published version of the manuscript.

Funding

This paper has been financially supported by the Chinese National Natural Science Foundation (No. 51874279), the Program of Ministry of Industry and Information Technology (TC220A04W-1,167), and the Shanxi Key Research and Development Program of China (2021020101003).

Data Availability Statement

The datasets generated and analyzed in the current study may be obtained from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Mitchell, G.W. Australasian Coal Mining Practice; Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2009. [Google Scholar]
  2. Peng, S. Longwall Mining; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  3. Huang, J.; Zhang, Y.; Zhu, L.; Wang, T. Numerical simulation of rock cutting in deep mining conditions. Int. J. Rock Mech. Min. Sci. 2016, 84, 80–86. [Google Scholar] [CrossRef]
  4. Kang, H.; Cho, J.; Park, J.; Jang, J.; Kim, J.; Kim, K.; Rostami, J.; Lee, J. A new linear cutting machine for assessing the rock-cutting performance of a pick cutter. Int. J. Rock Mech. Min. Sci. 2016, 88, 129–136. [Google Scholar] [CrossRef]
  5. Dai, X.W.; Huang, Z.W.; Shi, H.Z.; Wu, X.G.; Chao, X. Cutting force as an index to identify the ductile-brittle failure modes in rock cutting. Int. J. Rock Mech. Min. Sci. 2012, 146, 104834. [Google Scholar] [CrossRef]
  6. Li, X.F.; Wang, S.B. A theoretical method for estimating the torque of the shearer drum and the verification by full-scale shearer cutting experiments. Measurement 2022, 199, 111552. [Google Scholar] [CrossRef]
  7. Reid, D.C.; Hainsworth, D.W.; Ralston, J.C.; Mcphee, R.J. Shearer guidance: A major advance in longwall mining. Field Serv. Robot. Recent Adv. Res. Appl. 2006, 24, 469. [Google Scholar]
  8. Barham, D.K. Control and Automation of Underground Operations in the UK Mining Industry. IFAC Proc. Vol. 1985, 18, 139–143. [Google Scholar] [CrossRef]
  9. Kelly, M.; Hainsworth, D.; Reid, D. Longwall automation—A new approach. In Proceedings of the 3rd International Symposium “High Performance Mine Production”, CRISO Exploration & Mining, Aachen, Germany, 11–12 June 2003; pp. 5–16. [Google Scholar]
  10. Wang, S.B.; He, Y.; Wang, S.J.; Zhang, B.Y.; Ge, S.R. Study on the alignment method and experiment of scraper conveyor. J. China Coal Soc. 2017, 11, 3044–3050. (In Chinese) [Google Scholar]
  11. Kelly, M.; Hainsworth, D.; Reid, D. State of the art in longwall automation. In Proceedings of the 2005 Australian Mining Technology Conference “New Technologie”, Fremantle, Australia, 27–28 September 2005; The Australasian Institute of Mining and Metallurgy: Victoria, Australia, 2005. [Google Scholar]
  12. Peng, S.S.; Du, F.; Cheng, J.; Li, Y. Automation in U.S. longwall coal mining: A state-of-the-art review. Int. J. Min. Sci. Technol. 2019, 29, 151–159. [Google Scholar] [CrossRef]
  13. Li, S. Measurement & control and localisation for fully-mechanized working face alignment based on inertial navigation. Coal Sci. Technol. 2019, 8, 169–174. (In Chinese) [Google Scholar]
  14. Mundry, S.; Gajetzki, M.; Hoseinie, S.H. Longwall automation—Productivity and coal quality enhancement. Int. J. Min. Reclam. Environ. 2015, 29, 357–367.1. [Google Scholar]
  15. Reid, D.; Henderson, P.D.H.; Roberts, B.; Chacko, A.; Mcleod, C.R.M. Interconnection of Landmark Compliant Longwall Mining Equipment—Shearer Communication and Functional Specification for Enhanced Horizon Control; CSIRO Exploration and Mining: Floreat, Australia, 2005. [Google Scholar]
  16. Holm, M.; Beitler, S.; Arndt, T.; Mozar, A.; Junker, M.; Christian, B. Concept of shield-data-based horizon control for longwall coal mining automation. In Proceedings of the 16th IFAC Symposium on Automation in Mining, Mineral and Metal Processing August, San Diego, CA, USA, 25–28 August 2013. [Google Scholar]
  17. Wang, S.B.; Wang, S.J. Longwall mining automation horizon control: Coal seam gradient identification using piecewise linear fitting. Int. J. Min. Sci. Technol. 2022, 32, 821–829. [Google Scholar] [CrossRef]
  18. Beitler, S.; Holm, M.; Arndt, T.; Mozar, A.; Junker, M.; Bohn, C. A shield-data-based horizon control approach for thin seam coal mining utilizing plow technology. In Proceedings of the International Symposium on Automation and Robotics, Montréal, QC, Canada, 11–15 August 2013. [Google Scholar]
  19. Wang, G.F.; Xu, Y.X.; Ren, H.W. Intelligent and ecological coal mining as well as utilization technology in China: Review and prospects. Int. J. Min. Sci. Technol. 2019, 29, 161–169. [Google Scholar] [CrossRef]
  20. Yuan, L. Scientific conception of precision coal mining. J. China Coal Soc. 2017, 1, 1–7. (In Chinese) [Google Scholar]
  21. Wang, C.F.; Rong, Y. Concept, architecture and key technologies for transparent longwall face. Coal Sci. Technol. 2019, 7, 156–163. (In Chinese) [Google Scholar] [CrossRef]
  22. Yuan, L.; Zhang, P.S. Key technology and path thinking of dynamic reconstruction of mine transparent geological model. J. China Coal Soc. 2023, 48, 1–14. (In Chinese) [Google Scholar]
  23. Wang, G.F.; Ren, H.W.; Zhao, G.R.; Zhang, D.S.; Wen, Z.G.; Meng, L.Y.; Gong, S.X. Research and practice of intelligent coal mine technology systems in China. Int. J. Coal Sci. Technol. 2022, 9, 24. [Google Scholar] [CrossRef]
  24. Berry, M.W.; Mohamed, A.; Yap, B.W. Supervised and Unsupervised Learning for Data Science; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  25. Ruan, Y.; Chen, J.; Fan, Z.; Wang, T.; Mu, J.; Huo, R.; Huang, W.; Liu, W.; Li, Y.; Sun, Y. Application of K-PSO Clustering Algorithm and Game Theory in Rock Mass Quality Evaluation of Maji Hydropower Station. Appl. Sci. 2023, 13, 8467. [Google Scholar] [CrossRef]
  26. Shi, Z.; Zhao, H. Deep Multi-view clustering based on reconstructed self-expressive matrix. Appl. Sci. 2023, 13, 8791. [Google Scholar] [CrossRef]
  27. Ruan, Y.; Liu, W.; Wang, T.; Chen, J.; Zhou, X.; Sun, Y. Dominant Partitioning of Discontinuities of rock masses based on DBSCAN algorithm. Appl. Sci. 2023, 13, 8917. [Google Scholar] [CrossRef]
  28. Amiruzzaman, M.; Rahman, R.; Islam, M.R.; Nor, R.M. Evaluation of DBSCAN algorithm on different programming languages: An exploratory study. In Proceedings of the 2021 5th International Conference on Electrical Engineering and Information & Communication Technology, Dhaka, Bangladesh, 18–20 November 2021; Military Institute of Science and Technology: Dhaka, Bangladesh, 2021. [Google Scholar]
  29. Reddy, M.V.; Vivekananda, M.; Ruvn, S. Divisive hierarchical clustering with K-means and agglomerative hierarchical clustering. Int. J. Comput. Sci. Trends Technol. 2017, 5, 6–11. [Google Scholar]
  30. Nazari, Z.; Kang, D.; Asharif, M.R.; Sung, Y.; Ogawa, S. A New Hierarchical Clustering Algorithm. In Proceedings of the ICIIBMS 2015, Track2: Artificial Intelligence, Robotics, and Human-Computer Interaction, Okinawa, Japan, 28–30 November 2015; pp. 148–152. [Google Scholar]
  31. Madhulatha, T.S. An overview on clustering methods. IOSR J. Eng. 2012, 4, 719–725. [Google Scholar] [CrossRef]
  32. Zhong, L.; Yuan, J.S.; Zhang, W.H. Fuzzy C-mean algorithm with morophology similarity distance. In Proceedings of the 2009 Sixth International Conference on Fuzzy Systems and Knowledge Discovery, Tianjin, China, 14–16 August 2009; pp. 90–94. [Google Scholar]
  33. Kate, R.J. Using dynamic time warping distance as features for improved time series classification. Data Min. Knowl. Discov. 2016, 2, 283–312. [Google Scholar] [CrossRef]
  34. Dunn, J.C. A fuzzy relative of the ISODATA process and its use in detecting compact well-separated clusters. J. Cybern. 1974, 3, 32–57. [Google Scholar] [CrossRef]
  35. Dunn, J.C. Well separated clusters and optimal fuzzy partitions. J. Cybern. 1974, 1, 95–104. [Google Scholar] [CrossRef]
  36. Davies, D.L.; Bouldin, D.W. A Cluster Separation Measure. IEEE Trans. Pattern Anal. Mach. Intell. 1979, PAMI-1, 224–227. [Google Scholar] [CrossRef]
  37. Chou, C.H.; Su, M.C.; Lai, E. A new cluster validity measure and its application to image compression. Pattern Anal. Appl. 2004, 7, 205–220. [Google Scholar] [CrossRef]
  38. Smoliński, A.; Malashkevych, D.; Petlovanyi, M.; Rysbekov, K.; Lozynskyi, V.; Sai, K. Research into Impact of Leaving Waste Rocks in the Mined-Out Space on the Geomechanical State of the Rock Mass Surrounding the Longwall Face. Energies 2022, 15, 9522. [Google Scholar] [CrossRef]
  39. Malashkevych, D.; Petlovanyi, M.; Sai, K.; Zubko, S. Research into the coal quality with a new selective mining technology of the waste rock accumulation in the mined-out area. Min. Miner. Depos. 2022, 16, 103–114. [Google Scholar] [CrossRef]
Figure 1. Representation of a small portion of a longwall mining operation [7].
Figure 1. Representation of a small portion of a longwall mining operation [7].
Applsci 13 10242 g001
Figure 2. The coal seam geological model of the 18,201 longwall face (a), and difference distributions of the floor in face direction (b) and face-advancing direction (c).
Figure 2. The coal seam geological model of the 18,201 longwall face (a), and difference distributions of the floor in face direction (b) and face-advancing direction (c).
Applsci 13 10242 g002
Figure 3. Schematic of coal seam floor clustering.
Figure 3. Schematic of coal seam floor clustering.
Applsci 13 10242 g003
Figure 4. Schematic of morphology similarity measurement according to PLR of floor series.
Figure 4. Schematic of morphology similarity measurement according to PLR of floor series.
Applsci 13 10242 g004
Figure 5. Schematic of coplanarity measure between surface Mi and Mj.
Figure 5. Schematic of coplanarity measure between surface Mi and Mj.
Applsci 13 10242 g005
Figure 6. The flowchart of the agglomerative hierarchical clustering algorithm.
Figure 6. The flowchart of the agglomerative hierarchical clustering algorithm.
Applsci 13 10242 g006
Figure 7. The segmental floor located between 0 and 2100 m in face-advancing direction.
Figure 7. The segmental floor located between 0 and 2100 m in face-advancing direction.
Applsci 13 10242 g007
Figure 8. The segmental floor dendrogram (a) and variation of validity measure index (b) of clustering based on morphology similarity.
Figure 8. The segmental floor dendrogram (a) and variation of validity measure index (b) of clustering based on morphology similarity.
Applsci 13 10242 g008
Figure 9. The morphology of cluster centers based on morphology similarity for segmental floor.
Figure 9. The morphology of cluster centers based on morphology similarity for segmental floor.
Applsci 13 10242 g009
Figure 10. The segmental floor dendrogram (a) and variation of validity measure index (b) of clustering based on coplanarity.
Figure 10. The segmental floor dendrogram (a) and variation of validity measure index (b) of clustering based on coplanarity.
Applsci 13 10242 g010
Figure 11. The gradient surface of cluster centers based on coplanarity for segmental floor.
Figure 11. The gradient surface of cluster centers based on coplanarity for segmental floor.
Applsci 13 10242 g011
Figure 12. The alignment between morphology-based and coplanarity-based clusters.
Figure 12. The alignment between morphology-based and coplanarity-based clusters.
Applsci 13 10242 g012
Figure 13. The reconstructed floor (a) and the distribution of gradient in face and face-advancing directions (b).
Figure 13. The reconstructed floor (a) and the distribution of gradient in face and face-advancing directions (b).
Applsci 13 10242 g013aApplsci 13 10242 g013b
Figure 14. The whole floor dendrogram (a) and variation of validity measure index (b) of clustering based on morphology similarity.
Figure 14. The whole floor dendrogram (a) and variation of validity measure index (b) of clustering based on morphology similarity.
Applsci 13 10242 g014
Figure 15. The whole floor dendrogram (a) and variation of validity measure index (b) of clustering based on coplanarity.
Figure 15. The whole floor dendrogram (a) and variation of validity measure index (b) of clustering based on coplanarity.
Applsci 13 10242 g015
Figure 16. The reconstructed floor of the whole floor.
Figure 16. The reconstructed floor of the whole floor.
Applsci 13 10242 g016
Table 1. The position and cover zone of the cluster centers based on morphology similarity.
Table 1. The position and cover zone of the cluster centers based on morphology similarity.
No. of ClustersPosition of Cluster Center (X/m)Zone of Cluster (X/m)
124.80–45.6
265.646.4–80.8
3101.681.6–119.2
4136.8120–153.6
5197.6154.4–223.2
6238.4224–247.2
7267.2248–290.4
Table 2. The position, cover zone, and gradient in the face direction of the cluster center based on coplanarity.
Table 2. The position, cover zone, and gradient in the face direction of the cluster center based on coplanarity.
No. of ClustersPosition of Cluster Center Surface (X/m)Gradient in Face Direction of Cluster CenterZone of Cluster (X/m)
16.4~7.29.7516°0–12.8
230.4~31.29.7694°13.6–46.4
380~80.89.7251°47.2–111.2
4144.8~145.69.6532°112–180
5201.6~202.49.5967°180.8–231.2
6257.6~258.49.5149°232–290.4
Table 3. The cover zone, position, and stretching angle of the alignment clusters.
Table 3. The cover zone, position, and stretching angle of the alignment clusters.
No. of ClustersZone of Cluster (X/m)Position of Generating Line (X/m)Stretching Angle
10–12.824.89.7516°
213.6–46.49.7694°
347.2–80.865.69.7251°
481.6–111.2101.6
5112–119.29.6532°
6120–153.6136.8
7154.4–180197.6
8180.8–223.29.5967°
9224–231.2238.4
10232–247.29.5149°
11248–290.4267.2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guan, Z.; Wang, S.; Wang, J.; Ge, S. Longwall Face Automation: Coal Seam Floor Cutting Path Planning Based on Multiple Hierarchical Clustering. Appl. Sci. 2023, 13, 10242. https://doi.org/10.3390/app131810242

AMA Style

Guan Z, Wang S, Wang J, Ge S. Longwall Face Automation: Coal Seam Floor Cutting Path Planning Based on Multiple Hierarchical Clustering. Applied Sciences. 2023; 13(18):10242. https://doi.org/10.3390/app131810242

Chicago/Turabian Style

Guan, Zenglun, Shibo Wang, Jingqian Wang, and Shirong Ge. 2023. "Longwall Face Automation: Coal Seam Floor Cutting Path Planning Based on Multiple Hierarchical Clustering" Applied Sciences 13, no. 18: 10242. https://doi.org/10.3390/app131810242

APA Style

Guan, Z., Wang, S., Wang, J., & Ge, S. (2023). Longwall Face Automation: Coal Seam Floor Cutting Path Planning Based on Multiple Hierarchical Clustering. Applied Sciences, 13(18), 10242. https://doi.org/10.3390/app131810242

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop