Abstract

The extraction of underground mineral sources has a significant negative impact on the local environment, results in land surface subsidence. As far as subsidence monitoring technology is concerned, leveling is the most accurate. However, leveling can only obtain discrete point data but not the whole area information of the subsidence basin. In this study, Differential Interferometric Synthetic Aperture Radar (D-InSAR) combined with Unmanned Aerial Vehicle (UAV) technology is used to study the subsidence characteristics of the whole working panel. In this analysis, the Huainan mining area is tested as a research area; measured data are compared to the elevation accuracy of the Digital Surface Model (DSM) data, which can be considered for the subsequent works of the mining area. Based on the subsidence affected area, the ground object-type information is recorded to provide basic information for the ecological restoration work after mining so that the data before and after mining can be obtained synchronously. Finally, the differential interference results and Digital Orthophoto Map (DOM) data are combined to assess the spatiotemporal evolution of working panel subsidence and its influence on surface features. The main novelty of the proposed work is combining UAV and D-InSAR to get more accurate analysis of mining subsidence. It can be done using the proposed method.

1. Introduction

The depletion of land resources caused by coal mining is a well-known issue around the world. Surface buildings, water bodies, and roadways are all harmed when collapse pits, fissures, and subsidence steps develop. This has a significant impact on the residents of the mining area’s health and quality of life [1, 2]. Therefore, real-time and high-precision monitoring, prediction, and early warning of land surface disasters are critical for safe and efficient production as well as environmental protection in mining areas [35].

The conventional method of monitoring is to place observation stations in the form of discrete points along the working surface’s trend and inclination lines, establishing a monitoring mode with “points” as “lines” and “lines” as “polygons” [6, 7]. Despite its excellent precision, the limited monitoring scope, low efficiency, big workload, and expensive cost cannot be overlooked [810]. Simultaneously, this monitoring mode lowers the dependability of monitoring the entire mining face [11, 12]. The D-InSAR technique is used to synthesize elevation maps and get surface deformation differences over time, and it has the advantage of being available 24 hrs a day, on a large scale under any climatic conditions [1316]. It is widely employed in monitoring the subsidence evolution in mining areas such as coal mines [1720]. The subsidence of the land surface in the mining area, however, is large due to the radar wavelength limitation, and the land surface subsidence suppresses the monitoring threshold of the D-InSAR technique during the return visit time interval, generating phase unwrapping errors. Therefore, it is difficult to compute the land surface subsidence of two periods adequately [21, 22]. The InSAR method, on the contrary, mainly records the land surface elevation value, which does not reflect the properties and scope of surface objects such as agriculture, buildings, and highways. Moreover, the InSAR time-series approach’s potential is frequently limited by a lack of artificial objects and sparse vegetation cover, which makes acquiring deformation characteristics more difficult due to a lack of observation stations[23], particularly in the mine zone.

The UAV tilt photogrammetry technique can be well utilized to obtain the Digital Surface Model (DSM) and Digital Orthophoto Map (DOM) through optical photos imaging [24]. It is not limited by wavelength or other parameters, and it can monitor and observe the surface subsidence of the mining area. Mining subsidence monitoring has steadily been a research priority. For example, Ignjatović Stupar et al. [25] used the Velenje coal mine as a research area and collected observation data in 2017 by using a combination of the total station (TS) method, the Global Navigation Satellite System (GNSS), and Unmanned Aerial Vehicle (UAV) technologies. In this study, the height difference of reference points was evaluated using statistical analysis and test techniques of similar observations. Furthermore, utilizing the ultra-high-resolution normal perforation approach, Puniach et al. [26] provided a process for the automatic calculation of the horizontal displacement field in underground mining. The approach of using the UAV photogrammetry technique to monitor the dynamic land surface subsidence basin caused by underground coal mining and produce mining subsidence characteristics in a short time was presented by Dawei et al. [27].

