The quantitative classification of granites and their metallogenetic relations in West Qinling, Gansu Province, China

ABSTRACT The quantitative classification of granite and their metallogenetic relations have never been discussed. The Q-system clustering analysis and discriminant analysis methods were alternately used to quantitatively analyze the 11 oxide data in granite samples from the West Qinling area of Gansu Province, and then to construct the quantitative classification series models of granite (oxide). The granites samples are divided into three categories and eight sub-categories. The classification of granites is biased toward prospecting. According to the spatial correlation between eight types of granites and copper deposits, lead and zinc deposits, gold deposits, etc. (within 3 km of the intrusion) in the West Qinling area in Gansu Province, the “metallogenic related intrusions” are sought, and the prospecting target areas are defined. Furthermore, they provide reliable basis for regional geological prospecting.


Research background
There are more than 20 classification schemes for granite classification (Zhang, Wang, Xiong, & Li, 2008), such as type I, type S, type M, type A, collision type, volcanic arc type, mid-ocean ridge type, intraplate type, Haixi type, Caledon type, Andes type, modified type, remelting type, syntectic type and Adakite type, Himalayan type, Guangxi type, Zheming type, and Nanling type (Barbalin, 1999;Bonin, 1990;Chappell & White, 1974, 1992Collins, Beams, White, & Chappell, 1980;Eby, 1990;Ishihara, 1977;Maniar & Piccoli, 1989;Pitcher, 1983;Whalen, Currie, & Chappell, 1987;Zhang et al., 2008). These classification schemes mainly focus on the chemical components, some focus on the structural environments of granite formation, some focus on the classification of pressure relations, and others focus on the comprehensive factors (such as composition, source, and environment). However, there was little report on the quantitative classification of granite. Compared with the traditional qualitative geological research, quantitative classification is guided by the idea of big data (data-intensive), taking the existing rich geological information as the research object and using multivariate statistical analysis as the method (Liu, Li et al., 2018a;Liu, Li et al., 2018b;Zhang & Zhou, 2017;Zhou & Chen, 2018). Quantitatively identify the correlation between results and variables to achieve analysis and prediction (Luo, Zhang, & Peng, 2005;Luo et al., 2018;Wang et al., 2018). In addition, the quantitative classification of granite adopts the full data model, which extracts the valuable data from the whole data, pays close attention to the correlation between granite and mineralization, dilutes the causal relationship, and attaches great importance to the efficiency and result problems.
The West Qinling orogenic belt in Gansu Province is located between the Yangtze Craton and the Qaidam-Qianlian Orogenic Belt (Figure 1). The intrusive activity of granitic magma is very frequent and the intrusion of collisional granite magma is the dominant factor in the formation of gold deposits in this area (Han, Liu, Dong, & Ouyang, 2008;Jin, Zhang, He, & Jia, 2005;Liang, 2008;Qiu, Gong, Lu, & Liu, 2008;Yin & Zhao, 2006;Zhang et al., 2007). The geological data formed by the long-term accumulation of granite in the West Qinling orogenic belt are characterized by multidimension, temporality and spatiality, correlation and nonlinearity. In this paper, the quantitatively classification of granites in the West Qinling area are conducted to discuss their relationship with the spatial distribution of gold deposits, copper deposits, lead-zinc deposits, and other deposits, and further to determine the oreforming related intrusion, and mark the prospecting target areas, which has great significance for effective metallogenic prediction and optimization of the prospecting target area, based on ideas, methods of the geological big data and the data-driven model by big data applications.

