Identifying Geochemical Anomalies Associated with Gold Mineralization Using Factor Analysis and Spectrum–Area Multifractal Model in Laowan District, Qinling-Dabie Metallogenic Belt, Central China

: The Laowan deposit is a typical gold deposit in the Qinling-Dabie metallogenic belt, which produces the most gold resources in Central China. After being explored for decades, follow-up exploration requires additional theoretical support. In this study, the factor analysis (FA) and spectrum–area (S–A) multifractal model were used to process multi-element geochemical data from 369 samples collected in the study area for identifying the geochemical anomalies associated with gold mineralization. The results showed that: (1) the mean Au content in this region is up to 1000 times higher than the Au background values of the upper crust of the South Qinling unit; (2) the factor analysis revealed that Au, Ag, Cu, As, Sb, and S can be used as direct ore-prospecting criteria; (3) the observed elemental zonation is consistent with the zonation of metallic elements in the magmatic–hydrothermal system. This supports the magmatic–hydrothermal origin of the Laowan deposit; (4) the spectrum–area fractal model can help to decompose the geochemical patterns in a complex geological setting. The decomposed geochemical anomaly map obtained by the S–A multifractal model indicated that highly anomalous areas have a great relationship with the Au occurrence and can be a guidance for further exploration in the study area.


