As zonas de falha são compostas por um núcleo, onde a maior parte da deformação é acomodada, e uma zona de dano, com deformação menos intensa. A zona de dano pode atuar como caminho de fluxo preferencial devido à presença de fraturas, ou como barreira devido às bandas de deformação. Portanto, sua caracterização é essencial para a adoção de estratégias de produção adequadas em campos de petróleo. Os métodos geofísicos geralmente utilizados, porém, dificilmente permitem a identificação das zonas de dano devido à baixa resolução sísmica. Como alternativa, empregam-se observações em afloramentos superficiais. Contudo, há uma grande dispersão de dados, que pode estar relacionada a uma variedade de fatores, como as propriedades da rocha protólita e os mecanismos de deformação atuantes. Neste sentido, este trabalho apresenta duas metodologias baseadas no método dos elementos finitos (MEF) para analisar a formação e evolução das zonas de dano em escala de reservatório. Na primeira abordagem, a zona de falha é totalmente representada através de um meio contínuo, enquanto que na segunda, a falha é representada como um plano por meio de uma descontinuidade. Em ambas aproximações, a zona de dano é estabelecida através das regiões plastificadas. Os resultados numéricos obtidos se aproximaram das observações de campo e possibilitaram a identificação das vantagens e limitações das duas abordagens baseadas no MEF. Por fim, os resultados também permitiram identificar os principais parâmetros geomecânicos que influenciam o desenvolvimento das zonas de dano, bem como os diferentes mecanismos de deformação que ocorrem ao longo da zona de dano.