Research methods
According to the data-driven model, using multivariate statistical analysis methods, to make classification and simulation (modeling) and excavating the implicit correlation between data, the metallogenic prediction and optimization of the prospecting target area can be made. In order to ensure the reliability of the statistical results, we need to filter the data first.
(1) Data source: The data of major elements collected from 538 magmatic rocks, 333 gold deposits, of which 260 are above the ore site, and 73 are mineralization sites, 129 lead-zinc deposits, including 90 points above the mine site and 39 mineralization points, and 140 copper deposits in the West Qinling area of Gansu Province.
(2) Data selection: In order to eliminate the effects of strong alteration and carbonation, samples with SiO 2 <66%, H 2 O >7%, and CO 2 >7% are removed. Meantime, the samples with extra high values which are 1-2 orders of magnitude higher than the average value are rejected, so as not to affect the calculation of the average value. After that, a total of 146 qualified granite bodies and 196 sets of main elements were selected. (3) According to the spatial distribution relationship between granite bodies and deposits (points) in West Qinling area, it is stipulated that the deposit within 3 km from the edge of the granite intrusion is the most favorable area for metallogenesis, those within 3-5 km is the general area for metallogenesis, and those are far away from more than 5 km has less correlation with the granitic intrusions. In this study, the most favorable area was selected as the preferred area for prospecting target area. (4) Selecting the major oxides of SiO 2 , TiO 2 , A1 2 O 3 , Fe 2 O 3 , FeO, MnO, MgO, CaO, Na 2 O, K 2 O and P 2 O 5 , etc. Take the large contribution for classification as the main categorical variables for quantitative classification and serial modeling, and to dig the implicit correlation among data.

Granite quantitative model
More than 20 kinds of granite classification methods have been determined under the traditional logical thinking method. There is no classification method based on big data ideas, methods, and techniques so far. The statistics of granite oxide parameters in the West Qinling area of Gansu Province (Table 1) show that the larger the standard deviation, the more intense the change in granites. The large standard deviation is an important classification index. According to the classification of two-dimensional projections such as the predecessor "TAS classification diagram", they are no longer suited to the needs of contemporary research due to the one-sidedness of the results. To more realistically reflect the spatial correlation between granites and deposits in the region, it is necessary to establish a true and reasonable classification method for granites in this area. According to the existing geological data, 11 major oxides of SiO 2 , TiO 2 , A1 2 O 3 , Fe 2 O 3 , FeO, MnO, MgO, CaO, Na 2 O, K 2 O, and P 2 O 5 are selected as the variables to form a granite classification database (full sample, full variable). First, the Q-type system clustering analysis is performed on the database (the intra-group connection is used as the clustering method, and the correlation is the calculation model) (Zhao, Wei, & Jin, 1990). In the Q-type clustering pedigree map, the samples with obvious classification are selected as Class I samples, and the other samples are Class II samples (Luo et al., 2005;Luo, Hou, Zhang, Peng, & Niu, 2006;Luo et al., 2018;Wang et al., 2018). Two types of discriminant analysis are made for the new database of I and II samples. The quantitative classification model of Class I granite in the West Qinling area of Gansu province is constructed. With R 0 = 2.48 as the critical value, the sample with R 0 ≥ 2.48 would be identified as Class I granite, and the Q-type system cluster analysis would be performed again for the samples with R 0 < 2.48. The quantitative classification model of Class II and III granites would be used to construct a series of quantitative classification models of granites (oxide) in the West Qinling area of Gansu Province and the relationship between each classification lineage (Figure 2).
At present, the granites in the West Qinling area are only divided into three categories and eight sub-categories. The quantitative classification model can further divide the classification model of the next level according to its own needs, until the samples are not clearly classified in the cluster analysis of the Q-type system. Because the quantitative classification series model is based on the content of various oxides in granite samples in the West Qinling area, it has certain limitations. The classification of granite in other areas may be inaccurate. The quantitative classification model of granites in different areas (orogenic belts) should be constructed according to the demand. In order to more accurately study the correlation between granite and other geological events and products (ore deposits), such quantitative classification will be more practical.

Application of granite quantitative model
The R 0 of the granite sample is calculated according to the formula in the quantitative model, and the classification can only start from the first model, and then enter the next judgment model according to the judgment condition of R 0 , and cannot directly judge the quantitative model. If R 0 = 2.48 is the critical value, the sample with R 0 ≥ 2.48 is determined as Class I granite, and the sample with R 0 < 2.48 enters the second model. If the sample of R 0 ≥ 0.23 enters the right model for classification, the sample with R 0 < 2.48 is classified into the left model, and so on. The classification of the final sample can be determined according to the R 0 threshold of different models. The calculation formula of the R 0 for granite sample in each model: R 0 = constant + ∑ content of granite oxide × coefficient corresponding to each oxide content in the model

