Introduction

In recent years, with the rapid development of extensive cities in China, urbanization has provided people the convenience and efficiency of modern life, while problems such as a large population coupled with a small land area, a shortage of resources, traffic congestion, urban waterlogging, and environmental pollution have become increasingly prominent, seriously affecting the quality of life of people and the sustainable development of cities. The efficient and rational development and utilization of the underground space could resolve these urban diseases. Scientific assessment of the amount of available resources, development difficulty, potential value and comprehensive quality of the urban underground space is an important basic task to realize rational urban underground space planning1,2.

At present, many scholars have evaluated urban underground space resources from different perspectives. To evaluate the amount of resources available for the development and utilization of the underground space in Tianjin, Wang Yongli3 used a GIS-fuzzy comprehensive evaluation model. Through evaluation, a database of underground space resources was established, and a distribution map of underground space resources was generated. The resources available for the development and utilization of the underground space were characterized and measured. To evaluate the difficulty of underground space development, Li Pengyue et al.4 established a geological suitability evaluation model for the development and utilization of urban underground space resources based on their collaborative development and utilization and the comprehensive response to geological environmental problems. Analytical methods and fuzzy comprehensive evaluation methods were used to evaluate the difficulty of developing underground space resources in Chengdu. To evaluate the development potential of the urban underground space, Liu Jianan et al.5 adopted Wuhan city as an example, comprehensively considered three factors, namely, geological environment suitability, existing construction restrictions, and social and economic benefits, to construct an evaluation index system, and used the analytic hierarchy process for calculating a comprehensive underground space index of the development potential and conducting potential zoning by systematically considering ecological protection and cultural relic protection restrictions in Wuhan city. This index system and potential zones provided a scientific reference for the subsequent development and utilization of the urban underground space. Wang Xia et al.6 used the multilevel grey evaluation method to construct an evaluation index system for the development potential of urban underground space resources and evaluated the development potential of the underground space in Nanjing. This evaluation method required few samples and could be used to address the information inconsistency in the evaluation process. The problem of completeness and uncertainty provides new ideas for the evaluation of the development potential of urban underground space resources. To explore comprehensive quality assessment of underground space resources, Wu Xinzhen7 used GIS technology and a fuzzy comprehensive evaluation method to assess the suitability and potential development value of underground space resources in Wuhu city, conducted a comprehensive quality assessment of underground space resources, and delineated underground space zoning management and control strategies. The evaluation method and results guided the subsequent execution of urban underground space comprehensive utilization planning.

There are many ways to evaluate urban underground space resources. The more mature methods include the analytic hierarchy process8,9,10,11,12, fuzzy comprehensive evaluation method13,14,15, extension evaluation method16,17, negative list method18, variable fuzzy set theory19, multilevel grey evaluation method6, fuzzy c-means clustering algorithm20, random forest algorithm21, or a combination or improvement of the above methods4,7,22,23,24,25,26,27,28,29,30. The above methods have been widely employed to evaluate the amount of resources available for development, development difficulty, potential value and quality of underground space resources. However, there are certain limitations when used for the comprehensive evaluation of underground space resources from multiple perspectives, such as insufficient selection of indicator factors and weak hierarchy. Therefore, in this paper, we used multiple methods and combined models to comprehensively evaluate urban underground space resources. Since there are many factors that affect the development of underground space resources and each evaluation factor exhibits different attributes, indicators, and qualitative and quantitative standards, this task is complex and multilayered. Therefore, it is impractical to use a deterministic model to calculate the weights of evaluation factors.

Based on the above analysis, in this paper, we constructed a comprehensive evaluation index system for urban underground space resources, evaluated the underground space development resource quantity, development difficulty, and potential value and established comprehensive quality assessment models. Evaluation indicators, including urban environmental constraints, geological conditions, socioeconomic conditions and many other factors, were determined. The method of item-by-item elimination of restrictive factors was used to evaluate the amount of underground space resources available for development, and the analytic hierarchy process was used to determine the weight of each evaluation index. The evaluation results were visualized through GIS technology to comprehensively evaluate underground space resources in the study area.

With the rapid economic development and population growth in the main urban area of Xianyang city, which is an important development city in Shaanxi Province, urban land resources are in short supply, and the contradiction between humans and land is prominent. The large-scale development and utilization of the urban underground space are inevitable during sustainable urban development. Therefore, it is of practical significance to comprehensively evaluate the underground space resources in Xianyang city.

