#### Abstract

Due to the acceleration of urbanization in China, the land use and new building construction in mining subsidence areas are increasing. However, a large high-rise building (100 m × 80 m × 100 m) was planned to be built above gob areas in the Xiangcheng Mine subsidence area. According to the analysis of the borehole exploration data, borehole TV observations, and theoretical calculations, there were broken bedrock and developmental fractures in the fracture and caved zones above the gob areas. In addition, the depth affected by the building load is less than 9 m from the top of the fracture zone; thus, the stabilities of the gobs are low. Then, through the theoretical analysis of the space stress, compression settlement of foundation soil, and residual settlement of the gob and coal pillar under loads, combining the similar simulation study on settlement characteristics of foundation soil and gobs, the total settlement of the foundation can be regarded as the vector superposition of the compression settlement and residual settlement; then, the maximum settlement is calculated as 555 mm and the maximum tilt is 14.9 mm/m under the building load, which are greater than the permissible values of 200 mm and 2 mm/m, indicating that the foundation’s stability is inadequate. Therefore, the treatment measures of the deep pile group foundations in the soil layer and grouting reinforcement in the gobs were simultaneously put forward to effectively improve the stability of the building foundations.

#### 1. Introduction

An increasing number of mining subsidence areas are created due to the high intensity and large-scale mining of coal resources. According to the statistics from 2017, there has been more than 20,000 km^{2} of subsidence areas in China, involving the construction of 5000 km^{2} of urban and rural land, which seriously impacts the sustainable development of mining areas after resource depletion. With the acceleration of urbanization and modernization, land resource reuse and construction reclamation in subsidence areas are crucial [1, 2]. Some examples of such residential area reconstruction include the Beipiao mining area and photovoltaic power generation in the Datong mining area. However, destruction to the new buildings or structures constructed above gob areas has also occurred, such as settling, tilting, cracking, or collapse [3]. Some scholars [4–8] have studied the instability of gobs and foundations and put forward failure laws. They found that the residual deformation and instability of gobs are the crucial factors contributing to that destruction. Some scholars [9–12] have researched the movement and destruction of the overlying strata in longwall mining, which can generally be divided into the caved zone, fracture zone, and continuous deformation zone (shown in Figure 1). Research shows that, in the caved zone, the rocks are severely broken and fall to the coal floor with irregular shapes, and a large degree of collapse occurs in the middle of the gob, but some unfilled cavities exist near the boundary. And in the fracture zone, the rocks break into blocks, horizontal and vertical fractures are developed, and the separation is severe; however, near the top, the rock strata are relatively complete, and separation and fractures decrease. Some scholars [14–17] found that residual rock fractures and cavities are still present in the gob areas, as observed by drilling detection, caving rock mass is continuously compacted, and the surface is generating continuous settlement and deformation even after a long period through monitoring. References [18, 19] pointed out that when these fractures and cavities and residual pillars in the gob are influenced by the building load (shown in Figure 1); they would accelerate the movement of the unstable bedrock and the settlement of unstable coal pillars, thereby causing the destruction of the new building foundation.

The Huaibei Mining Group has produced more than 1 billion tons of coal resources, resulting in a large-scale mining subsidence area. Huaibei Mining Group planned to construct an office building and residential buildings above the gobs at the Xiangcheng Mine. Unlike low-rise residential buildings, the office building was uniquely built as a large building (designed size of 100 m × 80 m × 100 m). The range of the foundation would span three gob areas; moreover, the gobs were shallowly buried (from 90 to 115 m) with thick loose layers and thin bedrock, and the range affected by the building load was large. Thus, the additional stress transferred into the soil, fracture, and caved zones, causing instability of the gobs and building foundation. Therefore, studying the regularity of foundation settlement of large buildings is very important for analyzing foundation’s stability, selecting effective treatment measures, and ensuring construction safety.

