In this paper, the secular full-orbit simulations of runaway electrons with synchrotron radiation in tokamak fields are carried out using a relativistic volume-preserving algorithm. Detailed phase-space behaviors of runaway electrons are investigated in different dynamical timescales spanning 11 orders. In the small timescale, i.e., the characteristic timescale imposed by Lorentz force, the severely deformed helical trajectory of energetic runaway electron is witnessed. A qualitative analysis of the neoclassical scattering, a kind of collisionless pitch-angle scattering phenomena, is provided when considering the coupling between the rotation of momentum vector and the background magnetic field. In large timescale up to 1 s, it is found that the initial condition of runaway electrons in phase space globally influences the pitch-angle scattering, the momentum evolution, and the loss-gain ratio of runaway energy evidently. However, the initial value has little impact on the synchrotron energy limit. It is also discovered that the parameters of tokamak device, such as the toroidal magnetic field, the loop voltage, the safety factor profile, and the major radius, can modify the synchrotron energy limit and the strength of neoclassical scattering. The maximum runaway energy is also proved to be lower than the synchrotron limit when the magnetic field ripple is considered.