Description of the study area

The scope of the study area is the main urban area of Xianyang city. Xianyang city is located in the hinterland of the Guanzhong Plain in Shaanxi Province, bordering Xi'an to the east, Chang'an and Hu County to the south, Yangling, Fufeng, and Linyou to the west, and the Tongchuan and Qingyang and Pingliang areas in Gansu Province to the north. Moreover, the terrain is high in the north and low in the south, with low mountains in the north and plains in the south. From north to south, the landforms include medium and low mountains, loess hills, loess plateaus and Weihe River terraces, as shown in Fig. 1.

Figure 1
figure 1

Overview of the location of the study area.

The scope of the study area is the main urban area of Xianyang city. Xianyang city is located in the hinterland of the Guanzhong Plain in Shaanxi Province, bordering Xi’an in the east, Chang’an Hu County in the south, Yangling, Fufeng and Linyou in the west, and the Tongchuan and Qingyang and Pingliang areas in Gansu in the north. The terrain is high in the north and low in the south, with low mountains in the north and plains in the south. From north to south, the landforms include medium and low mountains, loess hills, loess plateaus and Weihe terraces, as shown in Fig. 1. The strata are Quaternary, Palaeogene and Neogene strata, and the lithology within 0–30 m mainly includes clayey soil, sand, gravel soil, loess, and gravel. The groundwater type in the area is karst groundwater. As of 2019, there were 245 underground spaces in the main urban area of Xianyang city, with a total development area of 1.549 million m2, and the per capita underground space area was approximately 2.2 m2. The type of underground space utilization mainly entailed underground garages, accounting for 82% of the underground space development ratio in Xianyang city. Moreover, underground space development was limited to shallow depths, with a scattered distribution, thus lacking horizontal connectivity, while a network could not be formed.

Study methods and evaluation process

Evaluation method

To comprehensively evaluate underground space resources, the quantity, quality, and favourable and restrictive factors of resource development of underground space resources should be assessed within a certain scope. In this paper, the underground space resources available for development in the study area were evaluated by the method of eliminating influencing factors one by one. Because there are many factors affecting the comprehensive quality of underground space resources and these factors are mutually restricted, interrelated, multilevel and complex, the analytic hierarchy process was used for evaluation, and the ArcGIS spatial analysis function was employed to visualize the evaluation results. A flow chart of the evaluation process is shown in Fig. 2.

Figure 2
figure 2

Evaluation method structure diagram.

Evaluation model

Estimation model for the underground space resources available for development

The degree of exploitability of underground space resources is affected by many factors. The most effective method for estimating the amount of underground space resources potentially available for development is item-by-item elimination of restrictive factors. First, we must determine the types of urban land use and apply remote sensing methods to obtain land use type data such as ground buildings, water bodies, green spaces, and farmland31,32. Then, we must divide the study area into difficult-to-develop areas, insufficient-to-develop areas, and fully developable areas according to the development and utilization of resources (Table 1). Unexploitable areas can be obtained by superimposing the underground space range under each constraint factor, and the amount of unexploitable resources can be deducted from the total amount of underground space resources to obtain the amount of resources available for development33. The calculation model of urban underground space resources under each land type is summarized in Table 2. The model can be expressed as:

$$ {\varvec{V}} = {\varvec{V}}_{{\varvec{n}}} - \bigcup\limits_{{{\varvec{i}} = 1}}^{k} {{\varvec{V}}_{{\varvec{i}}} } $$
(1)

where Vn denotes the natural underground space resources in the planning area, V denotes the potential exploitable space resources, and Vi (i = 1, 2… k) is the amount of space resources under restrictive factors such as geological hazards, water bodies and developed underground spaces.

Table 1 Factors influencing the exploitable degree of underground space resources.
Table 2 Amount of urban underground space resources under the various land use types.

Evaluation model for the difficulty of underground space resource development

Engineering geological conditions and urban environmental constraints are the main factors affecting the difficulty of underground space resource development, and the development difficulty can be evaluated according to shallow (0 ~ 15 m) and subshallow (15 ~ 30 m) areas on the basis of various factors affecting different spatial domains. Regarding the main urban area of Xianyang city, the influencing factors of the difficulty of underground space resource development are shown in Tables 3 and 4, and the evaluation models can be expressed as follows:

