The accurate prediction of high pressure phase equilibria is crucial for the development and the design of chemical engineering processes. Among them the modeling of complex systems, such as petroleum fluids with water, has become more and more important with the exploitation of reservoirs in extreme conditions. The aim of this work is to explore the capability of the NRTL-PR model to predict the solubility of solid polycyclic aromatic hydrocarbons in water. For this purpose, we first validate our methodology for fluid phase equilibria predictions of aromatic hydrocarbons and gas (CO2, C2H6) mixtures. Finally, we consider the prediction of the solid solubility of PAH in water, by fitting group parameters either only on SLE data or on both LLE and SLE data of aromatic hydrocarbon-water binary systems.