We numerically evaluated the performance of multi-coil shimming in human brain using an array of up to 1,000 coils with realistic off-resonance distributions of 37 healthy participants. The average and variation of shim current distributions were revealed. Singular Value Decomposition suggested orthogonal current modes to reduce off-resonance. The first 6 current modes accounted for about 90% of the variance of shim current distributions. They achieved shimming performance comparable to the 5th-order and 4th-order spherical harmonic in global and slice-selective shimming, respectively.