This correspondence presents a new second-order statistical approach to blind identification of single-input multiple-output (SIMO) autoregressive and moving average (ARMA) system models. The proposed approach exploits the dynamical autoregressive information of the model contained in the autocorrelation matrices of the system outputs but does not require the block Toeplitz structure of the channel convolution matrix used by classical subspace methods. For the multi-channel model with the same autoregressive (AR) polynomial, sufficient conditions and an efficient identification algorithm are given such that the multi-channel model can be uniquely identified up to a constant scaling factor. Furthermore, an extension of the result to blind identification of multi-channel models with different AR polynomials is presented. Simulation results are given to show the effectiveness of the proposed approach.