Application of Variable Weight Theory in the Suitability Evaluation of Regional Shallow Geothermal Energy Development

: Blindly developing shallow geothermal energy (SGE) may lead to di ﬃ culties in reinjection, a low coe ﬃ cient of performance (COP), and the waste of funds. In order to avoid these problems and improve development e ﬃ ciency, it is necessary to scienti ﬁ cally conduct a development suitability evaluation and classify the suitability levels before development. This article takes Chengde City, Hebei Province, China as an example and constructs an evaluation index system based on the analysis of regional basic conditions. The Analytic Hierarchy Process (AHP) was used to calculate the constant weight. A Fuzzy C-means Clustering (FCM) analysis was used to determine the critical value of evaluation index classi ﬁ cation and the interval threshold of variable weights. The parameters of the variable-weight model were calculated using the method of model backcal-culation. Based on the constructed variable-weight model, the suitability of SGE development for groundwater sources in the study area was evaluated, and it was divided into ﬁ ve levels: most suitable area, suitable area, relatively suitable area, less suitable area, and unsuitable area. Through a veri ﬁ cation analysis of engineering examples and a comparison with the evaluation results under traditional constant-weight models, it was found that the evaluation results based on variable weight theory have be tt er data discreteness and a higher accuracy. Research has shown that varia-ble-weight models can adjust the weight of each index based on its state value. Through this study, the accuracy of the suitability evaluation for regional SGE development can be improved. This can provide a certain reference for the suitability evaluation of SGE development in other regions.


