Slug flow is a two-phase flow pattern that requires large effort in its characterization and modeling, due to special characteristics of the phase’s spatial distribution, which causes flow intermittency. This flow occurs in a wide range of gas and liquid flow rates in pipes of medium and small diameters, with periodic variation of density, void fraction and pressures in pipe-cross-section. This work presents a numerical study of the slug regime through horizontal pipes using the Two-Fluid Model. A detailed statistical analysis of the flow was carried out with characterization of main slug parameters, such as slug length, front and tail velocities and slug frequency along the pipeline. The numerical results were analyzed through distribution histograms as well as probability density function (PDF) of the hydrodynamic variables, showing the stochastic characteristic of slug flow pattern. Further, comparisons with experimental data from the literature were performed, showing very good agreement.