In this paper, we present a nonlinear mathematical model, describing the spread of high-risk alcohol consumption behavior among college students in Colombia. We proved the existence and stability of the alcohol-free and drinking state equilibrium by means of Lyapunov function and LaSalle’s invariance principle. Also, we apply optimal control to study the impact of a preventive measure on the spread of drinking behavior among college students. Finally, we use numerical simulations and available data provided by the United Nations Office on Drugs and Crime (UNODC) and the Colombian Ministry of Justice to validate the obtained mathematical model.