The scattering of an arbitrarily incident Bessel beam by complex particles with inclusions was investigated by utilizing the method of moments and its fast algorithm based on surface integral equations. The mathematical description of the electromagnetic field components of the arbitrarily incident Bessel beams were derived by using the vector expressions of zero-order Bessel beam in combination with rotation Euler angles. Based on equivalence principle, the surface integral equations were established, which could solve the scattering problems involving complex particles with inclusions of arbitrary shape and structure. The scattering of an arbitrarily incident Bessel beam by several selected particles with inclusions was simulated. The numerical results indicate that the Differential Scattering Cross Section (DSCS) for a Bessel beam is smaller than that for a plane wave and increases with the decreasing of the half-cone angle. Moreover, the position offset of the beam center decreases the DSCS, and the DSCS will be quite sensitive to the incident angles. The obtained results are expected to provide useful help for techniques of laser detection on particle, diagnosis, and manipulation.