CoDA-Based Geo-Electrochemical Prospecting Prediction of Uranium Orebodies in Changjiang Area, Guangdong Province, China

: In order to achieve a breakthrough in the exploration of uranium orebodies in the


Introduction
Uranium mineral resources, which are important strategic resources and energy minerals, have always been crucial in ensuring the strategic development of various countries.China's breakthroughs in uranium exploration are currently under close scrutiny.The confirmed mineral resources are facing depletion, necessitating urgent exploration in deeper or peripheral areas for mineral resources.Geo-electrochemical methods are forms of exploration techniques that use deep penetration, and they exhibit characteristics such as greater detection depth and the ability to identify weak anomalies compared to traditional soil geochemical methods [1][2][3].Using this method allows for a more significant emphasis on potential anomalies present in deeper regions.
The geo-electrochemical technique was introduced to China from the former Soviet Union by Chinese scientists Fei Xiquan, Xu Bangliang, and Deng Kangle in the 1980s.Since then, it has stimulated significant research work in China [1,2].Luo [3][4][5][6][7][8] used the geo-electrochemical method to analyze over 30 different types of mineral deposits with overburden layers of varying thicknesses.The team repeatedly adjusted the technical parameters of the method and conducted comparative studies to summarize technical parameters suitable for different buried mineral resources, including independent power supply and fast positioning.This research has achieved good results in the exploration of precious metals, non-ferrous metals, and radioactive metal deposits in China, involving more than ten types of ores, such as gold, lead, zinc, tungsten, arsenic, silver, copper, antimony, and uranium.This method has also been widely used in the United States, Canada, the United Kingdom, Australia, India, and Southeast Asia, and is one of the most important methods for searching and exploring buried mineral resources in the world [9][10][11][12][13][14][15].
However, despite the numerous successful applications of the geo-electrochemical extraction method for prospecting hidden metallic ores, there remains significant room for further development.For instance, Liu [7] highlighted that, although there have been nearly 200 articles on the geo-electrochemical extraction method indexed in CNKI, only a mere 3% pertain to articles focusing on data processing.With the advancement of geochemical exploration, data processing methodologies have rapidly evolved from univariate to multivariate analyses and from linear to nonlinear analyses.The traditional approach of averaging values with standard deviations no longer holds a competitive edge and fails to capture the essence of the underlying data.As a deep penetration method, geo-electrochemical extraction warrants a comprehensive data analysis to unearth more insightful information, thus serving the prospecting of mineral resources more effectively.
As a type of compositional data, geo-electrochemical data is subject to the "closure effect", where the increase or decrease in the content of one element in a sample of the results changes the content of other elements, which may lead to a reverse change in the correlation between elements (positive correlation to negative correlation) when processing geo-electrochemical data.To address this problem, Aitchison [16] proposed the logarithmic ratio transformation method after studying compositional data.This method represents Aitchison's space data in a logarithmic ratio coordinate system and transforms it into Euclidean space in order to overcome the closure effect of compositional data, making it possible to analyze the data using multivariate statistical methods [17].This method has been proven effective in numerous research examples [18][19][20][21][22][23][24][25][26][27][28][29].Scholars generally believe that data transformed by this method can better reflect the true distribution patterns of data in space, accurately identifying the background and anomalous information of elements.
The geo-electrochemical methods have the capability to identify deep-seated anomaly information.However, due to various factors, such as geological activities and natural influences, and including the superimposition effect of multi-phase geological movements and the differential mobility of elements along fault migration, the geochemical patterns of elements remain ambiguous [3,7].As a result, identifying weak anomalies becomes challenging because they are often masked by high background values [30][31][32], leading to the misclassification of actual anomalies as background data.This misclassification can then affect the accurate delineation of the target areas.To address these issues, Cheng [33,34] has developed fractal and multi-fractal analysis methods to process geochemical data in order to differentiate geochemical elements as "anomaly" and "background".These methods can quickly and effectively extract information regarding mineral exploration, pollution, or other abnormal sources, while also facilitating a better understanding of the elements' geochemical distribution patterns.
Therefore, this study aims to analyze the compositional data of geo-electrochemical data in the study area through feasibility experiments in order to explore more potential meanings of the data through comparison with the raw data features.At the same time, in consideration of the multiple episodes of hydrothermal activity and the prolonged period of geological evolution in the study area, the overlapping of various geological processes has led to a complex and heterogeneous distribution of geochemical elements.In order to address this issue, this study intends to employ the S-A model to separate anomalies from the background, with the aim of better highlighting weak anomalies and determining the location of the target area.

Regional Geological Conditions
The Changjiang uranium ore field is situated south of the Zhuguang-Wanyang granitic rock (Figure 1).It occupies a geotectonic position at the junction of the post-Caledonian uplift in Fujian and Jiangxi provinces and the Hercynian-Indosinian depression in Hunan, Guangxi, and Guangdong provinces [35].The area is also located at the intersection of three uplift belts: Wanchangshan, Dayu-Jiufeng, and Zhuguang-Wanyang.During the Mesozoic era, the area experienced significant magmatic activity, leading to the formation of various intrusive rock bodies.The prominent ones include the intrusive rocks of the Zhuguangshan complex granite and the Guidong complex granite.
processes has led to a complex and heterogeneous distribution of geochemical elements.In order to address this issue, this study intends to employ the S-A model to separate anomalies from the background, with the aim of better highlighting weak anomalies and determining the location of the target area.

