To prevent gas hydrate formation and the plugging of the pipe, monoethylene glycol, MEG, is injected into the wellhead. To develop the thermodynamic model for the optimization of MEG injection, liquid liquid equilibrium, LLE, studies are necessary. In the present work, original data are presented of binary, ternary, and quaternary systems involving MEG, water, hexane, octane, toluene, and p-xylene under atmospheric pressure and in the temperature range between 278.15 and 333.15 K. The NRTL thermodynamic model was used for the solubility prediction of binary, ternary, and quaternary systems. The adjustable parameters of the model were determined through regression of binary mixtures data from literature. Experimental and predicted solubilities are in a quite good agreement, in view of the low solubility data involved in the different systems.