The X-ray emission from blazars has been widely investigated using several space telescopes. In this work, we explored statistical properties of the X-ray variability in the blazars S5 0716+714, OJ 287, Mkn 501, and RBS 2070 using the archival observations from the XMM-Newton telescope between the period 2002-2020. Several methods of timing and spectral analyses, including fractional variability, minimum variability time-scale, power spectral density analyses, and countrate distribution, were performed. In addition, we fitted various spectral models to the observations, as well as estimated hardness ratio. The results show that the sources are moderately variable within the intraday time-scale. Three of the four sources exhibited a clear bi-modal pattern in their countrate distribution, revealing possible indication of two distinct countrate states, that is, hard and soft countrate states. The slope indices of the power spectral density were found to be centred around 0.5. Furthermore, the spectra of the sources were fitted with single power law, broken power law, log-parabolic, and blackbody + log-parabolic models (the latter only for OJ 287). We conclude that for most of the observations log-parabolic model was the best fit. The power-spectral-density analysis revealed the variable nature of PSD slopes in the source light curves. The results of this analysis could indicate the non-stationary nature of the blazar processes on intraday time-scales. The observed features can be explained within the context of current blazar models, in which the non-thermal emission mostly arises from kilo-pc scale relativistic jets.