For EMC simulation, the vector fitting model is transformed into time domain state equation model. Then the system is balanced. Then the reduced model can be obtained by removing the states corresponding to the small HSV (Hankel Singular Values). The order of the reduced model is determined by the singular curvature spectrum. Finally, balanced truncation model reduction method is used on a sample chip PDN (Power Distribution Network) and compared with the performance of the AWE (Asymptotic Waveform Evaluation) method. Simulation results show that the proposed method can operate in a wide frequency range and has smaller error and faster speed.