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.

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.

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.

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.

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.

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.

Numerical modeling of the anomalous post-rift subsidence in the Baiyun Sag,Pearl River Mouth Basin

The Baiyun Sag is the deepest sag in the Pearl River Mouth Basin in northern continental margin of South China Sea,with the maximum sediment thickness over 12.5 km above the basement including >6.5 km sediments above the 30 Ma breakup unconformity.According to the theoretical models for the rifted basins,the post-rift subsidence is driven solely by the thermal contraction and can be calculated as the function of the lithospheric stretching factor.A method combining the forward modeling and reverse backstripping was designed to estimate lithospheric stretching factor.Using the 2D forward modeling based on the flexural cantilever model,we simulated the multi-rifting process of the Baiyun Sag with constrain of the backstripped profiles.By doing this the lithospheric stretching factor was obtained,and then the theoretical post-rift subsidence was calculated.The calculated theoretical subsidence was much smaller than the observed subsidence given by backstripping.Along the 1530 line in the Baiyun Sag,the anomalous post-rift subsidence is over 2 km in the sag center,and varies slightly to the north and south edges of the sag.This suggests that the anomalous post-rift subsidence continues beyond the sag both in the continental shelf to the north and in the continental slope to the south.The sensitivity tests in the forward modeling process indicate that only the use of low-angle faults (≤13°) can we simulate the shape of the backstripped profile.

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.

Comparison of L-band and X-band differential interferometric synthetic aperture radar for mine subsidence monitoring in central Utah

Differential interferometric synthetic aperture radar(DIn SAR), a satellite-based remote sensing technique, has potential application for measuring mine subsidence on a regional scale with high spatial and temporal resolutions. However, the characteristics of synthetic aperture radar(SAR) data and the effectiveness of DIn SAR for subsidence monitoring depend on the radar band(wavelength). This study evaluates the effectiveness of DIn SAR for monitoring subsidence due to longwall mining in central Utah using L-band(24 cm wavelength) SAR data from the advanced land observing satellite(ALOS)and X-band(3 cm wavelength) SAR data from the Terra SAR-X mission. In the Wasatch Plateau region of central Utah, which is characterized by steep terrain and variable ground cover conditions, areas affected by longwall mine subsidence are identifiable using both L-band and X-band DIn SAR.Generally, using L-band data, subsidence magnitudes are measurable. Compared to X-band, L-band data are less affected by signal saturation due to large deformation gradients and by temporal decorrelation due to changes in the surface conditions over time. The L-band data tend to be stable over relatively long periods(months). Short wavelength X-band data are strongly affected by signal saturation and temporal decorrelation, but regions of subsidence are typically identifiable over short periods(days). Additionally,though subsidence magnitudes are difficult to precisely measure in the central Utah region using X-band data, they can often be reasonably estimated.

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.

Lithospheric Stretching,Subsidence and Thermal History Modeling: Application to Yinggehai,Qiongdongnan and Songliao Basins in East China

he lithosphere stretching, subsidence and basement thermal history of the Songliao, Yinggehai and Qiongdongnan basins in East China have been investigated by using backstripping technique and stretching models. The subsidence curve of the Songliao basin characterizes a rift basin history evolving from synrift rapid subsidence to postrift thermal subsidence. The calculated stretching factor of the basin is about 1.8 to 2.0 and the horizontal extension is about 55-60 km. The reconstructed subsidence curves of the Qiongdongnan and Yinggehai basins show that the basins underwent several rapid subsidence episodes, which has been interpreted as the result of a multiple stretching of the lithosphere. The total mean stretching factors of the Qiongdongnan and Yinggehai basins, calculated by using a multiple instantaneous stretching model, are about 1.7-1.9 and 2.5-3.0 respectively. The stretching factors of these basins, estimated from subsidence, are consistent with those determined from crustal thinning. The reconstructed thermal histories of these basins are in agreement with the geological observation.

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.

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.

Development of a new algorithm for implementing the edge effect offset for subsidence calculations

