Um dos principais objetivos da engenharia de petróleo é desenvolver e aplicar técnicas capazes de aumentar a produtividade de poços de petróleo, incluindo a estimulação do poço por operações de fraturamento hidráulico. Estudos sobre a propagação de fraturas podem ser feitos analiticamente para algumas situações simplificadas envolvendo homogeneidade, isotropia e condições de contorno simples do meio geológico, ou pela aplicação de métodos numéricos, como o método dos elementos finitos, para casos mais complexos. A presente pesquisa apresenta análise numérica de fraturamento hidráulico utilizando o método estendido dos elementos finitos (XFEM), em conjunto com o modelo constitutivo de dano da Zona Coesiva (MZC). No método estendido dos elementos finitos a geometria da fratura se torna independente da malha, permitindo a propagação da fratura através do domínio, dispensando sucessivas gerações de malha necessárias na aplicação do método convencional dos elementos finitos. Os resultados numéricos obtidos foram comparados com soluções analíticas assintóticas no caso limite em que o regime da propagação é dominado pela rigidez da rocha, obtendo uma boa concordância. Adicionalmente, foram investigados os efeitos de diferentes parâmetros do fluido de injeção e as características de propagação da fratura quando a interface entre diferentes camadas geológicas é inclinada, mostrando dependência do ângulo de inclinação, das propriedades do material e das tensões in-situ.