Thermo-hydro-mechanical (THM) coupled simulation of the land subsidence due to aquifer thermal energy storage (ATES) system in soft soils

Aquifer thermal energy storage(ATES)system has received attention for heating or cooling buildings.However,it is well known that land subsidence becomes a major environmental concern for ATES projects.Yet,the effect of temperature on land subsidence has received practically no attention in the past.This paper presents a thermo-hydro-mechanical(THM)coupled numerical study on an ATES system in Shanghai,China.Four water wells were installed for seasonal heating and cooling of an agriculture greenhouse.The target aquifer at a depth of 74e104.5 m consisted of alternating layers of sand and silty sand and was covered with clay.Groundwater level,temperature,and land subsidence data from 2015 to 2017 were collected using field monitoring instruments.Constrained by data,we constructed a field scale three-dimensional(3D)model using TOUGH(Transport of Unsaturated Groundwater and Heat)and FLAC3D(Fast Lagrangian Analysis of Continua)equipped with a thermo-elastoplastic constitutive model.The effectiveness of the numerical model was validated by field data.The model was used to reproduce groundwater flow,heat transfer,and mechanical responses in porous media over three years and capture the thermo-and pressure-induced land subsidence.The results show that the maximum thermoinduced land subsidence accounts for about 60%of the total subsidence.The thermo-induced subsidence is slightly greater in winter than that in summer,and more pronounced near the cold well area than the hot well area.This study provides some valuable guidelines for controlling land subsidence caused by ATES systems installed in soft soils.

Two stages of subsidence and its formation mechanisms in Mid-Late Triassic Ordos Basin,NW China

Based on a large number of newly added deep well data in recent years,the subsidence of the Ordos Basin in the Mid-Late Triassic is systematically studied,and it is proposed that the Ordos Basin experienced two important subsidence events during this depositional period.Through contrastive analysis of the two stages of tectonic subsidence,including stratigraphic characteristics,lithology combination,location of catchment area and sedimentary evolution,it is proposed that both of them are responses to the Indosinian Qinling tectonic activity on the edge of the craton basin.The early subsidence occurred in the Chang 10 Member was featured by high amplitude,large debris supply and fast deposition rate,with coarse debris filling and rapid subsidence accompanied by rapid accumulation,resulting in strata thickness increasing from northeast to southwest in wedge-shape.The subsidence center was located in Huanxian–Zhenyuan–Qingyang–Zhengning areas of southwestern basin with the strata thickness of 800–1300 m.The subsidence center deviating from the depocenter developed multiple catchment areas,until then,unified lake basin has not been formed yet.Under the combined action of subsidence and Carnian heavy rainfall event during the deposition period of Chang 7 Member,a large deep-water depression was formed with slow deposition rate,and the subsidence center coincided with the depocenter basically in the Mahuangshan–Huachi–Huangling areas.The deep-water sediments were 120–320 m thick in the subsidence center,characterized by fine grain.There are differences in the mechanism between the two stages of subsidence.The early one was the response to the northward subduction of the MianLüe Ocean and intense depression under compression in Qinling during Mid-Triassic.The later subsidence is controlled by the weak extensional tectonic environment of the post-collision stage during Late Triassic.

Residual subsidence time series model in mountain area caused by underground mining based on GNSS online monitoring

The residual subsidence caused by underground mining in mountain area has a long subsidence duration time and great potential harm,which seriously threatens the safety of people's production and life in the mining area.Therefore,it is necessary to use appropriate monitoring methods and mathematical models to effectively monitor and predict the residual subsidence caused by underground mining.Compared with traditional level survey and InSAR(Interferometric Synthetic Aperture Radar)technology,GNSS(Global Navigation Satellite System)online monitoring technology has the advantages of long-term monitoring,high precision and more flexible monitoring methods.The empirical equation method of residual subsidence in mining subsidence is effectively combined with the rock creep equation,which can not only describe the residual subsidence process from the mechanism,but also predict the residual subsidence.Therefore,based on GNSS online monitoring technology,combined with the mining subsidence model of mountain area and adding the correlation coefficient of the compaction degree of caving broken rock and the Kelvin model of rock mechanics,this paper constructs the residual subsidence time series model of arbitrary point on the ground in mountain area.Through the example,the predicted results of the model in the inversion parameter phase and the dynamic prediction phase are compared with the measured data sequence.The results show that the model can carry out effective numerical calculation according to the GNSS monitoring data of any point on the ground,and the model prediction effect is good,which provides a new method for the prediction of residual subsidence in mountain mining.