Currently, many scholars [20–23] have researched the evaluation methods (such as mechanical analysis, similarity and numerical simulations, and nonlinear analysis) for the stability of foundations. However, related field studies on large buildings are insufficient in previous evaluations, especially for the study on the nonuniform settlement of foundation soil and gobs under the large-scale and high-strength load. In this paper, according to the field exploration, mechanical analysis, similarity simulation, and the establishment of settlement models, the effects of different settlement and deformation of the building foundation due to gob areas and coal pillars under the large building load are analyzed, a prediction formula for the settlement and deformation of the large building foundation is developed, and then the stabilities of the gobs and foundation are studied. Finally, effective governance measures are proposed. Different from the previous research on the grouting reinforcement of large foundation by the author himself [13], this paper enhances the theoretical research on the stability of large building foundations, improves the evaluation method of gob areas, and provides a good case for the similar constructions in mining subsidence areas.

#### 2. Exploration of the Gobs and Analysis of Overlying Rock Structure

##### 2.1. Distribution of Gobs

The construction site (delimited range of approximately 400 m × 380 m) and office building were located in seven mining areas of the Xiangcheng coal mine (shown in Figure 2). There were mainly seven gobs (shown in Figure 2); only one layer of coal seam was mined for more than 30 years, with strike longwall blast mining, dip angles of 2 to 8°, an average mining height of 2.5 m, a width of 20 to 100 m, and a buried depth of 90 to 115 m. And the width of the coal pillars between gob was from 2 to 43 m.

##### 2.2. Exploration and TV Imaging Test of the Boreholes

There were mainly six exploration boreholes constructed in gobs (shown in Figure 2). Four boreholes (K1, K2, K3, and K6) were located near the center of the gobs, others were close to the boundary, and all the boreholes incurred wind inhalation or drill felling in the process of drilling, deducing that the rocks were broken and residual cavities existed when being near the boundary.

Based on the drilling exposure of boreholes, the thickness of the loose layer is approximately 66 to 74 m, and the layers are mainly sandy clay, clay, and sand. The thickness of weathered bedrock is approximately 9 to 15 m, mainly for weathered mudstone and weathered sandstone. The thickness of unweathered bedrock is approximately 12 to 28 m, mainly for sandstone and mudstone. The buried depth of the coal seam is 89 to 107 m.

According to the drilling core, the weathered bedrock (depth from 70 to 82 m) was seriously weathered, and fractures were less developed. In the upper sandstone layer (depth from 82 to 88 m), the rock blocks were more complete and had fewer fractures. However, at the bottom of the sandstone layer and the mudstone below (depth from 88 to 97 m), the rock blocks were broken and incomplete, and fractures were developed. The bottom of the bedrock was mainly mudstone (depth from 97 to 110 m), the rock blocks were fragmented, and more cavities were present; thus, the bottom bedrock was located in the caved zone. Based on the mechanics’ test of the rock samples, the compressive strengths of the sandstone and mudstone were less than 20 MPa, indicating that the strength of the bedrock was weak. Thus, the bedrock was prone to occurring large settlement and deformation when being under the large building load.

Based on the observed results of the borehole TV imaging test, the fractures had more developed (depths from 83 to 94 m) and rock mass had collapsed (depths from 94 to 102 m) in the fracture zone and caved zone. Some images of the K6 borehole (depths from 87.1 to 104.5 m) are shown in Figure 3; the observation results of other borehole images are similar. Vertical fractures developed, and rock collapsed from 98.6 m to a location close to the gob in the caved zone. Therefore, there were still rock fractures and cavities in the gob areas.

##### 2.3. Analysis of the Bedrock Structure of a Gob

According to the observation results of rock characteristics and borehole TV imaging tests, the bedrock structure can be divided into weathered bedrock, sandstone layers, and collapsed mudstone (as immediate roofs) from top to bottom (shown in Figure 4). The weathered bedrock was broken, with low bearing capacity and in a state of compression, so it had few fractures. The fracture zone developed to the top of the sandstone based on observations; thus, the sandstone layer broke into blocks, but with relatively high strength and strong bearing capacity. According to the key layer theory [24], the sandstone layer can be used as the main key layer to directly affect the movement of the rock and soil above it. The strength of the immediate roof of the mudstone was low, and it collapsed during mining, so the rock mass broke more. The drilling detection results show that there are still a certain range of cavities and fractures in the rock mass of the caved zone in the middle of the gob, and the range increases near the boundary of the gob, which is the main space for the residual settlement of the gob.

