[pt] O presente trabalho foca em uma investigação numérica do escoamento através de meios porosos com macroporos, como os carbonatos. A presença de vugs e fraturas tem um forte efeito na característica do escoamento. O fluxo através da matriz porosa é geralmente descrito pela equação de Darcy e o fluxo através dos macroporos pela equação de Stokes. O acoplamento entre essas duas abordagens distintas traz grande complexidade para a modelagem de tais fluxos. Neste trabalho, utilizamos a formulação de Brinkman que é capaz de descrever ambos os fluxos através da matriz porosa e macroporos com uma única equação diferencial. Resolvemos o sistema de equações diferenciais usando o modelo de elementos finitos e implementamos o código com auxílio da plataforma FEniCS. Primeiro, resolvemos o fluxo 1-D através de placas paralelas com uma das paredes sendo um material poroso. O objetivo foi comparar as previsões obtidas pela formulação de Brinkman com aquela obtida usando a condição de contorno de Beavers-Joseph. Então, resolvemos um fluxo 2-D através de um meio poroso com macroporos. A geometria da estrutura porosa foi obtida a partir de cortes 2D de imagens tomográficas de carbonatos. O objetivo desta análise foi avaliar uma permeabilidade equivalente em função da área e estrutura dos macroporos. / [en] The present work focuses on a numerical investigation of the flow through porous media with macropores, such as carbonates. The presence of vugs and fractures have a strong effect on the flow characteristics. The flow through the porous matrix is usually described by the Darcy equation and the flow through the macropores by Stokes equation. The coupling between these two distinct approaches brings great complexity to the modeling of such flows. In this work, we use Brinkman formulation that is able to describe the flow, both in the porous matrix and macropores with a single differential equation. We solved the set of differential equations using the finite element model and implemented the code in the FEniCS platform. We first solved the 1-D flow through parallel plates with one of the walls being a porous material. The goal was to compare the predictions obtained with the Brinkman formulation to that obtained by using the Beavers-Joseph boundary condition. Then, we solved a 2-D flow through a porous medium with macropores. The geometry of the pore structure was obtained from 2D slices of tomographic images of carbonates. The goal of this analysis was to evaluate an equivalent permeability as a function of macropores area and structures.