[pt] Com o objetivo de aumentar a capacidade dos poços petrolíferos, métodos convencionais de recuperação são utilizados, os quais consistem na injeção de água ou gás para a manutenção da pressão do reservatório. A produção do óleo ocorre através do deslocamento do mesmo no espaço poroso, onde a água, fluido deslocante, é injetada para ocupar gradualmente o espaço do óleo, fluido deslocado. Devido aos efeitos capilares e às heterogeneidades do meio poroso, uma parcela de óleo residual acaba ficando retida no reservatório, apresentando baixo fator de recuperação de óleo devido a elevada viscosidade do óleo em relação à viscosidade do fluido injetado e altas tensões interfaciais entre os fluidos. A adição de polímeros à água garante um aumento na sua viscosidade, melhorando a razão de mobilidade água/óleo no meio poroso. Uniformizando a frente de avanço e melhorando a eficiência de varrido devido à melhora no deslocamento do óleo. O presente trabalho analisa o comportamento viscoelástico do polímero, isolando o efeito viscoso e elástico em função das taxas de cisalhamento e extensão, implementado em um modelo de simulação de injeção de polímeros na escala de reservatórios. O efeito das propriedades reológicas da solução polimérica mostram nos resultados de produção uma frente de avanço mais estável e consequentemente uma melhora na taxa de recuperação de óleo quando avaliou-se o comportamento puramente cisalhante. Entretanto uma melhora na taxa de recuperação e na estabilidade da frente de avanço para o comportamento puramente extensional só pode ser observado quando o número de capilaridade foi aumentado consideravelmente. / [en] Aiming to increase the capacity of oil fields, conventional recovery methods are used. These methods consist in the injection of water or gas to maintain the reservoir pressure. The oil production typically takes place by displacing this oil in the porous media, where the displacing fluid (water) is injected to gradually occupy the space of the displaced fluid (oil). Since due to capillary effects and the heterogeneity of the porous media, a residual oil portion ends up trapped in the reservoir. These methods lead to low values of oil recovery factor, which occurs mainly by two factors: high viscosity of the reservoir s oil in relation to the viscosity of the injected fluid and high interfacial tension between the fluids. The addition of polymers to the water ensures an increase in the viscosity of the injected fluid, improving mobility ratio between water and oil in the porous media. Thus, standardizing forward swept and improving the swept efficiency due to improved oil displacement, which reduces the formation of preferential paths in the reservoir, usually called fingers. This paper analyzes the viscoelastic behavior of the polymer, by isolating the viscous and elastic effect in function of its extension and shear rates, implemented in a polymer injection simulation model in a reservoir scale. The effect of the rheological properties of the polymer solution show in the production results a more stable injection front and consequently an oil recovery rate improvement when evaluated as a purely shear behavior. However an improvement in the recovery rate and stability of the injection front for pure extensional behavior can only be observed when the capillary number is increased considerably.

[pt] Os efeitos catastróficos de terremotos que frequentemente acontecem na costa ocidental da América do Sul motivaram a realização desta pesquisa. O Equador, em particular, encontra-se localizado em uma zona de alta atividade sísmica devido a processos de subducção de placas tectônicas que originam sismos de grande magnitude, como o terremoto de magnitude Mw igual 7.8 acontecido em 16 de abril de 2016 na cidade costeira de Pedernales, que produziu grandes danos à infraestrutura e perdas de vidas humanas. A presente pesquisa tem como objetivo estudar o comportamento dinâmico de dois tipos de estruturas de terra: um talude de encosta composto de solo coluvial de altura aproximada de 70 m e a barragem de enrocamento de Mazar construída com altura de 166 m para fins de aproveitamento hidrelétrico. A barragem de Mazar, em função da sua altura, é classificada como barragem de grande tamanho pela International Commission on Large Dams (ICOLD) e qualificada pela Norma Equatoriana da Construção (NEC, 2014) como estrutura de importância essencial, com obrigatoriedade de verificação do seu desempenho sísmico para dois terremotos de projeto: o terremoto operacional base (OBE) e o máximo terremoto provável (MCE) com períodos de retorno de 475 anos e 2500 anos, respectivamente. As análises efetuadas neste trabalho incluíram a geração de sismos de projeto, investigação da estabilidade estática e pseudoestática de taludes, determinação de deslocamentos permanentes mediante métodos analíticos simplificados e simulações numéricas. Diferenças significativas nos deslocamentos permanentes máximos foram observadas empregando as abordagens analítica e numérica. / [en] The catastrophic effects of earthquakes that frequently occur on the West coast of South America motivated this research. Ecuador, in particular, is located in a zone of high seismic activity due to processes of tectonic plate subduction that give rise to earthquakes of great magnitude, such as the seismic motion of magnitude Mw equal 7.8 that struck the coastal city of Pedernales on April 16 2016, producing great damages to the infrastructure and losses of human lives. The present research aims to study the dynamic behavior of two types of earth structures: a slope composed of colluvial soil of approximately 70 m height and the Mazar rockfill dam constructed for hydroelectric generation. The Mazar dam, according to its height (166 m), is classified as a large dam by the International Commission on Large Dams (ICOLD) and qualified by the Ecuadorian Construction Standard (NEC, 2014) as a structure of essential importance, with mandatory analysis of seismic performance under two seismic design ground motions: the operating basis earthquake (OBE) and the maximum credible earthquake (MCE) with return periods of 475 years and 2500 years, respectively. The analyses carried out in this work included the generation of the OBE and MCE earthquakes, investigation of static and pseudo static stability of slopes, determination of permanent displacements using simplified analytical methods and numerical simulations. Significant differences in the maximum permanent displacements were observed using analytical and numerical approaches.