#### 3. Nonlinear Settlement of the Gob and Residual Coal Pillar under Loads

##### 3.1. Distribution of the Additional Stress Field of the Large Building Foundation

###### 3.1.1. Additional Stress of the Foundation

The additional stress of the foundation is caused by the building load, which can cause compression settlement of the foundation. The office building uses rectangular deep foundation; thus, the additional stress of the foundation can be calculated by the Mindlin solution method [25]. The vertical additional stress (*σ*_{fz} for short) generated by the uniformly distributed load below the corner point of the rectangular foundation at depth *z* can be calculated using formula (1).where *a* is the length of the rectangle, *b* is the width of the rectangle, *z*_{0} is the depth of building foundation, *α*_{z} is the additional stress coefficient of the point (0, 0, *z*), *P* is the uniform load, *μ* is Poisson’s ratio of the soil, *z*_{1}, *z*_{2}, *R*_{1}, *R*_{2}, *r*_{1}, *r*_{2}, *r*_{3}, and *r*_{3} are the variable symbols.

The additional stress at any point in the foundation under the rectangular uniform load can be obtained by using the additional stress at the corner through the superposition principle. The magnitude and distribution of the additional stress are positively correlated with the additional stress coefficient, and the coefficient of any point can be initially and quickly calculated by programming with MATLAB software.

###### 3.1.2. Spatial Distribution of the Additional Stress of the Building Foundation

The length and width of the rectangular foundation of the office building are 100 m × 80 m, the foundation depth is 10 m, and Poisson’s ratio of the soil is 0.3. The contour map of the additional stress coefficient in the planes (for *z* equal to 15 m, 40 m, and 70 m) was drawn (shown in Figures 5(a)–5(c)), in which the rectangular frame signifies the foundation boundary.

**(a)**

**(b)**

**(c)**

The additional stress coefficient decreases from the center to the surroundings and outside of the rectangular boundary in the horizontal direction. The smaller the depth, the smaller the attenuation degree of the coefficient in the central region; thus, the higher the additional stress concentration, the more uniform the plane stress, but the greater the change degree at the positions near the two sides of the rectangular boundary, which is prone to uneven settlement. With increasing vertical depth, the additional stress coefficient decreases, and the effect of the load on the soil is reduced, but the influence range increases. At a depth of 70 m below the building foundation, the additional stress coefficient is 0.16–0.26.

##### 3.2. Settlement Calculation of the Caving Rock Mass in the Middle of the Gob

Salamon [26] researched the compaction characteristics of broken rock and determined that the stress-strain relationship of the collapsed rock mass can be calculated using formula (2).where *σ* is the axial stress, *ε* is the axial strain, *ε*_{m} is the maximum possible strain, and *E*_{0} is the initial tangent modulus.

The relationship [27] between *ε*_{m} and *K*_{pc} (residual expansion coefficient of the collapsed rock mass) and *E*_{0} and *σ*_{c} (strength of the collapsed rock mass) can be expressed by formula (3) and formula (4). Considering the effect of the additional stress caused by building loads and combining equations (2), (3), and (4), the compression settlement of the broken rock mass () can be calculated using formula (5).where *h*_{c} is the height of the collapsed rock mass and *σ*_{j} is the additional stresses of the building load.

According to the drilling exploration results, *h*_{c} ranges from 7.44 m to 8.30 m, with an average of 7.87 m and *σ* values of 1.75 MPa. The measured uniaxial compressive strength of mudstone (*σ*_{c}) ranges from 0.75 to 19.26 MPa, with an average of 10 MPa. The building load is 0.675 MPa, and the additional stress coefficient is 0.2, so *σ*_{j} (equal to 0.675 MPa multiplied by 0.2) is 0.135 MPa. According to the observed results of fracture distribution of borehole TV imaging test, the *K*_{pc} value is 1.1. Thus, the calculated value of is 196 mm.

