Quantitative Analysis and Evaluation of Coal Mine Geological Structures Based on Fractal Theory

: With the increasing depth of coal mining, the quantitative evaluation of the degree of geological structure development is becoming increasingly important for the control of mine water hazards in coal mining areas. Understanding the complexity of geological structure development can improve the safety and efﬁciency of coal production. At present, various evaluation indicators of the geological structure development cannot fully reﬂect the complexity of faults and folds, and the evaluation process is usually affected by subjective human factors. In this paper, the fractal dimension from fractal theory is used as the evaluation indicator to quantitatively analyze and evaluate the complexity of fault and fold structure in the mining area. To verify the evaluation results, the mathematical geology method is applied in an analysis of the trend surface of fault and fold networks. The results indicate that the fractal dimension can be applied for the quantitative analysis and evaluation of the complexity of fault and fold networks. In addition, the outcome of this work provides new insights into how to characterize the fault and fold structures of coal mining areas in northern China, and has some important implications to ensure the coal production safety. Gezhuang mine, features well-developed folds. A well-developed subsurface fold structure indicates a relatively high concentration of tectonic stress.


Introduction
Faults and folds are the main geological structures, and their existence in the coal mining area will seriously threaten the safety of coal mining [1][2][3]. Quantitative evaluation of mine fault and fold complexity is of great significance for the control of mine water hazards, the optimization of mining operations and the reduction of invalid drilling footage, which contributes to more safely and reasonably exploiting coal resources [3][4][5]. The geological structures are also the main factor controlling the karst development of coal mines in northern China. In order to prevent water inrush accidents caused by exposed faults and collapse columns, a comprehensive analysis on the complexity of fault and fold networks is necessary to determine the karst development and water-rich areas distribution in mining area [6][7][8][9]. To date, the evaluation indicators often used for characterizing fault and fold structures include the fold extension length, the fold direction, the fault throw, the average length of fault strike, and the uniformity of fault and fold distribution. However,

Study Area
The Kailuan mining area, located in Hebei Province, North China, is a typical waterrich mining area, and its geological structure and hydrogeological conditions are complex. The Ordovician limestone aquifers under the seam are not only water-rich with a strong flow but also exhibit a high hydraulic head pressure, and the prevention and control of water inrush from the Ordovician limestone are crucial for safe production in this mining area. At present, the main coal seams (the upper group of coals) are far from the Ordovician limestone in the Qian Jiaying and Fan Gezhuang mines, and the sum of the thicknesses of the separated and thin limestone layers exceeds 100 m. Normally, Ordovician groundwater is unlikely to cover a mining face with such a large vertical thickness, but if the mining encounters a complicated structure in the mining area, the structure and rock fractures can act as water channels, and the Ordovician groundwater may progressively encroach, abruptly increasing the risk of water inrush in the mine.
The structure in the study area has become complicated following Permian-Carboniferous coal formation. The area has experienced different stress fields in the Hercynian-Indochina period, the early and middle Yanshanian movement, the late Yanshanian period, the end of the Yanshanian-early Himalayan and the middle Himalayan periods, and the various multiperiod phases of compression and changes in the orientation of the stress field have produced the current structural assemblage [25]. We have found that the Lv Jiatuo anticline in the south, the Fan Gezhuang TaTuo syncline, the JinKou syncline area, and the northeastern region of the Qian Jiaying mine are high-compression areas with stress oriented Energies 2021, 14, 1925 3 of 14 in the SW and NE directions. The stress in the study area is relatively concentrated in these areas. A sketch map of the various structures in the research area is shown in Figure 1.
Energies 2021, 14, x FOR PEER REVIEW 3 of 15 stress oriented in the SW and NE directions. The stress in the study area is relatively concentrated in these areas. A sketch map of the various structures in the research area is shown in Figure 1.

