Submitted by Fabio Sobreira Campos da Costa (fabio.sobreira@ufpe.br) on 2017-07-14T12:28:59Z No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) Dissertacao Bautista Nilson.pdf: 11338318 bytes, checksum: ae22c70eb1719f066a5eeb3de436c953 (MD5) / Made available in DSpace on 2017-07-14T12:28:59Z (GMT). No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) Dissertacao Bautista Nilson.pdf: 11338318 bytes, checksum: ae22c70eb1719f066a5eeb3de436c953 (MD5) Previous issue date: 2016-01-26 / CNPQ / A modelagem e a simulação numérica do transporte de solutos, como por exemplo traçadores, em meios porosos heterogêneos e anisotrópicos, tais como aquíferos e reservatórios de petróleo constituem-se num grande desafio de natureza matemática e numérica. A modelagem de falhas selantes, canais, poços inclinados, pinchouts e outras características complexas demanda o uso de malhas não-estruturadas e não-ortogonais, capazes de se adaptar naturalmente ao domínio em estudo. Os pacotes computacionais utilizados comumente na indústria do petróleo, na sua grande maioria, se baseiam no Método das Diferenças Finitas com Aproximação de Fluxo por Dois Pontos (Two-Point Flux Approximation - TPFA) e no Método de Ponderação à Montante de Primeira Ordem (First Order Upwind Method - FOU), devido a sua facilidade de implementação e sua eficiência computacional. Infelizmente, os métodos TPFA são incapazes de produzir soluções convergentes em malhas não-ortogonais ou para tensores de dispersão ou permeabilidades completos e os métodos FOU produzem soluções com difusão numérica excessiva, exigindo malhas demasiadamente refinadas para obtermos soluções confiáveis. Uma alternativa ao TPFA, e que permite o uso de tensores completos e malhas não-ortogonais, é o Método dos Elementos Finitos de Galerkin (MEF), porém este método não produz soluções localmente conservativas, o que pode ser um problema sério para a modelagem de problemas envolvendo leis de conservação, como no escoamento em meios porosos. Outra alternativa são os Métodos de Volumes Finitos (MVF). Nas suas variantes mais robustas, estes métodos são capazes de lidar com malhas poligonais quaisquer e tensores de dispersão e permeabilidades completos e com razão de anisotropia arbitrária, além de produzir aproximações discretas de alta ordem e localmente conservativas. Neste contexto, no presente trabalho, apresentamos uma formulação MVF centrado na célula para a modelagem do transporte de um traçador não-reativo num escoamento monofásico em meios porosos heterogêneos e anisotrópicos. Para a discretização dos termos elípticos, tanto da equação de pressão quanto da equação de Advecção-Dispersão-Reação (ADRE), utilizou-se um MVF com aproximação de fluxo por múltiplos pontos que faz uso do estêncil diamante (MPFA-D) e para a discretização dos termos hiperbólicos, usamos o método FOU e um MVF do tipo MUSCL (Monotone Upstream Centered Scheme for Conservation Laws). A fim de testar nossa formulação, resolvemos alguns problemas benchmark encontrados na literatura. / Modeling and numerical simulation of solutes (e.g. Tracers) in heterogeneous and anisotropic porous media such as aquifers and oil reservoirs, constitute a bigger challenge of mathematics and numerical nature. Modeling sealants faults, channels, inclined wells, pinch outs and other complex features of these geological formations demand the use of unstructured and not orthogonal meshes, able to adapt naturally to the domain under study. The computational packages used commonly in the oil industry, mostly, are based on the Finite Difference Method with Two Point Flow Approximation (TPFA) and the Amount First Order Upwind method (FOU), due to its ease of implementation and its computational efficiency. Unfortunately, TPFA methods are unable to produce conver-gent solutions in non-orthogonal meshes or in permeability or dispersion full Tensor and FOU methods produce solutions with excessive numerical diffusion, requiring excessively refined mesh to obtain reliable solutions. An interesting alternative to TPFA, which allows the use of full tensor and not orthogonal meshes, is the Galerkin Finite Element Method (FEM), but this method does not produce solutions locally conservative, which can be a serious problem for modeling problems involving conservation laws as the flow in porous media. An interesting alternative is the Finite Volume Methods (MVF). In its most robust embodiments, these methods are able to cope with any polygonal mesh and full permeability or dispersion tensors and with an arbitrary anisotropy ratio, beyond producing discrete approximations of high order and locally conservative. In this context, the present study, we present one MVF formulation cell centered to modeling the transport of a non-reactive tracer in single-phase flow in heterogeneous and anisotropic porous media. For the elliptical discretization terms, both, the pressure equation as the equation advection-dispersion-reaction (ADRE), we used The FVMF multipoint flow approximation that uses the diamond stencil (MPPA-D) and for the discretization of hyperbolic terms, we use the FOU method and an MVF type MUSCL (Monotone Upstream Centered Scheme for Conservation Laws). In order to test our formulation, we solve some benchmark problems in the literature.

