Coupling Empowerment Based Application of Extension Method for Geothermal Potential Assessment

Plenty ofmathematics researches provide feasibility to calculate the weights of geothermal controlling factors and have been applied in geothermal potential assessment. In order to avoid the disadvantages of subjective and objective weighting calculation methods, an extension theory integrated weighting method was put forward, by combining with the process of AHP and mean variance method. The improved method can reach an agreement on subjective understanding of impact factors’ roles and data-based calculation weights.Then by replacing the point values with intervals, the extension theory was used in classification of geothermal assessment, according to extension judgment matrix.The evaluation results showed perfect performance in classification of impact factors, especially in Wudalianchi area, where 10 out of 11 selected impact factors agreed well with the actual evaluation grades. The study can provide a guidance for primary stage of geothermal investigation including the impact factor selection, weights calculation for impact factors, and the factors’ classification in geothermal assessment.


Introduction
Due to fossil fuel energy such as oil, coal, and natural gas which may account for >85% of entire energy consumption, the issues of energy and environment still need to be given deep consideration even if this adverse trend has been alleviated for now [1,2].In China, a great deal of geothermal energy remains undiscovered or fully exploited, thus fossil fuel energy can be deemed to be a competitive one undoubtedly [3,4].Since the 1990s, geothermal resources have been widely applied in various fields, improved people's lives significantly, and promoted social development drastically [5][6][7][8].According to a report proposed by the Chinese Academy of Engineering (2013), the entire utilized shallow geothermal energy (5.0 × 10 3 mwt) with medium and low temperatures is expected to increase threefold by 2050.In addition, the installed capacity for generating high temperature geothermal energy especially for enhanced geothermal system (EGS) has been greatly increased significantly in the last few decades.Thus, considering the important role of geothermal resources, the processing treatments including the geothermal identification, impact factor selection, geothermal potential regions assessment, and utilization conditions evaluation should be performed in early stage exploration [9].
In general, the commonly applied methods for geothermal potential evaluation usually are fuzzy comprehensive evaluation method [10], Monticello method [11], the back propagation (BP) neural network [12], the weights of evidence (W of E) method [13], logistic regression model, and so forth [14] (Williams et al., 2008).In these methods, the impact factors are single fixed and some selected databases are not quantitative.In addition, some methods are not easy to implement because the used impact factors not only have different influencing roles of the surrounding conditions but also possess inner nonlinear relationship between the factor pairs [15][16][17].Thus, the impact factors' corresponding weights should be taken into account before modeling application.The assignment of weights which can be obtained in these methods usually embraces two cases: judging on the basis of the subjective experience or calculating coefficients directly according to the objective database of factors.Both the two cases resulted in an unreliability of weight for each impact factor because of the artificial uncertainty or data deficiency [18,19].Therefore, a coupling empowerment extension theory was proposed as follows where the limitations on statistical methods or knowledge driven methods can be avoided [20].By making full use of the subjective selection ability in the analytic hierarchy process (AHP) and the objective ability of digging out source data in the mean variance method, the obtained index weights for impact factors were calculated, reflecting both subjectivity and objectivity.Then, the singular data were also removed reasonably in data processing.Finally, the extension method combining the weights was put forward in assessment of the geothermal potential regions.

Methodology
2.1.Analytic Hierarchy Process.AHP was first proposed by Saaty [21], and the system methods have been successfully applied to make optimization decisions by decomposing complex problems into several multiple indicator levels [22,23].Thus, the AHP method can also be used here in decision system to determine the weights of impact factors with characteristics of various objective distributions, multicriteria evaluation, and a situation that is difficult to quantify.Undoubtedly, the method was applied in a straightforward manner and performed easily.A schematic diagram of AHP for geothermal potential assessment was shown in Figure 1, where each subjective weight for impact factors was calculated by the analytic hierarchy process.

The Mean Variance Method.
Mean variance method is a method for solving decision weight coefficient with multiindex [24,25], where the weight for each impact factor can depend on the relative dispersion degree of sample database.
In general, as an uncertainty open system, geothermal occurrences are controlled by the various factors and vary significantly across different regions and periods.Thus, it is difficult to accurately express the geothermal distribution only according to the traditional objective weighting method due to the uncertainty and deficient database [26,27].Thus, with the help of the mean variance theory, a progress description for conversion of objective weights was put forward as follows, where the weight expressed a favorability degree of the uncertainty geothermal assessment appropriately: First, set a standardized evaluation matrix as follows: where  , is standardized index for the th class of the th impact factors and  is standardized evaluation matrix.Second, calculate the average standardized index with the equation as follows: where   represents the average standardized index for th impact factors.
Third, calculate the mean variance values for each impact factor; the formula was proposed as follows: where   is the mean variance index for th impact factor.Hence, an objective weight value   for various influencing factors is determined and the equation was put forward as follows: where   is the objective weight value for th impact factor.