Evolution mechanism and treatment timing of penetrating fissures

The Inner Mongolia mining area in western China are characterized by the development of numerous penetrating fissures,resulting in severe land damage.It is significant to reveal the underlying evolution mechanism and identify treatment timing for restoring the ecological environment.The Guanbanwusu mining subsidence area in Inner Mongolia,China was selected as the research case for this study.The evolution mechanism of different penetrating fissures was revealed by field measurement,physical simulation and theoretical analysis.The treatment timing prediction model for the mining subsidence area was established based on the enhanced Weibull time function.The results show that the ground fissures are mainly step-type and collapse-type fissures.The breaking form of overlying strata determines their vertical opening and horizontal dislocation.The high mining intensity in the western mining area results in a shortened period of dynamic fissure expansion and reduced closure degree.The damage extent of the overlying strata exhibits zoning characteristics both vertically and horizontally.The relative standard deviation of the prediction model is only 3.7%.Concurrently,the prediction model is employed to determine the optimal timing for treatment in the study area,estimated to be 259 days.Subsequently,once this threshold is reached,the study area undergoes treatment and restoration of its e cological environment.This study addresses the knowledge gap in this field by highlighting the interconnectedness between rock strata structure and evolution mechanism of penetrating fissures,thereby providing a method for determining the treatment timing in mining subsidence areas.

Using neural network modeling to improve the detection accuracy of land subsidence due to groundwater withdrawal

Despite the high efficiency of remote sensing methods for rapid and large-scale detection of subsidence phenomena,this technique has limitations such as atmospheric impact and temporal and spatial decorrelation that affect the accuracy of the results.This paper proposes a method based on an artificial neural network to improve the results of monitoring land subsidence due to groundwater overexploitation by radar interferometry in the Aliabad plain(Central Iran).In this regard,vertical ground deformations were monitored over 18 months using the Sentinel-1A SAR images.To model the land subsidence by a multilayer perceptron(MLP)artificial neural network,four parameters,including groundwater level,alluvial thickness,elastic modulus,and transmissivity have been applied.The model's generalizability was assessed using data derived for 144 days.According to the results,the neural network estimates the land subsidence at each ground point with an accuracy of 6.8 mm.A comparison between the predicted and actual values indicated a significant agreement.The MLP model can be used to improve the results of subsidence detection in the study area or other areas with similar characteristics.

InSAR-derived surface deformation characteristics and mining subsidence parameters in mountain coal mines

Mining-induced surface deformation disrupts ecological balance and impedes economic progress.This study employs SBAS-InSAR with 107-view of ascending and descending SAR data from Sentinel-1,spanning February 2017 to September 2020,to monitor surface deformation in the Fa’er Coal Mine,Guizhou Province.Analysis on the surface deformation time series reveals the relationship between underground mining and surface shifts.Considering geological conditions,mining activities,duration,and ranges,the study determines surface movement parameters for the coal mine.It asserts that mining depth significantly influences surface movement parameters in mountainous mining areas.Increasing mining depth elevates the strike movement angle on the deeper side of the burial depth by 22.84°,while decreasing by 7.74°on the shallower side.Uphill movement angles decrease by 4.06°,while downhill movement angles increase by 15.71°.This emphasizes the technology's suitability for local mining design,which lays the groundwork for resource development,disaster prevention,and ecological protection in analogous contexts.

A semi-infinite beam theoretical model on predicting rock slope subsidence induced by underground mining

When the mining goaf is close to the cliff,rock slope subsidence induced by underground mining is significantly affected by its boundary conditions.In this study,an analytical method is proposed by considering the key strata as a semi-infinite Euler-Bernoulli beam rested on a Winkler foundation with a local subsidence area.The analytical solutions of deflection are derived by analyzing the boundary and continuity conditions of the cliff.Then,the analytical solutions are verified by the results from experimental tests,FEM and InSAR,respectively.After that,the influence of changing parameters on deflections is studied with sensitivity analysis.The results show that the distance between goaf and cliff significantly affects the deflection of semi-infinite beam.The response of semi-infinite beam is obviously determined by the length of goaf and the bending stiffness of beam.The comparisons between semi-infinite beam and infinite beam illustrate the ascendancy of the improved model in such problems.