The influence of terrain variations on image resolution, photo overlap degree, and monitoring accuracy was examined [28]. Then, they attached a TOF sensor to the UAV and set a programmed PID controller code to keep it flying at the same height above the ground, effectively realizing the UAV following monitoring. The UAV technique achieves elevation accuracy of roughly 50 mm in mining areas with difficult topography and around 30 mm in flat terrain areas. The accuracy of the elevation value obtained by the UAV technique is lower than that of the initial stage due to the small amount of subsidence, whereas the accuracy of the small deformation value obtained by the D-InSAR technique is less than 10 mm due to the small amount of subsidence at the initial stage [29, 30]. Based on the advantages of the above two technologies, the study proposes a four-step procedure as below:(1)Because land surface deformation is minor early in the mining process, D-InSAR technology is used to monitor the subsidence of the mining area, record the surface deformation information, and assess the extent of subsidence influence.(2)The RTK and UAV techniques are performed to obtain the surface deformation information, respectively. The accuracy of the data obtained by the two technologies are compared with the measured leveling data. It has been established that the elevation data obtained by the UAV technique can meet the needs of subsidence monitoring parameter inversion in the mining area.(3)Overall mining demand planning before and after mining, as well as the collection of surface object type information and statistics during subsidence monitoring, are all important data sources for ecological restoration of the mining subsidence area after mining.(4)Combining the data of UAV and D-InSAR techniques, the vertical deformation image obtained from SAR data is combined with DOM to reveal the spatiotemporal subsidence evolution of the mine area.

An enhanced phase optimization approach integrating UAV and D-InSAR is provided in response to the aforementioned difficulties. The research region is the subsidence area in the Huainan mining area in China, and this method is used to process SAR and UAV data covering the study area. The evolution of surface subsidence in the mining area is then indicated by obtaining deformation in the mining area. The final results of the synthesized study point to the possibility of identifying potential dangers and developing additional strategies to reduce the impact of mining regions on the environment and resources.

2. Materials and Methods

2.1. Study Area

The Banji coal mine is located at the intersection of Lixin County, Yingshang County, and Yingdong District, China. In the context of administrative division, it belongs to Huji Town, Lixin County, and Bozhou City. The minefield is situated 25 km north of Lixin County and 30 km south of Yingshang County. The study area extends between longitudes 116° 09E′ and 116° 30′ E and latitudes 32º51′ N 32° 56′ N, covering an area of 3358.84 km2 (Figure 1). The land cover/land use of the mine area is characterized by dry land, arbor forest land, other forest lands, mining land, rural homestead, railway land, highway land, rural roads, river water surface, pit and pond water surface, ditches, idle land, and agricultural lands.

The Banji mine’s capital construction stage runs from 2019 to 2020. On April 5, 2021, mining in the study area commenced. The Dongyi mining region, with a strike line length of 1070 meters, a dip line length of 270 meters, and recoverable deposits of 940000 tonnes, is the first planned mining area. According to the information, there are 36 boreholes in the research area of the Banji mine. The loose layer thickness in the monitoring area is 500 m to 648 m with an average of 600 m. The overlying bedrock thickness ranges from 30 m to 330 m, with an average of 130 m and the coal seam floor mining depth ranges from 646 m to 876 m, with an average of 786 m. The working area has the characteristics of a thick loose layer, significant mining depth, and thin bedrock in general. The mine has no operating panel mining history and no mining subsidence area to date.

2.2. UAV Data Acquisition and Processing

The Saier 102s five-lens camera on the Dajiang M300 four-rotor UAV can collect and shoot image information from five different directions: vertical, front view, rear view, left view, and right view. As a result, obtaining the essential data are more convenient. Table 1 lists the major characteristics of the UAV and camera.

The study area’s intervisibility and light intensity were assessed to be good during the field assessment. Roads, wheat-growing fields, buildings, trees, and other land objects on the land surface. Leveling and RTK observation stations are mainly distributed on the side of the path or ditch system. Individual observation points near the wheat field are excluded based on the DOM and field observation. The other points can be used to confirm the tilt photogrammetry’s elevation correctness.

Based on an S-shaped total flight path, with respect to the study area, a 2.0 km2 region was included. Before flying, the UAV must uniformly distribute the image control points based on the route’s range and the actual geomorphological features and then utilize RTK to collect the image control points’ 3D coordinates. At the same time, the RTK and leveling data of monitoring points are collected before and after the flight. The number of photographs collected was 17650 and 17695, respectively.