##### 3.3. Settlement Calculation of the Hinged Structure of the Sandstone Layer at the Gob Boundary

###### 3.3.1. Mechanical Analysis of the Instability of the Hinged Structure

Reference [28] pointed out that the strata above the coal pillar and near the inclined boundary of the gob fractured, forming the hinged structure of rock block A, rock block B, and rock block C (shown in Figure 6). The fracture location is the elastic-plastic interface of the coal pillar. The coal pillar, immediate roof, and uncompacted caving rock mass are under rock block B; if the structure of rock block B loses stability under loading, sliding and rotary subsidence will occur, resulting in new subsidence. The stability of rock B has a great influence on the stability of the immediate roof and coal pillar.

The mechanical instability of a hinged structure is mainly rotary deformation instability or sliding instability of rock block B [29]. Reference [30] provided the “S-R” stability theory and instability criterion of hinged structures. Thus, considering the influence of additional stress due to building loads, it was determined that rock block B of the sandstone layer hinged structure below the office building would not experience rotary instability but was prone to sliding instability. After sliding instability, rock block B acts on the immediate roof, coal pillar, or caving rock below, the immediate roof is cut down, or the coal pillar ruptures when the stress exceeds the bearing limit, thus causing large settlement of rock block B.

###### 3.3.2. Settlement Analysis of Rock Block B

According to the mechanical analysis, the subsidence of fractured rock block B after sliding instability can form the following two modes.

Mode 1 is the collapse of the immediate roof under rock block B on one side of the coal pillar. After sliding instability of rock block B (Figure 7(a)), the load on the immediate roof and coal pillar below increase, the suspended part of immediate roof bends, making the stress of the coal pillar in the plastic zone increase when the stress exceeds the limit, resulting in the fracture of the immediate roof and the subsidence of rock block B (Figure 7(b)). Assuming that the destructive immediate roof and coal pillar are effectively broken into blocks and filled cavities, according to the principle of space conservation, the settlement of rock block B can be calculated using formula (6).where *M* is the mining thickness, *l*_{0} is the suspended distance of the immediate roof, *x*_{0} is the width of the plastic zone of the coal pillar, and Δ*S* is the subsidence of the rock block B.

The hinged structure above both sides of the coal pillar and coal pillar failure below rock block B are referred to as mode 2. When the load of rock block B of the hinged structure on both sides acts on the immediate roof and the coal pillar (Figure 8(a)), compared with the one side (Figure 7(a)), the stress of the immediate roof and coal pillar is higher, and the coal pillar is more prone to failure than the immediate pillar. If the width of the coal pillar (*B* for short) is small, the overall coal pillar is destroyed (Figure 8(b)), and the settlement of rock block B can be calculated using formula (7). However, if the width of the coal pillar is large, analogous to model 1, the coal pillar in the plastic zone of both sides would be destroyed, and the settlement can be calculated using formula (6).

##### 3.4. Calculation of the Instability Width of the Residual Coal Pillar

The load on the coal pillar can be regarded as the sum of the load of the uncollapsed rock mass and the loose layer above (shown in Figure 9). Considering the influence of the additional stress provided by the building load, the load (*P*_{m}) and stress (*σ*_{m}) of the coal pillar can be calculated using formula (8) and formula (9).where *H*_{j} is the thickness of the bedrock, *γ*_{j} is the average volume force of the bedrock, *h*_{s} is the thickness of the loose layer, *γ*_{s} is the average volume force of the loose layer, and *φ* is the rock collapse angle.

Referencing the calculation of the limit load of the strip mining coal pillar, the limit load (*σ*_{jx}) and stability coefficient (*K*_{m}) of the residual coal pillar can be calculated using formula (10).

