The effective integration of geochemical data with multisource geoscience data is a necessary condition for mapping mineral prospects. In the present study, based on the maximum entropy principle, a maximum entropy model (MaxEnt model) was established to predict the potential distribution of copper deposits by integrating 43 ore-controlling factors from geological, geochemical and geophysical data. The MaxEnt model was used to screen the ore-controlling factors, and eight ore-controlling factors (i.e., stratigraphic combination entropy, structural iso-density, Cu, Hg, Li, La, U, Na₂O) were selected to establish the MaxEnt model to determine the highest potential zone of copper deposits. The spatial correlation between each ore-controlling factor and the occurrence of a copper mine was studied using a response curve, and the relative importance of each ore-controlling factor was determined by jackknife analysis in the MaxEnt model. The results show that the occurrence of copper ore is positively correlated with the content of Cu, Hg, La, structural iso-density and stratigraphic combination entropy, and negatively correlated with the content of Na₂O, Li and U. The model’s performance was evaluated by the area under the receiver operating characteristic curve (AUC), Cohen’s maximized Kappa and true skill statistic (TSS) (training AUC = 0.84, test AUC = 0.8, maximum Kappa = 0.5 and maximum TSS = 0.6). The results indicate that the model can effectively integrate multi-source geospatial data to map mineral prospectivity.

The smoothing effect of data interpolation could cause useful information loss in geochemical mapping, and the uncertainty assessment of geochemical anomaly could help to extract reasonable anomalies. In this paper, multiple-point geostatistical simulation and local singularity analysis (LSA) are proposed to identify regional geochemical anomalies and potential mineral resources areas. Taking Cu geochemical data in the Mila Mountain Region, southern Tibet, as an example, several conclusions were obtained: (1) geochemical mapping based on the direct sampling (DS) algorithm of multiple-point geostatistics can avoid the smoothing effect through geochemical pattern simulation; (2) 200 realizations generated by the direct sampling simulation reflect the uncertainty of an unsampled value, and the geochemical anomaly of each realization can be extracted by local singularity analysis, which shows geochemical anomaly uncertainty; (3) the singularity-quantile (S-Q) analysis method was used to determine the separation thresholds of E-type α, and uncertainty analysis was carried out on the copper anomaly to obtain the anomaly probability map, which should be more reasonable than the interpolation-based geochemical map for geochemical anomaly identification. According to the anomaly probability and favorable geological conditions in the study area, several potential mineral resource targets were preliminarily delineated to provide direction for subsequent mineral exploration.