Introduction
Shallow geothermal energy (SGE) refers to the energy that exists within a depth range of 200 m below the Earth's surface, including in the rock, soil, groundwater, and surface water.It mainly comes from solar energy, which radiates from the Earth's surface.Surface soil and rocks absorb heat and conduct it downwards, forming SGE [1].This is significantly different from deep geothermal energy, which is renewable thermal energy from deep within the Earth, originating from molten magma and the decay of radioactive materials [2].The development and utilization of SGE can be divided into three modes: buried-pipe ground source heat pump system, groundwater ground source heat pump system, and surface-water ground source heat pump system [3].This article mainly discusses the development suitability of groundwater ground source heat pump systems.As a kind of clean and renewable energy, SGE has the characteristics of a wide distribution, large reservoirs, high efficiency, energy saving, and no pollution [4].It can serve as an important energy source for promoting clean heating or cooling worldwide.In addition, it can also contribute to China's achievement of the "30 carbon peak" and "60 carbon neutrality" goals [5].It is well-known that the initial investment in the development of SGE is high; the blind development of ground source heat pump technology to develop SGE can not only not achieve the energy-saving effect, but also cause high operating costs due to the lower COP [6].Taking Chengde City of Hebei Province as an example (Figure 1), according to the survey statistics, there are more than 70 SGE development and utilization projects in the urban planning area, and the groundwater ground source heat pump system is the main one.Due to the lack of scientific site selection based on a suitability evaluation in the early stage, there was no targeted development and utilization plan design based on different suitability levels.As a result, some projects are built in less suitable or unsuitable areas; the development and utilization design does not match the hydrogeological conditions.After the project was put into operation, there were problems such as the poor operational efficiency, difficulty in recharging, and even inability to operate.These seriously constrain the development and utilization of SGE.Therefore, in the early stage of developing SGE, it is necessary to establish a suitability evaluation index system and model according to regional characteristics, delimit suitability zones, and macroscopically grasp the suitability of regional SGE development [7].The suitability of SGE development is a complex nonlinear problem.It is related to basic conditions such as the aquifer water-richness, hydraulic conductivity, groundwater mineralization, average groundwater temperature, recharge capacity, and aquifer thickness [8].In addition, it is closely related to the economy.In order to scientifically evaluate the suitability of SGE development, it is necessary to consider these factors comprehensively and use systematic methods and theories for a comprehensive evaluation [9].At present, there are many studies on the suitability evaluation methods for the development of SGE in China.The evaluation method first involves constructing evaluation indices tailored to local conditions.Then, evaluation methods such as the AHP, Fuzzy Comprehensive Method, Entropy Weight Method, Grey Relational Degree Method, and Expert Questionnaire Method are used to calculate weights.On the basis of determining the weights, the spatial information fusion function of the GIS platform is utilized, and then the development suitability is determined by the weighted superposition method [11][12][13][14][15].A series of research achievements have also been made in the suitability evaluation of SGE development abroad.For example, developing a computational program in a GIS environment to derive a regional model that can combine physical and economic variables to analyze the potential for SGE development in the Marche region of central Italy [16].By developing a method called G.POT, the suitability of developing SGE in a region for building heating and cooling is evaluated [17].An evaluation was conducted on the appropriate areas for introducing closed-loop ground source heat pump (GSHP) systems in standard independent residences in Akita Plain, Japan, using the drilling depth required for the GSHP system as a suitability index [18].We then assess the installation potential of different types of GSHP systems in the Aizu Basin (Japan) based on groundwater conditions [19].Based on the specific actual situation in each region, experimental research and the establishment of index systems are carried out to determine the influencing factors and use various research methods to determine relevant indices.Suitable evaluations for SGE development are carried out [20,21].On the basis of meeting the requirements for the suitability evaluation of SGE and saving mining capital, the evaluation results are more focused on practical applications.Especially in the evaluation process, future urban development planning, regional development, and local funding and technical support were considered [22][23][24].It can be seen that a series of achievements have been made in the evaluation of the suitability of SGE development.However, the current evaluation method is mainly based on the constant-weight evaluation method.If the value of a certain index within a certain evaluation unit is too low, its final evaluation value may be neutralized by the high value of other indices, causing the evaluation result to be inconsistent with reality.On the contrary, if a certain index value is too high, its final evaluation value may be neutralized by the low values of other indices, which can also cause the evaluation results to be inconsistent with reality.It can be seen that the use of fixed and unchanging weights to reflect the different state values of evaluation indices has certain limitations on the suitability of SGE development.This is mainly manifested in the inability to highlight the impact of changes in evaluation index values on evaluation results, regardless of the use of any constant-weight model.Especially when there is a significant change or mutation in the status value of evaluation indices, the limitations of the constant-weight evaluation become more prominent [25,26].
In order to overcome the limitations caused by fixed weights, the weight values can be adjusted based on the state values of each evaluation index within each evaluation unit, thus highlighting the impact of excessively high or low evaluation index values on the evaluation results, achieving more realistic evaluation results and improving the evaluation accuracy.This article introduces the theory of variable weights for the first time.We then systematically study the suitability evaluation method for SGE development based on variable weight theory from the aspects of establishing evaluation indices, quantifying indices, calculating constant-weight weights, and constructing variable-weight models.Taking Chengde City (Figure 1), Hebei Province as the research object, this study compares the differences in results between the constant-weight evaluation and variableweight evaluation methods.We then analyze the advantages of variable-weight evaluation models in terms of data discretization and evaluation accuracy.I hope to provide some reference for the research on the suitability evaluation of SGE development and a scientific basis for its development and utilization through conducting this research.
The study area is located in the northeast of Hebei Province, with an area of approximately 39,500 km 2 .It is situated in the transitional zone between North China and Northeast China, adjacent to Beijing and Tianjin to the south, Inner Mongolia to the north, and Liaoning Province to the east (Figure 1).The study area belongs to the temperate semihumid and semi-arid continental monsoon mountain climate, which is a transitional area from the temperate zone to the warm temperate zone.It has the characteristics of four distinct seasons, synchronous rain and heat, and a large temperature difference between day and night.The average annual rainfall is 648 mm, the average temperature is 9.5 °C, and the annual frost-free period is about 160 days.
The study area has well-developed faults and fold structures, which have gone through three main stages: basement formation, cap rock development, and strong activity.The exposed strata in the study area are relatively complete, ranging from the Archean to the Proterozoic, and also from the Mesozoic to the Cenozoic.Only the Paleozoic lacks the Silurian and Devonian strata, as well as the adjacent Upper Ordovician and Lower Carboniferous strata (Figure 2).The groundwater in the study area is strictly controlled by the geotectonic and geomorphic units.From north to south, it is divided into two hydrogeological regions: the Inner Mongolian Plateau and the northern Hebei mountainous region.According to the type of groundwater and its storage medium, the aquifer in the study area can be divided into four different types: loose rock pore water, clastic rock fracture pore water, bedrock fracture water, and carbonate rock fracture cave water (Figures 3 and 4).Among them, the first three types are widely developed in this area, while the fourth type has a relatively small distribution range.Within a depth range of 200 m, the aquifer has a thickness of 16-149.17m, a hydraulic conductivity of 0.0132-13.03m/d, a groundwater mineralization degree of 150.07-1593 mg/L, and an average groundwater temperature of 4.6-15.1 °C.

