A Study on Historical Big Data Analysis of Surface Ecological Damage in the Coal Mining Area of Lvliang City Based on Two Mining Modes
Abstract
:1. Introduction
2. Data and Methods
2.1. Experimental Scheme
2.2. Overview of the Study Area
2.2.1. Overview of Lvliang City
2.2.2. Overview of Liulin County
2.2.3. Overview of the Huajin Mining Area
2.3. A “Satellite–UAV–Ground–Underground” Four-in-One Method for Monitoring Ecological Damage on the Surface of Mining Areas
2.4. Methods for Analyzing the Evolution of Ecological Damage on the Surface of Mining Areas
2.4.1. Lvliang Coal Mine Information Acquisition
2.4.2. Remote Sensing Image Selection and Processing Analysis
2.4.3. Analysis of the Trend of Area Evolution of Coal Mining Subsidence Area
2.4.4. Evolutionary Analysis of Land Use Types in Coal Mining Areas
2.4.5. Evolutionary Analysis of Vegetation Indices in Coal Mining Areas
2.5. Analysis of Differences in Ecological Damage between Two Mining Models
- (1)
- Based on the NDVI of remote sensing images in each period obtained from the above calculations, explore and analyze the differences in the mean values of surface NDVI in the 121-mining method coal mine area, the 110-mining method coal mine area, and all the coal mine areas in the whole county of Liulin County. Taking the 2404 working face (110 mining method) and the 5201 working face (121 mining method) of Shaqu No. 1 mine in the Huajin mining area, which were simultaneously mined and emptied in 2019, as an example, we explored and analyzed the trend of the mean surface NDVI and its difference between the 2404 working face (110 mining method) and the 5201 working face (121 mining method) and the Huajin mining area.
- (2)
- Taking the 4502 working face (110 mining method) and the 4306 working face (121 mining method) in the 4# coal seam of Shaqu No. 1 mine in the Huajin mining area as an example, consider the following: 4502 working face (110 mining method) mining time, November 2020–February 2022; plumb depth from the ground about 290–460 m; backing length 1102–1092 m; inclined length 220 m; thickness of the coal seam 3.8–4.6 m; average thickness 4.2 m. For the 4306 working face (121 mining method): mining time June 2015–May 2017; thickness of coal seam 3.8–4.8 m; average thickness 4.2 m; inclination length 220.8 m; back mining length 1401 m. Based on the on-site research data and historical observation data, the surface cracks near the cutting eye, the dynamic cracks in the advancing direction of the working face, and the characteristics of subsidence damage and multiple parameters such as crack shape, crack width, crack depth, crack step misalignment, crack length, crack spacing, etc. (Table 3), were monitored to compare and analyze the differences in the characteristics of the surface subsidence damage of the working face of the two types of working face mining.
3. Results
3.1. Data on Changes in Coal Mining Subsidence Areas on the Lvliang City Scale
3.2. Data on Changes in Coal Mining Areas Surface Ecological Damage on the Liulin County Scale
3.3. Data on Changes in Coal Mining Subsidence Areas on the Huajin Mining Area Scale
4. Discussion
4.1. Trend and Prediction of Changes in Subsidence Area
- In the first stage, from 1949 to 1981, the coal mining subsidence area in Lvliang City grew from 11,766.67 ha to 21,566.65 ha, and the slope of the curve of its segmented fitting function was 326.69, with an actual growth rate of 306.25 ha/a.
- In the second stage, from 1981 to 2008, the subsidence area grew to 46,306.67 ha, and the slope of the curve of its segmented fitting function was 896.31, with an actual average growth rate of 916.29 ha/a, which is 1.99 times more than the average annual growth rate of the previous stage.
- In the third stage, from 2008 to 2023, the subsidence area grows to 92,191.47 ha, and the slope of the curve of its segmented fitting function is 3212.34, and the actual average growth rate is 3058.99 ha/a, which is 2.34 times more than the average annual growth rate of the previous stage.
- In the first stage, from 1949 to 1981, the coal mining subsidence area in Liulin County grew from 2626.67 ha to 4820.01 ha, and the slope of the curve of its segmented fitting function was 72.98, with an actual growth rate of 68.54 ha/a.
- In the second stage, from 1981 to 2008, the subsidence area grew to 10,346.67 ha, and the slope of the curve of its segmented fitting function was 200.22, with an actual average growth rate of 204.69 ha/a, which is 1.98 times more than the average annual growth rate of the previous stage.
- In the third stage, from 2008 to 2023, the subsidence area grows to 20,594.19 ha, with a segmented fit function curve slope of 717.59 and an actual average growth rate of 683.17 ha/a, which is 2.33 times more than the average annual growth rate of the previous stage.
- In the first stage, from 2003 to 2016, the area of coal mining subsidence in the Huajin mining area grew from 144.35 ha to 980.01 ha, and the slope of the curve of its segmented fitting function was 67.71, with an actual growth rate of 64.28 ha/a.
- In the second phase, from 2016 to 2023, the subsidence area grows to 2014.8 ha, with a segmented fit function curve slope of 148.44 and an actual average growth rate of 147.83 ha/a, which is 1.3 times more than the average annual growth rate of the previous stage.
4.2. Trends in Land-Use Types and NDVI Changes
- Firstly, population change is an important influencing factor for the change of both, according to statistics, the population of Lvliang City has been growing continuously in the past 30 years. With a net increase of about 600,000 people, the growing population will generate huge demand for land, and the two are bound to undergo significant changes with the population growth. Theoretically, the pressure of increasing population on the demand for food will drive an increase in the area of arable land, but the data in this paper show that population growth has not caused an increase in the area of arable land but rather reflected in the steady increase in the area of land used for construction in order to meet the housing needs of the new population, a large area of agricultural land such as arable land has been zoned for residential land.
- Secondly, economic development is one of the main influencing factors driving the changes of the two. According to statistics, the gross domestic product (GDP) of Lvliang City reached CNY 236.61 billion in 2023. The vigorous development of the energy industry as a pillar has increased economic income and also accelerated the construction of city infrastructures, but the construction and development encroached upon and destroyed a large amount of agricultural and forestry land, resulting in a significant increase in the structure of the Lvliang City land use and the NDVI index changes.
- Thirdly, government planning and urbanization is another major factor driving the changes of the two. The policy issued by the Lvliang municipal government to actively carry out rational land planning and remediation is an important reason influencing the continuous balanced and stable development of the two and plays a key role in regulating them. Meanwhile, the rapid development of urbanization drives changes in the two mainly in three main features: firstly, the transfer of rural population; secondly, the gathering of non-agricultural industries to cities and towns; and thirdly, the transfer of agricultural labor force to non-agricultural labor force. According to statistics, the urbanization rate of Lvliang City reached 56.1% in 2023, and the land for urban construction has been increasing over the past 30 years, while the land for agriculture and forestry has been decreasing.
4.3. Differences in Ecological Damage between the Two Mining Models
- The width of the surface crack near the location of open-off cut of the 4502 working face (110 mining method) is about 5 mm, and the depth of crack development is about 4.5 m. There is a regular distribution of cracks on the surface above along the advancing direction of the working face, with the cracks length smaller than the inclined length of the working face. The starting points are about 50 m away from both sides of the working face. The width of the cracks is about 5 cm, and the spacing of the cracks ranges from 28 to 35 m. The crack in the advancing direction produces a step misalignment with a misalignment of about 12 cm.
- The width of the surface crack near the location of open-off cut of the 4306 working face (121 mining method) is about 10 cm, and the depth of the crack profile reaches 12 m. A large number of cracks along the advancing direction of the working face are irregularly distributed and produce obvious step misalignments, the cracks are staggered and connected, with a length of about 9.5 m beyond the sides of the working face. The width of the cracks reaches about 30 cm, and crack spacing is between 12 and 16 m. The apparent step misalignment produced by the cracks in the advancing direction amounted to about 48 cm.
- The height of the collapse zone on the 121 working face is 7.7 m, while that of the 110 working face is 12.8 m, increasing by 66.2%;
- The height of the crack zone in the 121 working face is 63.3 m, while that of the 110 working face is 22.3 m, which is reduced by 64.8%;
- The surface settlement height of the 121 working face is 1.5 m, while that of the 110 working face is 0.45 m, which is reduced by 70%.
5. Conclusions
- (1)
- As of the end of 2023, the coal mining subsidence areas in Lvliang City, Liulin County, and Huajin Mining District were 92,191.47 ha, 20,594.19 ha, and 2014.8 ha, respectively. Based on the current growth trend, it is expected that the coal mining subsidence area in Lvliang City, Liulin County, and Huajin Mining District will reach 130,739.55 ha, 29,205.27 ha, and 3796.08 ha, respectively in 2035, which is an increase of up to 41.812%, 41.813%, and 88.41%, respectively, compared to the existing coal mining subsidence areas in 2023.
- (2)
- Taking all mining areas within Liulin County in 2023 as an example, the average surface NDVI value within the 121 mining method mining area is lower than the average value of the entire county’s mining area, while the 110 mining method mining area is higher than the average value of the entire county’s mining area. The surface NDVI of the working face for both 121 mining method and 110 mining method will be lower than the average level of the mining area due to the impact of mining, but the surface NDVI of the working face of the 110 mining method will be restored to the average level of the mining area more quickly than that of the 121 mining method with the passage of time.
- (3)
- There are obvious differences in the damage characteristics of the working face between the 121 mining method and the 110 mining method, and the width of the surface cracks, the amount of misalignment of the crack steps, the length of the cracks, and the distribution density of the cracks in the working face of the 110 mining method are smaller than those in the working face of the 121 mining method. 110 mining method can effectively reduce the ecological damage caused by coal mining subsidence, and its widespread application can effectively realize the ecological environmental protection of coal mine areas and green mining of coal resources and contribute to the high-quality development of the coal industry.
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Hu, Z.Q.; Li, Y.Y.; Li, G.S.; Han, J.Z.; Liu, S.G. Opportunities and challenges of land reclamation and ecological restoration in mining areas under carbon neutral target. Coal Sci. Technol. 2023, 51, 474–483. [Google Scholar] [CrossRef]
- Stala-Szlugaj, K.; Grudzinski, Z. World steam coal management. Gospod. Surowcami Miner. Miner. Resour. Manag. 2021, 37, 5–26. [Google Scholar] [CrossRef]
- Laurence, D. Establishing a sustainable mining operation: An overview. J. Clean. Prod. 2011, 19, 278–284. [Google Scholar] [CrossRef]
- Hu, Z.Q. Reflections on the High Quality Development of Mining Land Reclamation and Ecological Restoration during the 14th Five Year Plan Period. J. Intell. Mine 2021, 2, 29–32. [Google Scholar]
- Hu, Z.Q.; Xiao, W. Some thoughts on green development strategy of coal industry: From aspects of ecological restoration. Coal Sci. Technol. 2020, 48, 35–42. [Google Scholar] [CrossRef]
- Dang, X.H.; Xu, L.J.; Gao, Y.; Feng, Y.Y.; Gao, Y. Effects of Coal Mining Subsidence on Soil Evaporation and Vegetation Transpiration. J. Basic Sci. Eng. 2022, 30, 1389–1401. [Google Scholar] [CrossRef]
- Li, G.S.; Hu, Z.Q.; Li, P.Y.; Yuan, D.Z.; Feng, Z.J.; Wang, W.J.; Fu, Y.K. Innovation for sustainable mining: Integrated planning of underground coal mining and mine reclamation. J. Clean. Prod. 2022, 351, 131522. [Google Scholar] [CrossRef]
- Wang, C.Y.; Gao, Z.; Li, J.P. Study on Influence of Underground Mining on Surface Subsidence. Coal Technol. 2018, 37, 10–12. [Google Scholar] [CrossRef]
- Hendrychová, M.; Kabrna, M. An analysis of 200-year-long changes in a landscape affected by large-scale surface coal mining: History, present and future. Appl. Geogr. 2016, 74, 151–159. [Google Scholar] [CrossRef]
- Wang, Y.J. Mining Area Research Progress and Prospect on Ecological Disturbance Monitoring in mining area. Acta Geod. Cartogr. Sin. 2017, 46, 1705–1716. [Google Scholar]
- Dulias, R. Impact of mining subsidence on the relief of the Rybnik Plateau, Poland. Z. Fur Geomorphol. 2011, 55, 25–36. [Google Scholar] [CrossRef]
- Wódka, M.; Kamieniarz, S.; Wojciechowski, T.; Przylucka, M.; Perski, Z.; Sikora, R.; Karwacki, K.; Jureczka, J.; Nadlonek, W.; Krieger, W.; et al. Post-mining deformations in the area affected by the former “Siersza” Hard Coal Mine in Trzebinia (southern Poland). Geol. Q. 2024, 68, 3–22. [Google Scholar] [CrossRef]
- Zhang, C.Y.; Li, J.; Lei, S.G.; Yang, J.Z.; Yang, N. Progress and Prospect of the Quantitative Remote Sensing for Monitoring the Eco-environment in Mining Area. Met. Mine 2022, 3, 1–27. [Google Scholar] [CrossRef]
- Zhao, Y.L.; Ding, B.L.; He, T.T.; Xiao, W.; Ren, H. Monitoring of subsidence water body change based on Google Earth Engine. J. China Coal Soc. 2022, 47, 2745–2755. [Google Scholar] [CrossRef]
- Xiao, W.; Chen, J.L.; Zhao, Y.L.; Hu, Z.Q.; Lv, X.J.; Zhang, S. Identify maize chlorophyll impacted by coal mining subsidence in high groundwater table area based on UAV remote sensing. J. China Coal Soc. 2019, 44, 295–306. [Google Scholar] [CrossRef]
- Xiao, W.; Hu, Z.Q.; Zhang, J.Y.; Zhao, Y.L.; Yang, K. The status and prospect of UAV remote sensing in mine monitoring and land reclamation. China Min. Mag. 2017, 26, 71–78. [Google Scholar]
- Yang, K.; Zhao, Y.L.; Zhang, J.Y.; Cheng, C.; Zhao, P.P. Tree height extraction using high-resolution imagery acquired from an unmanned aerial vehicle (UAV). J. Beijing For. Univ. 2017, 39, 17–23. [Google Scholar] [CrossRef]
- Lei, S.G.; Bian, Z.F. Research progress on the environment impacts from underground coal mining in arid western area of China. Acta Ecol. Sin. 2014, 34, 2837–2843. [Google Scholar]
- Ma, D.; Zhao, X.M. Dynamie monitoring of land damage in Xishan coalfield of Shanxi province. Sci. Surv. Mapp. 2022, 47, 107–114. [Google Scholar] [CrossRef]
- Li, J.; Peng, S.P.; Zhang, C.Y.; Yang, F.; Sang, X. Quantitative remote sensing-based monitoring and evaluation of the ecological environment in mining areas: Technology framework and application. J. Min. Sci. Technol. 2022, 7, 9–25+88. [Google Scholar] [CrossRef]
- Rahmonov, O.; Krzysztofik, R.; Środek, D.; Smolarek-Lach, J. Vegetation-and environmental changes on non-reclaimed spoil heaps in Southern Poland. Biology 2020, 9, 164. [Google Scholar] [CrossRef]
- Popovych, V.; Voloshchyshyn, A. Features of temperature and humidity conditions of extinguishing waste heaps of coal mines in spring. News Natl. Acad. Sci. Repub. Kazakhstan Ser. Geol. Tech. Sci. 2019, 4, 230–237. [Google Scholar] [CrossRef]
- Abramowicz, A.; Rahmonov, O.; Chybiorz, R. Environmental Management and Landscape Transformation on Self-Heating Coal-Waste Dumps in the Upper Silesian Coal Basin. Land 2021, 10, 23. [Google Scholar] [CrossRef]
- Dulias, R. The Impact of Mining on the Landscape; Springer: Cham, Switzerland, 2016. [Google Scholar] [CrossRef]
- He, M.C.; Song, Z.Q.; Wang, A.; Yang, H.H.; Qi, H.G.; Guo, Z.G. Theory of longwall mining by using roof cuting short wall team and 110 method—The third mining science and technology reform. Coal Sci. Technol. Mag. 2017, 1, 1–9+13. [Google Scholar] [CrossRef]
- He, M.C. Theory and engineering practice for non-pillars mining with automagical entry formation and 110 mining method. J. Min. Saf. Eng. 2023, 40, 869–881. [Google Scholar] [CrossRef]
- Wang, Q.; He, M.C.; Wang, Y.C.; Liu, J.T.; Xue, H.J. Research progress and prospect of automatically formed roadway by roof cutting and pressure relief without pillars. J. Min. Saf. Eng. 2023, 40, 429–447. [Google Scholar] [CrossRef]
- The People’s Government of Lvliang City. Overview of Lvliang City. Available online: http://www.lvliang.gov.cn/zjll/llgl/llgk/ (accessed on 20 April 2024).
- The People’s Government of Lvliang City. Physical Geography. Available online: http://www.lvliang.gov.cn/zjll/llgl/zrdl/ (accessed on 20 April 2024).
- Lvliang Investigation Team of National Bureau of Statistics, Lvliang City Bureau of Statistics. 2018 National Economic and Social Development Statistics Bulletin of Lvliang City. Lvliang Daily, 16 April 2019; 002. [Google Scholar]
- The People’s Government of Lvliang City. Energy Production and Utilization Situation. Available online: http://www.lvliang.gov.cn/zjll/zyys/ny/ (accessed on 20 April 2024).
- The People’s Government of Liulin Connty. Geological Introduction. Available online: http://www.liulin.gov.cn/zjll/rwdl_67897/zrdl/202112/t20211217_1601770.shtml (accessed on 20 April 2024).
- Hu, Z.Q.; Wang, L.L.; Xu, X.L. A Method for Obtaining Ecological Damage Information of Mining Land Based on Multiple Data Sources. CN103049655B, 23 July 2014. [Google Scholar]
- Nádudvari, Á.; Abramowicz, A.; Fabiańska, M.; Misz-Kennan, M.; Ciesielczuk, J. Classification of fires in coal waste dumps based on Landsat, Aster thermal bands and thermal camera in Polish and Ukrainian mining regions. Int. J. Coal Sci. Technol. 2021, 8, 441–456. [Google Scholar] [CrossRef]
- Li, Q.Z.; Hu, Z.Q.; Zhang, F.; Song, D.Y.; Liang, Y.S.; Yu, Y. Multispectral Remote Sensing Monitoring of Soil Particle-Size Distribution in Arid and Semi-Arid Mining Areas in the Middle and Upper Reaches of the Yellow River Basin: A Case Study of Wuhai City, Inner Mongolia Autonomous Region. Remote Sens. 2023, 15, 2137. [Google Scholar] [CrossRef]
- Song, D.Y.; Hu, Z.Q.; Zeng, J.Y.; Sun, H. Influence of mining on vegetation in semi-arid areas of western China based on the coupling of above ground and below ground—A case study of Daliuta coalfield. Ecol. Indic. 2024, 161, 111964. [Google Scholar] [CrossRef]
- Demirel, N.; Düzgün, Ş.; Emil, M.K. Landuse change detection in a surface coal mine area using multi-temporal high-resolution satellite images. Int. J. Min. Reclam. Environ. 2011, 25, 342–349. [Google Scholar] [CrossRef]
- Zhang, F.; Hu, Z.Q.; Liang, Y.S.; Li, Q.Z. Evaluation of Surface Crack Development and Soil Damage Based on UAV Images of Coal Mining Areas. Land 2023, 12, 774. [Google Scholar] [CrossRef]
- Wang, Q.Y.; Zhao, Y.L.; Fang, S.D.; Yang, X.; Zhou, H.; Liu, J.F. Inversion of soil moisture in bare soil based on multi-spectral remote sensing. J. Min. Sci. Technol. 2020, 5, 608–615. [Google Scholar] [CrossRef]
- Chen, W.; Xu, Z.J.; Guo, Q. Estimation of soil organic matter by UAV hyperspectral remote sensing in coal mining areas. Trans. Chin. Soc. Agric. Eng. 2022, 38, 98–106. [Google Scholar]
- Zhao, Y.X.; Sun, B.; Liu, S.M.; Zhang, C.; He, X.; Xu, D.; Tang, W. Identification of mining induced ground fissures using UAV and infrared thermal imager: Temperature variation and fissure evolution. ISPRS J. Photogramm. Remote Sens. 2021, 180, 45–64. [Google Scholar] [CrossRef]
- Qi, L.Z.; An, S.K.; Tian, C.; Zhou, D.W.; Dong, X. Mining Subsidence Monitoring Method and Parameter Inversion Based on Unmanned Aerial Vehicle Lidar Technology. Sci. Technol. Eng. 2022, 22, 4752–4761. [Google Scholar]
- Cui, F.; Wang, R.; Chen, C.Q.; Wang, Y.L.; Bai, Y.; Yin, C.Y.; Cao, Y.F.; Cao, R.; Liao, L.X. Integrated advanced prediction of small structures by ground penetrating radar and downhole seismic with complex geological conditions in East Yunnan. Coal Sci. Technol. 2023, 51, 222–230. [Google Scholar] [CrossRef]
- Wang, X.J. Monitoring, Evolution and Self-Healing Characteristic of Land Damage Due to High Tension Coal Mining in Windy and Sandy Area; China University of Mining and Technology: Beijing, China, 2014. [Google Scholar]
- Ju, C.Y.; Fu, Y.K.; Chen, C.; Feng, Z.W.; Bai, M.B.; Hu, Z.Q. Influence of coal mining subsidence cracks on soil surface watercontent in Shennan Mining Area. Coal Sci. Technol. 2022, 50, 309–316. [Google Scholar] [CrossRef]
- Cui, X.M.; Deng, K.Z. Research review of predicting theory and method for coal mining subsidence. Coal Sci. Technol. 2017, 45, 160–169. [Google Scholar] [CrossRef]
- Han, J.Z.; Hu, Z.Q.; Wang, P.J.; Yan, Z.G.; Li, G.S.; Zhang, Y.H.; Zhou, T. Spatio-temporal evolution and optimization analysis of ecosystem service value-A case study of coal resource-based city group in Shandong, China. J. Clean. Prod. 2022, 363, 132602. [Google Scholar] [CrossRef]
Category | Detailed List |
---|---|
Overall information | 1. City-wide coal resource distribution map; |
2. List of mining enterprises in the city (including closed coal mines in the past 30 years); | |
3. Introduction to the coal production and mining situation in the city. | |
Data of each coal mine | 1. Basic data such as distribution information of various coal mines, coalfield reservoir information, recoverable reserves, and output value; |
2. Historical and current data on coal mining in various coal mines (mining scope, time, construction method, working face size, advancing speed, etc.); | |
3. Plan view of mining engineering in various coal mines, comparison map between upper and lower shafts, and comprehensive geological bar chart; | |
4. Relocation of coal mining villages in various coal mines; | |
5. Management and reclamation of coal mining collapse in various coal mines; | |
6. Mining subsidence data such as water conducting fracture zone, collapse zone height, and subsidence coefficient generated by various coal mines. |
Selection Principle | Satellite Data Type | Years of the Selected Data |
---|---|---|
1 phase every 2 years in the past 10 years | Landsat-8 | 2023, 2021, 2019, 2017, 2015, 2013 |
1 phase every 5 years in the past 20 years | Landsat-7 | 2008, 2003 |
The rest are phased once every 10 years | Landsat-5 | 1993 |
Number | Evaluation Indicators | Measurement Methods |
---|---|---|
1 | Crack length | GNSS (Global Navigation Satellite System) receiver |
2 | Crack width | Steel ruler |
3 | Crack depth | Tape measure |
4 | Crack spacing | GNSS receiver |
5 | Step displacement | Tape measure |
6 | Crack shape | GNSS receiver |
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. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Li, Q.; Hu, Z.; Zhang, F.; Guo, Y.; Liang, Y. A Study on Historical Big Data Analysis of Surface Ecological Damage in the Coal Mining Area of Lvliang City Based on Two Mining Modes. Land 2024, 13, 1411. https://doi.org/10.3390/land13091411
Li Q, Hu Z, Zhang F, Guo Y, Liang Y. A Study on Historical Big Data Analysis of Surface Ecological Damage in the Coal Mining Area of Lvliang City Based on Two Mining Modes. Land. 2024; 13(9):1411. https://doi.org/10.3390/land13091411
Chicago/Turabian StyleLi, Quanzhi, Zhenqi Hu, Fan Zhang, Yanwen Guo, and Yusheng Liang. 2024. "A Study on Historical Big Data Analysis of Surface Ecological Damage in the Coal Mining Area of Lvliang City Based on Two Mining Modes" Land 13, no. 9: 1411. https://doi.org/10.3390/land13091411
APA StyleLi, Q., Hu, Z., Zhang, F., Guo, Y., & Liang, Y. (2024). A Study on Historical Big Data Analysis of Surface Ecological Damage in the Coal Mining Area of Lvliang City Based on Two Mining Modes. Land, 13(9), 1411. https://doi.org/10.3390/land13091411