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.

Studying Subsidence Coefficient of Mirages and Dynamic Storage Volume of Karstic Springs of Khorram Abad, West of Iran

The studied subject is about subsidence coefficient of mirages and dynamic storage volume of Karstic springs in Khorram Abad in West of Iran. Subsidence coefficient indicates ability of groundwater discharge and hydrologic properties of the environment;meaning effective porosity and transfer coefficient of springs. In general, in developed Karstic zones, each direct line of subsidence curve indicates a discharge regime. Obtained results from the study show that subsidence branch of Golestan, Motahari, Niloofar (changaei), Navekech, Dore Robat mirages have subsidence coefficient with mild slope and low value, which demonstrate passage of water through a seams system at the karst springs. Q and whirlpool stone mirages have two subsidence coefficients, which indicate passage of water through two seam systems in Karstic environment of springs. According to obtained results, process of changes in subsidence branch in these mirages has had at the first a mild slope and low discharge coefficient and in continue, its discharge would be declined with sharper slope and high discharge coefficient. In order to estimate dynamic storage volume of springs, MAILET general equation is applied, which is suitable for subsidence branch of hydrographs of centralized springs discharge. Following, dynamic storage volume of studied springs is analyzed and obtained results are presented in this study respectively.

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.

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.

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.

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.

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%）.

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.

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.

Surface mining subsidence control based on grouting-recovery ratio theory

From an analysis of the relationships of grouting coefficients under the conditions of grouting reducing subsidence, it is concluded that the reducing subsidence rate is in direct proportion to the grouting-recovery ratio, and its proportional coefficient depends on the coefficient of reducing subsidence and collapsing coefficient.The coefficient for reducing subsidence volume φ=ΔV/V_(ash) is related to overburden lithology.The coefficient for reducing subsidence volume φ is always larger than 1.This improves the reducing subsidence rate,the grouting-recovery ratio, and the coefficient for reducing subsidence volume, which can maximize the control of surface subsidence.

The Uplift of the Longmenshan Thrust Belt and Subsidence of the West Sichuan Foreland Basin

Based on fission track dating of apatite, and measurement of vitrinite reflectance of rock samples from the Longmenshan (Longmen Mountain)area and the West Sichuan foreland basin and computer modelling it is concluded that (l)the Songpan-Garze fold belt has uplifted at least by 3-4 km with an uplift rate of no less than 0.3-0.4 mm/a since 10 Ma B.P.; (2) the Longmenshan thrust nappe belt has uplifted at least by 5-6 km with an uplift rate of more than 0.5- 0.6 mm /a since 10 Ma B.P.; (3) the Longmenshan detachment belt has uplifted by 1 - 2 km at a rate of 0.016-0.032 mm/a since 60 Ma B.P.; (4) the West Sichuan foreland basin has uplifted by 1.7-3 km at a rate of 0.028-0.05 mm/a since 60 Ma B.P.; (5) the uplift rate of the area on the west side of the Beichuan-Yingxiu-Xiaoguanzi fault for the last 10 Ma is 40 times as much as that on its east side; (6) the uplifting of the the Songpan - Garze fold belt and the subsidence of the West Sichuan foreland basin 60 Ma ago exhibit a mirro-image correlation, i.e. the rapid uplifting of the the Songpan-Garze fold belt was corresponding to the rapid subsidence of the basin;the Songpan-Garze fold belt has uplifted at a much greater rate than the West Sichuan foeland basin in the last 60 Ma;and (7) the palaeogeothermal gradient was 25℃ /km in the West Sichuan foreland basin.

Innovative backfilling longwall panel layout for better subsidence control effect-separating adjacent subcritical panels with pillars

In recent years,field trials of non-pillar longwall mining using complete backfill have been implemented successively in the Chinese coal mining industry.The objective of this paper is to get a scientific understanding of surface subsidence control effect using such techniques.It begins with a brief overview on complete backfill methods primarily used in China,followed by an analysis of collected subsidence factors under mining with complete backfill.It is concluded that non-pillar longwall panel layout cannot protect surface structures against damages at a relatively large mining height,even though complete backfill is conducted.In such cases,separated longwall panel layout should be applied,i.e.,panel width should be subcritical and stable coal pillars should be left between the adjacent panels.The proposed method takes the principles of subcritical extraction and partial extraction;in conjunction with gob backfilling,surface subsidence can be effectively mitigated,thus protecting surface buildings against mining-induced damage.A general design principle and method of separated panel layout have also been proposed.

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.

Spatial-temporal Characteristics of Land Subsidence Corresponding to Dynamic Groundwater Funnel in Beijing Municipality,China