Fractal Theory
In 1975, mathematician Benoit B. Mandelbrot for the first time proposed the fractal concept, which is derived from the Latin word fractus. Its original meaning is scored, irregular, or broken [26]. He published his first book, "Fractal: Form, Chance and Dimension," which marked the official introduction of fractal theory. At present, there is no clear physical definition of a fractal. In space, fractals exhibit the form of a non-integer dimension indicating the theoretical existence from a thought perspective, which, in a certain sense, is a general term for the shape and structure of self-similarity without assigning a characteristic length [27]. Fractal theory emphasizes that in regard to substances or phenomena in regional and local areas, the structure, function, and space in these areas contain particular self-similarity characteristics [28]. The principles of self-similarity and iteration are important tenets of fractal theory, which states that fractals remain invariant under general geometric changes. Fractals can be described by the fractal dimension, which reflects the basic characteristics [29]. The common fractal dimensions include the Hausdorff dimension , similar dimension , capacity dimension , information dimension , and box-counting dimension [30][31][32]. Suppose that a whole is divided into small figures of exactly the same size and shape, and the linearity of each small figure is times that of the original figure, then the Hausdoff dimension is: where ( ) represents the number of small figures contained in the whole . Suppose that a fractal object (the whole) is divided into ( , ) subsets (local units) of the same size, and each subset is similar to the original set with a similarity ratio , then the similar dimension of the fractal set defined as:

Fractal Theory
In 1975, mathematician Benoit B. Mandelbrot for the first time proposed the fractal concept, which is derived from the Latin word fractus. Its original meaning is scored, irregular, or broken [26]. He published his first book, "Fractal: Form, Chance and Dimension," which marked the official introduction of fractal theory. At present, there is no clear physical definition of a fractal. In space, fractals exhibit the form of a non-integer dimension indicating the theoretical existence from a thought perspective, which, in a certain sense, is a general term for the shape and structure of self-similarity without assigning a characteristic length [27]. Fractal theory emphasizes that in regard to substances or phenomena in regional and local areas, the structure, function, and space in these areas contain particular self-similarity characteristics [28]. The principles of self-similarity and iteration are important tenets of fractal theory, which states that fractals remain invariant under general geometric changes. Fractals can be described by the fractal dimension, which reflects the basic characteristics [29]. The common fractal dimensions include the Hausdorff dimension D H , similar dimension D s , capacity dimension D c , information dimension D i , and box-counting dimension D b [30][31][32].
Suppose that a whole U is divided into N small figures of exactly the same size and shape, and the linearity of each small figure is r times that of the original figure, then the Hausdoff dimension is: where N(r) represents the number of small figures contained in the whole U.
Suppose that a fractal object A (the whole) is divided into N(A, r) subsets (local units) of the same size, and each subset is similar to the original set with a similarity ratio r, then the similar dimension D s of the fractal set A defined as: Energies 2021, 14, 1925 4 of 14 Suppose that the graph is a finite set in n-dimensional Euclidean space R n , and a sphere with radius δ is used to fill the graph. If N(δ) is the minimum number of spheres, then the capacity dimension is: Consider the number of fractals set elements contained in each covering U in the Hausdoff dimension, and let P i represent the probability that the elements in the fractal set belong to the covering U, then the information dimension is: Suppose A ∈ H(R m ), where R m is the Euclidean space, covering A with a closed square box with side length 1/2 n . If Nn(A) represents the number of boxes contained in A, then the box-counting dimension is: As can be seen from the above, fractal dimension has different definitions and calculation methods. The calculation method of fractal dimension is selected according to the characteristics of fractal objects. Here, the fractal dimension of faults and folds are determined by changing the observation scale. We apply two-dimensional and orthogonal grids with the same size, shape, and grid length to cover the topographic map of structure in the mining area. Counting the number of grids containing any traces of faults and folds, we change the grid size from a to a i , and the number of new statistical grids is N(a i ). The calculation results reveal that a i and N(a i ) satisfy an exponential relationship as follows: If the statistical data of a i and N(a i ) are fitted in the double logarithmic coordinate system, the fitting graph is almost a straight line. The regression equation is: In Equations (6) and (7), a i is the size of the two-dimensional and orthogonal grid, C 1 and C 2 are constants, N(a i ) is the number of grids for the different mesh sizes containing fault and fold traces, and D 1 and D 2 are the fractal dimensions of the fault and fold structures, respectively.
The fault network is associated with a widely developed fault in the mining area. In recent years, research results have revealed that the structure of faults and folds with fractal characteristics involves a complex system [33][34][35]. The complexity of the fault and fold networks in the mining area is expressed by the degree of damage to geological block. It is difficult to quantitatively characterize the complexity of faults and folds by conventional mathematical and statistical methods. If we apply the fractal dimension from fractal theory, we can obtain more accurate results than with conventional statistical methods [3,20]. The fractal dimension is employed to quantitatively evaluate the fold and fault network complexity in the geological block to comprehensively represent various changes in fault and fold indexes, including the extension length, scale, and distribution. According to the determined fractal dimension, we plot contour diagrams of the fault and fold complexity levels in the mining area.

