In this paper, we focus on the spectra numerical computation for the integral equation with the absolute oscillation and power-law or logarithmic singularity. Finite section method is applied to transform the integral equation to an algebraic eigenvalue problem. The entries of the coefficient matrix appearing in the bivariate highly oscillatory singular integrals can be represented explicitly in Gamma or the exponential integral functions. The decay rate of the entries is established to construct the truncation scheme. Then the infinite algebraic eigenvalue problem can be simplified to be the finite one. The corresponding error of the infinite and finite algebraic systems is also bounded. Finally, the numerical experiments are provided to illustrate the theoretical analysis.