Risk assessment of coastal flooding disaster by storm surge based on Elevation-Area method and hydrodynamic models:Taking Bohai Bay as an example

The future inundation by storm surge on coastal areas are currently ill-defined.With increasing global sealevel due to climate change,the coastal flooding by storm surge is more and more frequently,especially in coastal lowland with land subsidence.Therefore,the risk assessment of such inundation for these areas is of great significance for the sustainable socio-economic development.In this paper,the authors use Elevation-Area method and Regional Ocean Model System(ROMS)model to assess the risk of the inundation of Bohai Bay by storm surge.The simulation results of Elevation-Area method show that either a 50-year or 100-year storm surge can inundate coastal areas exceeding 8000 km^(2);the numerical simulation results based on hydrodynamics,considering ground friction and duration of the storm surge high water,show that a 50-year or 100-year storm surge can only inundate an area of over 2000 km^(2),which is far less than 8000 km^(2);while,when taking into account the land subsidence and sea level rise,the very inundation range will rapidly increase by 2050 and 2100.The storm surge will greatly impact the coastal area within about 10-30 km of the Bohai Bay,in where almost all major coastal projects are located.The prompt response to flood disaster due to storm surge is urgently needed,for which five suggestions have been proposed based on the geological background of Bohai Bay.This study may offer insight into the development of the response and adaptive plans for flooding disasters caused by storm surge.

Mechanical responses of underground carriageway structures due to construction of metro tunnels beneath the existing structure:A case study

To understand the mechanical response pattern of the existing structure and ground due to the construction of metro tunnels underneath,the finite difference method is adopted to study the torsional deformation and stress variation of the existing structure and the effect of underground carriageway structures on the surface subsidence.The curves of the maximum differential subsidence,torsion angle,and distortion of the cross-section of the existing structure show two peaks in succession during traversing of two metro tunnels beneath it.The torsion angle of the existing structure changes when the two tunnels traverse beneath it in opposite directions.The first traversing of the shield tunnel mainly induces the magnitude variation in torsional deformation of the existing structure,but the second traversing of the subsurface tunnel may cause a dynamic change in the magnitude and form of torsional deformation in the existing structure.The shielding effect can reduce the surface subsidence caused by metro tunnel excavation to a certain extent,and the development trend of subsidence becomes slower as the excavation continues.

Identified the hydrochemical and the sulfur cycle process in subsidence area of Pingyu mining area using multi-isotopes combined with hydrochemistry methods

Groundwater serves as an important water source for residents in and around mining areas.To achieve scientific planning and efficient utilization of water resources in mining areas,it is essential to figure out the chemical formation process and the ground water sulfur cycle that transpire after the coal mining activities.Based on studies of hydrochemistry and D,^(18)O-H_(2)O,^(34)S-SO_(4)isotopes,this study applied principal component analysis,ion ratio and other methods in its attempts to reveal the hydrogeochemical action and sulfur cycle in the subsidence area of Pingyu mining area.The study discovered that,in the studied area,precipitation provides the major supply of groundwater and the main water chemistry effects are dominated by oxidation dissolution of sulfide minerals as well as the dissolution of carbonate and silicate rocks.The sulfate in groundwater primarily originates from oxidation and dissolution of sulfide minerals in coal-bearing strata and human activities.The mixed sulfate formed by the oxidation of sulfide minerals and by human activities continuously recharges the groundwater,promoting the dissolution of carbonate rock and silicate rock in the process.

Study on the Safety and Prevention Technology of Coal Mining under the River in Xingyuan Coal Mine

Coal mining-induced surface subsidence poses significant ecological and infrastructural challenges, necessitating a comprehensive study to ensure safe mining practices, particularly in underwater conditions. This project aims to address the extensive impact of coal mining on the environment, infrastructure, and overall safety, focusing on the Shigong River area above the working face. The study employs qualitative and quantitative analyses, along with on-site engineering measurements, to gather data on crucial parameters such as coal seam characteristics, roof rock lithology, thickness, water resistance, and structural damage degree. The research encompasses a multidisciplinary approach, involving mining, geology, hydrogeology, geophysical exploration, rock mechanics, mine surveying, and computational mathematics. The importance of effective safety measures and prevention techniques is emphasized, laying the foundation for research focused on the Xingyun coal mine. The brief concludes by highlighting the potential economic and social benefits of this project and its contribution to valuable experience for future subsea coal mining.

