A minimal form of quasi-linear approximation (QLA), recently proposed with a stochastic forcing and POD modes (Hwang & Ekchardt, J. Fluid Mech., 2020, 894:A23), has been extended by employing a resolvent framework. A particular effort is made to reach an extremely high Reynolds number by carefully controlling the approximation without loss of the general scaling properties in the spectra, while setting out the main limitations and accuracy of the proposed QLA with possibility of further improvement. The QLA is subsequently applied to turbulent channel flow up to Reτ = 10 6 (Reτ is the friction Reynolds number). While confirming that the logarithmic wallnormal dependence in streamwise and spanwise turbulence intensities robustly appears, it reveals some non-trivial difference from the scaling of the classical attached eddy model based on inviscid flow assumption. Firstly, the spanwise wavenumber spectra do not show any well-visible inverselaw behaviour due to the viscous wall effect prevailing in a significant portion of the lower part of the logarithmic layer. Secondly, the near-wall peak streamwise and spanwise turbulence intensities are found to deviate from ln Reτ scaling for Reτ 10 4 . Importantly, the near-wall streamwise turbulence intensity is inversely proportional to 1/U + cl (U + cl is the inner-scaled channel centreline velocity), consistent with the scaling obtained from an asymptotic analysis of the Navier-Stokes equations (Monkewitz & Nagib, J. Fluid Mech. 2015, 783:474-503). The same behaviour was also observed for the streamwise turbulence intensity in the logarithmic region, as was predicted by the asymptotic analysis. Finally, the streamwise turbulence intensity in the logarithmic region is found to become greater than the near-wall one at Reτ O(10 5 ). It is shown that this behaviour originates from the near-wall spectra associated with large-scale inactive motions, the intensity of which gradually decays as Reτ → ∞.