En este documento, se estudia la solución del modelo de Darcy-Forchheimer en medios heterogéneos de alto contraste. Este problema se resuelve mediante un método de elementos finitos mixtos (MFEM) en una malla fina (la solución de referencia), donde la presión se aproxima mediante elementos constantes por partes; mientras tanto, la velocidad se discretiza mediante los elementos de menor orden de Raviart-Thomas. La solución en una malla gruesa se realiza utilizando el método de elementos finitos multiscale generalizado mixto (GMsFEM mixto). La ecuación no lineal puede resolverse mediante la conocida iteración de Picard. Se presentan varios experimentos numéricos en un dominio heterogéneo bidimensional para mostrar la buena aplicabilidad del método multiscale propuesto.