Determination of the Fold Fractal Dimension
The current data are based on the fault data of the mining engineering plan of coal seam No. 9 in the Fan Gezhuang mine (Figure 2), using the grid coverage method to delineate the structure of the mine and calculate the fractal dimension of each partitioned unit. We plot a contour map of the mine based on the calculation result and quantitatively evaluate the complexity of the geological structure in the Fan Gezhuang mine. The fold and fault fractal dimensions of the Fan Gezhuang mine are calculated as an example, and we describe the acquisition process of the fractal dimension of the faults and folds in the study area below.

Determination of the Fold Fractal Dimension
The current data are based on the fault data of the mining engineering plan of coal seam No. 9 in the Fan Gezhuang mine (Figure 2), using the grid coverage method to delineate the structure of the mine and calculate the fractal dimension of each partitioned unit. We plot a contour map of the mine based on the calculation result and quantitatively evaluate the complexity of the geological structure in the Fan Gezhuang mine. The fold and fault fractal dimensions of the Fan Gezhuang mine are calculated as an example, and we describe the acquisition process of the fractal dimension of the faults and folds in the study area below.

Calculation of the Fractal Dimension of Folds in the Fan Gezhuang Mine
First, the unit division is evaluated. Combined with the existing excavation engineering plan of the Fan Gezhuang mine, whereby the grid length is successively halved and divided into four grid levels, the dimensions of the various grids are 1000 m × 1000 m, 500 m × 500 m, 250 m × 250 m, and 125 m × 125 m, thus forming dense grids with an exponential increase in the number of grid blocks, as shown in Figure 2. Thereafter, we sum the number of different grid blocks containing fold traces. The grid classification method is shown in Figure 3. Hence, the complexity of the fault network in the area of the Fan Gezhuang mine can be quantified. m × 500 m, 250 m × 250 m, and 125 m × 125 m, thus forming dense grids with an exponential increase in the number of grid blocks, as shown in Figure 2. Thereafter, we sum the number of different grid blocks containing fold traces. The grid classification method is shown in Figure 3. Hence, the complexity of the fault network in the area of the Fan Gezhuang mine can be quantified.
Second, the number of grids containing a fold is counted, and the fractal dimension of the folds is calculated. All the above grids contain fold traces, and their individual fractal dimensions represent the value at the center point of the partitioned unit. Combined with the analysis results obtained via double logarithmic linear regression of the partitioned unit areas and the whole Fan Gezhuang area, in the Fan Gezhuang mine, the fractal dimension of the folds is 1.144 (Figure 4), and the fractal dimension of most units is higher than 0.90. This indicates that folds are well developed in the Fan Gezhuang mine.
The fractal dimension of each grid is mostly between 0 and 1.41. According to Table  1, the primary fractal dimension exhibits a wide range, and the fractal dimension of folds is between 0.97 and 1.23 for each partition. The fractal dimension is high, which indicates the structure of folds in the mine area is complicated.  Second, the number of grids containing a fold is counted, and the fractal dimension of the folds is calculated. All the above grids contain fold traces, and their individual fractal dimensions represent the value at the center point of the partitioned unit. Combined with the analysis results obtained via double logarithmic linear regression of the partitioned unit areas and the whole Fan Gezhuang area, in the Fan Gezhuang mine, the fractal dimension of the folds is 1.144 (Figure 4), and the fractal dimension of most units is higher than 0.90. This indicates that folds are well developed in the Fan Gezhuang mine.   The fractal dimension of each grid is mostly between 0 and 1.41. According to Table 1, the primary fractal dimension exhibits a wide range, and the fractal dimension of folds is between 0.97 and 1.23 for each partition. The fractal dimension is high, which indicates the structure of folds in the mine area is complicated. Fold development is complicated in the northwestern part of the TaTuo syncline and the northeastern and southern parts of the Bi Gezhuang syncline, where the fractal dimensions are 1.199 and 1.227, respectively. The fold development in the partitioned units is much higher than that in the other units. The fold fractal dimension of the other units is low or almost zero. The fractal dimension value is much higher than that in the area with little fold development. Fold development and complexity exhibit high heterogeneity in the mining area.