$$ S_{d} = \alpha \sum\limits_{j = 1}^{m} {w_{j} } \sum\limits_{k = 1}^{n} {w_{jk} u_{jk} } $$
(2)

where ujk is the value of the kth indicator in the jth indicator layer; wjk is the weight of the kth indicator in the jth indicator layer; wj is the weight of the jth topic layer; Sd is the development difficulty evaluation result; and α is the difficulty reduction factor. Usually, the deeper the excavation is, the more difficult the construction. According to past experience, shallow and subshallow layers are considered, and the reduction coefficients are 1 and 0.9, respectively.

Table 3 Weight of the resource development difficulty index of the shallow (0 ~ 15 m) underground space.
Table 4 Weight of the resource development difficulty index of the subshallow (15 ~ 30 m) underground space.

Evaluation model for the potential value of underground space resource development

The potential value of underground space resources can be determined via the comprehensive assessment of the probability of obtaining the expected levels of social, economic and environmental benefits at the urban location of underground space resources. According to the evaluation of the Xianyang area, through a large number of investigations, the potential value evaluation index of underground space resource development was finally determined (Table 5). With the use of the multifactor weighted average method, a potential value assessment model for underground space resources based on socioeconomic conditions was established as follows:

$$ S_{e} = \alpha \sum\limits_{j = 1}^{m} {w_{j} u_{j} } $$
(3)

where uj is the value of the jth indicator, wj is the weight of the jth index, Se is the potential value, and α is the value reduction coefficient. The deeper the excavation is, the lower the potential value of development.

Table 5 Potential value index weight of underground space resources.

Comprehensive quality evaluation model for underground space resources

The assessment and calculation of the comprehensive quality of underground space resources can reflect the comprehensive utilization value of underground space resources under the influence of geology, social economy and other conditions, which can be obtained by superimposing the difficulty and potential value of underground space development. The comprehensive quality can be expressed as:

$$ S = w_{d} S_{d} + w_{e} S_{e} $$
(4)

where wd and we are the weights of the natural and socioeconomic factors, respectively; Sd and Se are the resource development difficulty and potential value assessment value, respectively; and S is the overall quality.

Evaluation index weight

First, we compare the degrees of influence of the indicators to construct a judgement matrix and determine the relative importance of each pair of factors. This value is generally determined by the nine-point scaling method34. The weight of each indicator can be calculated by the arithmetic mean method in the analytical hierarchy process, as expressed in Eq. (5). The constructed comparison judgement matrix must be assessed for consistency. The criterion for determining whether the matrix meets the consistency requirement is CR < 0.1. For CR ≥ 0.1, the matrix does not meet the consistency requirements and must be rebuilt. CR can be calculated with Eq. (6). The weights of the indicators at all levels can be obtained according to the AHP flow chart, as shown in Fig. 3. The weights of all evaluation indicators are shown in Tables 3, 4, 5 and 6.

$$ \omega_{i} = \frac{1}{n}\sum\limits_{j = 1}^{n} {\frac{{a_{ij} }}{{\sum\nolimits_{k = 1}^{n} {a{}_{kj}} }}} $$
(5)
$$ CR = \frac{CI}{{RI}} $$
(6)
$$ \lambda_{\max } = \sum\limits_{i = 1}^{n} {\frac{{\sum\nolimits_{j = 1}^{n} {a_{ij} \omega_{j} } }}{{n\omega_{i} }}} $$
(7)

where \(\omega_{i}\) is the index weight, \(a_{ij}\) is the normalized value of the index in Row i and Column j, CI = (\(\lambda_{\max }\) − n)/(n − 1), n denotes the order of the matrix, \(\lambda_{\max }\) is the largest characteristic root of the matrix, which can be calculated according to Eq. (7), and RI is the average random consistency index.

Figure 3
figure 3

Flow chart of weight calculation using the AHP.

Table 6 Comprehensive quality weight of underground space resources.

Grading and quantification of the evaluation indicators

The evaluation indicators for the difficulty and potential value of underground space resource development are graded and quantified. The quantitative scores are selected from the [0, 1] interval according to the linear interpolation method. In this study, we divide each indicator into 4 levels, and the evaluation indicators are graded and scored. The first level is assigned a value of 1, the second level is assigned a value of 0.67, the third level is assigned a value of 0.33, and the fourth level is assigned a value of 0. According to the actual situation of Xianyang city, the scores of the different levels of evaluation indicators are shown in Tables 7 and 8.