Regional Geological Conditions
The Changjiang uranium ore field is situated south of the Zhuguang-Wanyang granitic rock (Figure 1).It occupies a geotectonic position at the junction of the post-Caledonian uplift in Fujian and Jiangxi provinces and the Hercynian-Indosinian depression in Hunan, Guangxi, and Guangdong provinces [35].The area is also located at the intersection of three uplift belts: Wanchangshan, Dayu-Jiufeng, and Zhuguang-Wanyang.During the Mesozoic era, the area experienced significant magmatic activity, leading to the formation of various intrusive rock bodies.The prominent ones include the intrusive rocks of the Zhuguangshan complex granite and the Guidong complex granite.The study area is primarily composed of granite, supplemented by well-developed strata, with the exception of the Silurian.The Zhuguangshan intrusive complex is a multiphase and multi-stage intrusive body.Magmatism during the Indosinian period constitutes the primary component of this complex, exhibiting a north-south distribution at the surface.The magmatite includes fine-grained, spotted biotite granite from the Indosinian period, medium-grained, black biotite granite from the Yanshanian period, and some later intermediate basaltic dykes which are also exposed in the region [36].
The overall magmatic activity in the area was relatively frequent, mainly dominated by acidic magmatic activity, with accompanying intermediate-basic magmatic activity [37].The rocks were initially formed in the Caledonian and were influenced by large-scale tectonic magmatic intrusion activity during the Indosinian-Yanshan periods, resulting in widespread occurrences of muscovite-biotite granite, black mica granite, and black mica monzogranite.During the Yanshan period, there were also intrusions of diabase dikes.In addition, there were also alkali metasomatic rocks formed by late-stage alkali fluid activities [36].The uranium deposits and occurrences in the area are influenced by multiple periods and stages of magmatic intrusion and tectonic activity, often exhibiting obvious multi-period characteristics.

Geological Features of Study Area
The study area is situated in the southern region of the Zhu-Guang area, specifically within the Changjiang Uranium ore field.In terms of tectonic position, it is located at the intersection of the Fujian-Jiangxi fold belt, which formed after the Caledonian orogeny, The study area is primarily composed of granite, supplemented by well-developed strata, with the exception of the Silurian.The Zhuguangshan intrusive complex is a multi-phase and multi-stage intrusive body.Magmatism during the Indosinian period constitutes the primary component of this complex, exhibiting a north-south distribution at the surface.The magmatite includes fine-grained, spotted biotite granite from the Indosinian period, medium-grained, black biotite granite from the Yanshanian period, and some later intermediate basaltic dykes which are also exposed in the region [36].
The overall magmatic activity in the area was relatively frequent, mainly dominated by acidic magmatic activity, with accompanying intermediate-basic magmatic activity [37].The rocks were initially formed in the Caledonian and were influenced by large-scale tectonic magmatic intrusion activity during the Indosinian-Yanshan periods, resulting in widespread occurrences of muscovite-biotite granite, black mica granite, and black mica monzogranite.During the Yanshan period, there were also intrusions of diabase dikes.In addition, there were also alkali metasomatic rocks formed by late-stage alkali fluid activities [36].The uranium deposits and occurrences in the area are influenced by multiple periods and stages of magmatic intrusion and tectonic activity, often exhibiting obvious multi-period characteristics.

