Objetivo: El modelado de flujos y contaminantes enfrenta grandes desafíos que van desde la definición de geometrías complejas hasta la convergencia del algoritmo y la determinación de los parámetros del modelo. El método presentado aquí aborda este problema a través de una aplicación conjunta de los métodos de Lattice-Boltzmann y diferencias finitas. La solución numérica de las ecuaciones de Navier-Stokes y advección-difusión permite simular los niveles de dispersión de la contaminación espaciotemporal en una fuente hídrica. Metodología: El método resolvió las ecuaciones de Navier-Stokes a través de la propuesta de Lattice Boltzmann, que proporciona el campo vectorial de velocidades del fluido. La solución de la ecuación de difusión-advección por diferencia finita requiere este campo de velocidad conocido. La solución codificada en C++, primero digitaliza el mapa, etiqueta cada píxel del mapa como fluido o sólido. El código usa la imagen digital para configurar la geometría de la zona de flujo, luego simula la dinámica del fluido y la concentración del contaminante. La comparación entre las soluciones analíticas de modelos ideales advectivo-difusivo, y las obtenidas a través de este método, facilita la medición de la confiabilidad del método; aplicado a una sección del río Sogamoso (Colombia), el algoritmo arrojó la velocidad, la presión y la concentración del agente. Resultados: Se encontró que el término advectivo predomina en gran parte de cada uno de los tramos del río Sogamoso. Sin embargo, hay una sección en el tramo superior, en donde predomina la componente difusiva. Esto ocurre porque allí la velocidad es muy baja. Conclusión: Los resultados señalan la confiabilidad, estabilidad y robustez del método para simular fenómenos de advección-difusión en cuerpos con geometrías irregulares.