Evaluation Index System Construction
The suitability of regional SGE development is a huge and complex system, which is composed of many relatively independent and interrelated parts; reasonable evaluation index is very important to the evaluation result [29,30].On the basis of the comprehensive research results of experts and scholars in the field of SGE development in recent years, the evaluation index system is determined through the analysis and research of the existing basic data in the study area and the characteristics of the regional geological environment.The evaluation index system is divided into two levels.The first-level indices are divided into geological and hydrogeological conditions, groundwater conditions and drilling conditions; the second-level indices included in each first-level index are shown in Table 1.

Primary Index
Secondary Index Quantitative Basis Geological 1 and hydrogeological conditions 2 The water-richness of aquifers Based on the degree of regional water-richness Hydraulic conductivity Based on the value of hydraulic conductivity Recharge capacity Based on the regional stratigraphic lithology Aquifer thickness Based on the value of aquifer thickness

Groundwater condition
Mineralization of groundwater Based on the value of groundwater mineralization Average temperature of groundwater Based on the value of regional groundwater temperature

Drilling condition Drilling difficulty
Based on the regional drilling difficulty Notes: 1 This refers to various factors that affect the geological phenomena on the Earth's surface, including stratigraphic lithology, geological structure, magmatic activity, etc.This article mainly involves stratigraphic lithology. 2This refers to the general term for conditions related to the formation, distribution, and variation of groundwater, including groundwater recharge, burial, runoff, discharge, water quality, and quantity.This article mainly involves the water-richness, hydraulic conductivity, and thickness of the aquifer.
For the definition or related description of geological and hydrogeological conditions and their secondary indices, please refer to Appendix A. For the definition or related description of groundwater conditions and their secondary indices, please refer to Appendix B. For the definition or related description of drilling conditions, please refer to Appendix C.

Quantification of Evaluation Index
After the establishment of the evaluation index system, it is necessary to conduct scientific and reasonable quantification of each evaluation index.Among all the indices, the hydraulic conductivity, aquifer thickness, mineralization of groundwater, average temperature of groundwater, and drilling difficulty are directly collected measured data or calculated data, without quantitative processing.For the water-richness of aquifers and recharge capacity, in order to facilitate later evaluation and calculation, the same quantified value interval is adopted for each of these two evaluation indices.Combined with the actual situation in the study area, the grading standard and quantified value of the evaluation indices are proposed in this article.The details are shown in Table 2. Notes: 1 This refers to a relatively impermeable rock layer that separates groundwater.It cannot provide water, nor can water pass through it. 2This refers to water that can pass through it, but its water holding capacity is poor and it cannot provide water. 3This refers to q ≤ 0.1 L (s.m); q is a unit water inflow, with a unit water inflow based on a diameter of 91 mm and a pumping depth of 10 m. 4 This refers to 0.1 < q ≤ 1 L (s.m). 5 This refers to q > 1 L (s.m).
After the evaluation index quantization, in order to eliminate the influence of each evaluation index dimension, it is necessary to normalize the quantization value of each evaluation index.The six evaluation indices of the water-richness of aquifers, recharge capacity, hydraulic conductivity, average temperature of groundwater, aquifer thickness, and drilling difficulty are positively correlated with the suitability of SGE development.In this article, the maximum value method (Formula (1)) was adopted for normalization.As for the mineralization of groundwater, it is negatively correlated with the suitability of SGE development; the minimum method (Formula (2)) is adopted for normalization treatment.
where Aij is the normalized quantized value; Xij is the index value before normalization; Xmaxj is the maximum value in the quantized data before normalization; and Xminj is the minimum value in the quantized data before normalization.
After the evaluation indices are quantified and normalized, thematic layers of each evaluation index are constructed by using ArcGIS spatial interpolation function (Figure 5).

