Desarrollamos un modelo matemático para las características de transmisión de la doble variante SARAS-CoV-2 con vacunación de la variante 1 para abordar este nuevo aspecto de la enfermedad. El modelo es examinado teóricamente y se derivan requisitos adecuados para la estabilidad de sus puntos de equilibrio. El modelo incluye los equilibrios endémicos de la variante 1 y la variante 2, además de los equilibrios endémicos y libres de enfermedad. Se utilizan varios enfoques para la estabilidad global y local del modelo. Para ambas cepas, determinamos los números reproductivos básicos y . Para investigar la ocurrencia de las capas (olas), expandimos el modelo para incluir un análisis basado en la segunda derivada. Luego, el modelo se expande a su forma estocástica y se calculan resultados numéricos. Para fines numéricos, utilizamos el método de diferencia finita no estándar. También se registra algún análisis de error.