The surface deformation prediction system(SDPS) program has been developed as an engineering tool for the calculation of subsidence deformation indices through the implementation of various prediction methods.From basic user-defined input parameters, SDPS can determine subsidence indices, such as mining induced displacements, strains, tilt, etc., at any elevation between the seam and the horizontal or varying surface topography.A fundamental parameter in obtaining reliable ground deformation results is the determination of the edge effect offset.The value assigned to the edge effect offset corresponds to a virtual offsetting of boundary lines delineating the extracted panel to allow for roof cantilevering over the mined-out area.The objective of this paper is to describe the methods implemented in updating the edge effect offset algorithm within SDPS.Using proven geometric equations, the newly developed algorithm provides a more robust calculation of the offset boundary of the extracted panel for simplistic as well as more complex mining geometries.Given that an extracted panel is represented by a closed polyline, the new edge offset algorithm calculates a polyline offset into the extracted panel with respect to the user defined edge effect offset distance.Surface deformations are then calculated using this adjusted panel geometry.The Matlab program was utilized for development and testing of the new edge effect offset algorithm.After completing rigorous testing regimes, the new offset algorithm will be integrated into SDPS further increasing the speed and reliability of the program resulting in a retrospective increase in capability and flexibility.

Changes of Precipitation Infiltration Recharge in the Circumstances of Coal Mining Subsidence in the Shen-Dong Coal Field,China

Coal mining subsidence is a universal environmental-geological problem in mining areas. By selecting the Shen-Dong coal mining subsidence area as the research field, this paper studies the changes in precipitation infiltration recharge in the circumstances of coal mining subsidence by means of field geological investigation and laboratory simulation experiments, which is expected to provide a scientific basis for eco-environmental restoration in the mining area. The results indicate that at the unstable stage of subsidence, three types of subsidence in the Shen-Dong mining area have positive effects on the precipitation infiltration recharge, and the type of full-thickness bedrock subsidence has the greatest influence. In the stable stage of subsidence, the precipitation infiltration process after long- term drought and the moisture migration in the aeration zone undergo three different stages： evaporation-infiltration before precipitation, infiltration-upward infiltration-infiltration during precipitation and evaporation-infiltration after precipitation. During a heavy rainfall infiltration process, the wetting front movement in fine sand, coarse sand and dualistic structure of fine-coarse sand consists of two stages： the stage of wetting front movement during precipitation, in which the wetting front movement distance has linear relationship with the depth, and the stage of wetting front movement after precipitation, in which the wetting front movement distance has the power function relationship with the depth. The wetting front movement velocity is influenced by the rainfall amount and the lithology in the aeration zone. However, as the depth increases, the movement velocity will decay exponentially.

Numerical analysis of surface subsidence in asymmetric parallel highway tunnels

Tunnelling related hazards are very common in the Himalayan terrain and a number of such instances have been reported. Several twin tunnels are being planned for transportation purposes which will require good understanding for prediction of tunnel deformation and surface settlement during the engineering life of the structure. The deformational behaviour, design of sequential excavation and support of any jointed rock mass are challenging during underground construction. We have raised several commonly assumed issues while performing stability analysis of underground opening at shallow depth. For this purpose, Kainchi-mod Nerchowck twin tunnels(Himachal Pradesh, India) are taken for in-depth analysis of the stability of two asymmetric tunnels to address the influence of topography, twin tunnel dimension and geometry. The host rock encountered during excavation is composed mainly of moderately to highly jointed grey sandstone, maroon sandstone and siltstones. In contrast to equidimensional tunnels where the maximum subsidence is observed vertically above the centreline of the tunnel, the result from the present study shows shifting of the maximum subsidence away from the tunnel centreline. The maximum subsidence of 0.99 mm is observed at 4.54 m left to the escape tunnel centreline whereas the maximum subsidence of 3.14 mm is observed at 8.89 m right to the main tunnel centreline. This shifting clearly indicates the influence of undulating topography and inequidimensional noncircular tunnel.