The timescale of thermal quench (TQ) remains a long-lasting issue in tokamak plasmas, which has not been fully understood yet. In this work, based on our previous thermal diffusion model and further considering the scattering caused by the electrostatic turbulence into account, a more generalized evolution equation of the electron temperature is derived from the electron drift-kinetic equation. This equation is applicable to all collisional regimes and a wide range of the stochastic magnetic fields (SMF) amplitude . On the one hand, the heat flux induced by the drift is comparable to that induced by the SMF under DIII-D parameters, but becomes negligible under International Tokamak Experimental Reactor (ITER)-like parameters. However, the overall impact of electrostatic turbulence appears to be insignificant in both scenarios. This is because the direct transport caused by drift is basically counterbalanced by the reduction of the transport induced by the SMF due to the modification effect of electrostatic turbulence. On the other hand, numerical results indicate the timescale of the TQ is strongly dependent on the SMF amplitude and can reach the order of 100 when . The scaling of TQ timescale being approximately proportional to device size remains invariant under varying SMF amplitudes, and the extrapolations from DIII-D to ITER-like case based on this scaling demonstrate alignment with numerical results. It is also shown that the collisionality-dependent and SMF amplitude-dependent of thermal diffusivity should be carefully considered to predict the TQ timescale. Moreover, it is also discussed that SMF amplitude exceeding 4 × 10−3 would lead to unacceptable thermal loading on divertor target during TQ in ITER.