2.3. Methods
2.3.1. SAR Data Processing

The European Space Agency (ESA) provided the Sentinel-1A (S1A) products, which were processed to execute the following corrections and enhancements: multilooking and terrain adjustments, as well as calculating phase differences and D-InSAR data. To generate the deformation image, the SLC subswaths were divided and debursted independently before completing the coherence phase and vertical deformation using two different date images (e.g., May 11, 2021–May 23, 2021; May 23, 2021–June 4, 2021; June 16, 2021–June 28, 2021; June 21, 2021–July 3, 2021). Finally, the obtained images were multilooked and projected using a digital elevation model (DEM) acquired from the Shuttle Radar Topography Mission (SRTM) and a Universal Transverse Mercator (UTM) coordinates system of data from the World Geodetic System 1984 (WG84) to perform terrain correction.

To acquire the LOS deformation value and decompose it to obtain the vertical subsidence value for periodic monitoring, the double-track approach was employed. Because of the area’s unique location, there are three types of SAR data coverage, two of which can be employed to reduce SAR data monitoring time.

2.3.2. UAV Photo Acquisition Process

The most common methods of data processing are data import, aero triangulation, prickles, and coordinate system transformation, among others. Osgb format 3D models, DOM, and DSM are the products.

2.3.3. Data Accuracy Verification

The accuracy of UAV tilt photogrammetry is verified by combining Real-Time Kinematic (RTK) data and leveling data collected in the field, namely, plane coordinates (x, y), and elevation.

2.3.4. Land Surface Statistics and Spatiotemporal Law Analysis

Land conditions can be gathered using the collaborative monitoring mode of SAR data and the UAV tilt photogrammetry technique to provide basic information for postmining ecological restoration. The process of subsidence at the initial mining stage is investigated, and the spatiotemporal evolution operation of subsidence for mining is obtained. Figure 2 depicts the major experimental procedure.

3. Results

3.1. Model Building

DOM, DSM, and 3D model of mining area are generated (Figure 3). The obtained DOM has the characteristics of geometric map accuracy and image features, high accuracy, rich information, intuitive truth, and a short production cycle that can reflect the surface object information of a mining area truly and effectively. The DSM includes elevation information of surface buildings, trees, roads, fields, and other ground objects, which most truly reflects surface fluctuation features. Compared with DEM, the table features are more detailed [31].

Figure 3 shows the 3D model of the mining area. The actual landform features can be seen in the 3D model of the mining area, with an image size of 91175 × 63554 and a resolution ratio of 22 cm.

3.2. Plane Coordinate Accuracy Verification

Figure 4(a) depicts the horizontal and RTK monitoring points distribution. The black box represents the working panel position, the red triangle represents the monitoring points arranged along the dip line, and the yellow five-pointed star represents the observation points along the strike line, all of which conforms to the mining area’s subsidence monitoring point distribution criterion. One by one, examine the matching effects of leveling observation stations and points on the DOM graph. The schematic representation of the matching degree of observation stations randomly picked from the dip line and strike line is shown in Figures 4(b) and 4(c), demonstrating that the plane coordinate of DOM has a high-precision intuitively.

To verify the plane accuracy of UAV tilt photography, 7 and 19 verification points were set on the road along the strike line and the dip line of the working face, respectively, with an average interval of 80 meters between each verification point (Figure 4(d)). The verification points of the strike line were represented by a yellow five-pointed star, while the verification points of the dip line were represented by a red triangle. The root mean square errors (RMSEs) of x and y coordinates of the verification points are shown in Table 2. Because the road condition in the trend direction is poor, the accuracy of the verification point is slightly lower than that of the verification point. The road in the trending direction, on the other hand, is a cement road in good condition.

3.3. Elevation Accuracy Verification

The elevation values of leveling observation stations in the two directions of the strike line and dip line are extracted from the first phase and compared to RTK and leveling data. Figure 5 depicts the accuracy analysis of RTK data collected in the field and matching DSM land surface point elevation values. The[[parms resize(1),pos(50,50),size(200,200),bgcol(156)]] verified with leveling data. The precision of the dip line’s pointing accuracy is distributed within 5 cm, as shown in Figure 6(b). MS29, MS50, and MS56 have absolute errors of more than 10 cm, with MS56 having a 26 cm anomaly. The surrounding wheat seedlings may be close to the point, affecting the point’s elevation value, which can be removed based on the DOM image.