Method of Monitoring Three-dimensional Mining Surface Deformation Based on InSAR

In order to solve the problems of small monitoring range,long time and high cost of existing sedimentation observation methods,based on two-view sentinel No.1 radar images of Guqiao mining area in Huainan City from November 4,2017 to November 28,2017,surface change information was obtained in combination with D-InSAR,and the three-dimensional surface deformation was monitored by two-pass method and single line of sight D-InSAR method.The results show that during the research period of 24 d,the maximum deformation of the mining area reached 71 mm,and the southern subsidence was the most obvious,which was in line with the mining subsidence law.The maximum displacement from the north to the south was about 250 mm,while the maximum displacement from the east to the west was about 80 mm,and the maximum subsidence in the center was 110 mm.It is concluded that D-InSAR technique has a good effect on the inversion of the mining subsidence,and this method is suitable for three-dimensional surface monitoring in areas with similar geological conditions.The monitoring results have certain reference value.

Ganges and Brahmaputra River System: Need for a Multi-Country Plan for Cleaning, Mitigation, Restoration and Protection

The Ganges and Brahmaputra River system is in the plains of the northern Indian subcontinent. The river is a wide sluggish stream flowing through densely populated and fertile agricultural regions of the world. The Ganges is known as the Hinduism holy river. In Bangladesh, the Brahmaputra is joined by the Teesta River. The western branch of the Brahmaputra confluences with the Ganges and contains most of the river flow. The eastern branch joins the Meghna River near Dhaka. The basin covers parts of four countries including India, Nepal, China, and Bangladesh. Of greater concern, however, has been the degradation in quality of the river water itself. The primary objective of this research is to encourage the development of a multi-country clean-up, mitigation, and protection plan for the Ganges-Brahmaputra rivers. This article constitutes a real tool for the restoration, enhancement and protection of the Ganges-Brahmaputra River system and its environment. The Ganges and Brahmaputra rivers are known for stream bank erosion, shifting channels, and sandbars that continually emerge in their course. The Ganges and Brahmaputra watershed is home to hundreds of millions of people, with the result that the river’s water over much of its course is highly polluted. Arsenic contamination of groundwater in Bangladesh continues to be the largest case of human poisoning in history. Catastrophic floods have prompted the World Bank to prepare a long-term flood-control plan for the region. Scores of cities and towns contribute to treated sewage into the river and its main tributaries, and dozens of manufacturing facilities contribute industrial waste. Also contributing to high pollution levels are agricultural runoff, the remnants of partially burned or unburned bodies from funeral pyres, and animal carcasses. High levels of disease-causing bacteria, as well as such toxic substances as chromium, cadmium, and arsenic, have been found in the Ganges and Brahmaputra. External research and funding of adsorptive media systems to help mitigate the high arsenic levels in drinking water (river and groundwater) is needed. The Ganges-Brahmaputra River system is of colossal importance to its entire environment. Restoration and protection measures must be adopted appropriately and at the scale of the concerned countries.

Structural Style and Hydrocarbon Prospectivity of the Growth Faults Related Structures in the Bengal Basin

The structural and tectonic evolution of the Bengal Basin is characterized by a complex interplay of factors, including sedimentation, the rise of the Himalayan Mountains, and the movements of Jurassic syn-rift faults. This study aims to comprehend the progression of growth faults inside the basin by examining fault geometry, basin development, and structural relief patterns. We used high-quality 2D seismic lines from the PK-MY-8403, classical seismic interpretation techniques and modeling were carried out to reveal the plate tectonic conditions, stratigraphy, and sedimentation history of the basin. The break-up unconformity, Paleocene and Eocene submerged conditions, and crucial geological formations including the Sylhet Limestone, Barail Group, and Surma Group were among the notable features recognized in seismic section. With an emphasis on growth strata and pre-growth strata, significant variations in layer thickness and relief were remarked in different stratigraphic levels. Basin development events like the evolution of the Miocene remnant ocean basin, sedimentation in Oligocene, Eocene Himalayan collision, and the Pliocene reverse fault development are analyzed. In the early the Pliocene compressional forces outpaced sedimentation rates and syn-depositional normal faults of Oligocene time began to move in opposite direction. Syn-depositional growth faults may have formed in the Bengal Basin as a result of this reversal. This research provides a detailed comprehensive knowledge of growth fault development in the Bengal Basin following the seismic interpretation, modelling, and thickness/relief analysis. The outcomes point to a substantial hydrocarbon potential, especially in regions like the Eocene Hinge Zone, where the prospectivity of the area is enhanced by carbonate reefs and Jalangi shale. However, the existence of petroleum four-way closure in the investigated region requires further investigation.

