The complexity of the spatial-temporal structure of metallogenic element distribution is one of the hot topics in metallogenic system research. Exploring the effective approach to quantitatively characterize the spatial-temporal structural complexity is of huge scientific significance to unravel the dispersion and aggregation mechanisms of metallogenic elements in the formation and evolution processes of the metallogenic system. In this paper, according to the characteristics of geological data, a multifractal detrended fluctuation analysis base on fractal interpolation (FIMFDFA) is proposed, and the effectiveness of this method is verified using de Wijs model. The multifractal detrended fluctuation analysis base on fractal interpolation (FIMFDFA) model is applied to describe the multifractal characteristics of metallogenic element (Cu-Pb-Zn-Mo) distributions and mineralization feature of the No.16 borehole, from Jiama porphyry copper deposit in Tibet, China. The results show that four metallogenic elements exhibit multifractal scaling, indicating that the distribution of metallogenic element concentrations presents a complex nonlinear structure, and the multifractal strength of elements with different mineralization grades is different. The multifractal intensity from large to small is medium mineralization, strong mineralization, and weak mineralization, which reflects the difference in multifactorial parameters and has a certain indication for mineralization. The quantitative analysis of element distribution complexity and the study of mineralization indication will help us to better understand the metallogenic process and provide a reference to further evaluation of regional metallogenic potential and the implement of mining activity in the industry.