According to the drilling exploration results and parameter values (*l*_{0} of 2 m, *M* of 2.5 m, *φ* of 70°, and *K*_{m} of 1), the minimum stable width of the coal pillar under or without building loads is 10 m and 9 m, respectively. The width of the residual coal pillar under the building ranges from 2 to 26 m, but it is primarily less than 10 m. The settlement of the sandstone layer at the boundary of the gob can be calculated using formula (11) for different width ranges.

#### 4. Similar Simulation Study on Cavity Distribution and Settlement of Gobs under Loads

##### 4.1. Design of the Experimental Schematic

The two-dimensional simulation experiment platform (1800 mm × 160 mm × 1300 mm) was adopted by taking the gobs of Xiangcheng Mine as the engineering background. According to the similarity principle and setting the geometric similarity ratio and rock strength ratio of the model to 200 and 320, respectively, the equivalent load of the building is applied by the counterweight block. In the model, the bedrock thickness is 15 cm, the loose layer is 35 cm, three gobs are excavated, the width of the coal pillar is 2 cm between two gobs, displacement measuring points (8 rows and 27 columns) are laid in different depths of the loose bedrock layer (3 rows in bedrock), and the measuring points under the central region of the load are intensively arranged. The layout and actual shapes of the model are shown in Figures 10(a) and 10(b).

**(a)**

**(b)**

##### 4.2. Analysis of the Results of the Simulation Experiment

###### 4.2.1. Distribution of Cavities in the Bedrock and Loose Layer above the Gobs after Mining

The results of the cavity distribution of the bedrock and loose layer after excavating the three gobs are shown in Figure 11. The range of the cavities is large, and the range decreases with rock mass turning and sliding for the hinged structure on the gob boundary. The closer to the middle of the gob is, the better the compaction is and the fewer cracks there are. There are more cavities and fractures in the caved zone, mainly the bed separation in the fracture zone, and the greater the height, the smaller the separation. With increasing mining width, the bending and fracture of rock increases, the cavity range decreases on the boundary, and the separation space decreases in the middle of the gob.

The loose layer bends and sinks with the subsidence of rock, and fractures develop upwards for a caving angle of nearly 90°. The loose layer has difficulty forming a structure due to the low strength; thus, continuous bending upward within the influence range of mining, fractures on both sides, and separation in the middle develop. With increasing mining time, cracks are continuously compacted.

###### 4.2.2. Distribution of the Cavities of the Bedrock and Loose Layer above the Gobs after Building Loads Are Applied

The building load is simulated by placing the counterweight block and setting 3 load levels (laying a counterweight block for 1 to 3 layers). Under the action of the first level of loading, the upper soil layer of the loose layer is forced to compress and settle, the load is transferred to the deep soil, the fractures and separation are compacted (shown in Figure 12(b)), and the depth affected by the additional stress of the load is 30 m. Under the second level of loading, the influence range expands downward and on both sides, the fractures and separation are continuously compacted, the affected depth is approximately 60 m, and the bending and subsidence range of the loose layer increases. Under the third level of loading, the load influence spreads to the bedrock of the 574 gob, the strata bend and sink in the middle of the gob, so the separation decreases, and the hinged fractured rock blocks on both sides rotate and slide, so the cavities decrease. Finally, the coal pillar is broken and collapses, resulting in rapid subsidence of the upper bedrock until it is stable after contact with the rocks on both sides (Figure 12(c)).

**(a)**

**(b)**

**(c)**

###### 4.2.3. Settlement of the Bedrock and Loose Layer above the Gob under Loading