Table 7 Difficulty index grading score of underground space resource development.
Table 8 Grading and scoring of the potential value indicators for underground space resources.

Grading of the evaluation results

With the use of the aforementioned evaluation model and method, based on the ArcGIS spatial analysis function, the single-factor layers are weighted and superimposed, and the evaluation results are graded by the principle of score interval evaluation and distribution (Table 9).

Table 9 Grading of the evaluation results.

Evaluation results and analysis

Estimation of the available resources for underground space development

With the use of ArcGIS statistics, excluding existing buildings, cultural relic reserves, ecological green space protection ranges, geological disaster-prone areas and water source protection areas in the main urban area of Xianyang city, the factors that control the underground space amount can be obtained, and the amount of development accounts for approximately 25.11% of the total development. According to Eq. (1) and Table 2, the underground space resources in shallow (0 ~ 15 m) and subshallow (15 ~ 30 m) areas in the main urban area of Xianyang city were estimated, and the underground space resources that could be reasonably developed under the different land types were obtained (Table 10).

Table 10 Amount of underground space resources available for rational development in the main urban area of Xianyang city (unit: 100 million m3).

The calculation results indicate that the current underground space development level in Xianyang city is still low, the exploitable area of the underground space in Xianyang city is approximately 82.3 km2, and the amount of shallow underground space in the main urban area of 30 m accounts for approximately 79.5% of the total underground space resources of the shallow and subshallow layers. However, underground buildings are different from surface buildings, with the characteristics of irreversible construction, and underground structures and the built environment often interact and influence each other. For example, underground structures at a certain development depth impose certain restrictions on the development of other structures at the same level. The difference in the degree of response of the different layers of different types of underground structures to seismic waves affects the interaction in development layer site selection. The current sufficient amount of development does not suggest that the underground space can be developed arbitrarily and subjectively, and unreasonable development can cause the proportion of the underground space that can be developed to be sharply reduced. Regarding underground transportation, the development of underground structures at unreasonable depths exerts a restrictive impact on the selection of underground transportation routes, reducing the amount of development along the traffic line.

Assessment of the difficulty of underground space resource development

Based on the spatial analysis function in ArcGIS, the single-factor layers are superimposed to evaluate the development difficulty of shallow (0 ~ 15 m) and subshallow (15 ~ 30 m) underground space resources in the main urban area of Xianyang city, and a difficulty grading map of shallow to subshallow underground space resource development is generated (Figs. 4 and 5). Tables 11 and 12 provide the corresponding development difficulty grading descriptions.

Figure 4
figure 4

Distribution map of the difficulty distribution of shallow (0 ~ 15 m) underground space resource development in the main urban area of Xianyang city.

Figure 5
figure 5

Distribution map of the difficulty of resource development of the subshallow (15 ~ 30 m) underground space in the main urban area of Xianyang city.

Table 11 Grading description of the difficulty of resource development of the shallow (0 ~ 15 m) underground space in the main urban area of Xianyang city.
Table 12 Grading description of the difficulty of resource development of the subshallow (15 ~ 30 m) underground space in the main urban area of Xianyang city.

According to Figs. 4 and 5, the prohibition area in the main urban area of Xianyang city is mainly located at the southern edge of the Loess Plateau north of the Wei River, which is distributed along a northeast–east direction, mainly in the cultural relic cemetery protection area. The main urban area of Xianyang city is difficult to develop, and the general area occurs as a northeast–east band, concentrated in the low-terrace areas on both sides of the Wei River in the southern part of the main urban area, which are more difficult to develop in terms of underground space construction due to the large thickness of local sand and gravel layers and the complex hydrogeological conditions. The areas with low development difficulties in the main urban area of Xianyang city are mainly distributed in the loess plateau area in the northern part of the city, the thick loess layer exhibits suitable engineering geological characteristics, and the difficulty of underground space development and construction is relatively low.

