A new method is proposed for the simulation of nonstationary stochastic earthquake ground motion based on orthogonal Hilbert-Huang transform (OHHT) spectra of its sample observations. Orthogonal EMD method is introduced to treat with IMF. OHHT avoids leakage of energy compared with conventional HHT, and it can be used for analysis of nonstationary signal. The average of the Hilbert spectra over the samples is defined as the Hilbert spectrum of the earthquake motion process and used as the target in the simulation of the process. Random Phase is also introduced to simulate nonstationary earthquake. The statistic characteristics function of stochastic process is given in the paper. Simulations of two earthquake ground motion are carried out. The results show that the samples can accurately reflect amplitude and frequency nonstationary of records. Samples and record of same earthquake ground motion process have same statistic characteristics.