Introduction
Geochemical exploration is a significant method for finding various types of ore deposits in mineral exploration [1][2][3][4][5][6][7][8]. Many deposits have been discovered by geochemical data processing and anomaly identification [9]. Geochemical exploration data includes numerous elements. The complexity and diversity of geological processes lead to different element associations [10]. Some element associations are closely related to specific types of mineralization, which can be used to guide ore-prospecting [11].
Factor analysis (FA) is one of the most popular multivariate analysis methods for reducing the dimensionality of datasets. It is a useful tool for combining several correlated variables into a single factor that represents a certain inherent feature in order to summarize the multivariate information [12]. The Laowan gold deposit is located in the East of the Qinling-Dabie orogenic belt and the Laowan ductile shear zone. The Northern and Southern sides are bounded by the Laowan and the Songpa ductile shear belt. The main host rock of the Laowan deposit is a metamorphic rock series from the Guishan Formation of the middle Mesoproterozoic period, which has undergone many ductile and brittle shear activities ( Figure 2). The lithology is mainly mica-quartz schist and plagioclase-hornblende schist, and mylonitic silicified rocks [53,54,58,67]. The original rocks are basic volcano-clastic sedimentary [53]. The Laowan ductile shear zone is the ore-controlling structure of the deposit, which has obvious deformation marks and various deformation fabrics. In addition, the fault is well-developed in the region, including the Laowan fault on the Northern side and the Songpa fault on the Southern sides, and several groups of brittle fault structures in the N-NW (North-northwest) and N-NE (Northhnortheast) directions [58]. The N-NW trending fault is the main ore-bearing structure in the area. The Laowan fault is a post-mineralization structure [67]. The Songpa fault was dominated by a ductile strike-slip in the early stages and superimposed brittle fractures in the later stages. The fault was formed early and revived after emplacement by granite porphyry dykes in the later stages. This caused the granite porphyry dykes to be cataclastic, with the local mineralization of gold and silver. The Laowan gold deposit is located in the East of the Qinling-Dabie orogenic belt and the Laowan ductile shear zone. The Northern and Southern sides are bounded by the Laowan and the Songpa ductile shear belt. The main host rock of the Laowan deposit is a metamorphic rock series from the Guishan Formation of the middle Mesoproterozoic period, which has undergone many ductile and brittle shear activities ( Figure 2). The lithology is mainly mica-quartz schist and plagioclase-hornblende schist, and mylonitic silicified rocks [53,54,58,67]. The original rocks are basic volcano-clastic sedimentary [53]. The Laowan gold deposit is located in the East of the Qinling-Dabie orogenic belt and the Laowan ductile shear zone. The Northern and Southern sides are bounded by the Laowan and the Songpa ductile shear belt. The main host rock of the Laowan deposit is a metamorphic rock series from the Guishan Formation of the middle Mesoproterozoic period, which has undergone many ductile and brittle shear activities ( Figure 2). The lithology is mainly mica-quartz schist and plagioclase-hornblende schist, and mylonitic silicified rocks [53,54,58,67]. The original rocks are basic volcano-clastic sedimentary [53]. The Laowan ductile shear zone is the ore-controlling structure of the deposit, which has obvious deformation marks and various deformation fabrics. In addition, the fault is well-developed in the region, including the Laowan fault on the Northern side and the Songpa fault on the Southern sides, and several groups of brittle fault structures in the N-NW (North-northwest) and N-NE (Northhnortheast) directions [58]. The N-NW trending fault is the main ore-bearing structure in the area. The Laowan fault is a post-mineralization structure [67]. The Songpa fault was dominated by a ductile strike-slip in the early stages and superimposed brittle fractures in the later stages. The fault was formed early and revived after emplacement by granite porphyry dykes in the later stages. This caused the granite porphyry dykes to be cataclastic, with the local mineralization of gold and silver. The Laowan ductile shear zone is the ore-controlling structure of the deposit, which has obvious deformation marks and various deformation fabrics. In addition, the fault is well-developed in the region, including the Laowan fault on the Northern side and the Songpa fault on the Southern sides, and several groups of brittle fault structures in the N-NW (North-northwest) and N-NE (Northh-northeast) directions [58]. The N-NW trending fault is the main ore-bearing structure in the area. The Laowan fault is a post-mineralization structure [67]. The Songpa fault was dominated by a ductile strike-slip in the early stages and superimposed brittle fractures in the later stages. The fault was Minerals 2020, 10, 229 4 of 18 formed early and revived after emplacement by granite porphyry dykes in the later stages. This caused the granite porphyry dykes to be cataclastic, with the local mineralization of gold and silver.
Magmatic rocks are relatively developed in the mining area, including numerous granite porphyry and variable gabbro. The Laowan granitoid intrusion with an age of 133 Ma [68], 127 Ma (test by author, unpublished) is distributed across the Laowan fault zone and was formed in the middle-late Yanshanian period with the gold deposit [53]. The lithology is feldspar biotite granite. The granite porphyry dykes, with an age of 139 Ma [58], are mainly distributed along the Songpa fault. The edge of the dykes is highly fractured and formed gold, silver, and other polymetallic mineralization. The lithology is feldspar biotite granite.
The mineralization of the Laowan deposit can be divided into three stages ( Figure 3): (1) the quartz oxide stage corresponds to the early stage of mineralization, and the mineral assemblage is quartz-magnetite-pyrite-gold; (2) the quartz polymetallic sulfide stage and the mineral assemblage is quartz-pyrite-sphalerite-gold, which is the main mineralization stage; (3) the quartz carbonate stage represents the late stage of mineralization, and the mineral assemblage is quartz-calcite-pyrite. The main ore types are altered amphibolic rocks, quartz-sericite-pyrite type, and quartz-feldspathic-pyrite type [56]. The time of the Laowan gold deposit is 138 ± 2 Ma obtained by 39Ar/40Ar dating of hydrothermal muscovite [55]. Figure 4 shows the mineral composition and symbiosis of ore as well as the occurrence of gold. Magmatic rocks are relatively developed in the mining area, including numerous granite porphyry and variable gabbro. The Laowan granitoid intrusion with an age of 133 Ma [68], 127 Ma (test by author, unpublished) is distributed across the Laowan fault zone and was formed in the middle-late Yanshanian period with the gold deposit [53]. The lithology is feldspar biotite granite. The granite porphyry dykes, with an age of 139 Ma [58], are mainly distributed along the Songpa fault. The edge of the dykes is highly fractured and formed gold, silver, and other polymetallic mineralization. The lithology is feldspar biotite granite.
The mineralization of the Laowan deposit can be divided into three stages ( Figure 3): (1) the quartz oxide stage corresponds to the early stage of mineralization, and the mineral assemblage is quartz-magnetite-pyrite-gold; (2) the quartz polymetallic sulfide stage and the mineral assemblage is quartz-pyrite-sphalerite-gold, which is the main mineralization stage; (3) the quartz carbonate stage represents the late stage of mineralization, and the mineral assemblage is quartz-calcite-pyrite. The main ore types are altered amphibolic rocks, quartz-sericite-pyrite type, and quartz-feldspathicpyrite type [56]. The time of the Laowan gold deposit is 138 ± 2 Ma obtained by 39Ar/40Ar dating of hydrothermal muscovite [55]. Figure 4 shows the mineral composition and symbiosis of ore as well as the occurrence of gold.