Main characteristics of quantitative classification of granite
The difference between Class I granite, Class II granite, and Class III granite is positively correlated with the content of MgO and Fe 2 O 3 , and negatively correlated with the content of CaO and P 2 O 5 . The content of MgO and Fe 2 O 3 has a great contribution to the classification, accounting for 90.23%. Class I granite has higher iron-magnesium material than other Classes, and is relatively basic in granite (Table 2), which is classified into plagioclase granite and granodiorite in mineralogy classification.
Note: The oxide content of the above classification is the average of each category. The Class II and Class III granites are also classified mainly based on the oxide contents of MgO, TiO 2 , FeO, Fe 2 O 3 , Al 2 O 3 , MnO, Na 2 O, etc. Among them, the contents of MgO, TiO 2 , FeO, and Fe 2 O 3 have the largest contribution rate to the classification, accounting for 91.97%. The Class II granites are negatively correlated with MgO, TiO 2 , and FeO, and positively correlated with Fe 2 O 3 , whereas the Class III granites are positively correlated with MgO, TiO 2 , and FeO, and negatively correlated with Fe 2 O 3 , indicating that the iron-magnesium materials in the Class III granites are lower than those in Class I granite, but higher than those in Class II granite.
Class II granite, the most acidic granite, can be divided into three types of granite (II1, II2, II3), and Class II3 granite can be further divided into two types (II31, II32) granite. The classification of Class II1 granite is positively correlated with CaO and Fe 2 O 3 , and negatively correlated with FeO and MnO. The classification of Class II2 and Class II3 is reversed, and negatively correlated with CaO and Fe 2 O 3 , and positively correlated with FeO and MnO. These granitoids are classified into plagioclase granite, granite porphyry, (plaque-like) plaque monzonitic granite, and tonalite. The other classifications of Class II2 and Class II3 granites are negatively correlated with Al 2 O 3 , P 2 O 5 , Na 2 O, and SiO 2 , and positively correlated with TiO 2 . The Class II3 is negatively correlated with TiO 2 and positively correlated with Al 2 O 3 , P 2 O 5 , Na 2 O, and SiO 2 . Correlation, the Class II2 granite mineralogy, is classified into monzonitic granite, but it is acidic. The other classifications of Class II31 and Class II32 granites are negatively correlated with TiO 2 , P 2 O 5 , MgO, MnO, and SiO 2 , and positively correlated with Fe 2 O 3 . Class II32 granites are inversely related to TiO 2 , P 2 O 5 , MgO, MnO, and SiO 2 and negatively correlated with Fe 2 O 3 . The mineralogical classification of Class II31 granite is classified into granite, and the mineralogical classification of Class II32 granite is classified into granite porphyry, monzonitic granite, and syenogranite. For Class III granites, the SiO 2 content is between Class I and Class II granites, which can be divided into two types (III1, III2) granites, and Class III2 can be further divided into two types (III21, III22) granite. The classification difference between Class III1 and Class III2 granites is the positive correlation with P 2 O 5 , MgO, Fe 2 O 3 , and CaO in Class III1 granite and the negative correlation with P 2 O 5 , MgO, Fe 2 O 3 , and CaO in Class III2 granite. The Class III1 granite is classified into quartz diorite, a small number of granodiorite, dioritic porphyrite. The classification difference between Class III21 and Class III22 granites is positively correlated with FeO, P 2 O 5 , TiO 2 , and K 2 O, and negatively correlated with Fe 2 O 3 . The Class III22 granites are negatively correlated with FeO, P 2 O 5 , TiO 2 , and K 2 O, and positively correlated with Fe 2 O 3 . The mineralogical classification is classified into the monzonitic granite, but the Class III21 granite is relatively acidic compared to the Class III22 granite.

