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.

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.

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.

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.

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.

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.

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.

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.

Monitoring of Coastal Subsidence Changes Based on GNSS Positioning and GNSS-IR

Coastal subsidence monitoring typically employs Global Navigation Satellite System(GNSS)positioning technology.This method provides information only about subsidence below the station base.Sediments in coastal areas tend to accumulate quickly,and subsidence can change significantly due to compaction and alluvium.Therefore,monitoring subsidence above the base is essential to obtain overall coastal subsidence.A new technology called GNSS-Interferometric Reflectometry(GNSS-IR)has been recently developed,which can utilize multipath effects to monitor reflector height.Since the base of the GNSS station is deep and the base length remains constant,the height changes measured by the GNSS-IR technology can reflect subsidence above the base.Accordingly,this paper employs GNSS-IR technology to measure subsidence changes above the base.Additionally,GNSS positioning technology is used to obtain subsidence changes below the base,and the overall subsidence change is then calculated using both GNSS-IR and GNSS positioning technology.The Mississippi River Delta,known for its significant sediment thickness,was selected as the study area,and data from FSHS,GRIS,and MSIN stations was analyzed.The results demonstrate that GNSS-IR can be used to measure the subsidence rate above the base,and the corrected overall subsidence rate is equivalent to the relative sea level rise rate.

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.

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.

Monitoring Study of Long-Term Land Subsidence during Subway Operation in High-Density Urban Areas Based on DInSAR-GPS-GIS Technology and Numerical Simulation

During subway operation,various factors will cause long-term land subsidence,such as the vibration subsidence of foundation soil caused by train vibration load,incomplete consolidation deformation of foundation soil during tunnel construction,dense buildings and structures in the vicinity of the tunnel,and changes in water level in the stratum where the tunnel is located.The monitoring of long-term land subsidence during subway operation in high-density urban areas differs from that in low-density urban construction areas.The former is the gathering point of the entire urban population.There are many complex buildings around the project,busy road traffic,high pedestrian flow,and less vegetation cover.Several existing items requiremonitoring.However,monitoring distance is long,and providing early warning is difficult.This study uses the 2.8 km operation line between Wulin Square station and Ding’an Road station of Hangzhou Subway Line 1 as an example to propose the integrated method of DInSAR-GPS-GIS technology and the key algorithm for long-term land subsidence deformation.Then,it selects multiscene image data to analyze long-termland subsidence of high-density urban areas during subway operation.Results show that long-term land subsidence caused by the operation of Wulin Square station to Ding’an Road station of Hangzhou Subway Line 1 is small,with maximumsubsidence of 30.64 mm,and minimumsubsidence of 11.45 mm,and average subsidence ranging from 19.27 to 21.33 mm.And FLAC3D software was used to verify the monitoring situation,using the geological conditions of the soil in the study area and the tunnel profile to simulate the settlement under vehicle load,and the simulation results tended to be consistent with the monitoring situation.

Subsidence prediction method based on equivalent mining height theory for solid backfilling mining

Based on the characteristics of strata movement of solid backfilling mining technology, the surface subsidence prediction method based on the equivalent mining height theory was proposed, and the parameters selection guideline of this method was also described. While comparing the parameters of caving mining with equivalent height, the subsidence efficient can be calculated according to the mining height and bulk factor of sagging zone and fracture zone, the tangent of main influence angle of solid backfilling mining is reduced by 0.2-0.5（while it cannot be less than 1.0）. For sake of safety, offset of the inflection point is set to zero, and other parameters, such as horizontal movement coefficient and main propagation angle are equal to the corresponding parameters of caving mining with equivalent height. In the last part, a case study of solid backfilling mining subsidence prediction was described. The results show the applicability of this method and the difference of the maximum subsidence point between the prediction and the observation is less than 5%.

A model for extracting large deformation mining subsidence using D-InSAR technique and probability integral method

Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar （D-InSAR） alone, a new algorithm was proposed to extract large deformation mining subsidence using D-InSAR technique and probability integral method. The details of the algorithm are as follows：the control points set, containing correct phase unwrapping points on the subsidence basin edge generated by D-InSAR and several observation points （near the maximum subsidence and inflection points）, was established at first; genetic algorithm （GA） was then used to optimize the parameters of probability integral method; at last, the surface subsidence was deduced according to the optimum parameters. The results of the experiment in Huaibei mining area, China, show that the presented method can generate the correct mining subsidence basin with a few surface observations, and the relative error of maximum subsidence point is about 8.3%, which is much better than that of conventional D-InSAR （relative error is 68.0%）.