Calculation of Constant Weight
Based on the evaluation index system of SGE development suitability, the AHP is used to calculate and determine the weight of each main control factor.The specific calculation process is as follows.Firstly, the evaluation system structure model is divided into three levels (Figure 6), namely, target layer A (suitability for SGE development), attribute layer B (including geological and hydrogeological conditions, groundwater conditions, and drilling conditions), and element index layer C (aquifer richness, recharge capacity, hydraulic conductivity, average groundwater temperature, groundwater mineralization, aquifer thickness, and drilling difficulty).Secondly, the importance of the indices in this layer compared pairwise with those in the previous layer was evaluated using expert scoring methods, and a judgment matrix was constructed (Tables 3-6).Again, we calculate the maximum eigenvalue of each judgment matrix (λmax), consistency index (CI), and consistency coefficient (CR).Here, the calculation process of these parameters is as follows.Firstly, we calculate the product (Mi) of each row element in the judgment matrix according to Formula (3).Secondly, we calculate the n-th root of Mi according to Formula (4).Thirdly, we construct the vector w according to Formula (5) and normalize the vector according to Formula (6).Fourthly, we calculate λmax according to Formula (7).On this basis, we calculate CI according to Formula (8) and CR according to Formula (9).If CR < 0.1, it is considered that the judgment matrix has passed the consistency test; otherwise, it does not have satisfactory consistency [31].Taking the judgment matrix A~Bi as an example, with a matrix order of 3, we calculate it according to the above formula λmax = 3.0735; the index weight are shown in Table 7, with a CI of 0.0368.When the matrix order is 3, the corresponding RI is 0.58, and the CR is 0.0707.According to this method, the final weight calculation results of each index are shown in Table 7.The CR is less than 0.1; that is, it passes the consistency test.
where Mi is the product of the i-th row elements in the matrix, bij is the j-th column element in the i-th row, and n is the matrix order.
where W -i is n-th root of Mi, Mi is the product of the i-th row elements in the matrix, and n is the matrix order.
[ ] where, w is a vector, W -i is n-th root of Mi, and n is the matrix order.
where wi is the weight value of the evaluation index, W -i is n-th root of Mi, and n is the matrix order.
where λmax is the maximum eigenvalue of the judgment matrix, A is the judgment matrix, and W is the eigenvector corresponding to λmax.wi is the weight value of the evaluation index, and n is the matrix order.
where CI is the consistency index, λmax is the maximum eigenvalue of the judgment matrix, and n is the matrix order.

RI CI
where CR is the consistency coefficient, and RI is the average random consistency index, which can be obtained by consulting relevant data tables.Note: λmax = 1; CR = 0 < 0.1; and CI = 0 < 0.1.

Construction of Variable-Weight Evaluation Model
The basic principle of variable weight theory is that the weight of evaluation index will be adjusted based on the evaluation index status value of each evaluation unit.Therefore, the weights can better serve the role of corresponding evaluation indices in the evaluation system.It can change the fixed weight in traditional constant-weight evaluation, making comprehensive decision-making more reasonable and scientific [32].Based on the theory of variable weight, it is necessary to first construct a state vector and determine the variable-weight vectors for the penalty interval, non-incentive and non-penalty interval, initial incentive interval, and strong incentive interval.The segmented function of the variable-weight vector used in this study is shown in Formula (10).Secondly, the critical value and variable-weight interval of the evaluation index classification are calculated.Finally, on the basis of meeting the requirements of comprehensive suitability evaluation, the appropriate parameter values in the state vector that meet the expected weights are calculated in reverse.For the determination of the parameters of the state variable-weight vector, this article cites previous research to determine the mathematical function of state variable-weight vectors [33].
where c, a1, a2, and a3 are the weight adjustment parameters of the state variable-weight vector; and dj1, dj2, and dj3 are the threshold of variable-weight interval for the j-th factor.