The accuracy of UAV tilt photogrammetry is verified by comparing DSM elevation and measured leveling value, as shown in Figures 6(c) and 6(d). The strike line elevation absolute errors are 0.04 cm and 8.62 cm, respectively. All other sites are within 5 cm except for MLA076, which is 8.62 cm, and MLA074, which is 6.78 cm. Absolute dip line elevation errors are 0.04 cm and 13.59 cm, respectively. MS29 has an absolute error of 13.59 cm, MS50 has an absolute error of 9.33 cm, MS40 has an absolute error of 5.77 cm, and MS46 has an absolute error of 5.66 cm.

After comparing DSM elevation and leveling values, preliminary verification shows that MLA015-MLA039 points have accepted accuracy.

As demonstrated in Table 3, it can be seen from the comparison of root mean square error among all data that the root mean square error of DSM, RTK, and leveling data is smaller on the trend line. On the strike line, the root mean square error between DSM elevation value and RTK data is considerable, while between DSM elevation value and leveling data, the root mean square error is the modest. As a result, it is possible that the surrounding environment has an impact on RTK data from the MLA015-MLA039 observation stations.

Overall, as compared to RTK data, the accuracy of DSM obtained by UAV tilt photogrammetry can approach within 5 cm in elevation, and its error is lower. The overall accuracy is fairly consistent, which is sufficient for mining area parameter inversion.

On the ground of the mining area, a real-time dynamic monitoring system with one reference station (BJCMP) and twelve real-time monitoring stations has been created. On the top of an office building in the industrial square, the reference station is positioned. Near the air, the shaft is a real-time monitoring station (MCORS01). At the intersection of the 110801 study area’s strike and dip observation lines, one real-time monitoring stations is set up (MCORS02). Figure 7 depicts the distribution of real-time monitoring stations.

The elevation data of MLA015 and MLA039 were collected in the first phase by connecting the Qianxun base station, which could be influenced by the surrounding environment and cause errors. In the second experiment, the RTK1 + 2 mode was used, with a reference station and two mobile data collection stations mounted on the mining area’s office roof. The reference station on the mining area’s office roof is shown by the red five-pointed star in Figure 7.

The error analysis diagram of the verification points on the road in the two strike and dip directions is taken as illustrated in Figure 8. The strike for rammed land roads is set at 7 points, and the dip for cement (asphaltic) roads is set at 19 points, with an average distance of 80 meters between the two points. In Figure 8(a), the elevation difference is mostly within 5 cm, with the exception of the JC9 point, which is 6.29 cm. In Figure 8(b), the difference in the JC20 point difference is 5.72 cm, while the other locations are all within 5 cm. The root mean square error of the cement road (dip) monitoring point is 2.15 cm, and the RMSEs of the tamped land road (strike) monitoring point is 3.71 cm, as shown in Table 4, essentially meeting the requirements of mining subsidence parameter inversion.

4. Discussion

4.1. Workface Influence Range and Land Category Statistics

The D-InSAR technique allowed monitoring of small deformations, and subsidence or collapse of the mining area due to minimal subsidence at the beginning of mining [32]. As a result, based on the surface deformation range and geological and mining conditions obtained from SAR data, the influence range of the mining research area was defined at 600 meters outward expansion on both sides of the strike direction and 400 meters outward expansion on both sides of the dip. Ground objects such as buildings, crops, roads, and water systems have been damaged.

Table 5 shows that 98 houses totaling 1.42 hectares were statistically affected. The total size of the farmland area is 186.03 hectares. The main road is County Road 046, which is 1.33 km long and covers 0.77 hectares. The water system consists mainly of fish ponds and ditches, which cover a total of 20.78 hectares. As demonstrated in Figure 5, the statistical results can provide basic information for postmining ecological restoration.

4.2. Temporal and Spatial Analysis