Sampling
In this study, the geochemical dataset was based on the analyses of 369 hand samples collected within an area of approximately 10 km 2 . The samples covered four ore sections: Pengjialaozhuang, Shangshanghe, Laowan, and Beiyangzhuang ( Figure 5). The geochemical profiles were roughly perpendicular to the main ore-controlling structures in the study area. A global positioning system (GPS) was used to locate the samples, and the error on location was less than 5 m. The lithologies contained cataclasite, schist, granitoid, quartz-feldspathic veins, and quartz veins. The continuous sample picking method was used to uniformly collect rocks of the same lithology to form one sample. The sampling density was increased where the lithologies were complex. The weight of each sample was more than 500 g. Prior to elemental analysis, the rocks were dried, broken, washed, mixed, and divided (the error in division was less than 1%) in the detection institute, according to the standard

Sampling
In this study, the geochemical dataset was based on the analyses of 369 hand samples collected within an area of approximately 10 km 2 . The samples covered four ore sections: Pengjialaozhuang, Shangshanghe, Laowan, and Beiyangzhuang ( Figure 5). The geochemical profiles were roughly perpendicular to the main ore-controlling structures in the study area. A global positioning system (GPS) was used to locate the samples, and the error on location was less than 5 m. The lithologies contained cataclasite, schist, granitoid, quartz-feldspathic veins, and quartz veins. The continuous sample picking method was used to uniformly collect rocks of the same lithology to form one sample. The sampling density was increased where the lithologies were complex. The weight of each sample was more than 500 g. Prior to elemental analysis, the rocks were dried, broken, washed, mixed, and divided (the error in division was less than 1%) in the detection institute, according to the standard (No.DZ/T0258-2014) issued by the Ministry of Land and Resources of the China, Beijing, China. For each sample, 15 elements, including Au, Ag, Cu, Pb, Zn, As, Sb, Bi, Hg, W, Mo, Sn, F, Co, and Ni, were determined at the Mineral Resources Supervision and Testing Center of Henan province, Zhengzhou, China. The analytical schemes for each element are shown in Table 1.

Factor Analysis
FA, which was outlined by Reimann et al. [13], is a powerful tool for processing geochemical data [69]. This method simplifies the interpretation of a complex dataset by reducing highdimensional data to a few representative factors that cannot be observed directly, but represent certain features inherent in the original dataset: e.g., [12,70,71]. FA concentrates on the potential processes controlling the data distribution rather than the responses [12].
In usual model assumptions [69], the FA model is defined as:

Factor Analysis
FA, which was outlined by Reimann et al. [13], is a powerful tool for processing geochemical data [69]. This method simplifies the interpretation of a complex dataset by reducing high-dimensional data to a few representative factors that cannot be observed directly, but represent certain features inherent in the original dataset: e.g., [12,70,71]. FA concentrates on the potential processes controlling the data distribution rather than the responses [12].
In usual model assumptions [69], the FA model is defined as: where y denotes the random vector, with the factor f of the dimension k < D. e denotes the error term, and Λ denotes the loading matrix. Ψ = Cov(e) has a diagonal form, and the diagonal elements represent uniqueness or unique variance [12]. This can also be used to solve complex questions of geological genesis, as well as the mineralization superposition problem [72], and determine the vectors towards concealed ore deposits [13]. By grouping related variables to form factors, the reduction of variables and simplification of data interpretation can be realized [68]. The results include scree plots, factor score maps, and loading plots that have been widely used to interpret geochemical data: e.g., [72][73][74][75]. Generally, the data were log-transformed and tested to determine whether all variables were normally distributed before FA [13]. The Kaiser-Meyer-Olkin index [76] and the test of sphericity [77] are often used to measure if FA is useful and appropriate for the data. The relevant calculations of FA in this study were calculated by SPSS software (25.0, IBM, Chicago, IL, USA).