The groundwater flow in aquifers contain evaporite rocks can cause problems such as geo-mechanical subsidence or collapse. In this work, we focus on the development of numerical models to simulate the flow in porous and non-porous domains in order to study the dissolution phenomenon and fractures evolution over time. The first part of this thesis is devoted to developing new solutions for the validation of numerical models to simulate density driven flow in porous media. The new procedure consist of solving simultaneously the flow and the transport equations using the Levenberg-Marquardt algorithm. The use of this technique allowed us to develop, for the first time, semi-analytical solutions of saltwater intrusion in the case of small diffusion and in the case of a large density contrast. In the second part of this work, we studied the flow in evaporitic rocks. A numerical code was developed to solve the nonlinear system using advanced numerical methods. To validate this new model, we have developed a semi-analytical solution for the density Stokes flow. The third part of this work is devoted to transport with dissolution of rock salt. As a first step, we studied the influence of various parameters on the dissolution of salt in Adlertunnel located at a depth of 160 m in the region of Basel in Switzerland. In a second step, we are interested in the simulation of the fracture's evolution as a result of the dissolution. The numerical model takes into account the Stokes flow and mass transport effects and dissolution of the fracture walls.

Submitted by Fernanda Rodrigues de Lima (fernanda.rlima@ufpe.br) on 2018-07-23T21:05:22Z No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO Braian Schneider Brum.pdf: 15683511 bytes, checksum: 80e1652e52ac088c04dacdef6294de42 (MD5) / Approved for entry into archive by Alice Araujo (alice.caraujo@ufpe.br) on 2018-07-24T18:17:55Z (GMT) No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO Braian Schneider Brum.pdf: 15683511 bytes, checksum: 80e1652e52ac088c04dacdef6294de42 (MD5) / Made available in DSpace on 2018-07-24T18:17:55Z (GMT). No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO Braian Schneider Brum.pdf: 15683511 bytes, checksum: 80e1652e52ac088c04dacdef6294de42 (MD5) Previous issue date: 2017-01-31 / CNPQ / O problema do escoamento bifásico de água e óleo em reservatórios de petróleo heterogêneos e anisotrópicos pode ser descrito utilizando certas hipóteses simplificadoras, por um sistema de equações diferenciais parciais não-lineares, composto por uma equação elíptica da pressão e uma equação hiperbólica de saturação. A modelagem deste problema representa um grande desafio, devido à complexidade dos ambientes deposicionais, incluindo camadas inclinadas, canais, fraturas e a possível existência de poços direcionais. Nesses casos, é particularmente complexo construir malhas estruturadas capazes de modelar o problema físico adequadamente. Esta dissertação tem como objetivo a modelagem numérica do escoamento bifásico em reser-vatórios naturalmente fraturados (RNF). Para isso, são apresentadas duas formulações pelo Método de Volumes Finitos (MVF) centrado na célula para a solução da equação de pressão: a primeira utiliza uma Aproximação de Fluxo por Dois Pontos (TPFA) em conjunto com multi-plicadores de transmissibilidade, e a segunda, capaz de lidar com tensores de permeabilidade completos e malhas poligonais, utiliza uma Aproximação do Fluxo por Múltiplos Pontos (MPFA-O) acoplada ao Modelo de Fraturas com Dimensão Reduzida (Lower-Dimensional Fracture Model - LDFM). O problema hiperbólico da saturação é resolvido pelo método de ponderação a montante de primeira ordem (First Order Upwind Method - FOUM) com duas formulações, uma explícita e outra implícita. A abordagem LDFM é eficiente quando compa-rada a outras estratégias apresentadas na literatura, e faz uso de uma equação adicional associ-ada à fratura que é tratada como uma entidade geométrica com dimensão inferior à do problema original, ou seja, para problemas em 3-D uma fratura é representada por superfícies (2-D) e para problemas em 2-D a fratura tem apenas uma dimensão no espaço. Isso reduz considera-velmente o número de graus de liberdade do sistema. Para isso, a malha deve ser ajustada à distribuição das fraturas no domínio computacional, de modo que essas fraturas são associadas às superfícies de controle, representadas pelas arestas dos volumes finitos em 2-D. É importante observar, que o campo de velocidades nas superfícies de controle que coincidem com as fraturas é dependente, tanto das pressões nas fraturas, quando das pressões nos volumes de controle que representam a rocha matriz e que compartilham os mesmos vértices da respectiva superfície de controle. A acurácia da formulação proposta foi verificada através da resolução de alguns pro-blemas envolvendo meio fraturado, incluindo barreiras e canais. / The problem of two phase flow of oil and water in heterogeneous and anisotropic oil reservoirs can be described using certain simplifying hypotheses by a system of nonlinear partial differential equations, composed by an elliptic equation of pressure and a hyperbolic equation of saturation. The modeling of this problem represents a major challenge due to the complexity of the depositional environments, including inclined layers, channels, fractures and the possible existence of directional wells. In such cases, it is particularly complex to construct structured meshes capable of properly modeling the physical problem. This dissertation aims to model one of these complex structures, naturally fractured reservoirs (NFR). For this, are presented two formulations with Finite Volume Method (FVM), cell-centered, for the solution of the pressure equation: the first one using a Two-Point Flux Aproximation (TPFA) in conjunction with transmissibility multipliers; and the second formulation using the Multi-Point Flux Aproximation (MPFA-O), which can deal with full permeability tensors and arbitrary polygonal meshes, coupled with a Lower-Dimensional Fracture Model (LDFM). The hyperbolic problem of saturation is solved with First Order Upwind Method (FOUM) with two formulations, one explicit and another implicit. The LDFM approach is effective when compared to other strategies in the literature, and makes use of an additional equation associated to the fracture, that is treated as a geometric entity with a smaller dimension than the original problem, in other words, for 3-D problems a fracture is represented by surfaces (2-D) and for 2-D problems the fracture has only one dimension in space. This considerably reduces the number of degrees of freedom of the system. For this, the mesh must be adjusted to the distribution of fractures in the computational domain, when these fractures are associated to the control surfaces, represented by the edges of the finite volumes in 2-D. It is important to note that the velocity field on the control surfaces that coincide with the fractures is dependent, both pressures in the fractures and in the control volumes, which represent the matrix rock and which share the same vertices of the respective control surface. The accuracy of the proposed formulation was verified by solving some problems involving a fractured matrix, including barriers and channels.

