In this paper, we establish Pareto $Z$-eigenvalue inclusion intervals of tensor eigenvalue complementarity problems based on the spectral radius of symmetric matrices deduced from the provided tensor. Numerical examples are suggested to demonstrate the effectiveness of the results. As an application we offer adequate criteria for the strict copositivity of symmetric tensors.