Calculation of the Fractal Dimension of Folds in the Qian Jiaying Mine
Based on the analysis of the fractal theory and the calculation method of the fractal dimension, the fold fractal dimension of the Qian Jiaying mine is 0.9391 (Figure 1). Most of the fractal dimensions of the partitioned units are less than 0.71. This explains why the fold distribution is fairly homogeneous and the fold complexity is relatively low. Only in grid numbers 95, 75, 11,25,36,48,24,92,37,27,50, and 73 is the fold fractal dimension higher than 1.00. This is consistent with the Qian Jiaying mine area being divided into the fold areas in the northeastern and southwestern regions and the middle of the monoclinic structure. The fold fractal dimension of the northeastern region is high, that of the middle region is low, and that of the southwestern region is high.

Calculation of the Fractal Dimension of Folds in the Lin Xi Mine
By calculating and comprehensively analyzing the fractal dimension, we concluded that the fold fractal dimension of the Lin Xi mine is 0.9826 (Figure 4c). The whole region is divided into 43 grids, and the fold fractal dimension is zero in most units. Only in grid numbers 6,7,26,32,33,34,35, and 36 is the fold fractal dimension higher than 1.00. This implies that the fold in the mine has a small distribution range and is unevenly distributed.

Calculation of the Fractal Dimension of Folds in the Lv Jiatuo Mine
By dividing and covering the geological structure map of the Lv Jiatuo mine, we determine that the fold fractal dimension of the Lv Jiatuo mine is 0.886 (Figure 4d). The folds in the Lv Jiatuo mine are unevenly distributed. Except for the Lv Jiatuo anticline, which has a large influence range, the rest of the folds only slightly influence the mine.
By denoting the mesh size as a i and the grid number of the fold as N(a i ), we determine that the correlation between the units in the study area is good. We entered all the coordinates and fractal dimensions of the center points of all units into MATLAB software to plot the fold fractal dimension in the study area using a computer.
As shown in Figure 5, all the fold structures in the Fan Gezhuang mine are larger than those in the Qian Jiaying mine. Only the northeastern Qian Jiaying mine area, which is located near the Fan Gezhuang mine, features well-developed folds. A well-developed subsurface fold structure indicates a relatively high concentration of tectonic stress.

Calculation of the Fractal Dimension of Faults in the Fan Gezhuang Mine
We establish grids to cover the geological structure chart and calculate the dimension of the fault, and the fractal dimension of the Fan Gezhuang mine fault i (Figure 6a), while the majority of the fractal dimensions of the partitioned units are than 0.8, namely, between 1 and 1.96. This indicates that many faults are unevenly uted in the Fan Gezhuang mine study area. The fractal dimensions in grid num

Calculation of the Fractal Dimension of Faults in the Fan Gezhuang Mine
We establish grids to cover the geological structure chart and calculate the fractal dimension of the fault, and the fractal dimension of the Fan Gezhuang mine fault is 1.5766 (Figure 6a), while the majority of the fractal dimensions of the partitioned units

Calculation of the Fractal Dimension of Faults in the Fan Gezhuang Mine
We establish grids to cover the geological structure chart and calculate the fractal dimension of the fault, and the fractal dimension of the Fan Gezhuang mine fault is 1.5766 (Figure 6a), while the majority of the fractal dimensions of the partitioned units are higher than 0.8, namely, between 1 and 1.96. This indicates that many faults are unevenly distributed in the Fan Gezhuang mine study area. The fractal dimensions in grid numbers 14, 26, 18, 10, and 23, which are concentrated in the Ta Tuo syncline area of the Fan Gezhuang mine, are higher than 1.90. Fault structures are relatively well developed in these areas, and most faults are normal faults.

Calculation of the Fractal Dimension of Faults in the Qian Jiaying Mine
Combining the results of double logarithmic linear regression of the grid data for each unit and the fault network analysis of the entire Qian Jiaying mine, we determine that the fault fractal dimension of the Qian Jiaying mine is 1.4544 (Figure 6b). Most of the fractal dimensions of the partitioned units are greater than 0.9, and the values are between 0.90 and 1.96, which reveals large differences. This implies that many faults are well developed and unevenly distributed in the Qian Jiaying mine. The fractal dimensions in grid numbers 22,39,21,38,29,30, and 69, which are concentrated in the northeastern Qian Jiaying mine, are higher than 1.9. This indicates that the fault structures are relatively well developed in this area, and the faults are extensional normal faults.