A review of monitoring,calculation,and simulation methods for ground subsidence induced by coal mining

Subsidence data acquisition methods are crucial to mining subsidence research and an essential component of achieving the goal of environmentally friendly coal mining.The origin and history of the existing methods of field monitoring,calcula-tion,and simulation were introduced.It summarized and analyzed the main applications,flaws and solutions,and improve-ments of these methods.Based on this analysis,the future developing directions of subsidence data acquisition methods were prospected and suggested.The subsidence monitoring methods have evolved from conventional ground monitoring to combined methods involving ground-based,space-based,and air-based measurements.While the conventional methods are mature in technology and reliable in accuracy,emerging remote sensing technologies have obvious advantages in terms of reducing field workload and increasing data coverage.However,these remote sensing methods require further technological development to be more suitable for monitoring mining subsidence.The existing subsidence calculation methods have been applied to various geological and mining conditions,and many improvements have already been made.In the future,more attention should be paid to unifying the studies of calculation methods and mechanical principles.The simulation methods are quite dependent on the similarity of the model to the site conditions and are generally used as an auxiliary data source for subsidence studies.The cross-disciplinary studies between subsidence data acquisition methods and other technologies should be given serious consideration,as they can be expected to lead to breakthroughs in areas such as theories,devices,software,and other aspects.

Responses of breeding waterbird communities to environmental changes in subsidence wetlands in the North China Plain

In the context of global degradation and loss of natural wetlands,waterbirds have been increasingly using artificial wetlands as alternative habitats.However,waterbirds are facing various threats in these artificial wetlands,due to dramatic environmental changes induced by anthropogenic activities.Exploring the effects of these changes on the temporal dynamics of the waterbird communities can help understand how waterbirds adapt to environmental changes and thus formulate effective management and conservation plans.In this study,we carried out field surveys on waterbirds and environmental factors across 20 subsidence wetlands created by underground coal mining in the Huainan coal mining area in the breeding seasons of 2016 and 2021.We predicted that the waterbird assemblages(i.e.,number of individuals,species richness,Shannon-Wiener diversity,Pielou evenness and species composition) differed between the two years,and that these differences were correlated with the temporal changes in environmental factors.Across the surveyed wetlands,we recorded 26 waterbird species in 2016 and 23 in 2021.For individual wetlands,the number of waterbird individuals and species richness increased by 71.6% and 20.1%,respectively,over the five years,with no changes in Shannon-Wiener diversity and Pielou evenness.The overall increase in the number of bird individuals was mainly caused by an increase in vegetation gleaners and gulls that adapt well to anthropogenic activities.The species composition was significantly different between the two years,which was mainly caused by changes in the number of individuals of dominant species under influence of changes in human activities.For most wetlands,the temporal pairwiseβ-diversities could be explained by species turnover rather than nestedness,probably due to high mobility of waterbird species and dramatic changes in local environments.Our study suggests that waterbird communities could respond to environmental changes in subsidence wetlands,providing important implications for waterbird conservation in human-dominated artificial wetlands.

Determination of convergence of underground gas storage caverns using non-invasive methodology based on land surface subsidence measurement

Undergroundgas storage caverns aremonitoredfor environmental safety in termsof equipmentandpotential emissions,particularly methane emissions from the underground and above-ground parts of the storage facility.Periodical measurements of land surface deformations and costly echometric measurements of convergence of individual storage facilities are carried out.The aims of environmental monitoring are:(1)to eliminate potential hazards in the shortest time,(2)assess the overall impact of intensive operation of storage facilities on the environment,(3)developmonitoringmethods relevant to environmental protection,and(4)take actions in case of failure.The paper presents a solution to the problem of determination of the convergence of underground caverns in a salt rock mass based on the results of land surface subsidence measurements carried out using the Gauss-Markov equalization algorithm.Themethod makes it possible for ongoing control of cavern volume convergence after each subsidence measurement on the ground surface and determining the actual impact of the use frequency(injection-mediumconsumption)on the convergence in time.The presentedmethodology is universal and verified on caverns located in a salt rockmass.The Gauss-Markov inversion model is the first used in this area,hence its application is significant.