The Determination of Variable-Weight Intervals
Before determining the variable-weight parameters, it is necessary to divide the weight adjustment range of each evaluation index, which is used to define the weight adjustment amplitude of the main control factor index value.There are three steps required here.Firstly, it is necessary to use statistical analysis methods to determine the classification threshold values for each index.Secondly, it is necessary to calculate the threshold of the variable-weight interval for each index based on the calculation of the classification threshold.Finally, the variable-weight interval of each index can be determined based on the threshold of the variable-weight interval [34].For the calculation of critical values for each index classification, the FCM was used in this study.Firstly, we extract clustering sample data from the ArcGIS database.Secondly, we determine the fuzziness C (number of clustering layers) as 4, the weighted index m as 2, and the iteration error e as 1 × 10 −4 .Thirdly, FCM is implemented by input Python, and each evaluation index value is classified by calculation.Finally, determine the classification threshold fi for index values, as shown in Table 8.In this table, "f1, … f6" represents the classification threshold values of each index after normalization, rather than the actual values of each index.Taking the hydraulic conductivity as an example, its actual value ranges from 0.0132 to 13.03 m/d, and, after normalization, its range is 0 to 1.The FCM analysis was used to determine its classification critical values as 0.0833, 0.1042, 0.2083, 0.2292, 0.3958, and 0.4375, respectively, which are the data shown in Table 8.
On the basis of determining the classification critical value, we calculate the weight threshold of variable weights according to Formula (11): where fj is the classification critical value of the index value of the j-th factor; and dj is the variable-weight interval threshold for the j-th index.Therefore, as shown in Table 9, variable-weight threshold intervals for each master control factor are obtained [35].After determining the variable-weight intervals for different indices, it is necessary to calculate the adjustment parameters c, a1, a2, and a3 of the constructed variable-weight vector.The calculation of these parameters mainly adopts the method of model backcalculation, and the specific process includes three steps.Firstly, selecting an evaluation unit requires that the four index values should be located in four different variable-weight intervals.In addition, one of the remaining index values needs to be within the penalty interval.Secondly, based on objective facts and the preferences of decision-makers, the AHP is used to calculate the ideal variable weights of each evaluation index within the evaluation unit.Finally, we substitute the values of various indices, constant-weight values, and ideal variable-weight values within the evaluation unit into the variable-weight vector.Based on this, we establish an equation system about the parameters to be solved, and then calculate the parameter values [36].The parameter values are shown in Table 10.On the basis of determining the variable-weight interval and variable-weight parameters, the constructed variable-weight evaluation model is as follows [37]:  (13) where j = 1, 2, 3, 4 ….. 7; dj1, dj2, and dj3 are the thresholds for the variable-weight interval of the j-th main control factor; and Wi o is the constant weight of the evaluation index.

The Evaluation Results
Based on the obtained variable-weight parameter values, we substitute them into the MATLAB program, and, finally, obtain the suitability index values within each evaluation unit according to Formulas (13) and (14).The suitability index calculation data are linked to the layer attribute table through the attribute library.Based on the ArcGIS platform, the classification threshold is determined according to the natural classification method.The suitability of SGE development for groundwater sources in the study area is divided into five levels: most suitable area, suitable area, relatively suitable area, less suitable area, and unsuitable area.Based on the above evaluation methods, the suitability evaluation of SGE development in the study area was completed, and an evaluation subarea layer was formed (Figure 7).The statistical results of the evaluation are shown in Table 11.From the evaluation results, it can be seen that, within the study area, the most suitable area for the development of SGE with underground water sources is 5848.1 km 2 , accounting for 14.82%; the suitable area is 9794.5 km 2 , accounting for 24.82%; the relatively suitable area is 16,775.9km 2 , accounting for 42.50%; the low-suitability area is 6809.3km 2 , accounting for 17.25%; and the unsuitable area is 241.3 km 2 , accounting for 0.61%.The proportion of less suitable and unsuitable areas is about 17.86%, and the proportion of the most suitable and suitable areas is about 39.64%.