Geological Features of Study Area
The study area is situated in the southern region of the Zhu-Guang area, specifically within the Changjiang Uranium ore field.In terms of tectonic position, it is located at the intersection of the Fujian-Jiangxi fold belt, which formed after the Caledonian orogeny, and the Hercynian-Indosinian depression area in Hunan, Guangxi, and Guangdong provinces.
The geological history of the study area is characterized by significant tectonic and magmatic intrusions during the Indochina and Yanshan periods.Granite rocks, formed during these two periods, play a prominent role in the geological composition of the region.The overall magmatic activity within the area exhibits relatively frequent occurrences, predominantly characterized by acidic magmatic activity, concomitant with episodes of intermediate-basic magmatic events.The formation of the rock masses initially occurred during the Caledonian period, influenced by extensive tectonic-magmatic intrusion activities during the Indochina-Yanshan periods.This resulted in the widespread occurrence of diorite granite, biotite granite, and biotite-diorite granite, with intrusive occurrences of Yanshan period diabase veins.Additionally, there were later formations of alkali-altered rocks due to late-stage alkali metasomatism [36].The study area is situated at the intersection of the early Yanshan Changjiang rock mass and the Indochina Youdong granite body.The major rock types in the study area include Late Triassic two-mica granite and Middle Triassic two-mica granite.Additionally, there are several minor intrusion rocks and veins present.During the Yanshan period, intrusive diabase veins were also formed in the area (Figure 2).and the Hercynian-Indosinian depression area in Hunan, Guangxi, and Guangdong provinces.
The geological history of the study area is characterized by significant tectonic and magmatic intrusions during the Indochina and Yanshan periods.Granite rocks, formed during these two periods, play a prominent role in the geological composition of the region.The overall magmatic activity within the area exhibits relatively frequent occurrences, predominantly characterized by acidic magmatic activity, concomitant with episodes of intermediate-basic magmatic events.The formation of the rock masses initially occurred during the Caledonian period, influenced by extensive tectonic-magmatic intrusion activities during the Indochina-Yanshan periods.This resulted in the widespread occurrence of diorite granite, biotite granite, and biotite-diorite granite, with intrusive occurrences of Yanshan period diabase veins.Additionally, there were later formations of alkali-altered rocks due to late-stage alkali metasomatism [36].The study area is situated at the intersection of the early Yanshan Changjiang rock mass and the Indochina Youdong granite body.The major rock types in the study area include Late Triassic two-mica granite and Middle Triassic two-mica granite.Additionally, there are several minor intrusion rocks and veins present.During the Yanshan period, intrusive diabase veins were also formed in the area (Figure 2).The region has well-developed faults, including the NNE-trending Mianhuakeng Fault, the NW-trending Youdong Fault, the NNW-trending Ajiling Fault, the near NStrending Jiulongyan Fault, and the Shuishi Fault.In addition, there are various secondary faults in different directions.Among them, the NNE-trending Mianhuakeng Fault and the NW-trending Youdong Fault are two larger-scale faults in the region.The former is the main ore-controlling structure in the area, located on the northwest side of the survey area, and is often associated with mylonites, lamprophyres, and siliceous veins, commonly accompanied by silicification, sericitization, and chloritization; the latter is in the southern part of the study area, mainly composed of lamprophyre and silicified cataclastic rock.The junctions and compounds of the two faults have important control significance for uranium mineralization in the southern Zhuguang granitic rock area.The region has well-developed faults, including the NNE-trending Mianhuakeng Fault, the NW-trending Youdong Fault, the NNW-trending Ajiling Fault, the near NStrending Jiulongyan Fault, and the Shuishi Fault.In addition, there are various secondary faults in different directions.Among them, the NNE-trending Mianhuakeng Fault and the NW-trending Youdong Fault are two larger-scale faults in the region.The former is the main ore-controlling structure in the area, located on the northwest side of the survey area, and is often associated with mylonites, lamprophyres, and siliceous veins, commonly accompanied by silicification, sericitization, and chloritization; the latter is in the southern part of the study area, mainly composed of lamprophyre and silicified cataclastic rock.The junctions and compounds of the two faults have important control significance for uranium mineralization in the southern Zhuguang granitic rock area.
The Shulouqiu uranium deposit is a typical granite-type uranium deposit, which is located in the northern part of the Changjiang area and belongs to the southern part of the Zhuguangshan complex rock mass.This deposit has magmatic and tectonic environments similar to the Shulouqiu uranium deposit.The ore type of the Shulouqiu deposit is relatively uniform, primarily composed of porphyritic microcrystalline quartz (or varicolored microcrystalline quartz) and fluorite-uraninite-pyrite-type ores.Uraninite stands as the primary ore mineral in the Shulouqiu deposit, often presenting spherical or annular characteristics [35,36].The surrounding rock alterations have experienced multiple episodes of hydrothermal activity.Alterations associated with mineralization manifest as illitization, potassic alteration, sodic plagioclase alteration, hematization, chloritization, and carbonatization, among others.Among alterations related to mineralization, the most prominent include fluoritization, silicification, hematization, and clay mineralization.

Method and Technology
3.1.Geo-Electrochemical Technology and Test Analysis 3.1.1.The Geo-Electrochemical Method The geo-electrochemical method, also known as the "Partial Metal Extraction Method" [3], originated in the former Soviet Union.It was initially proposed by Ryss in the 1960s and corresponding field equipment was developed for its application in mineral exploration.In the book "Geo-electrochemical Exploration Method", published in 1983 [38], Ryss comprehensively discusses the theory, method, and application of this method and refers to a series of related methods, such as the geo-electrochemical method, elemental occurrence form method, and diffusion extraction method, collectively as the geo-electrochemical method.
The geo-electrochemical method breaks the dynamic balance of ions in underground rocks under the influence of an artificial electric field, causing anions and cations associated with mineralization to migrate towards the extraction electrode and eventually adsorb onto specific geo-electrochemical devices (foams) (Figure 3).By quantitatively analyzing the element content in the foam, the distribution pattern and geo-electrochemical characteristics can be determined, and the location of buried ore bodies can be inferred based on a comprehensive analysis of geological conditions, thereby achieving the goal of prospecting for mineral deposits [3,4,13].
similar to the Shulouqiu uranium deposit.The ore type of the Shulouqiu deposit is relatively uniform, primarily composed of porphyritic microcrystalline quartz (or varicolored microcrystalline quartz) and fluorite-uraninite-pyrite-type ores.Uraninite stands as the primary ore mineral in the Shulouqiu deposit, often presenting spherical or annular characteristics [35,36].The surrounding rock alterations have experienced multiple episodes of hydrothermal activity.Alterations associated with mineralization manifest as illitization, potassic alteration, sodic plagioclase alteration, hematization, chloritization, and carbonatization, among others.Among alterations related to mineralization, the most prominent include fluoritization, silicification, hematization, and clay mineralization.

The Geo-Electrochemical Method
The geo-electrochemical method, also known as the "Partial Metal Extraction Method" [3], originated in the former Soviet Union.It was initially proposed by Ryss in the 1960s and corresponding field equipment was developed for its application in mineral exploration.In the book "Geo-electrochemical Exploration Method", published in 1983 [38], Ryss comprehensively discusses the theory, method, and application of this method and refers to a series of related methods, such as the geo-electrochemical method, elemental occurrence form method, and diffusion extraction method, collectively as the geoelectrochemical method.
The geo-electrochemical method breaks the dynamic balance of ions in underground rocks under the influence of an artificial electric field, causing anions and cations associated with mineralization to migrate towards the extraction electrode and eventually adsorb onto specific geo-electrochemical devices (foams) (Figure 3).By quantitatively analyzing the element content in the foam, the distribution pattern and geo-electrochemical characteristics can be determined, and the location of buried ore bodies can be inferred based on a comprehensive analysis of geological conditions, thereby achieving the goal of prospecting for mineral deposits [3,4,13].