The settlement change trends of the bedrock layer and loose layer in the adjacent measuring line are close, so line 1 (location in sandstone layer) and line 6 (near central depth of loose layer) are selected to analyze the settlement change in different stages. Figure 13(a) shows that the settlement of the gob and coal pillar increase under loading, but the change degree is different. Stratal separation in the bedrock decreases, and subsidence increases in the middle of the gob but decreases gradually to both sides. As the load increases, the settlement increases on the gob boundary for fractured rock block rotation. When the coal pillar (location from 97 to 99 cm) is broken, the settlement increases from 0.41 m to 0.85 m, the rock strata on both sides in the 8 cm range rotate and sink, and the settlement increases from 0.43 m to 0.53 m. Figure 13(b) shows that the settlement values of each measuring point in the loose layer are close when the load is small. With the increase in the load, the settlement under the load increases apparently and decreases gradually to both sides, and the bedrock subsidence increases the settlement of the loose layer. When the coal pillar is broken, the settlement of the roof and loose layer above it increases rapidly. The settlement range of the loose layer is larger than that of the gob below, and the loose layer bends and sinks within the influence range of mining and loads; thus, the settlement of the loose layer under loading can be regarded as the superposition of the compaction settlement itself and the residual settlement of the gob.

**(a)**

**(b)**

#### 5. Prediction of the Settlement and Deformation of the Foundation under the Building Load

##### 5.1. Calculation Model of the Compression Settlement of the Foundation Soil

###### 5.1.1. Spatial Stress Analysis and Layered Total Settlement of the Foundation Soil

To simplify the calculation, the influence range of the rectangular load (*a* × *b*) on the foundation soil below can be regarded as a cuboid of *l*_{a} × *l*_{b} × *H*_{z} (shown in Figure 14(a)). According to the design specification [31] and considering the influence of the shallow quicksand layer and gob, the regarded depth is where the additional stress is 10% of the gravity stress (*σ*_{z} for short) of the soil under the foundation, which is the depth affected by the building load (*H*_{z} for short, where the top is the ground). The connection line of the influence boundary on the soil surface with a load on the centerline of the long and short sides of the rectangle is calculated as the length of *l*_{a} and *l*_{b}, *H*_{z} is the maximum depth below the center point of the foundation, the spatial relationship is shown in Figure 14(b), and the vertical additional stress at any point in the affected spatial range of soil can be quickly calculated based on the analysis in Section 3.1.1.

**(a)**

**(b)**

Considering that the soil has the characteristics of stratification and compressibility, the settlement can be calculated using the layered summation method. The settlement at a point (*x*, *y,* 0) on the surface of the foundation can be regarded as the accumulation of the compression settlement of each unit of soil on the straight line from the point to the depth of *H*_{z}, so the settlement value of the unit at point (*x*, *y*, *z*) (d*S*_{(x,y,z)} for short) and point (*x*, *y*, 0) (S_{(x,y,0)} for short) can be calculated using formula (12). It can also be divided into multiple soil layers, and each soil layer is integrated and then summing.where *α*_{z(x, y, z)} is the vertical additional stress coefficient at point (*x*, *y*, *z*) and *E*_{s(x, y, z)} is the compression modulus of the unit soil at point (*x*, *y*, *z*).

###### 5.1.2. Compression Settlement under the Office Building

The thickness of the loose layer is 70 m, and the average bulk density is 20 kN/m^{3}. The maximum load of the office building is 675 kPa, taking the additional stress from the building load equal to 10% of the gravity stress of the soil as the criterion for the affected depth. The curves of the additional stress under different loads and gravity stress are shown in Figure 15. It can be seen that the greater the load, the greater the affected depth. When the building load is 675 kPa, the maximum affected depth is 78.1 m, indicating that the additional stress affects the bedrock. Similar to the calculation of affected depth, the planerange of the influenced soil is 50 m extended outward from the rectangular boundary, so the space range of the influenced foundation soil is 200 m × 180 m × 70 m.

Referencing the experimental data of adjacent soil samples from the Haizi mine, combined with the layered lithology, thickness, and burial depth data of the K6 borehole, the compression modulus can be calculated and seen in Table 1.