The Result Analysis
According to the evaluation results, the most suitable areas for the development of SGE in the study area are mainly distributed in the central and western parts of Xinglong County, the western part of Kuancheng County, the southwestern part of Pingquan City, the southern part of Shuangqiao District, the southwestern part of Weichang County, and the northwestern part of Fengning County.This area has a good water-richness, large aquifer thickness, and strong recharge capacity, making it an ideal area for the development of SGE based on groundwater sources.
According to the evaluation results, the suitable areas for the development of SGE based on groundwater sources in the study area are distributed in the southern and eastern parts of Xinglong County, the central and southern parts of Luanping County, Chengde County, Kuancheng County, the northeastern part of Shuangqiao District, the southern part of Shuangluan District, the northeastern and northwestern parts of Longhua County, the southeastern part of Pingquan City, and the central part of Weichang County.These areas have a good recharge capacity, high average groundwater temperature, and easy drilling, making them suitable for the development of SGE based on groundwater sources.
According to the evaluation results, the relatively suitable areas for the development of SGE based on groundwater sources in the study area are distributed in the northwest and southeast of Weichang County, central and northern parts of Fengning County, southwest and eastern parts of Longhua County, northern parts of Shuangluan District, northern parts of Shuangqiao District, eastern and northern parts of Chengde County, western and northern parts of Pingquan City, southwestern and southeastern parts of Luanping County, Yingzi District, and northern parts of Xinglong County.These areas have a moderate difficulty in drilling, average groundwater temperature, and water-richness, and belong to the more suitable areas for the development of SGE based on groundwater sources.Feasibility studies need to be conducted before development, and development plans should be designed according to local conditions.
According to the evaluation results, the less suitable areas for the development of SGE based on groundwater sources in the study area are distributed in the northeast of Weichang County, the central and southern parts of Fengning County, the western and northern parts of Luanping County, the western parts of Shuangluan District, the southern parts of Longhua County, and the central and eastern parts of Pingquan City.These areas have a high groundwater mineralization, poor recharge capacity, and difficult drilling, and belong to the less suitable areas for the development of SGE based on groundwater sources.Before development, a detailed hydrogeological exploration is required in the engineering construction area to fully demonstrate the feasibility of development and minimize development risks.
According to the evaluation results, the unsuitable areas for the development of SGE based on groundwater sources in the study area are mainly distributed in the central part of Pingquan City and Weichang County, while other areas have smaller distribution areas.These areas are not suitable for the development of SGE based on groundwater sources, and should be avoided during the site selection stage of the project.

The Comparison of Evaluation Results
Based on the calculated constant weight, a constant-weight evaluation model (Formula ( 14)) is used to evaluate the suitability of SGE development in the study area.We determine the classification threshold based on the natural breakpoint classification method on the ArcGIS platform.The suitability of SGE development is divided into five levels: most suitable area, suitable area, relatively suitable area, less suitable area, and unsuitable area (Figure 8).The evaluation results are shown in Table 12 and compared with the evaluation results of the variable-weight model in Table 13 and Figure 9.  ) where V is the suitability index for SGE development (the larger the value is, the more suitable it is for development), which is the weighted total score of all evaluation factors.Wj is the weight of the evaluation indices, Xj is the standardized value of the indices, and m is the number of evaluation indices.By comparison, it can be concluded that the area of the low-suitability area has the largest change, followed by the relatively suitable area, while the area of the suitable area has the smallest change, followed by the most suitable area.The evaluation results based on the variable-weight model show that the most suitable area, suitable area, and relatively suitable area are all larger than the evaluation results of the constant-weight model, while the less suitable area and unsuitable area are smaller than the evaluation results of the constant-weight model.It can be seen that the evaluation results based on the variableweight model are more inclined towards suitability.
By comparing the main distribution areas of the results based on the two evaluation models, it can be seen that there is a significant difference in their distribution areas (Table 14).Taking the most suitable area as an example, the evaluation results based on the constant-weight model show that it is mainly distributed in Xinglong County, Kuancheng County, and Chengde County.The evaluation based on the variable-weight model shows that it is mainly distributed in Weichang County, Xinglong County, and Fengning County.It can be seen that, for the spatial distribution of the most suitable area, the results of the two models are completely different.Similarly, the same applies to the main spatial distribution of unsuitable areas.The evaluation results based on the constant-weight model show that it is mainly distributed in Fengnian County, Weichang County, and Longhua County.The evaluation based on the variable-weight model shows that it is mainly distributed in Weichang County, Longhua County, and Chengde County.This is also a big difference.