As an urban environmental geotechnical engineering activity, the development of urban underground space resources is difficult due to regional engineering geological conditions and urban environmental constraints. According to Tables 11 and 12, the difficulty of resource development of the shallow underground space in the main urban area of Xianyang city is affected by rock and soil problems, as well as the constraints of ecological protection requirements and existing building facilities, and the difficulty of resource development of the subshallow underground space is relatively high, mainly reflected in the obvious impact of groundwater on underground space development. In general, with increasing development level, the development difficulty correspondingly increases, and it is mainly controlled by the geological conditions.

Evaluation of the potential value of underground space resources

By weighting and integrating the distribution of social factors such as location, population, economy, transportation, and land use type in the main urban area of Xianyang city, a potential value distribution map of the underground space in the main urban area of Xianyang city is obtained (Fig. 6).

Figure 6
figure 6

Classification map of the potential value of the underground space resources in the main urban area of Xianyang city.

As shown in Fig. 6, the area with a high potential value of the underground space in the main urban area of Xianyang city is the expansion area around the core area and the core area on the south bank of the Wei River, which exhibits the highest degree of development in terms of the population density, economic level, transportation and location conditions. The area involved covers 17.99 km2, of which the area north of the Wei River is the expansion of the core area, mainly involving residential and commercial areas around the urban core area, while the area south of the Wei River is mainly distributed at the centre of Century Avenue. The area was developed late, the planning concept is relatively new, and it provides a higher development potential by virtue of the advantages of its geographical location. The area with a high potential value mainly involves residential and industrial land around the northern urban area of Weihe, covering a total of 32.31 km2, affected by the distribution of industrial plants and railway barriers. The area has a low population density and lacks large-scale commercial office facilities. Considering the parking and civil defence needs of the residential area and storage use in the industrial zone, the area still exhibits high development potential. The area with a general potential value is mainly located in the built-up area north of the Wei River, the southern extension area of the Wei River and the surrounding area of the airport, covering a total of 95.05 km2. The area has not yet been built up on a large scale, and many areas are vacant. Considering that the future development plan Xianyang mainly aims to expand the city to these areas, there is also a certain development value of the underground space. The area with a low potential value involves the 437.75 km2 area outside the built-up area of the city, which has not yet been developed on a large scale except for the airport area. The land mostly encompasses farmland and cultural relic sites, and the development value is low.

Comprehensive quality evaluation of underground space resources

With the use of the evaluation model and method described above, based on the spatial analysis function in ArcGIS, the development difficulty and potential value of underground space resources in the main urban area of Xianyang city were superimposed to generate a comprehensive quality grading map of the shallow to subshallow underground space resources (Figs. 7 and 8).

Figure 7
figure 7

Comprehensive quality grading map of the shallow (0 ~ 15 m) underground space resources in the main urban area of Xianyang city.

Figure 8
figure 8

Comprehensive quality grading map of the subshallow (15 ~ 30 m) underground space resources in the main urban area of Xianyang city.

The comprehensive quality grading map of the shallow to subshallow underground space resources in Xianyang city clearly reflects the comprehensive resource endowment in the main urban area, among which the core areas such as Renmin Road, Weiyang Road and Century Avenue cover a total area of 21.52 km2, and the comprehensive quality grade of the underground resources in the shallow area (0 ~ 15 m) is the highest, while the area with the highest grade in the subshallow area (15 ~ 30 m) is reduced to 4.37 km2 due to the increase in construction difficulty. The comprehensive quality grade of the underground resources in the surrounding extension area of the urban core area and the airport area is high, and the shallow and subshallow areas are 55.98 km2 and 65.49 km2, respectively. Affected by the adverse geological conditions and location distribution, the comprehensive quality grade of the core area and the extended area around the Weihe River is low, and the areas in the shallow and subshallow areas under this grade are 110.62 km2 and 112.32 km2, respectively. The comprehensive quality grade of the underground space resources in the Beiyuan area except the airport and other surrounding areas is the lowest, and the areas in the shallow and subshallow areas are 327.22 km2 and 346.71 km2, respectively. The cultural relic sites, water source protection areas and ecological reserves in the planning area are prohibited development areas due to the requirements of relevant norms and protection, with a shallow restricted area of 67.76 km2 and a subshallow restricted area of 54.21 km2.

Discussion

