The Guilaizhuang gold deposit in Pingyi County, Shandong Province, is the first discovered tellurium-gold type epithermal deposit in China. Its deep and peripheral areas still hold significant prospecting potential, and there is an urgent need to establish an effective exploration model. In this study, 672 primary halo samples collected from 58 exploration drill holes at the Guilaizhuang gold deposit were analyzed. A three-dimensional primary halo analysis method was adopted to overcomes the limitations of traditional two-dimensional approaches. Through multivariate statistical analysis, primary halo indicator elements were identified. A 3D geochemical anomaly model was constructed to reveal the characteristics of primary halos and their correlation with mineralization while evaluating the prospecting for concealed deep ore bodies. The results show that: (1) The Guilaizhuang gold deposit exhibits an axial zoning sequence of As-Sb-Se (frontal halo), Ag-Au-Te-Pb-Zn (proximal halo), and Mo-V (rear halo). Notably, rear halo elements (Mo, V) are enriched at the top of the orebody clusters, showing a distinctive "reverse zoning" phenomenon. (2) Gold-indicating element ratios (Ag/Au, As/Sb, Zn/Pb) display a NE-SW directional zoning, consistent with the strike of the structural-cryptoexplosive breccia zone controlled by the Guilaizhuang F1 fault. (3) The frontal-rear halo multiplicative ratio[w(As×Sb×Se)/w(Mo×V)] exhibits anomalous inflection points at elevations of -750 m and -1700 m, indicating the presence of concealed orebodies at depth. The findings suggest that the mineralization at Guilaizhuang involved multi-stage hydrothermal fluids superposition. Ore-forming fluids originating from the deep southwest part of the deposit migrate upward along the structural conduits formed by the Guilaizhuang F1 fault to the near-surface, forming structural-cryptoexplosive breccia type and quasi-layered carbonate-hosted type gold orebodies. It is inferred thar the deep prospecting potential of the deposit remains promising, with the deep southern region below -750 m elevation along exploration line 35 identified as a priority target zone for further exploration.