The k-space unfolding matrix for parallel imaging in MRI was examined and was found to have the structure of block circulant quasi-band, providing a rigorous mathematical justification for the interpolation kernel in the commonly used GRAPPA method. The optimal GRAPPA kernel size and dimension were closely related to sensitivity spectrum.