The Discrepancy Analysis of Evaluation Results
This article uses the standard deviation ellipse tool in ArcGIS to analyze the discreteness of evaluation results.This tool can analyze the directionality and discreteness of data, where the long axis of the ellipse represents the directionality of the data, and the short axis represents the discreteness of the data.The longer the short axis, the better the discreteness of the data, and vice versa.From Figure 10 and Table 15, it can be seen that among the five evaluation levels, the elliptical short axis lengths of the most suitable, suitable, and less suitable areas based on the variable-weight model are all longer than those based on the constant-weight model.The elliptical short axis lengths of the more suitable areas are relatively close, while the elliptical short axis lengths of the unsuitable areas differ greatly.However, due to the small unsuitable area of the variable-weight evaluation results (only 0.61%).Therefore, from an overall analysis, it is found that the evaluation results based on the variable-weight model have better discrete data.This is because the variable-weight model starts from the whole to the local, fully considering the changes in various parameters within a certain evaluation unit, and allocates different evaluation index weights for each evaluation unit based on the evaluation nature and purpose, which is more in line with reality.

Conclusions
This article takes Chengde City, Hebei Province as the research object, analyzes the basic conditions, and constructs an evaluation index system and quantitative methods.AHP, FCM, and model inverse algorithms were used to calculate relevant parameters and establish a variable-weight evaluation model to evaluate and partition the suitability of SGE development.The main conclusions were as follows: (1) Compared with the evaluation results of the constant-weight model, the evaluation graph based on the variable-weight model has relatively good discreteness, and the variable-weight evaluation model can adjust the weights of each evaluation index based on the index state values.It can meet the preferences of decision-makers for evaluation indices in different combination states, thus overcoming the limitations of fixed weights in traditional constant-weight evaluation models.(2) Through the verification and analysis of existing projects, the evaluation results based on the variable-weight model have a higher accuracy.The evaluation method based on the variable-weight theory can accurately reflect the suitability of SGE development in different regions of the study area, and provide a reference for its scientific development and utilization.
reverse case, the less suitable it is [38,39].The recharge capacity is greatly affected by geological conditions, and the recharge capacity is different in different rock strata.On the basis of a 1:200,000 hydrogeological survey in the early stage, this article classifies the recharge capacity of different rock formations by consulting experts and the literature.

Aquifer thickness
The thickness of the aquifer refers to the total thickness of the aquifer in the soil layer or rock layer.It mainly affects the water yield of a single well.The greater the aquifer thickness is, the more suitable it is for exploiting the SGE of an underground water source, and, when it is the reverse case, the less suitable it is.The aquifer thickness data in this article are collected from the hydrologic exploration boreholes in the study area.

Figure 2 .
Figure 2. Brief geological map of the study area [27].

Figure 5 .
Figure 5. Thematic map of suitability evaluation indices for shallow geothermal energy (SGE) development.

Figure 6 .
Figure 6.The AHP structural model based on evaluation index system.

Figure 7 .
Figure 7. Suitability evaluation map for SGE development in the study area based on variableweight model.

Figure 9 .
Figure 9.Comparison chart of evaluation results based on constant-weight evaluation model and variable-weight evaluation model.

Figure 10 .
Figure 10.Comparison of standard deviation ellipses of evaluation results based on constant-weight evaluation model and variable-weight evaluation model.

Figure 11 .
Figure 11.Local comparison of evaluation results based on constant-weight evaluation model and variable-weight evaluation model.

Table 1 .
The main controlling factors affecting the suitability of shallow geothermal energy (SGE) development.

Table 2 .
Evaluation index grading and evaluation criteria.

Table 7 .
Evaluation index weight.

Table 8 .
Evaluation index classification critical value (normalized processing).

Table 9 .
Variable-weight interval based on FCM clustering.

Table 10 .
Adjustment parameters of state variable-weight vectors.

Table 11 .
Statistical table for suitability evaluation of SGE development in the study area based on variable-weight model.

Table 12 .
Statistical table for suitability evaluation of SGE development in the study area based on constant-weight model.

Table 14 .
Comparison table of main distribution regions of evaluation results for different models in the study area.

Table 15 .
Comparison table of short axis length based on different evaluation models.

Table 18 .
Comparison table of existing engineering evaluation results based on different evaluation models.