La circulation de l’eau souterraine dans les aquifères contenants des roches évaporitiques peut provoquer des problèmes géo-mécaniques tels que l'affaissement du sol ou l'effondrement. Dans ce travail, nous nous intéressons au développement de modèles numériques permettant de simuler les écoulements dans les milieux poreux et non poreux ainsi que les phénomènes de dissolution et d’évolution des fractures dans le temps. La première partie de cette thèse est consacrée au développement de nouvelles solutions pour la validation des modèles numériques simulant les écoulements densitaires en milieux poreux. La nouvelle procédure consiste à résoudre simultanément les deux systèmes d’écoulement et de transport en utilisant l’algorithme de Levenberg-Marquardt. L’utilisation de cette technique nous a permis de développer, pour la première fois, des solutions semi-analytiques d’intrusion d’eau salée dans le cas de faible diffusion ainsi que dans le cas d’un grand contraste de densité. Dans la deuxième partie de ce travail, nous nous sommes intéressés aux écoulements dans les fractures des roches évaporitiques. Un code de calcul a été développé pour résoudre ce système non linéaire en utilisant des méthodes numériques adaptées. Pour valider ce nouveau modèle, nous avons développé une solution semi-analytique pour les écoulements densitaires de Stokes. La troisième partie de ce travail est consacrée au transport avec dissolution de la roche salée. Dans un premier temps, nous avons étudié l’influence de différents paramètres sur la dissolution du sel dans l’Adler tunnel situé à une profondeur de 160 m dans la région de Bâle en Suisse. Dans un second temps, nous nous sommes intéressés à la simulation de l’évolution dune fracture sous l’effet de la dissolution. Le modèle numérique développé prend en compte les écoulements de Stokes ainsi que le transport de masse avec effets densitaires et la dissolution des parois de la fracture. / The groundwater flow in aquifers contain evaporite rocks can cause problems such as geo-mechanical subsidence or collapse. In this work, we focus on the development of numerical models to simulate the flow in porous and non-porous domains in order to study the dissolution phenomenon and fractures evolution over time. The first part of this thesis is devoted to developing new solutions for the validation of numerical models to simulate density driven flow in porous media. The new procedure consist of solving simultaneously the flow and the transport equations using the Levenberg-Marquardt algorithm. The use of this technique allowed us to develop, for the first time, semi-analytical solutions of saltwater intrusion in the case of small diffusion and in the case of a large density contrast. In the second part of this work, we studied the flow in evaporitic rocks. A numerical code was developed to solve the nonlinear system using advanced numerical methods. To validate this new model, we have developed a semi-analytical solution for the density Stokes flow. The third part of this work is devoted to transport with dissolution of rock salt. As a first step, we studied the influence of various parameters on the dissolution of salt in Adlertunnel located at a depth of 160 m in the region of Basel in Switzerland. In a second step, we are interested in the simulation of the fracture’s evolution as a result of the dissolution. The numerical model takes into account the Stokes flow and mass transport effects and dissolution of the fracture walls.