[pt] O entendimento dos processos que levam à instabilidade de taludes rochosos é de grande importância em análises de risco. Normalmente, desplacamentos e quedas de blocos de rocha ocorrem em épocas sem chuva, levando a crer que existem outros fatores deflagradores que devem ser analisados. Estudos relativamente recentes tem levantado a importância das tensões térmicas na propagação de fissuras e noacúmulo de tensões de tração. De forma a avaliar o impacto da ciclagem térmica sobre a rocha, foram realizadas análises termomecânicas utilizando o software de elementos finitos Abaqus 6.14. Nas simulações efetuadas foram utilizadas geometrias simplificadas, representativas de placas rochosas reais de áreas costeiras do Brasil. O maciço rochoso foi considerado isotrópico e homogêneo e os parâmetros térmicos utilizados (condutividade, calor específico e coeficiente de expansão térmica) foram adotados de acordo com valores obtidos na literatura para rochas de origem e constituição similar às rochas presentes na costa da região sudeste do Brasil (Gnaisses e Granitos). As condições de contorno de temperatura foram obtidas através de campanhas experimentais extraídas de estudo anterior. Admitiu-se ainda, apenas a ocorrência do modo I de propagação de fraturas, sendo a análise da propagação restrita a comparação entre K1c e os valores K1 desenvolvidos. Os resultados das simulações numéricas mostraram que a ocorrência de baixas temperaturas e de altos gradientes térmicos entre as faces externa e interna da placa ocasionam os maiores valores de K1, cuja magnitude é comparável aos valores de K1c típicos para rochas encontradas na região sudeste do Brasil. / [en] The understanding of the processes that lead to instability of rock slopes is of great importance in risk analysis. Normally, detachments and rockfalls occur in times without rain, leading to belief that there are other triggering factors that need to be evaluated. Relatively recent studies have raised the importance of thermal stresses in the propagation of cracks and in the accumulation of tensile stresses. In order to evaluate the impact of thermal cycling on the rock, thermomechanical analyzes were performed using Abaqus 6.14 finite element software. In the simulations, simplified geometries were used, representative of real rock plates of coastal areas of Brazil. The rock mass was considered isotropic and homogeneous and the thermal parameters (conductivity, specific heat and coefficient of thermal expansion) were adopted according to values obtained in the literature for rocks of origin and constitution similar to the rocks present in the coast of the southeastern region of Brazil (Gneisses and Granites). The temperature contour conditions were obtained through experimental campaigns extracted from a previous study. Only the occurrence of mode I of fracture propagation was allowed, with propagation analysis being restricted to the comparison between K1c and the K1 values developed. The results of the numerical simulations showed that the occurrence of low temperatures and high thermal gradients between the outer and inner faces of the plate cause the highest values of K1, whose magnitude is comparable to the values of K1c typical for rocks found in the southeastern region of Brazil.

[pt] Reservatórios carbonáticos contêm entre 50 e 60 por cento de petróleo e gás do mundo. No entanto, estas rochas têm apresentado problemas operacionais consideráveis durante as etapas de perfuração e produção. Tais problemas são originados pelas dificuldades na caracterização adequada destas rochas em função da complexa distribuição espacial de sua micro e macro estrutura. Enquanto ambas as escalas possuem importância no entendimento do comportamento de carbonatos, o presente trabalho trata da análise da microestrutura. Nesse âmbito, a proposta deste trabalho integra o processamento e análises de imagens digitais adquiridas mediante microtomografia de raios-X, elaboração de malhas de elementos finitos e simulação numérica de forma a prever propriedades elásticas, com o objetivo de correlacionar a microestrutura e o módulo de Young. Foram utilizadas imagens digitais de amostras de afloramento (travertinos) considerados como possíveis análogas a rochas carbonáticas de reservatório. A metodologia implementada permitiu obter sub-amostras que incluem uma gama de microestruturas e porosidades numa única amostra para assim simular numericamente o módulo de Young. Os resultados dos ensaios virtuais foram comparados com os resultados de ensaios reais executados em amostras da mesma rocha e dados encontrados na literatura. A aplicação da metodologia e os resultados obtidos indicam o potencial e as limitações atuais desta técnica. Conclui-se neste estudo que o valor do módulo de Young das amostras simuladas numericamente é afetado pela distribuição espacial, conectividade dos poros e pela microporosidade da rocha. Os valores encontrados no ensaio virtual tendem a serem superiores aos obtidos em ensaios reais devido à dificuldade de inserir, no modelo de micro estrutura, os contatos entre os grãos assim como regiões mais compressíveis não determinadas na aquisição das imagens. / [en] Carbonate reservoirs contain between 50 and 60 per cent of oil and gas in the world. However, these rocks have presented considerable operational problems during drilling and production steps. Such problems are caused by the difficulties in proper characterization of these rocks due the complex spatial distribution of their micro and macrostructure. Although both scales have importance in understanding the behavior of carbonates, the present work deals with the analysis of the microstructure. In this context, the proposal of this work integrates the processing and analysis of digital images acquired by x-ray microtomography, finite element mesh generation and numerical simulations to predict elastic properties in order to correlate the microstructure and the Young s modulus. We used digital images of samples of outcrop (travertine) considered as a possible analogous to carbonate rocks of reservoir. The methodology implemented allowed get sub-samples that include a range of microstructures and porosities in a single sample to determine numerically the Young s modulus. The virtual test results were compared with laboratory test results performed on samples from the same rock and data found in the literature. The application of the methodology and the results obtained indicate the potential and current limitations of this technique. This study concluded that the value of Young s modulus of numerically simulated samples is affected by spatial distribution, pore connectivity and microporosity of these rocks. The values found in the virtual test tend to be higher than those obtained in laboratory tests due to the difficulty of inserting, in the model of micro structure, contacts between the grains as well as more compressible regions not certain on the acquisition of the images.

[pt] Berger e Podolak apresentaram uma interpretação geométrica do resultado seminal de Ambrosetti e Prodi sobre o comportamento das soluções de certas equações diferenciais parciais elípticas semi-lineares. Consideram-se extensões deste ponto de vista, a partir das quais se desenvolve um algoritmo numérico para resolver as equações. / [en] Berger and Podolak obtained a geometric interpretation of the seminal result of Ambrosetti and Prodi regarding the behavior of solutions of certain semilinear elliptic partial differential equations. We consider extensions of such interpretation to develop a stable numerical algorithm that solves the equations.

[pt] 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. / [en] One of the main objectives of petroleum engineering is to develop and apply techniques capable of increasing the productivity of oil wells, including the stimulation of well by hydraulic fracturing operations. Studies on the propagation of fractures can be done analytically for some simplified situations involving homogeneity, isotropy and simple boundary conditions of the geological medium, or by the application of numerical methods, such as the finite element method, for more complex cases. The present research presents a numerical analysis of hydraulic fracturing using the extended finite element method (XFEM), in conjunction with the damage constitutive model of Cohesive Zone (MZC). In the extended finite element method the fracture geometry becomes independent of the mesh, allowing the propagation of the fracture through the domain without successive mesh generations as necessary in the conventional finite element method. The computed numerical results were compared with asymptotic analytical solutions in the limit case in which the propagation regime is dominated by the rigidity of the rock with good compatibility. In addition, this study investigates the effects of different parameters of the injection fluid and the fracture propagation characteristics when the interface between different geological layers is inclined, shows dependency between the angle of inclination with the properties of the material and the in-situ stresses.