Measuring Line, Measuring Point Layout, and Device
According to the research objectives of the National Key Research and Development Program of China (2017YFC0602600), a total of seven northeast to southwest geoelectrochemical survey lines were completed in the Changjiang area, with a total area of 13.57km 2 , using a sampling interval of 800 × 100 m (Cj2, 4, 6) and 400 × 100 m (Cj7, 8, 9, 10), and a total of 241 samples were collected (Figure 2).
The piece of equipment used in this study is an independent power supply dipole geo-electrochemical device developed and improved by the Buried Ore Deposit Prediction Research Institute of the Guilin University of Technology.The device consists of a constant voltage and current-controlled power supply, refined carbon rods, polyurethane foam, and filter paper bands.After being connected by wires, the device can be used in the field.The extraction electrode is composed of highly conductive, refined carbon rods and pre-treated Minerals 2024, 14, 15 6 of 18 foam to eliminate the testing error that may be caused by the background value of the foam.All foams were pre-treated for normalization [3,4] with a mixture of 20% HCl and 0.5% (NH 2 ) 2 CS and a soaking time of 12 h.The geo-electrochemical technical parameters for this study are as follows: the power supply voltage was 9 V, the extraction pole spacing was 1 m, the extraction liquid was 10% diluted nitric acid (500 mL poured into each borehole), the power supply time was 24 h, and the depth of the extraction electrode buried in the soil horizon B was 30 cm.
Following the principle of "from known to unknown", this study conducted a feasibility test of the method and technique by setting sampling points on the known profile of the No. 2 exploration line of the Shulouqiu uranium deposit and conducted a comparative analysis by plotting the element content extracted using the geo-electrochemical method.

Fieldwork Process
The sampling work in this study is strictly in accordance with the technical specifications of geo-electrochemical extraction [3].Two cylindrical pits with a spacing of 1 m and a depth of approximately 30 cm on either side of a vertical measurement line in the field were excavated.Subsequently, 500 mL of extraction liquid was measured twice using a measuring cup, and the extraction electrodes were immersed in the liquid to ensure full contact with moisture.The electrodes were then buried in the excavated pits, parallel with each other, and the pits were sequentially backfilled in the order of soil extraction.The external power supply was connected, and the device was left to supply power continuously for 24 h.After retrieving the device, the foam samples were moved to the laboratory, individually packaged at each sampling point, and then sent for testing and analysis.

Sample Analysis and Testing
A total of 241 samples were collected for geo-electrochemical extraction in this study, including 241 low-voltage dipole extraction anode and cathode samples.Sample analysis was completed by the Nonferrous Metals Guilin Mineral Geological Testing Center of the Guilin Mineral Geological Research Institute, using the ICP-MS (Inductively Coupled Plasma Mass Spectrometry) analysis method.Ten elements were analyzed, including Ti, V, Co, Ni, As, Mo, Sb, Pb, Th, and U. Previous studies have shown that acid treatment can effectively reduce the background value of solid adsorbents in geoelectric extraction [39], which is much lower than the adsorbed element content in the sample.Therefore, the background value does not cause significant interference in the analysis of the elements.The samples are pretreated by wet digestion for analysis.Specifically, 0.2500 g of foam sample collected in the field was weighed and placed in a polyethylene beaker, and an appropriate amount of HNO 3 -H 2 O 2 was added.After being left overnight, it was heated on a 120 • C constant-temperature heating plate until it was dried, and then the residue was dissolved in aqua regia after cooling in the laboratory.The detailed testing procedures can be found in the related literature by Shi [39].

Geo-Electrochemical Technology and Test Analysis
Geochemical data, like other typical compositional data, exhibit a distinct characteristic known as the "closure effect" [40][41][42][43].Generally, the sum of different element contents in geochemical samples is a constant value, and the closure effect results in pseudo-correlations among geochemical elements, making the analysis results uncertain.To eliminate the negative effects caused by this effect, appropriate methods are usually applied to transform the data involved in the experiment [44,45].The transformed data are usually more interpretable in the application of classical statistical methods than the untransformed data [44,46,47].
Additive log-ratio transformation (alr), centered log-ratio transformation (clr), and isometric log-ratio transformation (ilr) are the three commonly used transformations for compositional data [16,17].However, each of the three transformation methods has its advantages and disadvantages.After the raw data are transformed by the alr or ilr method, one variable is reduced.The clr-transformed data have the same number of variables as the raw data, but the sum of the variables is zero.
Considering the advantages and disadvantages of the three transformation methods mentioned above, Filzmoser [44] proposed using the ilr transformation to open the closure effect of the data and then expressing the loadings of the ilr transformation in the clr space using a standard orthogonal basis.This overcomes the closure effect of the geochemical data and ensures the consistency of the number of variables before and after the transformation, which is more conducive to geochemists making reasonable interpretations.In addition, Filzmoser [44] also proposed the Robust Principal Component Analysis (RPCA) method to deal with outliers in the data, which can effectively avoid biased covariance matrix estimates due to the existence of outliers when using Principal Component Analysis (PCA).The theoretical formulas for each transformation method are presented below.
(1) Taking an original data matrix X ij (m × n, where m represents the number of samples and n denotes the number of types of elements) as an example, two distinct transformations will be presented.clr transformation: x is the elemental content at that point; j = 1,2,. ..,n; g(x i ) is the geometric mean value of the "i"-th sample.
(2) ilr transformation: The ilr transformation is derived based on the clr transformation and the standard orthonormal basis vectors "v j " [44].