Submitted by Pedro Barros (pedro.silvabarros@ufpe.br) on 2018-09-21T19:43:58Z No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) TESE Fernando Raul Licapa Contreras.pdf: 6260685 bytes, checksum: d6d66a2bde6c2b142b3b9f666d333f4e (MD5) / Approved for entry into archive by Alice Araujo (alice.caraujo@ufpe.br) on 2018-09-24T23:13:56Z (GMT) No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) TESE Fernando Raul Licapa Contreras.pdf: 6260685 bytes, checksum: d6d66a2bde6c2b142b3b9f666d333f4e (MD5) / Made available in DSpace on 2018-09-24T23:13:56Z (GMT). No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) TESE Fernando Raul Licapa Contreras.pdf: 6260685 bytes, checksum: d6d66a2bde6c2b142b3b9f666d333f4e (MD5) Previous issue date: 2017-01-10 / FACEPE / No presente trabalho apresentamos formulações numéricas para a discretização de um modelo matemático que descreve o escoamento bifásico de água e óleo em reservatórios de petróleo heterogêneos e anisotrópicos. Estas equações são resolvidas através da metodologia IMPES (Implicit Pressure Explicit Saturation) que é, em particular, um método segregado para tratar escoamento bifásico (água-óleo) em que, o campo de pressão é obtida implicitamente e o campo de saturação é obtida explicitamente. Para resolver a equação de pressão propormos formulações numéricas baseadas no método de volumes finitos lineares (a saber: MPFA-HD e MPFA-H) e no método de volumes finitos não lineares (a saber: NLFV-PP e NLFV-DMP), estas formulações em geral lidam com meios altamente heterogêneos e anisotrópicos sobre malhas poligonais em geral. De outra parte, neste trabalho são estudadas as varias propriedades numéricas que satisfazem as formulações numéricas propostas. Além disso, no caso do método não linear que satisfaz o principio do máximo discreto (DMP), mostramos que este possui um estêncil pequeno desde que os pontos de interpolação sejam calculados a partir dos pontos harmônicos. Para resolver a equação de saturação, utilizamos um método multidimensional nomeada MOOD (Multi-dimensional Optimal Order Detection), esta método é baseada numa discretização espacial de alta ordem com procedimento de limitação a posteriori das variáveis de estado (i. e. saturação) não físicos. Em cada volume de controle o grau polinomial ótimo é determinado para construir uma aproximação das soluções respeitando o principio do máximo discreto. Mediante um processo iterativo o grau polinomial é reduzido quando o volume de controle em questão não satisfaz as condições de estabilidade. Neste método, uma estratégia baseada na Técnica de Mínimos Quadrados é usada para a reconstrução polinomial. Todos métodos propostos são localmente conservativos, e todas as incógnitas são completamente centradas nos volumes de controle, e possuem a habilidade de representar reservatórios com propriedades geológicas complexas utilizando malhas poligonais quaisquer. O desempenho dessas foi avaliado a partir da solução de problemas relevantes encontrados na literatura. / In the present work we present numerical formulations for the discretization of a mathematical models that describes the two-phase flow of water and oil in heterogenous and anisotropic oil reservoirs. These equations are solved through Implicit Pressure Explicit Saturation (IMPES) methodology which is in particular a segregated method to treat two-phase (water-oil) flow in which the pressure field is implicitly obtained and the saturation field is obtained explicitly. In order to solve the pressure equation we propose numerical formulations based on the finite linear volume method (MPFA-HD and MPFA-H) and non-linear finite volume method (NLFV-PP and NLFV-DMP). Formulations generally deal with highly heterogeneous and anisotropic media over polygonal meshes in general. On the other hand, in this work the numerical properties that satisfy the proposed numerical formulations are studied. In addition, in the case of the non-linear method that satisfies the discrete maximum principle (DMP), we show that it has a small stencil since the interpolation points are calculated from the harmonic points. To solve the saturation equation, we use a multidimensional method named MOOD (Multidimensional Optimal Order Detection), this method is based on a high order spatial discretization with a posteriori limitation procedure of the non-physical state variables (saturation). In each control volume the optimum polynomial degree is determined to construct an approximation of the solutions respecting the discrete maximum principle. By means of an iterative process the polynomial degree is reduced when the control volume in question does not satisfy the stability conditions. In this method, a strategy based on the Least Squares technique is used for polynomial reconstruction. All proposed methods are locally conservative, and all unknowns are completely centered on control volumes, and have the ability to represent reservoirs with complex geological properties using any polygonal meshes. The performance of these was evaluated from the solution of relevant problems found in the literature.

