With meshfree and fully Lagrangian features of particle methods, smoothed particle hydrodynamics (SPH) is suitable to achieve high-accurate simulations of multiphase flows with large interfacial deformations, discontinuities, and multi-physics. Multiphase flow simulations with SPH methods have been reported abundantly and the specific implementations are much different. In this review, the basic SPH method and issues about fluid pressure, surface tension and solid boundary are discussed. And various implementations of SPH for multiphase flow simulation are mainly summarized as: (1) Lagrangian solver for the two-fluid model (TFM): The two phases are discreterized into two independent groups of SPH particles and coupled by the explicit interphase interaction; (2) multiphase SPH: The multiphase SPH method is considered as the natural extension of SPH method on multiphase flow simulation, and the interphase interaction is implicitly described by SPH parameters; (3) coupling of SPH and other discrete methods: The two phases with large differences each adopt different discrete methods to give play to the advantages of different Lagrangian methods; and (4) coupling of SPH and grid-based methods: The grid method handles the simple main-flow to obtain the balance between accuracy and efficiency. Also, some issues associated with SPH simulations of multiphase flows, such as the physicalization of simulation parameters and the improvement of accuracy and efficiency, are suggested as requiring attention.