El problema de programación no lineal multiobjetivo de nivel bi-áspero de múltiples opciones (MR-BLMNPP) ha aumentado notablemente en diversas aplicaciones reales. En el modelo actual, las funciones objetivo tienen un parámetro de múltiples opciones, y las restricciones se representan como un conjunto áspero. En la primera fase, se utilizan las diferencias divididas de Newton (NDDs) para formular un polinomio de las funciones objetivo. Luego, basándose en la teoría de conjuntos ásperos, el modelo se convierte en un Modelo de Aproximación Superior (UAM) y un Modelo de Aproximación Inferior (LAM). En la segunda fase, se presentan dos modelos basados en la Técnica de Preferencias de Orden por Similitud a la Solución Ideal (TOPSIS) para resolver el MR-BLMNPP. Se aplica un modelo basado en TOPSIS de max-min difuso y programación de metas difusas (FGP) para abordar el conflicto entre las funciones de distancia biobjetivo modificadas. También se presenta un algoritmo para resolver MR-BLNPP. La aplicabilidad y eficiencia de los dos modelos basados en TOPSIS sugeridos en este estudio se presentan a través de un algoritmo y una ilustración numérica. Finalmente, el estudio presenta un modelo de planificación de producción de nivel bi-áspero (BL-PPM) como una aplicación ilustrativa.