Human serum albumin (HSA) has seven common fatty acid (FA) binding sites. In this study, we used the molecular mechanics Poisson-Boltzmann surface area method to identify high affinity FA binding sites on HSA in terms of binding free energy. Using multiple HSA-FA (myristate, palmitate) complex models constructed by molecular dynamics simulations, two methods were performed in molecular mechanics Poisson-Boltzmann surface area, the "three-trajectory method" and the "single-trajectory method". The former, which is less precise than the latter but may be more accurate as it includes the effects of conformational change upon binding, was used to classify high and low affinity sites. As a result, Sites 2, 4, and 5 were identified as high affinity sites for both FAs. The latter method, which is precise because energies are calculated from snapshots of the same trajectory for HSA-FA complex, was performed to compare the magnitude of binding free energy for these sites. The order of magnitude was 5 > 4 > 2, identical to that of a previous publication by others. In this way, a combination of the two methods was effectively used to identify high affinity sites. This study therefore provides an insight into the quantitative identification of high affinity FA binding sites on HSA.