The settlement of any point on the foundation soil surface is calculated by formula (12), and the settlement of the three-dimensional map and contours of the soil surface under the building load is obtained (shown in Figure 16). The maximum settlement of soil occurs below the load of 359 mm and gradually decreases from the center to the surroundings. Near the load boundary, the settlement value decreases rapidly; thus, the uneven settlement is large. The tilt curve of the soil surface in the midline of the short side of the rectangular foundation is produced (shown in Figure 16(c)). The tilt from the center to both sides first increases and then decreases, and the maximum tilt at the boundary is 14 mm/m. In the range of 20 m near both sides of the boundary, the change degree of tilt is the largest, indicating that the foundation is prone to uneven settlement.

**(a)**

**(b)**

**(c)**

##### 5.2. Calculation of the Residual Settlement of the Gobs under the Building Load

The loose layer will be accompanied by the settlement of the sandstone layer, so the settlement of the sandstone layer can be regarded as the equivalent mining thickness (*M*_{d} for short). The residual settlement can be predicted by the probability integral method.

###### 5.2.1. Prediction by the Probability Integral Method

According to the probability integral method [32], taking the equivalent mining thickness at a point (*s*, *t*) in the affected area of the gob as *M*_{d} (*s*, *t*), so the settlement at any point (*x*, *y*) on the ground surface can be calculated using formula (13). Then, the tilt (*i*_{(x,y,φ)}), horizontal movement (*U*_{(x,y,φ)}), curvature (*K*_{(x,y,φ)}), and horizontal deformation (*ε*_{(x,y,φ)}) in a given *φ* direction can be calculated using formula (14).where *q*_{c} is the residual sinking coefficient, *H* is the buried depth, *α* is the dip angle of coal seam, tan*β* is the main effecting angle tangent, *r* is the radius of influence, and *b* is the displacement factor.

###### 5.2.2. Affected Area of the Gob and Equivalent Mining Thickness

The affected area of the gob, delimited by the vertical section method [33], was larger than that of the building (shown in Figure 17). According to the analysis of observation results from adjacent mines, the parameters (*φ*, *δ*, *γ*, *β*) were selected; thus, the length and width of the affected area were calculated as 291 m and 254 m.

According to the distribution of gobs and coal pillars in the affected area (shown in Figure 18), combining the analysis of settlement in the middle of the gob, settlement of the fractured structure on the boundary, and settlement of the broken coal pillars, the affected area can be divided into four types, including the midarea (*M*_{d} values of 196 mm), coal pillar instability area (*M*_{d} values of 196 mm), coal pillar stability area (*M*_{d} values of 196 mm), and potential subsidence area (*M*_{d} values of 1250 mm).

###### 5.2.3. Calculation Results of the Residual Settlement and Deformation of the Gobs

The parameters of the probability integral method are selected based on the similarities with the data of adjacent mines (shown in Table 2). The residual settlement, tilt, curvature, and horizontal deformation of the gob were calculated according to the range of different regions and the equivalent mining thickness in the affected area of the gob, and the contour maps of the settlement and tilt are shown in Figures 19 and 20 . According to the contour, the maximum settlement, tilt, curvature, and horizontal deformation below the building foundation were 196 mm, 1.45 mm/m, 0.06 (10^{−3}/m), and 1.44 mm/m, respectively.

##### 5.3. Prediction of the Total Settlement and Deformation of the Foundation

The total settlement of the foundation can be regarded as the vector superposition of the compressive settlement of the soil and residual settlement of the gob. Using the calculation coordinate system of the residual settlement, the total settlement (*S*_{ft} (*x, y*)) can be calculated using formula (15). Then, the tilt (*i*_{(x,y,φ)}), curvature (*K*_{(x,y,φ)}), and horizontal deformation (*ε*_{(x,y,φ)}) in a given *φ* direction can be calculated using formula (16). Thus, the maximum total settlement, tilt, curvature, and horizontal deformation below the building foundation were 555 mm, 14.9 mm/m, 0.06 (10^{−3}/m), and 1.44 mm/m, respectively.

#### 6. Stability Analysis and Treatment Measures of the Large Building Foundation

##### 6.1. Stability Evaluation of the Foundation

###### 6.1.1. Relationship between the Influence of the Building Load and Fracture Zone of the Gob

