The relativistic corrections to the color-octet J/ψ hadroproduction at the Tevatron and LHC are calculated up to O(v 2 ) in nonrelativistic QCD factorization frame. The short distance coefficients are obtained by matching full QCD with NRQCD results for the partonic subprocessThe short distance coefficient ratios of relativistic correction to leading order for color-octet1 , and 3 P [8] J at large p T are approximately -5/6, -11/6, and -31/30, respectively, for each subprocess, and it is 1/6 for color-singlet state 3 S[1]1 . If the higher order long distance matrix elements are estimated through velocity scaling rule with adopting v 2 = 0.23 and the lower order long distance matrix elements are fixed, the leading order cross sections of color-octet states are reduced by about a factor of 20 ∼ 40% at large p T at both the Tevatron and the LHC. Comparing with QCD radiative corrections to color-octet states, relativistic correction is ignored along with p T increasing. Using long distance matrix elements extracted from the fit to J/ψ production at the Tevatron, we can find the unpolarization cross sections of J/ψ production at the LHC taking into account both QCD and relativistic corrections are changed by about 20 ∼ 50% of that considering only QCD corrections. These results indicate that relativistic corrections may play an important role in J/ψ production at the Tevatron and the LHC. PACS numbers: 12.38.Bx,12.39.St,13.85.Ni * Electronic address: