Em geral, a análise da estabilidade e a análise das respostas de poços de petróleo são realizadas de forma determinística em relação às propriedades mecânicas e hidráulicas do meio rochoso. No entanto, sabe- se que os meios rochosos e em particular rochas sedimentares, podem mostrar um considerável grau de heterogeneidades, em micro, meso e macro-escala. Essas heterogeneidades produzem variabilidade espacial nas propriedades mecânicas e hidráulicas dos meios rochosos. Essa variabilidade mostra em geral um caráter espacial pronunciado. O presente estudo propõe o desenvolvimento de procedimentos de análise numérica, utilizando elementos finitos, de processos fluido mecânicos acoplados, monofásicos e bifásicos, que levem em conta a variabilidade espacial de propriedades hidráulicas e mecânicas e a variabilidade das condições iniciais de tensões e poro pressões. Nesse estudo, empregam-se os procedimentos numéricos desenvolvidos em duas fases distintas. Na análise probabilística da estabilidade de poços e na análise probabilística das respostas dos poços durante a produção, considerando o acoplamento fluido mecânico com fluxo bifásico.