Multi-Fractal S-A Model
The spatial distribution of geochemical elements is affected by changes in geological history, making the geochemical field an unsteady field [33].The complexity of the geochemical field is the result of the superposition of various geological processes over a long period of geological evolution.A large amount of geochemical data has shown that most geochemical anomalies exhibit fractal characteristics [48][49][50][51].In recent years, geologists have proposed various methods using fractal techniques to extract information on geochemical anomalies.Among these methods, the three methods proposed by Cheng, namely the element content-area fractal model (C-A model), energy spectrum density-area fractal model (S-A model), and the local singularity analysis method [31,[52][53][54][55][56][57][58], have been extensively employed and have produced substantial outcomes.
The S-A method can measure the general self-similarity of the anisotropy corresponding to the geochemical anomaly in the energy spectrum space, achieving the goal of separating anomalies [31,59].Currently, the S-A method has become one of the standard methods for anomaly decomposition.Extracting anomaly information using the S-A method exhibits a diversity of forms and demonstrates self-similarity characteristics in the frequency domain.This self-similarity can be characterized by an exponential model.

A(≥ S) ∝ S −β
(5) In Equation ( 5), S represents the spectrum energy density, A (≥S) represents the unit area greater than a certain threshold (S), β is the multi-dimensional fractal index of anisotropy, also known as the fractal dimension, which can be obtained on a log-log plot, and ∝ indicates proportionality.With the increase of the S value, the corresponding A will decrease, and the change rule mainly depends on the β fractal dimension.
To identify the background values and anomalies in the study area, this study employs ilr data transformation and compares the performance of elements on PCA and RPCA.This requires the selection of suitable principal components as the geochemical associations and using the S-A model to extract geochemical association anomalies and backgrounds in the study area.The software used in this study is the Arcfractal plugin within ArcGIS [60], developed by Zuo's team at the GPMR State Key Laboratory of China University of Geosciences (Wuhan, China).

Feasibility Analysis
The Shulouqiu deposit is situated in the upper left area of the study region (Figure 2).The deposit's host rocks consist of the Youdong and Changjiang rock mass.The Youdong rock mass is of the Indosinian granite variety, while the Changjiang rock body belongs to the Yanshan granite, mostly intruding the Youdong rock mass and primarily composed of medium-grained biotite granite.Ore-bearing faults within and surrounding the area are highly developed, mainly dominated by the NNW-trending ductile shear zones.The central section of the ore-bearing fault zone comprises silicified breccias, quartz veins, and other siliceous frameworks, while the sides exhibit successive hydrothermal alterations, including silicification, sericitization, and chloritization.
The experimental research carried out on the known profile can be found in Figure 4.
As can be seen from Figure 4, the anomalies of the main ore-forming elements U, Th, and Pb, which are highly correlated with U, extracted by geo-electrochemical are mostly reflected between the 6th and 12th measuring points (dashed line part in Figure 4) within the range of the verified drilling exploration work.The peak points are at the 8th and 11th measuring points.Meanwhile, multi-elements such as Co, Ni, V, and Ti also show different degrees of anomalies, which have a high degree of coincidence with the range of the concealed ore body, indicating good ore prospectivity.(Ore bodies and mineralized bodies are distinguished based on the minimum industrial grade and the minimum mining thickness.Anything meeting or exceeding the minimum industrial utilization criteria for each ore type is regarded as an ore body, while anything below is considered a mineralized body.) According to previous studies, the geochemical indicator w(Th)/w(U) value can be used as a reference in the exploration of Granite uranium deposits [61].Scholars [62] tend to believe that the w(Th)/w(U) value could reflect the degree of uranium activation.Granite with w(U)/w(Th) > 1/3 is in favor of uranium mineralization.Figure 4 shows that the anomaly range of the w(U)/w(Th) value corresponds well to the range of the known orebodies at depth, as well as element U.In summary, the geo-electrochemical method is feasible and effective for exploring granite-type uranium deposits and can be used in the exploration of concealed uranium orebodies in the Changjiang area.used as a reference in the exploration of Granite uranium deposits [61].Scholars [62] tend to believe that the w(Th)/w(U) value could reflect the degree of uranium activation.Granite with w(U)/w(Th) > 1/3 is in favor of uranium mineralization.Figure 4 shows that the anomaly range of the w(U)/w(Th) value corresponds well to the range of the known orebodies at depth, as well as element U.In summary, the geo-electrochemical method is feasible and effective for exploring granite-type uranium deposits and can be used in the exploration of concealed uranium orebodies in the Changjiang area.

Descriptive Statistical Analysis
In order to further explore the intrinsic connections among 10 elements in the study area, this study analyzed the statistical parameters of the raw data, logarithmically transformed data, and ilr-transformed data (Table 1).Obviously, from the coefficient of variation, the 10 elements all have large dispersion.U is the largest, followed by Co, Pb, As, and Mo.Indicating these elements may be enriched or depleted in local areas.
Based on the kurtosis and skewness of the raw data, it can be seen that all elements do not follow the normal distribution, causing the classical statistical analyses to produce erroneous results [44].By comparing the data characteristics after logarithmic and ilr transformation, it can be clearly seen that each element no longer has significantly high kurtosis or skewness, and the transformed data approximates the normal distribution.
Meanwhile, the Box plots and density curves (Figures 5 and 6) indicate the Pb, Ti, and U elements have a large number of high-value outlier points, which may suppress the anomalies of these elements in other areas of the study.The density curves of the data after logarithmic and ilr transformation are approximates of the normal distribution.

Principal Component Biplot
Figure 7 presents the biplots of the first two principal components (PC1 and PC2) obtained from PCA and RPCA analyses of the raw data, logarithmically transformed data and ilr-transformed data.As shown in Figure 7a,b, the biplots of the raw data and logarithmically transformed data exhibit a strong "one-sided" trend, indicating the presence of closure effects [58].However, after applying the ilr transformation, the biplot of the In summary, from the statistical parameters and the statistical figures, the logarithmically and ilr-transformed data are approximates of the normal distribution.Comparatively speaking, the ilr transformation has the more powerful effect in this case, because of the lower values of skewness and kurtosis, and a more concentrated tendency in the frequency distribution.

Principal Component Biplot
Figure 7 presents the biplots of the first two principal components (PC1 and PC2) obtained from PCA and RPCA analyses of the raw data, logarithmically transformed data, and ilr-transformed data.As shown in Figure 7a,b, the biplots of the raw data and logarithmically transformed data exhibit a strong "one-sided" trend, indicating the presence of closure effects [58].However, after applying the ilr transformation, the biplot of the transformed data shows a radial pattern, and the closure effect is eliminated.
Minerals 2023, 13, x FOR PEER REVIEW 12 geochemical association [61], and the latter is a geological background element ass tion, which are the primary ore-forming elements.U is the dominant ore-forming elem in the study area, while Th is typically a companion element to uranium mineralizat

Geochemical Association Mapping
For further comparison, the geochemical association mapping of PC1 from PCA RPCA are carried out.The used interpolation method is Inverse Distance Weigh (IDW).Obviously, the geochemical association maps of PC1 from PCA using raw and logarithmically transformed data are relatively similar (Figure 8a,b).They both hibit high-value anomalies that occupy almost half of the study area, with the anom distributed in two main blocks in the east and west.In contrast, after applying th transformation, the PC1 for classical and robust methods are relatively similar (Fi 8c,d).This transformation better highlights the anomalies in the north and south that w previously suppressed.
Based on the geological map (Figure 2), the large anomalous area on the west o transformed PC1 (Figure 8c,d) may be due to the intersection of the Late Triassic twogranite, Middle Triassic two-mica granite, and Early Jurassic biotite granite.Intru Specifically, the biplots of the raw data and logarithmically transformed data (Figure 7a,b) show similar distributions, with all elements having positive loadings on PC1 and Ti, Th, and V having negative loadings on PC2.In contrast, the biplots of the ilr-transformed data (Figure 7c,d) exhibit a completely different distribution pattern, with a more dispersed variable distribution, indicating that the closure effect has been eliminated.
The biplot of PCA and RPCA is quite similar, with all ten elements evenly distributed in the four quadrants.However, there are still subtle differences; for example, the Th, U, Pb, As, and Co have positive loadings on PC1 in classical PCA, while Th, U, and Pb have positive loadings on PC1 in RPCA.Based on the geochemistry, Pb is a decay product of U, and their geochemical behavior should be more similar.The primary ore mineral in the study area is pitchblende uranium ore, comprising mainly UO 2 , ThO 2 + REE 2 O 3 , PbO, SiO 2 , and CaO [35].Therefore, the elements U and Pb exhibit a closer relationship in this region.
Geochemical association of Th-U-Pb is a very important indicator for uranium deposit exploration [35], pointing to the fact that the RPCA can gain more reliable results in this case (Figure 7).
Accordingly, based on the biplot of RPCA, 10 elements can be divided into two distinct element groups: Th-U-Pb and Ti-V-Co-Ni-Mo-Sb-As; the former is a mineralization geochemical association [61], and the latter is a geological background element association, which are the primary ore-forming elements.U is the dominant ore-forming element in the study area, while Th is typically a companion element to uranium mineralization.

Geochemical Association Mapping
For further comparison, the geochemical association mapping of PC1 from PCA and RPCA are carried out.The used interpolation method is Inverse Distance Weighting (IDW).Obviously, the geochemical association maps of PC1 from PCA using raw data and logarithmically transformed data are relatively similar (Figure 8a,b).They both exhibit highvalue anomalies that occupy almost half of the study area, with the anomalies distributed in two main blocks in the east and west.In contrast, after applying the ilr transformation, the PC1 for classical and robust methods are relatively similar (Figure 8c,d).This transformation better highlights the anomalies in the north and south that were previously suppressed.shape is distributed mainly along the contact boundary.The small anomalous area on the north is located at the contact point between Middle Triassic two-mica granite and Early Jurassic biotite granite.The small anomalous area on the east is also related to the intrusion of the Cretaceous two-mica granite.Additionally, the anomalous area on the south side is associated with the faults.Overall, the anomalous distribution areas of ilr-transformed PC1 are controlled by the complex secondary faults, the Yanshan-Indosinian intrusive granite, and the contact zone of different granite (Figure 8c,d).

Multifractal S-A Model Results
Based on the results above, the PC1 from RPCA is used to further determine the geochemical association anomalies in this study area.The S-A model is used to separate the geochemical background and anomalies.The filtering values are determined based on the spectrum energy density (S) and area (A) plot (Figure 9).
The inflection point values were determined by the curve segment fitting (specifi- Based on the geological map (Figure 2), the large anomalous area on the west of ilrtransformed PC1 (Figure 8c,d) may be due to the intersection of the Late Triassic two-mica granite, Middle Triassic two-mica granite, and Early Jurassic biotite granite.Intrusive rocks of the Yanshan and Indosinian periods converge in this area, and the anomalous shape is distributed mainly along the contact boundary.The small anomalous area on the north is located at the contact point between Middle Triassic two-mica granite and Early Jurassic biotite granite.The small anomalous area on the east is also related to the intrusion of the Cretaceous two-mica granite.Additionally, the anomalous area on the south side is associated with the faults.Overall, the anomalous distribution areas of ilr-transformed PC1 are controlled by the complex secondary faults, the Yanshan-Indosinian intrusive granite, and the contact zone of different granite (Figure 8c,d).

Multifractal S-A Model Results
Based on the results above, the PC1 from RPCA is used to further determine the geochemical association anomalies in this study area.The S-A model is used to separate the geochemical background and anomalies.The filtering values are determined based on the spectrum energy density (S) and area (A) plot (Figure 9).

. Mineral Prospectivity Analysis
The uranium mineralization in the southern Zhuguang region is closely associated with mantle-crust magmatic processes, fault tectonics, and hydrothermal activities.Under the influence of mantle fluids, pre-Cambrian uranium-bearing metamorphic basement rocks underwent re-melting, resulting in the Indosinian and early Yanshan uranium-rich granite intrusion.Various metamorphic events subsequently activated uranium within these rock bodies.Intense late Yanshan tectonic activities facilitated the upward movement of mantle alkaline fluids, causing the intrusion of acidic magmas from The inflection point values were determined by the curve segment fitting (specifically, refer to Figure 9, where "lg" denotes the logarithm with a base of 10).Three straight lines are derived through the lgS-lgA plot and are separated with the inflection points S 1 = 4.60 and S 2 = 7.55 ("lg" represents the logarithm with a base of 10).The background and anomaly maps are obtained using the inflection value and inverse Fourier transformation.In these maps, S < S 1 represents noise, S > S 2 represents backgrounds, and S 1 < S < S 2 represents anomalies (Figures 9 and 10).
By comparing Figure 10a with Figure 8d, it can be observed that the anomalies in the western part are weakened and the anomalies in the eastern part are strengthened after the S-A anomaly decomposition.At the same time, the background area can be divided into one low background and one high background area, where the blue low-background part is mainly distributed in the Middle Triassic two-mica granite, and the red high-background part is almost distributed around the contact surface of different periods of the granite.

. Mineral Prospectivity Analysis
The uranium mineralization in the southern Zhuguang region is closely associated with mantle-crust magmatic processes, fault tectonics, and hydrothermal activities.Under the influence of mantle fluids, pre-Cambrian uranium-bearing metamorphic basement rocks underwent re-melting, resulting in the Indosinian and early Yanshan uranium-rich granite intrusion.Various metamorphic events subsequently activated uranium within these rock bodies.Intense late Yanshan tectonic activities facilitated the upward movement of mantle alkaline fluids, causing the intrusion of acidic magmas from crustal or crust-mantle hybrid sources and mafic magmas from the mantle.This further transformed and activated uranium within the rock bodies, providing abundant thermal energy, hydrothermal fluids, and mineralizers for uranium mineralization.Residual magma or other fluids during this period, influenced by fault tectonics, extracted uranium from within and surrounding granitic rocks, preferentially precipitating and enriching ores under favorable conditions [36].
The area exhibits highly developed ore-controlling structures.The occurrence and scale of uranium deposits are often controlled by tectonic alteration zones.They are primarily distributed within the fracture zones of the granitic intrusions, with some located

. Mineral Prospectivity Analysis
The uranium mineralization in the southern Zhuguang region is closely associated with mantle-crust magmatic processes, fault tectonics, and hydrothermal activities.Under the influence of mantle fluids, pre-Cambrian uranium-bearing metamorphic basement rocks underwent re-melting, resulting in the Indosinian and early Yanshan uranium-rich granite intrusion.Various metamorphic events subsequently activated uranium within these rock bodies.Intense late Yanshan tectonic activities facilitated the upward movement of mantle alkaline fluids, causing the intrusion of acidic magmas from crustal or crust-mantle hybrid sources and mafic magmas from the mantle.This further transformed and activated uranium within the rock bodies, providing abundant thermal energy, hydrothermal fluids, and mineralizers for uranium mineralization.Residual magma or other fluids during this period, influenced by fault tectonics, extracted uranium from within and surrounding granitic rocks, preferentially precipitating and enriching ores under favorable conditions [36].
The area exhibits highly developed ore-controlling structures.The occurrence and scale of uranium deposits are often controlled by tectonic alteration zones.They are primarily distributed within the fracture zones of the granitic intrusions, with some located in the outer zones of the granitic intrusions, under the influence of fault structures and the intrusion boundaries of the granites.
The main uranium-enriched rock mass in the study area is the Late Triassic biotite granite, particularly its conjunction with other rock masses, which favors mineralization.The NNW and NE trending faults serve as primary channels for hydrothermal migration, thereby facilitating mineralization.The positions of intersection and combination of faults in various directions are of particular interest.Additionally, the Yanshan period intrusion of the Cretaceous hydrothermal veins also plays a proactive role in uranium enrichment and mineralization.
Based on the processed RPCA1 anomaly map obtained through the S-A model, the high anomaly areas align well with the aforementioned favorable geological characteristics for mineralization.It is considered that this method effectively identifies the enrichment of uranium and associated elements, providing indicative significance in delineating target areas.

Target Division
The favorable prospecting geological features are closely related to the distribution of faults, contact surfaces of granite, and intrusions of basic veins in the research area.Combining the geochemical indicators with geological background, the distribution of the PC1 from RPCA are significantly important indicators for mineral prospectivity mapping.Finally, based on this understanding, four prospecting prediction targets were identified in the study area (Figure 11).

Target Division
The favorable prospecting geological features are closely related to the distribution of faults, contact surfaces of granite, and intrusions of basic veins in the research area.Combining the geochemical indicators with geological background, the distribution of the PC1 from RPCA are significantly important indicators for mineral prospectivity mapping.Finally, based on this understanding, four prospecting prediction targets were identified in the study area (Figure 11).Target I This target spans two survey lines, with an area of approximately 1.06 km 2 (Figure 11).PC1 displays a high-value anomaly in this area (Figure 10).The NW-oriented Youdong fault and two faults-oriented NNW jointly control this target.Meanwhile, a diabase vein passes through the center of the anomaly along the Youdong fault, and the NE side of the anomaly is the boundary of Early Jurassic two-mica granite and Late Triassic Target I This target spans two survey lines, with an area of approximately 1.06 km 2 (Figure 11).PC1 displays a high-value anomaly in this area (Figure 10).The NW-oriented Youdong fault and two faults-oriented NNW jointly control this target.Meanwhile, a diabase vein passes through the center of the anomaly along the Youdong fault, and the NE side of the anomaly is the boundary of Early Jurassic two-mica granite and Late Triassic two-mica granite, with a complex boundary shape, which is conducive to the enrichment of U mineralization.Therefore, it is classified as the target I.
Target II This target area is located near the Ajiling faults, spanning two survey lines (Cj4 and Cj6), with an area of approximately 1.148 km 2 .After S-A decomposition, the anomaly is significantly strengthened, the center of the anomaly is located at the intersection of the NNW-oriented and near-NS-oriented faults.Meanwhile, a basalt vein passes through the target.Therefore, it is classified as the target II.

Target III
This target lies in the northern part of the study area, near the Mianhuakeng fault, and covers an area of approximately 0.266 km 2 .The III target lies at the intersection of the Mianhuakeng fault and the NNW-oriented fault, a basalt vein of passes through the southern side, and the contact zone of Late Triassic two-mica granite and Early Jurassic biotite granite appear in this target.

Target IV
This target area is situated close to the Tuobeiling, with a smaller area of approximately 0.271 km 2 .The IV target is mainly located on the Youdong fault, and a basaltic dyke also passes through this area.The intersection of the NNW-oriented fault and the Youdong fault is located on the east side.

Figure 2 .
Figure 2. Geological map of the study area.

Figure 2 .
Figure 2. Geological map of the study area.

Figure 4 .
Figure 4. Known exploration line of the Shulouqiu uranium deposit.Figure 4. Known exploration line of the Shulouqiu uranium deposit.

Figure 4 .
Figure 4. Known exploration line of the Shulouqiu uranium deposit.Figure 4. Known exploration line of the Shulouqiu uranium deposit.

Figure 6 .
Figure 6.Density curves of (a) logarithmically transformed and (b) ilr-transformed datasets of metallogenetic elements.4.3.CoDA Results 4.3.1.Principal Component Biplot Figure 7 presents the biplots of the first two principal components (PC1 and PC2) obtained from PCA and RPCA analyses of the raw data, logarithmically transformed data, and ilr-transformed data.As shown in Figure 7a,b, the biplots of the raw data and

Figure 7 .
Figure 7. Biplots of the PC1 and PC2 were obtained by the raw dataset with PCA (a) and loga mically transformed dataset with PCA (b) and ILR-transformed dataset with PCA (c) and ILR-t formed dataset with RPCA (d).

Figure 7 .
Figure 7. Biplots of the PC1 and PC2 were obtained by the raw dataset with PCA (a) and logarithmically transformed dataset with PCA (b) and ILR-transformed dataset with PCA (c) and ILR-transformed dataset with RPCA (d).

Figure 8 .
Figure 8. Geochemical association maps of the PC1 (a) using the raw data, (b) using logarithmically transformed data, (c) using ilr-transformed data PCA, and (d) RPCA using ilr-transformed data.

Figure 8 .
Figure 8. Geochemical association maps of the PC1 (a) using the raw data, (b) using logarithmically transformed data, (c) using ilr-transformed data PCA, and (d) RPCA using ilr-transformed data.

Minerals 2023 ,
13, x FOR PEER REVIEW 14 of 19 background part is almost distributed around the contact surface of different periods of the granite.

Figure 10 .
Figure 10.Geochemical association anomaly (a) and background (b) obtained by S-A model.

Figure 10 .
Figure 10.Geochemical association anomaly (a) and background (b) obtained by S-A model.

Figure 11 .
Figure 11.Proposed prospecting prediction targets in the study area.

Figure 11 .
Figure 11.Proposed prospecting prediction targets in the study area.

Table 1 .
Statistical parameters of geo-electric chemistry element content (num = 241).The coefficient of variation is the ratio of the standard deviation to the mean value of the raw data; the mass fraction unit for all elements is 10 −6 .