Based on Hubbard model with the Hartree-Fock approximation, we study the properties of quantum discord (QD) between the nearest-neighbor sites A and B in zigzag graphene nanoribbons thermalized with a reservoir at temperature T. Several influences of the site position, on-site Coulomb repulsion U, temperature, and ribbon width on QD are discussed in detail. The results show that QD is robust against thermal fluctuations, and QD for the leg pairs along the zigzag chain near ribbon edges is always larger than that for the rung pairs linking two adjacent zigzag chains. QD for the rung pairs increases and then approaches to saturation as the ribbon width increases, where the velocity of saturation is strongly correlated to U. Moreover, for rung pairs the values of U at the QD peaks perform the scaling behaviors with increasing ribbon width.