Surface Processes Driving Intracontinental Basin Subsidence in the Context of India–Eurasia Collision:Evidence from Flexural Subsidence Modeling of the Cenozoic Southern Tarim Basin along the West Kunlun Foreland,NW Tibetan Plateau

The India–Eurasia collision has produced a number of Cenozoic deep intracontinental basins,which bear important information for revealing the far-afield responses to the remote collision.Despite their significance,their subsiding mechanism remains the subject of debate,with end-member models attributing it to either orogenic or sedimentary load.In this study,we conduct flexural subsidence modeling with a two-dimensional finite elastic plate model on the Hotan-Mazatagh section along the southern Tarim Basin,which defines a key region in the foreland of the West Kunlun Orogen,along the NW margin of the Tibetan Plateau.The modeling results indicate that the orogenic load of West Kunlun triggers the southern Tarim Basin to subside by up to less than ~6 km,with its impact weakening towards the basin interiors until ~230 km north from the Karakax fault.The sedimentary load,consisting of Cenozoic strata,forces the basin to subside by ~2 to~7 km.In combination with the retreat of the proto-Paratethys Sea and the paleogeographic reorganization of the Tarim Basin,we propose that surface processes,in particular a shift from an exorheic to an endorheic drainage system associated with the consequent thick sedimentary load,played a decisive role in forming deep intracontinental basins in the context of the India-Eurasia collision.

Formation of the Zengmu and Beikang Basins,and West Baram Line in the southwestern South China Sea margin

The Zengmu and Beikang basins,separated by the West Baram Line(WBL)in the southwestern South China Sea margin,display distinct geological and geophysical features.However,the nature of the basins and the WBL are debated.Here we explore this issue by conducting the stratigraphic and structural interpretation,faults and subsidence analysis,and lithospheric finite extension modelling using seismic data.Results show that the WBL is a trans-extensional fault zone comprising normal faults and flower structures mainly active in the Late Eocene to Early Miocene.The Zengmu Basin,to the southwest of the WBL,shows an overall synformal geometry,thick folded strata in the Late Eocene to Late Miocene(40.4-5.2 Ma),and pretty small normal faults at the basin edge,which imply that the Zengmu Basin is a foreland basin under the Luconia and Borneo collision in the Sarawak since the Eocene.Furthermore,the basin exhibits two stages of subsidence(fast in 40.4-30 Ma and slow in 30-0 Ma);but the amount of observed subsidence and heat flow are both greater than that predicted by crustal thinning.The Beikang Basin,to the NE of the WBL,consists of the syn-rift faulted sub-basins(45-16.4 Ma)and the post-rift less deformed sequences(16.4-0 Ma).The heat flow(~60 mW/m2)is also consistent with that predicted based on crustal thinning,inferring that it is a rifted basin.However,the basin shows three stages of subsidence(fast in 45-30 Ma,uplift in 30-16.4 Ma,and fast in 16.4-0 Ma).In the uplift stage,the strata were partly folded in the Late Oligocene and partly eroded in the Early Miocene,which is probably caused by the flexural bulging in response to the paleo-South China Sea subduction and the subsequent Dangerous Grounds and Borneo collision in the Sabah to the east of the WBL.

Ground subsidence mechanism of a filling mine with a steeply inclined ore body

Long-term field monitoring finds that serious surface subsidence can still occur even if the high strength cemented fill method is adopted.Combining the results of numerical simulations with global position system(GPS)monitoring,we took a typical filling mining mine with a steeply inclined ore body as an example,and explored its ground subsidence mechanism.The results show that the ground subsidence caused by the mining of steep ore body is characterized by two settlement centers and a significantly uneven spatial distribution,which is visibly different from ground subsidence characteristic of the coal mine.The subsidence on the hanging wall is much larger than that on the footwall,and the settlement center tends to move to the hanging wall with the increase of mining depth.The backfill improves the strength and surrounding rock bearing capacity,which leads to a lag of about 3 years of the subsidence.However,under the actions of continuous and repeated mining disturbances,the supporting effect of the backfill can only reduce the amplitude of the deformation,but it cannot prevent the occurrence of settlement.