Monitoring large-area mining subsidence by GNSS based on IGS stations

In order to monitor large-area mining subsidence accurately, a high-precision global navigation satellite system （GNSS） monitoring network was established based on the nearby international GNSS service （IGS） stations taken as reference points. Given the non-linear motions of IGS stations, the robust Kalman filtering （RKF） model was presented to determine the datum of multi-period monitoring network considering the velocity and weekly solution of IGS stations. The theory proposed was applied to monitoring mining subsidence in northern Anhui coal mine in China. According to the case study, the RKF model to establish monitoring datum is better than the prediction method and the weekly solution from IGS analysis centers （ACs）, and the corresponding precision of deformation can reach up to millimeter level with 4 h observation. The research provides an efficient and accurate approach for monitoring large-area mining subsidence.

Evaluation of PS-DInSAR technology for subsidence monitoring caused by repeated mining in mountainous area

The high resolution Terra SAR-X dataset was employed with DIn SAR and persistent scatterer interferometry（PSI） technique for subsidence monitoring in a mountainous area. For DInS AR technique, the generally used SRTM and relief-DEM, which was derived from aerial topographic map, were used to evaluate the influence of external DEM. The results show that SRTM could not fully compensate the complex topography of the research area. The corner reflectors installed during the acquisition of SAR dataset were used to estimate the accuracy of geocoding. The terrain corrected geocoding results based on relief-DEM were much better than using SRTM, with the root mean square error（RMSE） being 6.35 m in X direction and 11.65 m in Y direction（both in UTM projection）, around one pixel of the multilooked intensity image to be geocoded. For PSI technique, the results from time-series analysis of multi-baseline differential interferograms were integrated to restrict only persistent scatterer candidates near the boundary of subsiding area for regression analysis. The results demonstrate that PSI can refine the boundary of subsidence, which could then be used to derive some angular parameters to help people to learn the law of subsidence caused by repeated excavation in this area.

The Characteristics and Causes of Land Subsidence in Tanggu Based on the GPS Survey System and Numerical Simulation

Land subsidence is a severe hazard threatening Tanggu, a flat lowland area, and evidences of land subsidence can be seen throughout the city. A new reasonable GPS network was set up in this area from 2008 to 2010. The monitoring data show that land subsidence was serious and two main subsidence cones were obviously formed in the region. One emerged at Hujiayuan, with the maximum rate reaching 60 ram/a, and the influence region enlarged prominently from 2005 to 2010. The other one occurred at Kaifaqu, which became obvious only after 2005, and it showed a decreasing tendency with time. To analyze the causes of ground settlement, a correlation between groundwater withdrawal and land subsidence was firstly made. The results confirmed that over-exploitation of groundwater was the major cause for the severe settlement in Hujiayuan. Meanwhile, the subsidence of Kaifaqu was also related to groundwater withdrawal before 2005. However, the relationship became unconspicuous after 2005. To find the cause of this abnormity, a three-dimensional finite element numerical model, coupled with groundwater flow and subsidence, was built. The simulation results indicate that the subsidence induced by high-rise buildings is serious, but the affected range is limited and it also shows a decreasing trend with time, corresponding to the subsidence characteristics at Kaifaqu. Therefore, more attention should be paid to this hazard induced by engineering construction besides groundwater withdrawal, as more high-rise buildings are under construction in Tanggu.

Longwall surface subsidence control by technology of isolated overburden grout injection

Surface subsidence is a typical ground movement due to longwall mining, which causes a series of environmental problems and hazards. In China, intensive coal extractions are commonly operated under dense-populated coalfields, which exacerbates the negative subsequences resulted from surface settlement. Therefore, effective approaches to control the ground subsidence are in urgent need for the Chinese coal mining industry. This paper presents a newly developed subsidence control technology: isolated overburden grout injection, including the theory, technique and applications. Relevant procedures such as injection system design, grouting material selection, borehole layout, grout take estimation and injection process design are proposed. The applicability of this technology has been demonstrated through physical modelling, field measurements, and case studies. Since 2009, the technology has been successfully applied to 14 longwall areas in 9 Chinese coal mines. The ultimate surface subsidence factors vary from 0.10 to 0.15. This method has a great potential to be popularized and performed where longwall mining are implemented under villages and ground infrastructures.