Statistical analysis of correlation between quantitative classification and main ore-forming spatial distribution
According to the spatial distribution relationship between granite bodies and various deposits in the West Qinling area, different types of granite bodies with gold deposits, copper deposits, lead-zinc deposits distributed and within or within 3 km of the rock mass are collectively referred to as "metallogenesis related rock mass". The statistical analysis of the correlation between the quantitative classification results of granites and the spatial distribution of deposits (Table 3) shows that different minerals and granite classifications have unique mineralization correlations (Figure 3).
(1) There are disagreements with different sub-categories in five intrusions (the samples are classified into different sub-categories) based on the classification methods. According to the number of samples, the intrusions with classification disagreements are classified to the type with more samples. The statistical results of quantitative classification of granites show that the West Qinling granites are mainly distributed in Table 3. Quantitative classification results of granites with ore deposits spatial correlation statistical analysis.
Quantitative classification Unit Ⅰ Ⅱ Ⅲ Total Ⅱ1 Ⅱ2 Ⅱ31 Ⅱ32 Ⅲ1 Ⅲ21 Ⅲ22 (2) Some of the same ore deposits are located within 3 km of the boundary of two or more rock masses, so the numbers of statistical deposits are greater than the actual number of deposits, but it does not affect the determination and statistics for ore-forming related rocks masses. The gold-bearing ore-related intrusions are mainly distributed in the Class II1, II2, and III1, and the ore deposits distribution are densely distributed. There is a great prospecting potential within 3 km of the granite body. There are little spatial correlations between the Class II3, III22 granites, and gold deposits, between the Class II32, II31, III21, III22 granites, and copper deposits, and between the Class III21, III22 granites, and lead-zinc deposits. The rock mass is mainly distributed in the Class III1 granites for copper deposits and the Class II1, III1 for lead-zinc deposits. In general, the gold, copper, and lead-zinc deposits have intimately spatial correlation with Class III1 granite, and it is looking for Class III1 granite and its surrounding area has a large prospecting space ( Figure 3). Therefore, quantitative classification can more accurately study the classification of granites and their spatial correlation with the deposit. (3) Although the above statistical results are not perfect due to the fact that some of the samples are not collected from granites in the West Qinling, it can generally provide a certain prospecting direction for further prospecting in the West Qinling area. (4) Taking 22 attributes such as rock mass tectonic belt, quantitative classification, qualitative classification, formation age, rock alteration and mineralization characteristics as variables, combined with geology, geophysical exploration, geochemical exploration, remote sensing and other information, according to Q-type hierarchical cluster analysis (using clustering within the group as a clustering method and correlation as a computational model), ideas and methods can be used to construct a quantitative prediction model for the Au, Cu, Pb, and Zn target areas in the West Qinling area of Gansu Province, which can be quickly and accurately delineated prospecting target areas.

Conclusions
(1) The modeling ideas and methods of the "granular (oxide) quantitative classification series model in the West Qinling area of "Gansu Province" are not perfect and need further improvement. The sub-classes of granite quantitative classification can be further modeled to divide the smaller classification, and the model itself is established only based on the granite data in the West Qinling, which has certain limitations. Different regions should build different quantitative classification models of granites according to the practical demand, and strive to be more accurate study of the correlation between granite and other geological events with their products (ore deposits), which will make such quantitative classification more practical and meaningful.
(2) Quantitative classification series models are controlled by multiple variables of different combinations. The quantitative classification results are preferably controlled by the test data, so it is crucial to ensure the authenticity of the original data and improve the representativeness of the samples. (3) Quantitative classification is more accurate in studying the classification of granite and their spatial correlation with the deposit. The classification results of granite in West Qinling are biased toward the classification of ore deposits, mainly concentrated in the three types of Class II1, II2, and III1 granites. In general, gold deposits with Class II1, II2, and III1 granites, copper deposits with Class III1 granites, and lead-zinc deposits with Class II1 and III1 granites have close spatial correlation.