The land surface movement law of thick alluvium has been the subject of numerous research studies [33] with the majority of them focused on numerical simulation, similar material simulation, and observed strike and dip line data. The mobile communication has powerful data processing capabilities for accurately revealing land surface dynamic changes in a short time [3437]. In terms of numerical simulation and comparable material simulation, these two methods are subject to some uncertainty, and there will be a gap between them and reality. The profile line formed by the discrete data from strike and dip line observation stations is insufficient for analyzing the temporal and spatial evolution law of mining subsidence in the form of area. Conventional methods cannot actually and intuitively illustrate the influence range of surface subsidence and the temporal and spatial evolution law of mining subsidence [3842].

The differential interferogram during the two periods of Sentinel-1 SAR data is shown in the left side of Figure 9. The three-dimensional depiction of the subsidence influence range is displayed on the right side. It is depicted as expanding the subsidence to 25 times due to the modest subsidence effects. In Figure 9(d) the solid black line indicates the range of subsidence of less than 20 mm, and the solid yellow line shows the range of subsidence of less than 40 mm. On April 7, 2021, the working face begins mining, and on May 24, 2021, June 4, 2021, and June 23, 2021, the mining progress is 115.5 m, 157 m, and 232.5 m, respectively. This reveals that the D-InSAR technique has the capability to reveal the minimal and large amount of land surface deformation observation data [43, 44].

From May 11, 2021, to May 23, 2021, the central area does not reach County Road 046; during the period of May 23, 2021 to June 4, 2021, the main influence area is extended to County Road 046. In addition, during the period from June 16, 2021 to June 28, 2021 and from June 21, 2021 to July 3, 2021, the main influence area is extended to the east side of County Road 046. It can be observed that the mining subsidence area expands in the mining direction between two dates. In the early mining stage Figure 9(a), the image of land surface changes expands to reveal the minimal changes. In this image, the convergence at the boundary is very slow and reveals a long-term slight subsidence zone. The surface deformation value is small and displays the subsidence basin in an elliptical to circular form, with symmetry on both sides of the dip and asymmetry on both sides of the strike. The subsidence area reaches its mamima in the center, and the elliptical subsidence basin becomes nearly rounded. Figures 9(c) and 9(d) revealed that the, the main subsidence area is nearly rounded form, with symmetrical development on both sides of strike and dip that reveals subsidence positively related to mining extraction., The overall subsidence of the land surface Figures 9(a)9(c) is dominant, and the subsidence in the middle is not cleary apparent, but the solid yellow line area has large surface deformation showing the characteristics of a concentrated subsidence Figure 9(d).

Thus, its monitoring is essential for the safety and economics of coal mining [45, 46]. Therefore, risk assessment and management of coal mining using spatiotemporal mobile imaging systems is extremely important for sustainable development [47].

5. Conclusion

In this study, the DOM, DSM, and 3D model of the surface of the mining area working face panel are obtained based on the UAV tilt photogrammetry technique to better reveal the subsidence law in the mining area of Banji coal mine, China. The elevation values extracted by DSM, RTK elevation data, and leveling data are compared in pairs based on the evening process. The accuracy of RTK data is significantly different due to the influence of the surrounding environment. UAV tilt photogrammetry data provides a higher level of stability than RTK data, making it suitable for parameter inversion and subsidence prediction.

The processed SAR data derived from Sentinel-1 data provided information on the subsidence of the working panel on a regular basis; however, it cannot reflect the land surface characteristics and geomorphic features. The influence range and degree of the retreating progress on the surface features are analyzed on the unified scale of space and time by combining the SAR data with the DOM data from the mining area.

Integrate leveling data, RTK, SAR, DOM, DSM, and 3D model data generated by UAV tilt photogrammetry allowed us to reveal the subsidence law of the working panel in the whole basin on a unified scale of space and time. It has been realized that monitoring discrete points forming a linear form to a 3D + time scale allows providing a basis for the comprehensive study of the mining subsidence law.

Overall, this study verified the accuracy of DSM. In the future, Data will be upgraded in the future to allow for parameter inversion of the mining area, mining forecast of the working panel, and a more adequate database for postmining ecological restoration work.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project was supported by the Key Research and Development Program of Anhui Province (Grant no. 202104a07020014) and the Major Science and Technology Project of Anhui Province (Grant no. 202103a05020026).