Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / As técnicas de injeção de traçadores têm sido amplamente utilizadas na investigação de escoamentos em meios porosos, principalmente em problemas envolvendo a simulação numérica de escoamentos miscíveis em reservatórios de petróleo e o transporte de contaminantes em aquíferos. Reservatórios subterrâneos são em geral heterogêneos e podem apresentar variações significativas das suas propriedades em várias escalas de comprimento. Estas variações espaciais são incorporadas às equações que governam o escoamento no interior do meio poroso por meio de campos aleatórios. Estes campos podem prover uma descrição das heterogeneidades da formação subterrânea nos casos onde o conhecimento geológico não fornece o detalhamento necessário para a predição determinística do escoamento através do meio poroso. Nesta tese é empregado um modelo lognormal para o campo de permeabilidades a fim de reproduzir-se a distribuição de permeabilidades do meio real, e a geração numérica destes campos aleatórios é feita pelo método da Soma Sucessiva de Campos Gaussianos Independentes (SSCGI). O objetivo principal deste trabalho é o estudo da quantificação de incertezas para o problema inverso do transporte de um traçador em um meio poroso heterogêneo empregando uma abordagem Bayesiana para a atualização dos campos de permeabilidades, baseada na medição dos valores da concentração espacial do traçador em tempos específicos. Um método do tipo Markov Chain Monte Carlo a dois estágios é utilizado na amostragem da distribuição de probabilidade a posteriori e a cadeia de Markov é construída a partir da reconstrução aleatória dos campos de permeabilidades. Na resolução do problema de pressão-velocidade que governa o escoamento empregase um método do tipo Elementos Finitos Mistos adequado para o cálculo acurado dos fluxos em campos de permeabilidades heterogêneos e uma abordagem Lagrangiana, o método Forward Integral Tracking (FIT), é utilizada na simulação numérica do problema do transporte do traçador. Resultados numéricos são obtidos e apresentados para um conjunto de realizações amostrais dos campos de permeabilidades. / Tracer injection techniques have been widely used to investigate flows in heterogeneous porous media, especially in problems related to numerical simulation of miscible flows in oil reservoirs and to contaminant transport in aquifers. Oil reservoirs are generally heterogeneous and may possess spatially significant variations in their properties on several length scales. These spatial variations are incorporated into the governing equations for flow problems in porous media on the basis of random fields. Random fields provide a natural description of rock heterogeneities in the typical case in which the geological knowledge of rock is much less detailed than is necessary to predict flow properties through it deterministically. In this thesis we adopt a scalar log-normal permeability field k(x) to reproduce the statistical distribution of the permeability values of a real medium, and the numerical generation of these random fields is based on a Successive Sum of Independent Gaussian Fields defined on multiple length scales. The aim of this work is to study the uncertainty quantification in inverse problems for tracer transport in heterogeneous porous media in a Bayesian framework and propose the permeability update based on observed measurements of spatially sparse tracer concentration at certain times. A two-stage Markov chain Monte Carlo (MCMC) method is used to sample posterior probability distribution with hierarchical priors and the Markov chain is constructed from random reconstruction of the permeability fields. To solve the Darcys law we use a mixed finite elements method which are suitable to compute accurately the relevant fluxes in heterogeneous permeability fields and a Lagrangian strategy, the Forward Integral Tracking (FIT) method, for the numerical simulation of tracer transport problem. Numerical results are presented for a set of sampled realizations of the permeability fields.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-03-02T18:09:12Z No. of bitstreams: 1 deniseschimitzdecarvalhotristao.pdf: 734334 bytes, checksum: 9fda9bda660d5bfec3204e328fe66d1c (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-03-06T19:58:58Z (GMT) No. of bitstreams: 1 deniseschimitzdecarvalhotristao.pdf: 734334 bytes, checksum: 9fda9bda660d5bfec3204e328fe66d1c (MD5) / Made available in DSpace on 2017-03-06T19:58:58Z (GMT). No. of bitstreams: 1 deniseschimitzdecarvalhotristao.pdf: 734334 bytes, checksum: 9fda9bda660d5bfec3204e328fe66d1c (MD5) Previous issue date: 2013-08-30 / O desenvolvimento de modelos matemáticos e métodos computacionais para a simulação de escoamentos em meios porosos é de grande interesse, devido à sua aplicação em diversas áreas da engenharia e ciências aplicadas. Em geral, na simulação numérica de um modelo de escoamento em meios porosos, são adotadas estratégias de desacoplamento dos sistemas de equações diferenciais parciais que o compõem. Este estudo recai sobre esquemas numéricos para leis de conservação hiperbólicas, cuja aproximação é não-trivial. Os esquemas de volumes finitos de alta resolução baseados no algoritmo REA (Reconstruct, Evolve, Average) têm sido empregados com considerável sucesso para a aproximação de leis de conservação. Recentemente, esquemas centrais de alta ordem, baseados nos métodos de Lax-Friedrichs e de Rusanov (Local Lax-Friedrichs) têm sido apresentados de forma a reduzir a excessiva difusão numérica característica destes esquemas de primeira ordem. Nesta dissertação apresentamos o estudo e a aplicação de esquemas de volumes finitos centrais de alta ordem para equações hiperbólicas que aparecem na modelagem de escoamentos em meios porosos. / The development of mathematical models and computational methods for the simulation of flow in porous media has a great interest because of its applications in engineering and other sciences. In general, in order to solve numerically the flow model in porous media the system of partial differential equations are decoupled. This study focus on the numerical schemes for the hyperbolic conservation laws, which solution is non-trivial. The finite volume schemes based on high order algorithm REA (Reconstruct, Evolve, Average) have been used with considerable success for the numerical solution of the conservation laws. Recently, high-order central schemes, based on the methods of Lax-Friedrichs and Rusanov (Local Lax-Friedrichs) have been presented, they reduce the excessive numerical diffusion presented in the first order schemes. In this dissertation we present the study and application of the high-order finite volume central schemes for hyperbolic equations as appear in the porous media flow modeling.

Orientador: Eduardo Cardoso de Abreu / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matemática Estatística e Computação Científica / Made available in DSpace on 2018-08-26T22:28:13Z (GMT). No. of bitstreams: 1 Ferraz_PaolaCunha_M.pdf: 6535346 bytes, checksum: 5f9c9ba53cd3e63fc60c09c90ad2c625 (MD5) Previous issue date: 2015 / Resumo: O foco deste trabalho é aproximação numérica de problemas envolvendo equações diferenciais parciais (EDPs), de natureza elíptica, no contexto de aplicações em dinâmica de fluidos em meios porosos. Especificamente, a dissertação pretende contribuir com uma implementação de um algoritmo multiescala e multigrid, recentemente introduzido na literatura, para resolução aproximada de sistemas de equações lineares de grande porte e mal condicionados, proveniente dessa classe de EDPs, tipicamente associada a problemas de Poisson de pressão-velocidade com condições de contornos típicas de fluxo em meios porosos. O problema concreto de Poisson discutido neste trabalho será desacoplado do sistema de transporte de EDPs de convecção-difusão, com convecção dominante, e linearizado por meio do emprego de uma técnica de decomposição de operadores. A metodologia para a discretização do problema elíptico de Poisson é elementos finitos mistos híbridos. A resolução numérica do sistema linear resultante deste procedimento será realizado via um método do tipo Gradientes Conjugados com Pré-condicionamento (PCG) multiescala e multigrid. Combinamos as metodologias multi-escala e multigrid de modo a capturar os distintos comprimentos de onda associados aos diferentes comprimentos de onda do operador diferencial auto-adjunto de Poisson, fortemente influenciado pela heterogeneidade das propriedades geológicas do meio poroso, em particular da permeabilidade absoluta, que pode exibir flutuações em várias ordens de grandeza. Experimentos computacionais em aplicações de problemas de dinâmica de fluidos em meios porosos são apresentados e discutidos para verificação dos resultados obtidos / Abstract: The focus of this work is the numerical approximation of differential problems involving partial differential equations (PDE's) of elliptic nature, in the context of modelling and simulation in fluid dynamics in porous media. The dissertation aims to contribute with an implementation of a multiscale multigrid algorithm, recently introduced in the literature, designed for solving ill-conditioned large linear systems of equations derived from those classes of PDE's, typically associated with Poisson problems of pressure-velocity with boundary conditions typical of flow in porous media. The Poisson problem discussed here is identified from the coupled convection-diffusion transport system counterpart of PDE's, with dominated convection, and by a linearization by means the use of an operator splitting approach. The methodology used for the discretization of the Poisson elliptic problem is by mixed hybrid finite elements. The numerical solution of the resulting linear system will be addressed by a multiscale multigrid preconditioned conjugate gradient (PCG) method. We combine both methodologies in order to capture the distinct wavelengths associated with the different wavelengths from the assosiated self-adjoint Poisson operator, strongly influenced by the heterogeneity of the geological properties of the porous media, in particular to the absolute permeability tensor, which in turn might exhibit very large fluctuations of orders of magnitude. Numerical experiments in applications of fluid dynamics problems in porous media are presented and discussed for a verification of the results obtained by direct numerical simulations with the multiscale multigrid algorithm under consideration / Mestrado / Matematica Aplicada / Mestra em Matemática Aplicada

Orientador: Eduardo Cardoso de Abreu / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matemática Estatística e Computação Científica / Made available in DSpace on 2018-08-27T06:42:03Z (GMT). No. of bitstreams: 1 Vieira_Jardel_M.pdf: 2696940 bytes, checksum: 9517d68c44824f91bd411caf141d2ef1 (MD5) Previous issue date: 2015 / Resumo: O foco desta dissertação de mestrado consiste em um estudo computacional de equações pseudo-parabólicas em mecânica dos fluidos e fenômenos de transporte de fluidos em meios porosos. Serão considerados problemas de valor de contorno e inicial associados a duas classes de modelos de equações de evolução pseudo-parabólicas: um modelo de advecção-difusão com termo pseudo-parabólico que exibe um certo caráter dispersivo e um outro modelo pseudo-parabólico "puro", i.e., sem a presença do termo de advecção. O primeiro modelo se relaciona com a modelagem física do fluxo de duas fases incompressíveis em dinâmica de fluidos em meios porosos, onde são considerados modelos de pressão capilar dinâmica, ou seja, em que os efeitos dinâmicos são também incluídos na diferença de pressão entre as fases fluidas. Uma discussão sobre a relevância física em aplicações e da importância matemática do sistema governante de equações para pressão capilar dinâmica em fenômenos de transporte de fluidos em meios porosos é também feita de modo a indicar algum suporte à escolha dos métodos estudados para aproximação numérica dos modelos consideradores. Além disso, um conjunto de experimentos numéricos é apresentado e discutido para avaliar a qualidade das soluções obtidas do estudo proposto, bem como para justificar variações dos métodos numéricos estudados. Especificamente, para o modelo pseudo-parabólico puro, os resultados são comparados com soluções analíticas para o caso linear. Para o modelo pseudo-parabólico com o termo de advecção, é avaliado se os resultados dos métodos numéricos empregados concordam qualitativamente com resultados da literatura / Abstract: The focus of this work consists of a computational study of pseudo-parabolic equations in fluid mechanics and transport phenomena in porous media. For concreteness, we consider initial-boundary value problems related to two classes of systems of evolution pseudo-parabolic equations: a advection-diffusion model, which in turn the pseudo-parabolic term exhibits a certain dispersive character, and a second of "purely" pseudo-parabolic nature, i.e., without the presence of advection term. The first model relates to the modeling of incompressible two-phase flow in porous media, which in turn takes into account the nonlinear dynamic capillary pressure effects, where the dynamic effects are also included into the pressure difference between the fluid phases. Further, a discussion of the physical and mathematical relevance of the governing system of equations for dynamic capillary pressure in porous media fluid transport phenomena is also made in order to drive the choice of the numerical approximations for the differential models under investigation. Moreover, a set of numerical experiments are presented and discussed to address the quality of the obtained solutions proposed study, as well as to justify variations of the numerical methods studied. Specifically, to the purely pseudo-parabolic model, the results are compared along with analytical solutions with respect to a linear case. On the other hand, to the nonlinear pseudo-parabolic model with advection term, it is performed numerical experiments in order to account the correct qualitative behavior of the computed solutions against the available results discussed in the recent literature / Mestrado / Matematica Aplicada / Mestre em Matemática Aplicada

Submitted by Renata Lopes (renatasil82@gmail.com) on 2016-01-13T12:34:27Z No. of bitstreams: 1 weslleydasilvapereira.pdf: 6406997 bytes, checksum: 17a5ee95515b6ad53c2e29db478c6e81 (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2016-01-25T17:29:23Z (GMT) No. of bitstreams: 1 weslleydasilvapereira.pdf: 6406997 bytes, checksum: 17a5ee95515b6ad53c2e29db478c6e81 (MD5) / Made available in DSpace on 2016-01-25T17:29:23Z (GMT). No. of bitstreams: 1 weslleydasilvapereira.pdf: 6406997 bytes, checksum: 17a5ee95515b6ad53c2e29db478c6e81 (MD5) Previous issue date: 2015-03-23 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / É crescente o interesse na utilização de métodos térmicos para recuperação de óleo de média e alta viscosidade. Um desses métodos é a combustão in situ, que consiste na liberação de calor no interior do reservatório através da combustão do ar injetado. As componentes mais pesadas do óleo atuam como combustível para as reações exotérmicas e o calor gerado reduz a viscosidade do óleo, estimulando o fluxo em direção aos poços de produção. Os modelos matemáticos para este método de recuperação em geral são complexos. Portanto, a obtenção de soluções analíticas para tais modelos é inviável, sendo necessária a utilização de simulações computacionais. Diversos trabalhos apresentam estudos analíticos e numéricos de modelos unidimensionais para a combustão em meios porosos. Em trabalhos anteriores, estimativas analíticas para modelos unidimensionais foram obtidas. Neste trabalho, tais estimativas são ligeiramente generalizadas através da inclusão da pressão prevalecente. É proposto um modelo bidimensional para o processo de combustão in situ em meios porosos heterogêneos que considera pressão variável. Soluções numéricas são obtidas utilizando o método de elementos finitos para a discretização espacial, o esquema de diferenças finitas de Crank-Nicolson para discretização no tempo e o método de Newton para resolução das equações não lineares resultantes. Estimativas analíticas para a temperatura e velocidade da onda de combustão são obtidas através de um modelo unidimensional simplificado. Tais estimativas são validadas com sucesso para o modelo geral através das simulações. Uma outra simplificação unidimensional do modelo geral é simulada numericamente através de duas abordagens: a primeira é similar à utilizada para a solução do modelo geral; e a segunda é escrita como um problema de complementaridade. Os problemas de complementaridade não-linear são resolvidos pelo algoritmo FDA-NCP. As duas abordagens numéricas utilizadas são comparadas com uma estimativa analítica para a onda térmica e mostram bons resultados. / There is a growing interest in using thermal methods for the recovery of medium and high viscosity oil. One of these methods is the in-situ combustion, which consists in release heat within the reservoir through combustion of the injected air. The heavier oil components are used as fuel for exothermic reactions and the generated heat reduces the oil viscosity, stimulating the flow towards the production well. In general, the mathematical models for this recovery method are complex. Therefore, the analytical solutions for such models are impossible, requiring numerical simulations. Several works present analytical and numerical studies of one-dimensional models for combustion in porous media. In previous works analytical estimates for one dimensional models were obtained. Here these estimates are slightly generalized by including the prevailing pressure. We propose a two-dimensional model for the in-situ combustion process in heterogeneous porous media, considering variable pressure. Numerical results are obtained using the finite element method for spatial discretization, Crank-Nicolson finite difference scheme for time discretization and Newton’s method for the arising nonlinear equations. Analytical estimates for combustion wave speed and combustion wave temperature are obtained using one-dimensional simplified model. These estimates are successfully validated in the general model through the simulation results. Another one-dimensional simplification of the general model is numerically simulated by two approaches: the first is similar to the one previously described; and the second one is written as a complementarity problem. The arising nonlinear complementarity problems are solved by the FDA-NCP algorithm. Both numerical approaches are compared to the analytical estimate for the thermal wave, showing good agreement.

[pt] Estima-se que 50 por cento das reservas mundiais de petróleo e gás sejam mantidas em reservatórios carbonáticos naturalmente fraturados. Um dos maiores desafios neste tipo de formações é a sua heterogeneidade. Além da presença de fraturas que conectam longitudinalmente o meio poroso, vugs em diferentes escalas e distribuições estão espalhados por toda a matriz porosa. Essas cavidades tornam as características do escoamento de fluidos significativamente diferentes daquelas dos reservatórios convencionais de estrutura porosa homogênea, e trazem a necessidade de avaliar propriedades petrofísicas equivalentes para o meio heterogêneo. Neste estudo, uma abordagem microfluídica é usada para determinar as curvas de permeabilidade relativa de água e óleo e os perfis de distribuição das fases em micromodelos 2D de meios porosos vugulares. Experimentos de injeção simultânea de água-óleo em estado estacionário foram realizados a diferentes fluxos fracionários de água, monitorando a dinâmica da queda de pressão e visualizando o deslocamento de fluidos na escala de poros. A aquisição de imagens em tempo real por microscopia de fluorescência permitiu examinar a evolução da saturação das fases. A comparação direta entre as curvas de permeabilidade relativa dos meios porosos vugulares com aquela da matriz porosa mostrou que a incorporação de vugs leva a (i) maior permeabilidade absoluta equivalente, especialmente com cavidades mais longas e em maior número, (ii) aumento da ocupação de óleo na matriz porosa, devido à invasão de água menos eficiente, e (iii) maior permeabilidade relativa à água, que flui preferencialmente pelo espaço vugular. Esses resultados são consistentes com a natureza molhável ao óleo dos micromodelos, uma vez que os vugs oferecem menor resistência capilar ao fluxo da fase não molhante. Nossa abordagem microfluídica de baixo custo provavelmente nos permitirá estudar sistematicamente configurações mais complexas de meios porosos heterogêneos / [en] It is estimated that 50 percent of world s oil and gas reserves are held in naturally fractured carbonate reservoirs. One of the biggest challenges in this type of formation is its heterogeneous nature. Besides the presence of fractures that longitudinally connect the porous medium, vugs at different scales and distributions are scattered throughout the porous matrix. These cavities cause fluid flow characteristics to significantly differ from those of conventional homogeneous pore structure reservoirs and bring the need to evaluate equivalent petrophysical properties of the heterogeneous medium. In this study, a microfluidic approach is used to determine the water and oil relative permeability curves and phase distribution profiles in 2D micromodels of vugular porous media. Steady-state water-oil injection experiments were performed in these devices at different fractional flows, while monitoring the dynamics of the pressure drop and visualizing the fluid displacement at the pore scale. Live-image acquisition through fluorescence microscopy made it possible to examine the evolution of the saturation of water and oil phases. The direct comparison between the relative permeability curves of well-characterized vugular porous media and their porous matrix showed that the incorporation of vugs leads to (i) higher equivalent absolute permeability, especially with longer cavities and higher vug density, (ii) increased oil occupancy in the porous matrix, due to less efficient water invasion into the porous matrix, and (iii) higher relative permeability to water, which flows preferentially through the vugular space. These results are consistent with the oil-wet nature of micromodels, since the vugs are offering less capillary resistance to the flow of the non-wetting phase. Our low-cost microfluidic approach will likely allow us to systematically study more complex vugular-fractured systems.

[pt] Tradicionalmente, os testes de formação em poços de petróleo buscam caracterizar o campo de permeabilidades a partir da interpretação dos transientes de pressão (PTA) nos períodos de fluxo e estática baseados em modelos isotérmicos de escoamento em meios porosos. Com o avanço da instrumentação dos testes, registros mais precisos de temperatura passaram a estar disponíveis e fomentaram a pesquisa baseada em modelos não isotérmicos que possibilitaram a análise a partir dos transientes de temperatura (TTA). Além da caracterização de parâmetros do reservatório como permeabilidade e porosidade com a interpretação dos transientes de temperatura, os dados de pressão obtidos a partir de um modelo não isotérmico representa de forma mais fidedigna o fenômeno físico sobretudo quando os testes são submetidos a maiores diferenciais de pressão. Este trabalho consiste no desenvolvimento de um simulador para teste de formação que considera a modelagem não isotérmica de reservatório unidimensional radial acoplado a um poço produtor e na utilização deste simulador, associado a métodos de otimização multivariável, para resolução do problema inverso da caracterização de parâmetros do reservatório. Alguns métodos de otimização foram testados e o algoritmo do Simplex de Nelder-Mead apresentou melhor eficácia. Foram estabelecidos três tipos de problemas e utilizados em três casos hipotéticos considerando inclusive a imposição artificial de ruídos nos sinais de pressão e temperatura utilizados para resolução do problema inverso. / [en] Traditionally, oil well formation tests aim to characterize the reservoir permeability field from pressure transient analysis (PTA) of drawdown and build up based on isothermal flow models in porous media. With the advancement of well test instrumentation, more accurate temperature records became available and have encouraged researches based on non-isothermal models that made possible the temperature transient analysis (TTA). In addition to the characterization of reservoir parameters such as permeability and porosity by TTA, the pressure data obtained from a non-isothermal model represent better the physical phenomenon, especially when the tests are subjected to greater drawdowns. This work consists in the development of a simulator for formation test that considers non-isothermal modeling of a unidimensional radial reservoir coupled to a production well and in the use of this simulator, associated with multivariable optimization methods, to solve the inverse problem of reservoir parameters characterization. Some optimization methods were tested and the Nelder-Mead Simplex algorithm presented better efficiency. Three types of problems were established and used in three hypothetical cases, including artificially imposed noise in pressure and temperature signals used to solve the inverse problem.

[pt] Nas últimas décadas, grandes investimentos foram feitos no desenvolvimento de modelos e métodos numéricos para prever e analisar os diferentes aspectos do processo de recuperação de óleo e gás. Neste contexto, os simuladores modernos devem ser capazes de incorporar uma grande variedade de opções para responder questões relacionadas ao gerenciamento de reservatórios de forma rápida e precisa. Neste trabalho, nós apresentamos um simulador de reservatórios baseado em uma arquitetura de plugins, onde diferentes formulações, solvers e modelos podem ser desenvolvidos, estendidos e aprimorados. A partir desta abordagem, utilizamos o modelo black-oil para implementar técnicas tradicionais e do estado da arte, como os métodos totalmente e adaptativamente implícito, os métodos de Newton-Raphson e Newton Inexato, controladores heurístico e PID para passo de tempo adaptativo e aproximações de fluxo de um ponto baseados no potencial de fase tradicional e C1-contínuo. Diversas configurações de plugins foram testadas e validadas com simuladores comerciais e seus desempenhos foram utilizados para determinar quais as mais adequadas para resolver problemas de escoamento multifásico. / [en] During the last decades, large investments were made towards the development of numerical models and methods to forecast and analyze the different aspects of oil and gas recovery. In this context, modern simulators must be able to incorporate a wide variety of options to answer questions related to reservoir management accurately and effectively. In this work, we present a reservoir simulator based on a plugin architecture, where different formulations, solvers, and models can be developed, extended, and enhanced. With this approach, we use the black-oil model to implement traditional and state-of-the-art techniques, including fully- and adaptiveimplicit methods, heuristic and PID time-step controllers, Newton-Raphson and Inexact Newton, and C1-continuous and conventional phase-potential single-point upstream weighting. Several plugin configurations were tested and validated with commercial simulators, and their performances were used to determine which are the most suitable to solve multiphase flow problems.

[pt] A microtomografia computadorizada de uma amostra de rocha possibilita uma caracterização do meio poroso e pode ser utilizada para estimar propriedades da rocha em macroescala, isto é, em escala de reservatório. Métodos baseados em mapas de distâncias e em algoritmos de erosão são as principais abordagens utilizadas para extração de uma rede de poros e gargantas a partir de imagens microtomográficas de rocha. Este trabalho propõe um método híbrido para a construção da rede, de modo que, durante o processo de modelagem na escala de poros, obtemos um esqueleto do espaço poroso por meio de um algoritmo de erosão e utilizamos um mapa de distâncias para construir uma rede de poros e gargantas. A determinação dos poros e gargantas a partir do esqueleto adota uma abordagem não-determinística possibilitando a geração de múltiplas redes com configurações distintas a partir de um mesmo esqueleto. Avaliamos a variabilidade dos cenários gerados e comparamos as estimativas para as propriedades petrofísicas com as obtidas pelo método de Bolas Máximas por meio dos resultados de uma simulação de fluxo monofásica na rede. / [en] Computerized microtomography of a rock sample enables a characterization of the porous medium and can be used to estimate rock properties at the macro-scale, i.e., reservoir-scale. Methods based on distance maps and thinning algorithms are the main approaches used for extracting a pore and throats network from microtomographic rock images. This paper proposes a hybrid method for constructing the network. So that during the pore-scale modeling process, we obtain a skeleton of the pore space by using a thinning algorithm and a distance map to build a network of pores and throats. The determination of pores and throats from the skeleton assumes a non-deterministic approach enabling the generation of multiple networks with distinct configurations from the same skeleton. We evaluate the variability of the generated scenarios and compare the estimates for the petrophysical properties with those obtained by the Maximum Ball Method through the results of a single-phase flow simulation on the network.

[pt] A teoria da homogeneização de equações diferenciais tornou-se um campo aberto de pesquisa em diversas áreas das ciências exatas e mostrou-se ser uma poderosa ferramenta para a compreensão do comportamento global de materiais heterogêneos. Apesar de ser conhecido que a dedução da lei de Darcy através das equações de Navier-Stokes já é tema debatido há décadas muitas questões continuam em aberto, principalmente a respeito de condições de contorno mais complexas, casos envolvendo fluxos multifásicos e técnicas de homogeneização numérica. Sabe-se que a lei de Darcy se apresenta sob forma de uma relação linear apenas para um intervalo de gradiente hidráulico e que este intervalo se sobrepõe ao intervalo de fluxo laminar do fluido através dos vazios do solo.Se propõe neste trabalho, então, compreender a perda de linearidade na lei de Darcy, a partir da teoria da homogeneização, modificando e explorando os resultados obtidos anteriormente na literatura. / [en] The theory of homogenization of differential equations has become an open field of research in several areas of the exact sciences and has proved to be a powerful tool for understanding the global behavior of heterogeneous materials. Despite knowing that the deduction of Darcy s law through the Navier-Stokes equations has been debated for decades, many questions remain open, mainly regarding more complex boundary cenditions, cases involving multiphase flows and the numerical homogenization techniques. It is known that Darcy s law is presented in the form of a linear relationship only for a range of hydraulic gradient that overlaps the range of laminar flow of fluids through soil voids. Therefore, it is proposed in this work to understand the loss of linearity in Darcy s law, based on the theory of homogenization, modifying and exploring the limit results obtained in the literature.