Spectrum-Area Multifractal Model
Based on the fractal theory, Cheng et al. [23] proposed the C-A fractal model, concerning the relationship between the concentration value (ρ) and the area (A) enclosed by the contour value : where A (ρ) represents the area enclosed by the contour value ρ; ν is the threshold defining contour ρ. a 1 and a 2 (greater than 0) are fractal dimensions that can be estimated from slopes of straight lines in the log-log plot of A (ρ) versus ρ. ∝ denotes proportionality [78]. Based on the argument that geochemical patterns in the spatial domain could be considered as superimposed signals of different frequencies, Cheng et al. [50] proposed the S-A model based on the C-A model to separate geochemical anomalies from the background levels and characterize the spectral energy density-area relationship in the frequency domain: where S represents the spectral energy density, and A (≥S) symbolizes the area of the wavenumber with a threshold above S. β denotes the anisotropic scaling exponent, d represents the degree of overall concentration, while ∝ denotes proportionality [78]. According to the anisotropic scale characteristics in the frequency domain, the S-A multifractal model decomposes the spatial model into several components. After lots of practice, it has proved to be useful for breaking down complex geochemical patterns: e.g., [7,40,[79][80][81].

Statistical Analysis
In this study, a descriptive statistical analysis was performed on the concentrations of 15 significant elements from 369 samples. The statistical parameters (mean and standard deviation) and mean background values of these elements in the upper continental crust of the South Qinling unit [82] are given in Table 2. The coefficient of variation (CV) is used to represent the dispersion of data. The mean of Au content is 1232.14 ng/g, which is much higher than the background values in the upper continental crust of the South Qingling unit (1.10 ng/g) [82]. The CV values for Au, Ag, Sb, Bi, Mo, and Hg are over five, indicating the inhomogeneous distribution of these elements, meaning that they were more likely to be concentrated proximally to ore.
The Q-Q plot is a powerful tool for visualizing element distributions [83]. The Q-Q plots ( Figure 6) indicated that the geochemical data processed by log-transformation did not obey normal distribution. The data parallel to the line represents the dominant background value in the geochemical data. In contrast, the data at both ends represent a relatively small number of anomalies, which are of great significance.
This study also took the contents of selected elements in different lithologies into account to consider the effect of hosting rocks. Due to the limited number of samples of metagabbro and granite porphyry dykes, and the complexity and diversity of the genesis of quartz veins, their statistical significances were limited, so are not discussed in this study.
As shown in Table 2, Au content was highest in the fracture zone, followed by mica-quartz schist in the Guishan Formation. Ag and Sb had the highest contents in marble, followed by the fracture zone. As content in mica-quartz schist was the highest, followed by the fracture zone. The contents of siderophile elements such as Cu, Co, and Ni were highest in the fracture zone, followed by the Guishan Formation, which reflects that the basement rock of the strata was basic rock. In summary, the fracture zone had good potential for gold mineralization, and the Guishan Formation was mean ore-bearing stratum. This study also took the contents of selected elements in different lithologies into account to consider the effect of hosting rocks. Due to the limited number of samples of metagabbro and granite porphyry dykes, and the complexity and diversity of the genesis of quartz veins, their statistical significances were limited, so are not discussed in this study.
As shown in Table 2, Au content was highest in the fracture zone, followed by mica-quartz schist in the Guishan Formation. Ag and Sb had the highest contents in marble, followed by the fracture zone. As content in mica-quartz schist was the highest, followed by the fracture zone. The contents of siderophile elements such as Cu, Co, and Ni were highest in the fracture zone, followed by the Guishan Formation, which reflects that the basement rock of the strata was basic rock. In summary, the fracture zone had good potential for gold mineralization, and the Guishan Formation was mean ore-bearing stratum.

Elemental Association
Some elements can be a guide to discover mineral deposits because they are spatially associated with specific types of mineralization. In contrast, the remaining elements might reflect other geological processes or events unrelated to mineralization [9]. An R-type factor analysis was used to determine the ore-forming elemental association. The principal component extraction [84] was used to extract uncorrelated components, which represented most of the variability carried on the data variables. According to the varimax rotation technique [85], the factor loadings matrix was rotated orthogonally to minimize the number of high-load variables on each factor. All variables in the factor analysis were used to construct a scatter plot (Figure 7). Four variables with eigenvalues greater than 1.0 were selected, providing a cumulative explained variance of 75.91%.

Elemental Association
Some elements can be a guide to discover mineral deposits because they are spatially associated with specific types of mineralization. In contrast, the remaining elements might reflect other geological processes or events unrelated to mineralization [9]. An R-type factor analysis was used to determine the ore-forming elemental association. The principal component extraction [84] was used to extract uncorrelated components, which represented most of the variability carried on the data variables. According to the varimax rotation technique [85], the factor loadings matrix was rotated orthogonally to minimize the number of high-load variables on each factor. All variables in the factor analysis were used to construct a scatter plot (Figure 7). Four variables with eigenvalues greater than 1.0 were selected, providing a cumulative explained variance of 75.91%.
The results of the FA are shown in Table 3. Communalities represent the common factor variance extracted from each element, with a higher value indicating a better explanation of variables. According to the results, Factor 1 (F1) was dominated by Cu, Co, Ni, V, and Fe; Factor 2 (F2) was dominated by Au, Ag, Cu, As, Sb, and S; Factor 3 (F3) was dominated by Bi, F, Mo, and Sn; and Factor 4 was dominated by Pb and Zn.   The results of the FA are shown in Table 3. Communalities represent the common factor variance extracted from each element, with a higher value indicating a better explanation of variables. According to the results, Factor 1 (F1) was dominated by Cu, Co, Ni, V, and Fe; Factor 2 (F2) was dominated by Au, Ag, Cu, As, Sb, and S; Factor 3 (F3) was dominated by Bi, F, Mo, and Sn; and Factor 4 was dominated by Pb and Zn. Group F1 consisted of siderophile elements that occur mainly in minerals, such as olivine and pyroxene, and are related to basic and ultrabasic rocks representing the Guishan Formation. Group F2 represented the main metallogenic element association. These elements mainly precipitate in mid-low temperature environments and are a sulphophile element, which corresponds to the mineral assemblage in the quartz-polymetallic sulfide stage. Therefore, this factor represents the main metallogenic stage. Cu was present both in group F1 and group F2 because it has both siderophile and sulphophile properties. Group F3 indicated high-temperature hydrothermal metasomatism, which is closely related to crust source granite. Group F4 may represent Pb-Zn mineralization or the primary halo near the gold ore bodies.
As illustrated in the factor score maps (Figure 8), the areas with the highest scores in group F1 were widely distributed between the Songpa and Laowan faults, and highly consistent with the distribution areas of the Guishan Formation. The factor scores of group F1 were also high in the local area, where the Guishan Formation was in contact with the Laowan granite. Areas with high factor scores in group F2, which was the main ore-forming element association, were mainly distributed in the middle of the Shangshanghe and Laowan ore sections, with other high anomalies distributed around granite porphyry dykes. On the other hand, the factor scores in group F2 in the area near the Laowan granite and the Southwest and Northeast of the study area were lower. The highest scores from group F3 were mainly around granite porphyry dykes. There were also low anomalies around the Laowan granite. The covered area of the Guishan Formation in the central part of the study area had low anomalies with small local areas, which may represent undiscovered concealed intrusions. Areas with high factor scores in group F4 were mainly distributed in areas around granite porphyry dykes and the Central area of the junction between the Shangshanghe and Laowan ore sections, coinciding with the anomalous areas in group F2 and F3. Group F1 consisted of siderophile elements that occur mainly in minerals, such as olivine and pyroxene, and are related to basic and ultrabasic rocks representing the Guishan Formation. Group F2 represented the main metallogenic element association. These elements mainly precipitate in midlow temperature environments and are a sulphophile element, which corresponds to the mineral assemblage in the quartz-polymetallic sulfide stage. Therefore, this factor represents the main metallogenic stage. Cu was present both in group F1 and group F2 because it has both siderophile and sulphophile properties. Group F3 indicated high-temperature hydrothermal metasomatism, which is closely related to crust source granite. Group F4 may represent Pb-Zn mineralization or the primary halo near the gold ore bodies.
As illustrated in the factor score maps (Figure 8), the areas with the highest scores in group F1 were widely distributed between the Songpa and Laowan faults, and highly consistent with the distribution areas of the Guishan Formation. The factor scores of group F1 were also high in the local area, where the Guishan Formation was in contact with the Laowan granite. Areas with high factor scores in group F2, which was the main ore-forming element association, were mainly distributed in the middle of the Shangshanghe and Laowan ore sections, with other high anomalies distributed around granite porphyry dykes. On the other hand, the factor scores in group F2 in the area near the Laowan granite and the Southwest and Northeast of the study area were lower. The highest scores from group F3 were mainly around granite porphyry dykes. There were also low anomalies around the Laowan granite. The covered area of the Guishan Formation in the central part of the study area had low anomalies with small local areas, which may represent undiscovered concealed intrusions. Areas with high factor scores in group F4 were mainly distributed in areas around granite porphyry dykes and the Central area of the junction between the Shangshanghe and Laowan ore sections, coinciding with the anomalous areas in group F2 and F3. The distribution correlations of the factor score anomalies were overlapping or neighboring. The zonation of ore elements was around granite porphyry dykes, suggesting that ore-forming fluids may have emanated from these intrusions, which also provided thermodynamic. The distribution of anomalies in group F2 and the main metallogenic element showed an obvious overlap, indicating that the Guishan Formation was the main ore-bearing formation. The areas where anomalies from group F4 overlapped with those of the main ore-forming elements may have been the primary halo around the gold ore bodies, while the areas where the anomalies of group F4 overlapped with granite porphyry dykes may represent another polymetallic mineralization.
Overall, the distribution of anomalies across the four factors was generally controlled by regional ductile shear structures. The neighboring characteristics of the distributions displayed a zonation pattern that approximated Bi-Mo and Pb-Zn in the middle and distal Au-Ag. The gold The distribution correlations of the factor score anomalies were overlapping or neighboring. The zonation of ore elements was around granite porphyry dykes, suggesting that ore-forming fluids may have emanated from these intrusions, which also provided thermodynamic. The distribution of anomalies in group F2 and the main metallogenic element showed an obvious overlap, indicating that the Guishan Formation was the main ore-bearing formation. The areas where anomalies from group F4 overlapped with those of the main ore-forming elements may have been the primary halo around the gold ore bodies, while the areas where the anomalies of group F4 overlapped with granite porphyry dykes may represent another polymetallic mineralization.
Overall, the distribution of anomalies across the four factors was generally controlled by regional ductile shear structures. The neighboring characteristics of the distributions displayed a zonation pattern that approximated Bi-Mo and Pb-Zn in the middle and distal Au-Ag. The gold mineralization system in a low-temperature environment was furthest from the magmatic activity center, while the Bi-Mo association in a high-temperature environment was close to the magmatic activity center, and the Pb-Zn mineralization system was in the middle of the two, showing superposition characteristics. This zonation feature was consistent with the zonation rule of metallic elements in a magmatic-hydrothermal liquid system related to intrusion [86]-that is, in a high-temperature system, Bi-Mo was concentrated closest to the center of the magmatic-hydrothermal metallogenic system, while gold, silver, and other elements were relatively far away from the magmatic-hydrothermal metallogenic system.
Based on the distribution of each factor and the metallogenic model, we inferred that ore-forming fluid might emanate from granite porphyry dykes, migrate along with the ductile shear structure in a NW-SE direction under the activation of Yanshanian magmatic activity, and, finally, precipitate in the secondary fault. From the center of the magmatic activity to the periphery, the element associations in the study area showed zonation patterns of Mo → Pb-Zn → Au-Ag. Combined with the age of the deposit (138 Ma) [55] versus the intrusions (127-139 Ma) [58,67] with the spatial distribution of element associations, the gold and polymetallic metallogenesis of the Laowan deposit had close spatial and temporal relationships with Yanshanian magmatism (especially the granite porphyry), which supported the magmatic-hydrothermal origin of the Laowan deposit. This study not only verified the previous views on the source of ore-formation and the genetic mechanism of the Laowan deposit, but also could be used as an independent statistical validation of the element zonation model proposed by the predecessors [58].

S-A Multifractal Analysis
The S-A fractal method was applied to decompose the mixed geochemical populations and identify geochemical anomalies in the first factor scores. Before multifractal analysis, the interpolation of point data to raster maps was used. ArcGIS software developed by Esri was used to set the search radius to 20 points, then, the inverse distance weighted (IDW) method was used for interpolation, and the raster data after interpolation was obtained. GeoFractal software (2.0, State Key Laboratory of Geological Processes and Mineral Resources, Wuhan, China) [87] was used for S-A fractal analysis.
Fourier transform was used to transform the spatial distribution pattern of the F1 factor scores into the power spectrum density and area, and the least-squares method was used to fit three lines ( Figure 9). Two thresholds (0.365 and 6.897) were selected, and the three lines produced represented noise, giving values of y = −0.09x + 5.22 with R2 = 0.81, an anomaly of y = −0.44x + 5.36 with R2 = 0.99, and a background of y = −0.63x + 6.04 with R2 = 1.00, from left to right (Figure 9). Minerals 2020, 10, 229 12 of 18 mineralization system in a low-temperature environment was furthest from the magmatic activity center, while the Bi-Mo association in a high-temperature environment was close to the magmatic activity center, and the Pb-Zn mineralization system was in the middle of the two, showing superposition characteristics. This zonation feature was consistent with the zonation rule of metallic elements in a magmatic-hydrothermal liquid system related to intrusion [86]-that is, in a hightemperature system, Bi-Mo was concentrated closest to the center of the magmatic-hydrothermal metallogenic system, while gold, silver, and other elements were relatively far away from the magmatic-hydrothermal metallogenic system. Based on the distribution of each factor and the metallogenic model, we inferred that oreforming fluid might emanate from granite porphyry dykes, migrate along with the ductile shear structure in a NW-SE direction under the activation of Yanshanian magmatic activity, and, finally, precipitate in the secondary fault. From the center of the magmatic activity to the periphery, the element associations in the study area showed zonation patterns of Mo → Pb-Zn → Au-Ag. Combined with the age of the deposit (138 Ma) [55] versus the intrusions (127-139 Ma) [58,67] with the spatial distribution of element associations, the gold and polymetallic metallogenesis of the Laowan deposit had close spatial and temporal relationships with Yanshanian magmatism (especially the granite porphyry), which supported the magmatic-hydrothermal origin of the Laowan deposit. This study not only verified the previous views on the source of ore-formation and the genetic mechanism of the Laowan deposit, but also could be used as an independent statistical validation of the element zonation model proposed by the predecessors [58].

S-A Multifractal Analysis
The S-A fractal method was applied to decompose the mixed geochemical populations and identify geochemical anomalies in the first factor scores. Before multifractal analysis, the interpolation of point data to raster maps was used. ArcGIS software developed by Esri was used to set the search radius to 20 points, then, the inverse distance weighted (IDW) method was used for interpolation, and the raster data after interpolation was obtained. GeoFractal software (2.0, State Key Laboratory of Geological Processes and Mineral Resources, Wuhan, China) [87] was used for S-A fractal analysis.
Fourier transform was used to transform the spatial distribution pattern of the F1 factor scores into the power spectrum density and area, and the least-squares method was used to fit three lines ( Figure 9). Two thresholds (0.365 and 6.897) were selected, and the three lines produced represented noise, giving values of y = −0.09x + 5.22 with R2 = 0.81, an anomaly of y = −0.44x + 5.36 with R2 = 0.99, and a background of y = −0.63x + 6.04 with R2 = 1.00, from left to right (Figure 9).  The anomaly and background separated by the above method are shown in Figure 10, and the drill holes that reached the ore bodies are projected on the figure (shown as black triangles). As seen in Figure 10, the background value was not evenly distributed in the study area. These drill holes were not only distributed in the areas with high background values, but also in areas with low background values, indicating that the traditional methods of identifying geochemical anomalies, such as calculating the mean and standard deviation, are not effective when the background value distribution is not uniform. The areas with high background values were mostly related to the Guishan Formation and granite porphyry dykes, suggesting that they may provide the metal source for the formation of the ore deposits. It can also be seen from the anomaly map in Figure 10 that the morphology of the anomalies was controlled by the ductile shear zone. The anomaly map of the S-A analysis provided a tool for heightening weak geochemical anomalies. Most drill holes with encountered ore were distributed in or near the areas with high anomaly values, and a little deviation may have been caused by the angles of the drill holes, indicating that this S-A fractal model had an excellent predictive effect in this study area and presented valid targets for further exploration. Combined with the viewpoints above, few prospecting targets were divided into the periphery of the Shangshanghe and Laowan ore sections, which are the main mining areas ( Figure 11). The anomaly and background separated by the above method are shown in Figure 10, and the drill holes that reached the ore bodies are projected on the figure (shown as black triangles). As seen in Figure 10, the background value was not evenly distributed in the study area. These drill holes were not only distributed in the areas with high background values, but also in areas with low background values, indicating that the traditional methods of identifying geochemical anomalies, such as calculating the mean and standard deviation, are not effective when the background value distribution is not uniform. The areas with high background values were mostly related to the Guishan Formation and granite porphyry dykes, suggesting that they may provide the metal source for the formation of the ore deposits. It can also be seen from the anomaly map in Figure 10 that the morphology of the anomalies was controlled by the ductile shear zone. The anomaly map of the S-A analysis provided a tool for heightening weak geochemical anomalies. Most drill holes with encountered ore were distributed in or near the areas with high anomaly values, and a little deviation may have been caused by the angles of the drill holes, indicating that this S-A fractal model had an excellent predictive effect in this study area and presented valid targets for further exploration. Combined with the viewpoints above, few prospecting targets were divided into the periphery of the Shangshanghe and Laowan ore sections, which are the main mining areas ( Figure 11).

Conclusions
In this study, the Laowan gold deposit in the Qinling-Dabie metallogenic belt illustrated how the geochemical exploration method could be applied to identify the geochemical anomalies associated with gold mineralization. Factor analysis was used to reveal the Au polymetallic association, which can be used as direct ore-prospecting criteria. We established a resource to potentially evaluate the model at the periphery of the mining area using S-A multifractal analysis. The following conclusions were obtained: (1) The factor analysis model in this study revealed that the gold polymetallic elements of the Laowan gold deposit was Au, Ag, Cu, As, Sb, and S. The distribution of factors also portrayed the metallogenic patterns that the ore-forming fluids were mainly from granite porphyry dykes and the Guishan Formation, and migrated along the NW-SE striking faults under the activation and remelting of Yanshanian magma; (2) From the center of magmatic activity, which is granite porphyry dykes, to the periphery, the element associations in the study area displayed Bi-Mo → Pb-Zn → Au-Ag. This zonation characteristic suggests that the Laowan gold deposit is magmatic-hydrothermal deposit and more likely related to Yanshanian magmatism; (3) The hybrid method combining FA and S-A fractal modelling is a useful tool for identifying geochemical anomalies. The former is used to integrate multi-element concentration values, and the latter is applied to decompose mixed geochemical patterns in a complex geological setting.

Conclusions
In this study, the Laowan gold deposit in the Qinling-Dabie metallogenic belt illustrated how the geochemical exploration method could be applied to identify the geochemical anomalies associated with gold mineralization. Factor analysis was used to reveal the Au polymetallic association, which can be used as direct ore-prospecting criteria. We established a resource to potentially evaluate the model at the periphery of the mining area using S-A multifractal analysis. The following conclusions were obtained: (1) The factor analysis model in this study revealed that the gold polymetallic elements of the Laowan gold deposit was Au, Ag, Cu, As, Sb, and S. The distribution of factors also portrayed the metallogenic patterns that the ore-forming fluids were mainly from granite porphyry dykes and the Guishan Formation, and migrated along the NW-SE striking faults under the activation and re-melting of Yanshanian magma; (2) From the center of magmatic activity, which is granite porphyry dykes, to the periphery, the element associations in the study area displayed Bi-Mo → Pb-Zn → Au-Ag. This zonation characteristic suggests that the Laowan gold deposit is magmatic-hydrothermal deposit and more likely related to Yanshanian magmatism; (3) The hybrid method combining FA and S-A fractal modelling is a useful tool for identifying geochemical anomalies. The former is used to integrate multi-element concentration values, and the latter is applied to decompose mixed geochemical patterns in a complex geological setting.