Coupling Empowerment Rule.
The abovementioned judgment matrix in AHP and the mean variance method were considered as subjective method and objective method, according to the subjective experience and entirely the law of objective data, respectively.Thus a coupling empowerment rule was proposed by coupling AHP with mean variance method, in order to reach an agreement on the intuitive understanding for the potential geothermal phenomena in the field and the law of objective survey data [28].
To generate coupling weights for impact factors, a detailed procedure was proposed as follows: First, set a distance function with two parameters as follows: where   is the obtained weight for th impact factor by the analytic hierarchy process,   is the calculated weight for th impact factor by the mean variance method, and (  ,   ) is deemed to be distance value th impact factor between weights   and   .Second, set the coupling weights as follows: where   is the coupling weight for th impact factor and combines a linear relationship between AHP and the mean variance method;  or  is the coupling coefficient of the two kinds of weights.Third, in order to reach an agreement on the subjective process and objective method, (5) was reset as follows to obtain the corresponding coupling coefficients: where " +  = 1" means the constraint condition of the coupling coefficient.The coupling coefficients  and  were calculated by combining ( 5) with (7), and then the integrated weight   for th impact factor was obtained by substituting (6) with the its corresponding coefficients  and .The coupling weight method can make full use of the attribute optimization ability of analytic hierarchy process and the objective ability of data mining of the mean variance method.Hence by using the method, the interference of human factors can be removed effectively and the singular data in data processing can be reasonably eliminated.The abovementioned calculation process of the coupling weights was performed in a schematic diagram Figure 2.

Extension Theory.
As a cross discipline of originality, the extension theory was first proposed by Wen [29], and the quantitative tool has been used to study the transformation process that included system's characteristics, connotation, denotation, feasibility, and logicality.Hence, a matterelement was considered as the logical cell of extension theory and described as a formula  = (  , , ), where the three fundamental elements of the matter-element , , and , respectively, were the names, characteristics, and attribute values.Matter-element can reflect the relationship between the quality and quantity with the equation  = ().Thus, a real axis was constructed to reflect this correlation properly.The domain ( 0 ,  0 ) between the point  0 and interval  0 = (, ) on the real axis was expressed as follows: The relationship between the distance ( 0 ,  0 ) and domain ( 0 ,  0 ) was proposed as follows, with the point  0 and the interval  0 = (, ): Correlation function is based on distance.It expands the qualitative description with "certain nature " to a quantitative description with "the degree of the nature ."The elementary dependent function was established which can make the correlation function be described with the formula as below: The procedure of the extenics classification algorithm for geothermal potential assessment was proposed as follows: First, define the matter-element for the geothermal potential assessment as follows: where   is the th rank of evaluated geothermal potential regions,   is the th character within the th rank, and   is the attribute value for th character of th rank.Second, define classical domain and the limited domain as follows: where   is the susceptibility degree.The characters here referred just to the evaluation indexes.  is a value scale th degree of th character, which is the classical domain for normalized evaluation index.Similarly, the limited domain was proposed as follows: where  denotes the overall partition of suitability grades and   is the limited domain with respect to character   .Third, define the matter-element to be evaluated.The evaluation information was collected by the kind of geothermal energy level and a matter-element was proposed as follows: where   is the matter-element for th impact factor,   is th the impact factors,   is the th interval value of th impact factor, the values are corresponding to the magnitudes, and   is the specific data collected from the geothermal evaluation.
Fourth, calculate the correlation degree of evaluation values for geothermal potential assessment.
The correlation   (  ) is the th factor of the th geothermal potential regions with respect to the geothermal category: where Fifth, calculate the correlation degree of the evaluation rank in geothermal evaluation.
The correlation of the geothermal potential   (  ) related to evaluation level is where   is the weighting distribution coefficient of the evaluation factor, and it was calculated by the formula below: Sixth, determine geothermal potential rank.After calculating the correlation   (  ) for impact factor   related to evaluation level with the formula   (  ) = max{  (  )}, the potential evaluation level for categories was obtained qualitatively.The value and the relationship of   (  ) can reflect the potential of the th impact factor and the degree of category.

