Seismic load is distributed to modes of a structure through the modal participation factor(MPF). The modal participation factor is essential to analyze structural response under earthquake load. MPF of a real structure differs from that of analytical mathematical model due to the error induced from analytical assumptions and during the construction. In this study, an identification method is proposed to calculate the 1st MPF of real structure based on $H^{infty}$ optimal model reduction. The MPF is obtained from the relationship between observability and controllability matrices realized from system identification and those of a prototype 2-degree state space model. The proposed method is verified thorough numerical examples.