Calculation of the Fractal Dimension of Faults in the Lin Xi Mine
Using the same method, we calculate the fractal dimension of the Lin Xi mine fault to be 1.0319 (Figure 6c). The majority of the fractal dimensions of the partitioned units approach 0, and the values are between 0.80 and 1.50. Many faults occur in a certain part of the Fan Gezhuang mine, and their distribution is very uneven. We observe that the fault structures are well developed in the northeastern and southeastern corners of the Lin Xi mine, as shown in Figure 7.
Using the same method, we calculate the fractal dimension of the Lin Xi mine fault to be 1.0319 (Figure 6c). The majority of the fractal dimensions of the partitioned units approach 0, and the values are between 0.80 and 1.50. Many faults occur in a certain part of the Fan Gezhuang mine, and their distribution is very uneven. We observe that the fault structures are well developed in the northeastern and southeastern corners of the Lin Xi mine, as shown in Figure 7.

Calculation of the Fractal Dimension of Faults in the Lv Jiatuo Mine
The fractal dimension of faults in the Lv Jiatuo mine is calculated, and we obtain a fault fractal dimension of 0.9964, as shown in Figure 6d. Most partitioned units exhibit a fractal dimension higher than 0, namely, between 0.80 and 1.00. The fault structures are not widely distributed, and fault development is concentrated and complicated in the Lv Jiatuo mine.
According to the obtained fractal dimension data, we applied MATLAB software to generate a synthesized contour map of the fault and fold fractal dimensions in the Qian Jiaying and Fan Gezhuang mines, as shown in Figure 8. The figure demonstrates that the fault and fold fractal dimensions of the Ta Tuo syncline area and Jin Kou syncline are high and that the developed structure is complex.

Calculation of the Fractal Dimension of Faults in the Lv Jiatuo Mine
The fractal dimension of faults in the Lv Jiatuo mine is calculated, and we obtain a fault fractal dimension of 0.9964, as shown in Figure 6d. Most partitioned units exhibit a fractal dimension higher than 0, namely, between 0.80 and 1.00. The fault structures are not widely distributed, and fault development is concentrated and complicated in the Lv Jiatuo mine.
According to the obtained fractal dimension data, we applied MATLAB software to generate a synthesized contour map of the fault and fold fractal dimensions in the Qian Jiaying and Fan Gezhuang mines, as shown in Figure 8. The figure demonstrates that the fault and fold fractal dimensions of the Ta Tuo syncline area and Jin Kou syncline are high and that the developed structure is complex.

Analysis of the Mathematical Geology
In addition to using the fractal dimension, we can also apply the mathematical geology method to assess the development level of geological structures, thereby verifying the analytical fractal calculation results. Here, the trend surface analysis of mathematical geology method is conducted based on the distribution of the structures, which is to use simple power polynomials to fit complex curved surfaces to flatten and fill real curved surfaces. Selecting the No. 9 coal seam with an extremely well-developed fault as an ex-

Analysis of the Mathematical Geology
In addition to using the fractal dimension, we can also apply the mathematical geology method to assess the development level of geological structures, thereby verifying the analytical fractal calculation results. Here, the trend surface analysis of mathematical geology method is conducted based on the distribution of the structures, which is to use simple power polynomials to fit complex curved surfaces to flatten and fill real curved surfaces. Selecting the No. 9 coal seam with an extremely well-developed fault as an example, we describe the pattern of fault and fold development. Figure 9 shows a contour map of the floor of the No. 9 coal seam, and the distribution of the geological structure in the Qian Jiaying mine is depicted. The geological structures in the No. 9 coal seam are well developed, and the central shape is a monoclinal structure dipping to NW. The fold structure expands in the northeast and southwest of the mining area, and the fold is very well developed in the northeast of the mining area. The fold axis is oriented approximately NW, and the fault is very well developed in the Qian Jiaying mine.
Energies 2021, 14, x FOR PEER REVIEW 12 of 15 curvature of the bend. This phenomenon illustrates that the fold structure in the Qian Jiaying mine is a small fold and is not notable.  The cubic residual diagram of the No. 9 coal seam floor in the Qian Jiaying mine shows that the distribution of the residual contours in the northeast of the No. 9 coal seam floor is clear, and that there are three extremely thick residual contours. The geological structure of this area is extremely complex. The fold structure in the region is represented by the distribution of the zero line, but the extension direction of the zero line in the southwest of the fold region accurately reflects the position of the fold axes in the division of the structure. This is consistent with the distribution of the fractal dimensions of folds and faults in the northeastern Qian Jiaying mine. This approach could be applied to quantify various indexes, including the structural complexity, extension length, size, and distribution of the structure. The quantitative evaluation of the complexity of the faults and folds by the fractal dimension is assessed with the conventional mathematical statistics method, With the use of the mathematical geology method, we analyzed the trend surface corresponding to the floor of the No. 9 coal seam in the Qian Jiaying mine. The one-variable multiple linear equation is selected as the trend surface model to fit the spatial surface of the geological structure. The results show that the fitting degree of the initial trend surface analysis is higher than 0.64, and the quadratic fitting degree reaches 0.77, which is a high fitting level. Moreover, the cubic fitting degree of the trend surface is very high, namely, up to 0.91. Therefore, this research applies the gridding method with MATLAB software to analyze the results of the cubic trend surface and to obtain the cubic trend surface chart and cubic residual diagram of the No. 9 coal seam floor (Figures 9 and 10, respectively). Based on the cubic trend surface diagram of the No. 9 coal seam floor in the Qian Jiaying mine, we deduce that only the monoclinal structure in the middle mining area is consistent with the contour map of the coal floor, and the overall structure is anomalous in the NE direction. The trend surface in the northeastern and southwestern regions is not consistent with the contour map of the No. 9 coal seam floor, and the trend surface indicates a higher curvature of the bend. This phenomenon illustrates that the fold structure in the Qian Jiaying mine is a small fold and is not notable.  The cubic residual diagram of the No. 9 coal seam floor in the Qian Jiaying mine shows that the distribution of the residual contours in the northeast of the No. 9 coal seam floor is clear, and that there are three extremely thick residual contours. The geological structure of this area is extremely complex. The fold structure in the region is represented by the distribution of the zero line, but the extension direction of the zero line in the southwest of the fold region accurately reflects the position of the fold axes in the division of the structure. This is consistent with the distribution of the fractal dimensions of folds and faults in the northeastern Qian Jiaying mine. This approach could be applied to quantify various indexes, including the structural complexity, extension length, size, and distribution of the structure. The quantitative evaluation of the complexity of the faults and folds by the fractal dimension is assessed with the conventional mathematical statistics method, and the fractal dimension method yields better results.