[pt] O estudo de iniciação e propagação de fraturas em rochas devido à energia gerada por explosões é uma tarefa desafiadora em mecânica computacional dada a natureza multifísica e multi-escala do fenômeno. Uma das técnicas utilizadas para simulação deste processo pelo método dos elementos finitos consiste em acompanhar a evolução de fraturas no tempo, com atualizações frequentes da malha de elementos, o que torna as análises demoradas e complicadas, com perda de precisão numérica no processo de atualização dos valores calculados em pontos da malha antiga para os pontos correspondentes da malha nova. O método estendido dos elementos finitos (XFEM) permite a incorporação de enriquecimentos locais, i.e. de um conjunto de funções de interpolação enriquecidas que fornecem valores das variáveis de interesse (deslocamentos, tensões) com maior precisão e eficiência computacional. Além disso, é importante ressaltar, que a presença da fratura, e sua propagação no tempo através da rocha, não é geometricamente modelada e a malha de elementos não precisa ser constantemente atualizada. Quatro diferentes abordagens são examinadas para simular o processo de fraturamento na rocha, com a comparação entre os respectivos resultados: o método XFEM, os elementos de interface coesivas, os elementos finitos com singularidade e a técnica de eliminação de elementos que remove elementos da malha, simulando o fraturamento, quando os mesmos atingem a ruptura de acordo com algum critério. Nesta pesquisa, o método XFEM é aplicado para investigar o desmonte de rocha com base no método dos nós fantasma onde as descontinuidades nos campos de deslocamentos são introduzidas através de novos graus de liberdade em elementos sobrepostos. O maciço rochoso considerado é um granito admitido isotrópico no meio homogêneo ou heterogêneo que tem comportamento elástico linear até o início da quebra, onde a propagação de fraturas utiliza o modelo de zona coesiva. Alguns exemplos numéricos são apresentados aspectos relacionados com o fraturamento de um maciço rochoso sujeito a explosão, a fim de discutir as vantagens e limitações. Além disso, os resultados numéricos são comparados com os obtidos por outros autores utilizando diferentes abordagens numéricas. / [en] The study of propagation of fractures in rocks due to the energy generated by explosions is a challenging task in computational mechanics given the multiphysics and multiscale nature of the phenomenon. One of the most widely used methods for simulation of this process is the finite element method, which follows the time evolution of fractures, with frequent updates of mesh elements to represent the new geometry of the newly fractured material. This approach, besides being computationally time consuming and difficult for the necessity of constant rebuilding meshes, also results in the loss of numerical accuracy when the variables of interest are mapped and interpolated from the old mesh to the Gauss points and nodal points the new mesh. The Extended Finite Element Method (XFEM) local enrichment functions to be easily incorporated into a finite element approximation. The presence of fracture is ensured by the special enriched functions in conjunction with additional degrees of freedom with greater accuracy and computational efficiency. Furthermore, it is important to note that does not require the mesh to match the geometry of the fracture. It is a very attractive and effective way to simulate initiation and propagation of a crack along an arbitrary, solution-dependent path without the requirement of remeshing. Four different approaches are examined to simulate the rock fracturing process, with comparison between respective results: the XFEM, the interelement crack method, the conventional finite element method (FEM) using a remeshing technique and based on the linear fracture mechanics and the element deletion method with Rankine failure-type material model to simulate discrete rock fracture. In this research, XFEM is applied to investigate rock blasting based on the phantom node method where discontinuities in the displacement fields are introduced through new degrees of freedom in overlapping elements. The rock mass considered is a sound granite admitted as an isotropic, homogeneous or heterogeneous medium that remains linear elastic right up the moment of breakage, and then the propagation of cracks using the cohesive zone model. Several numerical examples are presented aspects related to the fracturing of a rock mass under the effect of blast-induced dynamic pressure pulse, in order to discuss the advantages and limitations of each of the aforementioned approaches. Furthermore, the numerical results are compared with those obtained by other authors using different numerical approaches.