The development of the underground space in the main urban area of Xianyang city is in its infancy, and the amount of underground space resources available for development and utilization is considerable. According to the difficulty and potential value at each level of development, the comprehensive quality at the different levels is considered in an overall manner, and the rational selection and planning of development locations can greatly improve the efficiency of underground space utilization. From the perspective of the development difficulty, it is necessary to consider the comprehensive geological conditions of the region, which can affect the rationality, safety and economy of underground space development, construction, operation and management. From the perspective of the potential value, the population density, economic conditions, urban functional layout and location conditions undoubtedly have different demands for underground space development and construction, which affects the social value and operational benefits of underground space development. On the basis of considering the comprehensive quality at the different levels, development sequence determination and location selection of regional underground space locations can be achieved.

In the main urban area of Xianyang city, the comprehensive quality level of shallow underground space development is general and above, reaching approximately 188.12 km2, accounting for 32.26% of the total quantity, and the comprehensive quality of the underground space resources of Weiyang West Road, Renmin Road Street and Weiyang Street in Xianyang city is very high, which exhibits both the urgency of development and favourable geological conditions. There are slightly better geological conditions and some areas with a higher demand in the northeastern area of Maquan Street, part of Gudu Street and Xixian International Airport in Dizhang Street, so the comprehensive quality of resources in this area is also high, and the above areas can be considered during shallow development.

The comprehensive quality level of subshallow underground space development in the main urban area of Xianyang city is general and above, reaching approximately 182.18 km2, accounting for 31.24% of the total quantity, and the comprehensive quality of the underground space resources of Weiyang West Road, Renmin Road Street and Weiyang Street in Xianyang city remains high, but due to geological environment problems, subshallow development is impacted. There are favourable geological conditions and areas with certain needs in the northeast of Xiwu Street, the northern area of Maquan Street, some areas of Gudu Street and Xixian International Airport in Dizhang Street, so the comprehensive quality of the resources in this area is also high, and subshallow development should be considered and planned in these areas.

Conclusion

There are many factors that affect the comprehensive evaluation of urban underground space resources, and these factors are multilevel and complex. In this paper, we used the analytic hierarchy process to comprehensively evaluate underground space resources as a whole, organically combining qualitative and quantitative methods. We decomposed complex systems and transformed multicriteria decision-making problems that are difficult to fully quantify into multilevel single-objective problems for obtaining reasonable and clear evaluation process levels, and the results were simple and clear. Regarding comprehensive evaluation of the underground space resources in the main urban area of Xianyang city, the following conclusions can be drawn.

In terms of quantity, at present, the underground space resources available for development in the main urban area of Xianyang city account for approximately 25.11% of the total development quantity, the underground developable area is approximately 82.3 km2, and the shallow underground space resources (30 m) available for development in the main urban area reach 2.465 billion m3, accounting for approximately 79.5% of the total resources of the shallow underground space. The potential for development and utilization is enormous, underground space development mainly involves shallow depths, the degree of subshallow development is low, and the underground space development level is still relatively low.

In terms of the comprehensive quality, first, the main urban area of Xianyang city was evaluated in regard to the difficulty of underground space resource development, and the notable influences of the topography and geomorphology, rock and soil mass, structural stability and groundwater occurrence characteristics on the development of the underground space at different depths in Xianyang were studied and summarized. The specific role and spatial characteristics of the impact distribution were determined. The development of subshallow underground space resources is relatively difficult, mainly reflected in the obvious impact of groundwater on underground space development. Second, based on socioeconomic influencing factors such as population, traffic conditions, location conditions, land price and land use type, the potential value of underground space resources was evaluated, and the potential value distribution of the underground space in the study area was obtained. The areas with a high potential value were the areas along Renmin Middle Road and the middle of Weiyang Road, the core area of the south bank of the Wei River, and the residential and industrial land areas around the northern urban area of the Wei River, while the other regional areas exhibited a low degree of development and low potential value. Finally, according to the evaluation results for the difficulty of underground space resource development and the evaluation results for the potential value of underground space resources, the comprehensive quality of shallow and subshallow underground space resources in the main urban area of Xianyang city was graded. The comprehensive quality level of the underground space resources in core areas such as Renmin Road, Weiyang Road and Century Avenue was the highest, the comprehensive quality of the underground space resources in the expansion area around the urban core area and the airport area was the second highest, and the comprehensive quality level in the area outside the core area and the extension area around the Weihe River was average. The comprehensive quality level of the underground space resources in the Beiyuan area, except for the airport area and other surrounding areas, was low.