Este trabajo presenta los resultados de simulacion de flujo y transporte de agua subterranea en medio poroso heterogeneo. Esta simulacion se realizo a partir de una aproximacion de un medio fracturado por medio de un medio poroso equivalente utilizando el software de geoestadistica T-PROGS. Este software utiliza la probabilidad de transicion y el modelo de cadena de Markov como aproximacion estocastica. La representacion del medio poroso equivalente estocastico se definio segun las facies de roca, dependiendo de su conductividad hidraulica y densidad de fractura. Se definieron cuatro facies de roca: roca de roca escasamente fracturada (SFB), roca de roca escasamente a moderadamente fracturada (SMFB), roca madre moderadamente a muy fracturada (MHFB) y roca madre altamente fracturada (HFB). La ecuacion de flujo de agua subterranea se resolvio numericamente a traves de diferencias finitas, en las que se considero el tensor completo de conductividades hidraulicas y la simulacion de transporte se realizo mediante la implementacion del metodo Rastreo Aleatorio de Particulas (RAP), en las que la interpolacion de velocidad bilineal y la dispersividad tensorial son consideradas. Estos resultados se compararon con una simulacion de flujo y transporte implementada con el esquema de diferencias finitas (DF). Los escenarios de validacion presentan que las concentraciones obtenidas por medio del RAP pueden diferir del esquema DF dependiendo de como se distribuyen las particulas en el tiempo de ejecucion inicial. Los escenarios de simulacion del comportamiento de los contaminantes presentan diferencias cuando se simula en un medio poroso heterogeneo. En este caso, se evidencia un transporte mas rapido del contaminante mostrando direcciones preferenciales definidas que coinciden con zonas de conductividades hidraulicas mas elevadas.