[pt] Este trabalho tem como objetivo apresentar um método para previsão de recalque de grupos de estacas, solicitadas axialmente em maciços de solo estratificados, através do emprego de técnicas numéricas. O método de análise utilizado é baseado em uma versão simplificada do método dos elementos de contorno, sugerida por Poulos e Davis (1968) para o caso de solo homogêneo, isotrópico e linear elástico, desenvolvida a partir das soluções de Mindlin (1936). Neste caso, a modificação proposta constou da substituição de meio estratificado por meios homogêneos equivalentes, através da utilização de procedimentos sugeridos por Poulos (1979) e Yamashita et al (1987). Os valores de recalque obtidos com os programas computacionais desenvolvidos nesta pesquisa foram comparados com resultados correspondentes obtidos pelo métodos dos elementos finitos ou publicados na literatura por diferentes autores e através de diversas metodologias. Estas comparações permitem concluir que a utilização da proposição aproximada sugerida neste trabalho fornece, sob certas condições, valores bastante satisfatórios do recalque de grupos de estacas em maciços de solo estratificados. As vantagens de utilização deste procedimento estão diretamente associadas a uma significativa simplificação na preparação dos dados para as necessárias descrições geométricas e geotécnicas do problema bem como a redução do tempo de processamento. Fatores de interação para grupos de estacas foram produzidos para alguns casos específicos de grupos em maciços de solo formados por camadas duplas. / [en] This thesis presents a numerical analysis of pile groups deflections, under axial load in layered soils. The method uses a simplified version of the boundary element technique, as proposed by Poulos and Davis (1968) for a homogeneous, isotropic and linear elastic half-space, based on Mindlin s equation. In this case, the major step consisted in replacing the multi-layered system for an equivalent homogeneous half-space using the well known procedures suggested by Poulos (1979)and Yamashita et al (1987). Results of the deflection of pile groups obtained with the proposed method were compared to results obtained using finite element analysis and to results published in the technical literature obtained by several authors. It has been shown that the proposed method provides results which compare favorably with existing solution and basically presents the advantage of the simplification in handling the soil and geometrical datas necessary for the problem description. Interaction factors for groups of pile were produced for some specific cases consisting of a double layer system. / [es] Este trabajo tiene como objetivo presentar un método para previsión de recalque de grupos de estacas, solicitadas axialmente en macizos de suelo extratificados, a través del empleo de técnicas numéricas. EL método de análisis utilizado tiene como base en una versión simplificada del método de los elementos de contorno, sugerida por Poulos y Davis (1968) para el caso de suelo homogéneo, isotrópico y lineal elástico, desarrollada a partir de las soluciones de Mindlin (1936). En este caso, la modificación propuesta incluye la substituición de medio extratificado por medios homogéneos equivalentes, a través de la utilización de procedimentos sugeridos por Poulos (1979) y Yamashita et al (1987). Los valores de recalque obtenidos con los programas computacionales desarrollados en esta investigación fueron comparados con resultados correspondentes obtenidos por el métodos de los elementos finitos o publicados en la literatura por diferentes autores y a través de diversas metodologías. Estas comparaciones permiten concluir que la utilización de la proposición aproximada sugerida en este trabajo brinda, bajo certas condiciones, valores bastante satisfactorios del recalque de grupos de estacas en macizos de suelo extratificados. Las ventajas de utilización de este procedimiento están directamente asociadas la una significativa simplificación en la preparación de los datos para las necesarias descripciones geométricas y geotécnicas del problema así como la reducción del tiempo de procesamiento. Fueron produzidos factores de interación para grupos de estacas para algunos casos específicos de grupos en macizos de suelo formados por estratos duplas.

