A thermodynamic model has been developed to predict the distribution behavior of Ni, Cu, Co, Fe, S, As, Sb, and Bi in nickel smelting and direct high-grade nickel matte smelting processes. The model has been validated by numerous experimental data and industrial data with a wide range of operating conditions. The effect of operating conditions on the distributions of Ni, Cu, Co, As, Sb, and Bi among the gas, matte, and slag phases has been investigated. It was found that the distribution behavior of Ni, Co, Cu, As, Sb, and Bi in the nickel smelting furnace depends on process parameters such as the smelting temperature, matte grade, oxygen enrichment, Fe/SiO 2 ratio in the slag, Cu/Ni ratio in charge, and oil/air ratio. The parameters also have an influence on the behavior of Fe 3 O 4 in the slag.