Conclusions
In this work, the fractal theory is used to quantitatively evaluate the complexity of geological structures in the coal mine area. The fractal dimensions of folds and faults in four coal mines are calculated respectively and used as evaluation indexes to evaluate the development degree of faults and folds, and the mathematical geology analysis is conducted to verify the accuracy of the evaluation. The results of this method reflect certain The cubic residual diagram of the No. 9 coal seam floor in the Qian Jiaying mine shows that the distribution of the residual contours in the northeast of the No. 9 coal seam floor is clear, and that there are three extremely thick residual contours. The geological structure of this area is extremely complex. The fold structure in the region is represented by the distribution of the zero line, but the extension direction of the zero line in the southwest of the fold region accurately reflects the position of the fold axes in the division of the structure. This is consistent with the distribution of the fractal dimensions of folds and faults in the northeastern Qian Jiaying mine. This approach could be applied to quantify various indexes, including the structural complexity, extension length, size, and distribution of the structure. The quantitative evaluation of the complexity of the faults and folds by the fractal dimension is assessed with the conventional mathematical statistics method, and the fractal dimension method yields better results.

Conclusions
In this work, the fractal theory is used to quantitatively evaluate the complexity of geological structures in the coal mine area. The fractal dimensions of folds and faults in four coal mines are calculated respectively and used as evaluation indexes to evaluate the development degree of faults and folds, and the mathematical geology analysis is conducted to verify the accuracy of the evaluation. The results of this method reflect certain advantages, and are of great significance to coal production safety. The main conclusions are as follows: (1) The value of the fold fractal dimension of the Fan Gezhuang mine is 1.44, and this value is very high, which indicates the geological structures in the mine area is complicated. The fractal dimension of each partitioned unit is mainly between 0.97 and 1.23, and the obtained results on the structural complexity of the region reveal a high heterogeneity. The value of the fault fractal dimension is 1.5766, and the fractal dimension is between 1.00 and 1.96 in most partitioned units and is higher than 0.8 everywhere. Many uniformly distributed faults are concentrated in the Ta Tuo syncline area, and most faults are well-developed extensional normal faults. (2) The fold fractal dimension value in the Qian Jiaying mine is 0.9391. The fractal dimension of the majority of the partitioned units is less than 0.71, which shows the fractal dimension of the many uniformly distributed folds is high in the northeast, low in the middle, and high in the southwest. The fault fractal dimension is 1.4544, and the fractal dimension of most partitioned units is between 0.9 and 1.96, with most being higher than 0.9, which reveals a wide range. Many faults are uniformly distributed in the northeastern region, and most faults are well-developed extensional normal faults.