[pt] O presente trabalho apresenta a aplicação das metodologias de análise de estabilidade para taludes rochosos de grande altura comumente aplicados dentro do mundo da mineração. Mostra os principais fatores que influem dentro de uma análise de estabilidade como são a determinação das propriedades geomecânicas do maciço rochoso (rocha intacta e descontinuidades) e a caracterização geológica estrutural da região estudada. Esses parâmetros geralmente são obtidos a partir de ensaios de laboratório e de campo, além da informação que fornecem os mapeamentos geológicos. Outra forma de determinar estas propriedades de resistência é através de retro-análises em regiões onde tenha ocorrido algum tipo de colapso ou pela recopilação de dados que foram utilizados em análises anteriores e possam ser extrapolados. Também é apresentado um passo a passo das distintas etapas de um estudo de estabilidade e a obtenção de dados que finalizam na aplicação de softwares especializados na área de geotécnia e geomecânica, os quais permitirão determinar e pré-visualizar os possíveis problemas de instabilidade dentro dos taludes de um pit mineiro. Serão apresentados dois tipos de análises de estabilidade, o primeiro baseado na teoria do método de equilíbrio limite, o qual procura a possível superfície crítica de deslizamento dentro de um talude e que, condiciona a estabilidade da parede aos valores do fator de segurança. A segundo análise está relacionada com um modelo numérico, o qual aplica o método de elementos distintos ou de blocos, que permite a aplicação de um modelo constitutivo que descreva o comportamento dos materiais e das descontinuidades dentro do maciço rochoso, fornecendo como resultados a pré-visualização de variação de deslocamentos e vetores de velocidade para cada etapa de escavação. / [en] This thesis presents the application of procedures for the stability analysis of high slopes commonly applied in the mining world. It presents the main parameters that influence the stability analysis: determination of the rock mass geomechanical properties (intact rock and discontinuities) and the geological structures characterization of the studied region. These parameters generally are the result of laboratory and field tests, in addition, the information about the geological mapping. Another procedure for determining the properties strength may be through a back analysis of a collapsed zone or extrapolation of a data compilation from previous analyses. Apart from strength parameter determination and concept applying of rock mechanics, this research describes step by step the different stages of a stability analysis and its obtained data, this work then finishes with a geotechnical and geomechanical software application, which will determines and previews slope stability problems for a mining pit. This thesis includes two stability analysis procedures: the first method applies the limit equilibrium theory, which looks for the critical failure surface that depends on depends on a factor of safety. The second one applies a numerical model that uses distinct element method, which through a constitutive model, describes how the materials and discontinuities behave in a rock mass, obtaining as a result, the displacement and velocity vectors for each excavation stage.

[pt] Reservatórios carbonáticos contêm entre 50 e 60 porcento de petróleo e gás do mundo. No entanto, estas rochas têm apresentado problemas operacionais consideráveis durante as etapas de perfuração e produção. Tais problemas são originados pelas dificuldades na caracterização adequada destas rochas em função da complexa distribuição espacial de suas micro e macro estrutura. Enquanto ambas as escalas possuem importância no entendimento do comportamento de carbonatos, o presente trabalho trata da análise da microestrutura. Nesse âmbito, a proposta deste trabalho integra o processamento e análises de imagens digitais, elaboração de malhas de elementos finitos e simulação numérica de forma a prever propriedades elásticas, com o objetivo de correlacionar a microestrutura e o módulo de Young. Foram utilizadas imagens digitais de amostras de afloramento (travertinos) considerados como possíveis análogas a rochas carbonáticas de reservatório. A metodologia implementada permitiu obter subamostras que incluem uma gama de microestruturas e porosidades numa única amostra para assim simular numericamente o módulo de Young. Conclui-se neste estudo que os poros maiores (vuggys) são os principais responsáveis pelo comportamento do material quando o mesmo é carregado, e também a influência da estrutura porosa nos resultados, foi possível observar que a posição dos poros afeta o módulo de elasticidade do corpo de prova. / [en] Carbonate reservoirs contain between 50 and 60 percent of oil and gas reserves in the world. However, these rocks have presented considerable operating problems during the phases of drilling and production. Such problems are originated by the difficulties in the proper characterization of these rocks due to the complex spatial distribution of its micro and macro structure. While both scales have importance on the carbonate behavior understanding, this paper deals with the analysis of the microstructure. In this context, the aim of this work is part of the processing and digital image analysis, mesh preparation and finite element numerical simulation to predict elastic properties, in order to correlate the microstructure and the Young s modulus. Using digital images of outcrop samples (travertine) considered as possible analogous to carbonate rocks reservoir. The methodology implemented afforded sub-samples that include a range of porosities and microstructures for a single sample well numerically simulate the Young s modulus. It is concluded in this study that the larger pores (vuggys) are mainly responsible for the behavior of the material when it is loaded, and also the influence of the porous structure in the results, it was observed that the position of pores affects the modulus of elasticity of sample.