Due to long-term over-exploitation of groundwater in Beijing Municipality, regional groundwater funnels have formed and land subsidence has been induced. By combining a groundwater monitoring network, GPS monitor- ing network data, radar satellite SAR data, GIS and other new technologies, a coupled process model based on the dy- namic variation of groundwater and the deformation response of land subsidence has been established. The dynamic variation of groundwater fimnels and the land subsidence response process were analyzed systematically in Beijing. Study results indicate that current groundwater funnel areas are distributed mainly in the southwest of Shunyi District, the northeast of Chaoyang District and the northwest of Tongzhou District, with an average decline rate of groundwa- ter level of 2.66 rn/yr and a maximum of 3.82 m/yr in the center of the funnels. Seasonal and interannual differences exist in the response model of land subsidence to groundwater funnels with uneven spatial and temporal distribution, where the maximum land subsidence rate was about --41.08 mm/yr and the area with a subsidence rate greater than 30 mm/yr was about 1637.29 km2. Although a consistency was revealed to exist between a groundwater funnel and the spatial distribution characteristics of the corresponding land subsidence funnel, this consistency was not perfect. The results showed that the response model of land subsidence to the dynamic variation of groundwater was more revealing when combining conventional technologies with InSAR, GIS, GPS, providing a new strategy for environmental and hydrogeological research and a scientific basis for regional land subsidence control.

Concurrent mining and reclamation for underground coal mining subsidence impacts in China

Large scale underground mining of coal resources in China using longwall mining has resulted in ecological and environment problems, including surface subsidence that is considered serious due to competing interests of prime agricultural lands, food security, and regional economic development. The subsided lands must be rehabilitated soon after mining to be agriculturally productive to minimize loss of farmland. Similarly, precious water resources must also be managed during and after mining to protect this natural resource. Toward these goals, the concept of ＂Concurrent mining and subsidence reclamation （CMR）＂ was proposed by Professor Hu of the China University of Mining and Technology, Beijing （CUMTB）. Over the last two decades CMR concepts have evolved and successfully applied in the field in different parts of China. This innovative technology has increased available farmland during the mining process, and provided better land protection and food security in mining areas even with high groundwater table. The technology has been used in 5 of the 14 large coal bases in China. This paper describes the technology concepts, design and guiding principles for planning with two case studies from different regions to enhance its application both in China and in other countries.

Migration of Depocenters and Accumulation Centers and its Indication of Subsidence Centers in the Mesozoic Ordos Basin

Based on the integrated study of structure attributions and characteristics of the original basin in combination with lithology and lithofacies, sedimentary provenance analysis and thickness distribution of the Mesozoic Ordos Basin, it is demonstrated that the depocenters migrated counterclockwise from southeast to the north and then to the southwest from the Middle-Late Triassic to the Early Cretaceous. There were no unified and larger-scale accumulation centers except several small isolated accumulation centers before the Early Cretaceous. The reasons why belts of relatively thick strata were well developed in the western basin in several stages are that this area is near the west boundary of the original Ordos Basin, there was abundant sediment supply and the hydrodynamic effect was strong. Therefore, they stand for local accumulation centers. Until the Early Cretaceous, depocenters, accumulation centers and subsidence centers were superposed as an entity in the southwest part of the Ordos Basin. Up to the end of the Middle Jurassic, there still appeared a paleogeographic and paleostructural higher-in-west and lower-in-east framework in the residual basin to the west of the Yellow River. The depocenters of the Ordos Basin from the Middle-Late Triassic to the Middle Jurassic were superposed consistently. The relatively high thermal maturation of Mesozoic and Paleozoic strata in the depocenters and their neighborhood suggest active deep effects in these areas. Generally, superposition of depocenters in several periods and their consistency with high thermal evolution areas reveal the control of subsidence processes. Therefore, depocenters may represent the positions of the subsidence centers. The subsidence centers （or depocenters） are located in the south of the large-scale cratonic Ordos Basin. This is associated with flexural subsidence of the foreland, resulting from the strong convergence and orogenic activity contemporaneous with the Qinling orogeny.

Prediction of subsidence risk by FMEA using artificial neural network and fuzzy inference system

Construction of metro tunnels in dense and crowded urban areas is faced with many risks, such as sub- sidence. The purpose of this paper was the prediction of subsidence risk by failure mode and effect anal- ysis （FMEA） and fuzzy inference system （FIS）. Fuzzy theory will be able to model uncertainties. Fuzzy FMEA provides a tool that can work in a better way with vague concepts and without sufficient informa- tion than conventional FMEA. In this paper, S and D are obtained from fuzzy rules and 0 is obtained from artificial neural network （ANN）. FMEA is performed by developing a fuzzy risk priority number （FRPN）. The FRPN for two stations in Tehran No.4 subway line is 3.1 and 5.5, respectively. To investigate the suit- ability of this approach, the predictions by FMEA have been compared with actual data. The results show that this method can be useful in the prediction of subsidence risk in urban tunnels.