Application with Coupling Weight-Based Extenics Method
3.1.Database Preparation.Plenty of hydrothermal activities in China including volcanos, hot springs, and hot water wells have been well studied, and the basic geological and heat flow databases for geothermal emergencies regions were also accurately collected and sorted out [30,31].In this paper, there are 11 famous volcanic areas and seven specific geological conditions were selected to evaluate the geothermal potential regions.The selected impact factors, respectively, were stratum and lithology named  1 , structure characteristics ( 2 ), volcanic activities ( 3 ), water chemical characteristics ( 4 ), geothermal gradient ( ∘ C/100 m;  5 ), earthquake activities ( 6 ), and hot spring temperature ( ∘ C;  7 ).The index of each impact factor was shown in Table 1.
The selected impact factors were put forward as follows: Stratum and Lithology ( 1 ).A stratum is the place to store the geothermal energy where the thermal storage and conductivity vary significantly in different lithologies, and only the good stratum and lithology can store rich geothermal resources.
Structure Characteristics ( 2 ).Faults have been well developed in areas with plenty of active structures, which provided channels for the underground hot water migration.Thus, the active structure is closely related to geothermal energy.
Volcanic Activities ( 3 ).Magmatic activities in the volcano regions can provide the heat source for the geothermal field; hence they are directly related to the existence of geothermal energy.
Water Chemical Characteristics ( 4 ).Hydrochemical characteristics of the geothermal fluid (ionic composition and mineralization degree) have inner correlations with the geological tectonic activities and the corresponding thermal reservoir.
Geothermal Gradient ( 5 ).Geothermal gradients were composed of lines contracted from the higher temperature; thus a higher temperature for land surface can indicate the more abundant geothermal resources.
Earthquake Activities ( 6 ).Seismic activities are closely related to the active structure; hence a large seismic level in history usually indicates the development of the fault structure and volcanic activities.
Hot Spring Temperature ( 7 ).The land surface temperature anomalies are deemed to be a reference indicator for the geothermal emergencies, where higher temperature means more possibility of geothermal occurrences.
All the selected impact factors were determined by qualitative methods, and the influence of each factor was divided into ranges.According to a previous experience and a large number of statistical data [31], the variation range of influencing factors is shown in Table 1.

Coupling Weight Calculation.
According to the statistics and analysis of a large number of data, the selected impact factors were divided into 3 grades by single factor method and the results were presented in Table 2.In order to calculate conveniently, it is necessary to standardize the index factors mentioned in Table 1 and the equation was proposed as follows: where   is the attribute value for th impact factor of th selected geothermal regions,   max is the maximum attribute value,   min is the minimum attribute value of  factor, and    is the standardized value for th impact factor of th study areas.The calculated results were put forward as shown in Table 3.
The abovementioned parameter  0 for geothermal potential regions was determined with a reference in Table 2, and then the coupling weights were obtained according to process of Figure 2. The results of comprehensive weights were presented as shown in Table 4.

Matter-Element Assessment.
According to the standardized data in Table 3, the classical domain matrix for matterelement and joint domain matter-element were structured as follows: where  01 ,  02 , and  03 are the classified domain matterelement matrix for high, medium, and low level geothermal potential regions, respectively;   is the joint domain matterelement matrix.Thus, take Wudalianchi geothermal potential region as an example; the calculated results about comparison of evaluation grade and actual results were presented as shown in Table 5.
In general, the coupling weight usually indicates an influential degree between geothermal emergencies and the controlling factor objectively.The extension correlation degree and the comprehensive weights were presented in Table 5.With the classification in the matter-element assessment ((21); Table 4), the results reflected the nature of the geothermal potential evaluation well, where areas with volcanic activities, tectonic activities, higher temperature anomalies, and other influential factors agree with geothermal potential regions.According to results in Wudalianchi (Table 5), geothermal occurrences have inner connections with a higher temperature and the water with high ion content.However, there are disadvantage conditions including Quaternary sediment whose thickness varies significantly, in which volcanic activities with a low frequency and simple structure without enough surrounding impact factors showed medium or low evaluation grade in geothermal assessment, where hot springs hardly were discovered.Finally, the comparison of calculated and actual results in Wudalianchi area was still considered as successful, since 10 out of 11 classifications for impact factors were the same.

Conclusions
In this paper, to overcome the subjective and objective disadvantages in weighting calculation, an extension theorybased method was proposed by coupling the process of AHP and mean variance method, in order to reach an agreement on the intuitive understanding for the selected geothermal impact factors and their objective database.The improved coupling weights were used in classification of geothermal potential assessment, by replacing the point values with intervals according to extension judgment matrix.The evaluation results reflected the actual situation well and can play an indicative role in the partition of the geological conditions that may have inner connections with geothermal emergencies.An example put forward in Wudalianchi region was used in analysis and verified, and the results showed perfect performance where 10 out of 11 impact factors agree with the actual evaluation grades.Therefore, the study can objectively provide a guidance for initial phase of geothermal investigation.However, the application of improved classification method in geothermal potential evaluation, especially in the aspect of parameter selection, should be given further improvement because the selected impact factors may not be considered completely.

Figure 1 :
Figure 1: Hierarchical structure model of geothermal potential.

Figure 2 :
Figure 2: Schematic diagram for calculating the coupling weights.

Table 1 :
Impact factor indicators and variation range.

Table 2 :
Dimensionless quantization and standardized index level for impact factors.

Table 3 :
Standardized attribute values for geothermal potential regions.

Table 4 :
The coupling weights for all used impact factors.

Table 5 :
Comparison of calculated and actual results in Wudalianchi.