According to the analysis of the affected depth based on the building load and observation result of the fracture zone height (*H*_{li} for short) of the ZK6 borehole under the foundation, the maximum affected depth is 78.1 m, the height of the fractured zone is 12.6 m, and the affected depth has not reached the fractured zone, but the distance is only 8.7 m between the two which is less than the empirical value of 10 m; therefore the gob will cause the instability of the foundation.

###### 6.1.2. Analysis of the Settlement and Deformation of the Foundation

The maximum settlement and deformation of the foundation below the office building are shown in Table 3. According to the specification [31], the permissible settlement is 200 mm, and the permissible tilt is 2 mm/m for a high-rise building. The maximum settlement (359 mm) and tilt (14 mm/m) of the soil are greater than the permissible value, and the maximum residual settlement (196 mm) and tilt (1.45 mm/m) of the gob are close to the permissible value, indicating that the compression of loose soil has a great influence on the foundation, and the risk is greater after superposition; therefore, the stability of the building foundation is poor. The foundation is prone to large settlement and uneven tilt damage without treatment measures.

##### 6.2. Treatment Measures

Due to the large thickness of the loose layer below the construction site, the compression settlement of the soil under the building load is large. Additionally, there are several flow sand layers (approximate overall thickness of 11 m) in the shallow soil from 15 to 30 m, which are prone to shear slip due to seepage, thus increasing the uneven settlement of the soil. Considering that there are clay and clay gravel layers in the buried depth from 40 to 70 m of the loose layer with a high degree of consolidation and low compressibility, the pile group foundation is adopted, and the depth of the pile is 50 m, which passes through the flowing sand layers to transfer the load to the deep soil layer, thus effectively reducing soil settlement.

Grouting reinforcement in the gob is an effective measure to control residual settlement. The cement-fly ash slurry is injected into the gob through drilling to fill the cavities and fractures in the caved zone and fracture zone, thus reducing the residual settlement of the gob.

#### 7. Conclusions

Through the detection results of borehole drilling and color drilling TV observations of gobs, the bedrock structure over the gob can be divided into the weathered bedrock layer, fractured sandstone layer (key layer), and collapsed mudstone (immediate roof) layer. There is still some range of cavities and fractures in the caved zone in the middle of the gob, and the space increases near the boundary, which is the main space of residual settlement in the gob.

The maximum compressive settlement of the foundation soil below the office building is 359 mm, and the settlement is relatively uniform in the center area but the settlement and the tilt change greatly when being close to the load boundary. The maximum tilt is 14 mm/m, and the foundation is prone to generate uneven settlement near the boundary. Based on the random medium theory and equivalent mining height theory, the maximum residual settlement and tilt of the gob under the building load are 196 mm and 1.45 mm/m, respectively.

The total settlement of the foundation can be regarded as the vector superposition of the compression settlement and residual settlement. The maximum total settlement is 555 mm, and the maximum tilt is 14.9 mm/m of the foundation under the building load, which are greater than the permissible values (200 mm and 2 mm/m, resp.). In addition, the depth affected by the building load is less than 9 m from the top of the fracture zone, so the stability of the gob areas is poor. This indicates that the stability of the building foundation at the construction site was adverse; therefore, governance measures are necessary.

A deep pile group foundation (pile tip depth 50 m) is used to go through the flow sand layer in the shallow soil to transmit the load to the deep clay layer with small compressibility, thus effectively reducing the uneven settlement of the soil. In addition, grouting slurry fills cavities and cracks in the caved zone and fracture zone of gob areas, thus effectively reducing the residual settlement. The two treatment measures can effectively improve the stability of the building foundation.

#### Data Availability

The data used to support the findings of this study are included within the article partly or available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

The authors acknowledge Tiandi Science and Technology Co., Ltd. for providing the project data. Support from the Special Project of Science and Technology Innovation Fund of China Coal Technology and Engineering Group (no. 2019-ZD004), National Science and Technology Major Project (no. 2016ZX05068002-003), and National Natural Science Foundation of China (no. 51804162) is also gratefully acknowledged.