Submitted by Fabio Sobreira Campos da Costa (fabio.sobreira@ufpe.br) on 2016-07-01T15:05:14Z No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Souza_Tese_2015_09_22.pdf: 8187999 bytes, checksum: 664629aed28d692dce410fefbfe793dc (MD5) / Made available in DSpace on 2016-07-01T15:05:14Z (GMT). No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Souza_Tese_2015_09_22.pdf: 8187999 bytes, checksum: 664629aed28d692dce410fefbfe793dc (MD5) Previous issue date: 2015-09-22 / Anp / Sob certas hipóteses simplificadoras, o modelo matemático que descreve o escoamento de água e óleo em reservatórios de petróleo pode ser representado por um sistema não linear de Equações Diferenciais Parciais composto por uma equação elíptica de pressão (fluxo) e uma equação hiperbólica de saturação (transporte). Devido a complexidades na modelagem de ambientes deposicionais, nos quais são incluídos camadas inclinadas, canais, falhas e poços inclinados, há uma dificuldade de se construir um modelo que represente adequadamente certas características dos reservatórios, especialmente quando malhas estruturadas são usadas (cartesianas ou corner point). Além disso, a modelagem do escoamento multifásico nessas estruturas geológicas incluem descontinuidades na variável e instabilidades no escoamento, associadas à elevadas razões de mobilidade e efeitos de orientação de malha. Isso representa um grande desafio do ponto de vista numérico. No presente trabalho, uma formulação fundamentada no Método de Volumes Finitos é estudada e proposta para discretizar as equações elíptica de pressão e hiperbólica de saturação. Para resolver a equação de pressão três formulações robustas, com aproximação dos fluxos por múltiplos pontos são estudadas. Essas formulações são abeis para lidar com tensores de permeabilidade completos e malhas poligonais arbitrárias, sendo portanto uma generalização de métodos mais tradicionais com aproximação do fluxo por apenas dois pontos. A discretização da equação de saturação é feita com duas abordagens com característica multidimensional. Em uma abordagem mais convencional, os fluxos numéricos são extrapolados diretamente nas superfícies de controle por uma aproximação de alta resolução no espaço (2ª a 4ª ordem) usando uma estratégia do tipo MUSCL. Uma estratégia baseada na Técnica de Mínimos Quadrados é usada para a reconstrução polinomial. Em uma segunda abordagem, uma variação de uma esquema numérico Verdadeiramente Multidimensional é proposto. Esse esquema diminui o efeito de orientação de malha, especialmente para malhas ortogonais, mesmo embora alguma falta de robustez possa ser observada pra malhas excessivamente distorcidas. Nesse tipo de formulação, os fluxos numéricos são calculados de uma forma multidimensional. Consiste em uma combinação convexa de valores de saturação ou fluxo fracionário, seguindo a orientação do escoamento através do domínio computacional. No entanto, a maioria dos esquemas numéricos achados na literatura tem aproximação apenas de primeira ordem no espaço e requer uma solução implícita de sistemas algébricos locais. Adicionalmente, no presente texto, uma forma modificada desses esquemas “Verdadeiramente” Multidimensionais é proposta em um contexto centrado na célula. Nesse caso, os fluxos numéricos multidimensionais são calculados explicitamente usando aproximações de alta ordem no espaço. Para o esquema proposto, a robustez e o caráter multidimensional também leva em conta a distorção da malha por meio de uma ponderação adaptativa. Essa ponderação regula a característica multidimensional da formulação de acordo com a distorção da malha. Claramente, os efeitos de orientação de malha são reduzidos. A supressão de oscilações espúrias, típicas de aproximações de alta ordem, são obtidas usando, pela primeira vez no contexto de simulação de reservatórios, uma estratégia de limitação multidimensional ou Multidimensional Limiting Process (MLP). Essa estratégia garante soluções monótonas e podem ser usadas em qualquer malha poligonal, sendo naturalmente aplicada em aproximações de ordem arbitrária. Por fim, de modo a garantir soluções convergentes, mesmo para problemas tipicamente não convexos, associados ao modelo de Buckley-Leverett, uma estratégia robusta de correção de entropia é empregada. O desempenho dessas formulações é verificado com a solução de problemas relevantes achados na literatura. / Under certain simplifying assumptions, the problem that describes the fluid flow of oil and water in heterogeneous and anisotropic petroleum reservoir can be described by a system of non-linear partial differential equations that comprises an elliptic pressure equation (flow) and a hyperbolic saturation equation (transport). Due to the modeling of complex depositional environments, including inclined laminated layers, channels, fractures, faults and the geometrical modeling of deviated wells, it is difficult to properly build and handle the Reservoir Characterization Process (RCM), particularly by using structured meshes (cartesian or corner point), which is the current standard in petroleum reservoir simulators. Besides, the multiphase flow in such geological structures includes the proper modeling of water saturation shocks and flow instabilities associated to high mobility ratios and Grid Orientation Effects (GOE), posing a great challenge from a numerical point of view. In this work, a Full Finite Volume Formulation is studied and proposed to discretize both, the elliptic pressure and the hyperbolic saturation equations. To solve the pressure equation, we study and use three robust Multipoint Flux Approximation Methods (MPFA) that are able to deal with full permeability tensors and arbitrary polygonal meshes, making it relatively easy to handle complex geological structures, inclined wells and mesh adaptivity in a natural way. To discretize the saturation equation, two different multidimensional approaches are employed. In a more conventional approach, the numerical fluxes are extrapolated directly on the control surfaces for a higher resolution approximation in space (2nd to 4th order) by a MUSCL (Monotone Upstream Centered Scheme for Conservation Laws) procedure. A least squares based strategy is employed for the polynomial reconstruction. In a second approach, a variation of a “Truly” Multidimensional Finite Volume method is proposed. This scheme diminishes GOE, especially for orthogonal grids, even though some lack of robustness can be observed for extremely distorted meshes. In this type of scheme, the numerical flux is computed in each control surface in a multidimensional way, by a convex combination of the saturation or the fractional flow values, following the approximate wave orientation throughout the computational domain. However, the majority of the schemes found in literature is only first order accurate in space and demand the implicit solution of local conservation problems. In the present text, a Modified Truly Multidimensional Finite Volume Method (MTM-FVM) is proposed in a cell centered context. The truly multidimensional numerical fluxes are explicitly computed using higher order accuracy in space. For the proposed scheme, the robustness and the multidimensional character of the aforementioned MTM-FVM explicitly takes into account the angular distortion of the computational mesh by means of an adaptive weight, that tunes the multidimensional character of the formulation according to the grid distortion, clearly diminishing GOE. The suppression of the spurious oscillations, typical from higher order schemes, is achieved by using for the first time in the context of reservoir simulation a Multidimensional Limiting Process (MLP). The MLP strategy formally guarantees monotone solutions and can be used with any polygonal mesh and arbitrary orders of approximation. Finally, in order to guarantee physically meaningful solutions, a robust “entropy fix” strategy is employed. This produces convergent solutions even for the typical non-convex flux functions that are associated to the Buckley-Leverett problem. The performance of the proposed full finite volume formulation is verified by solving some relevant benchmark problems.