Spelling suggestions: "subject:"multipoint flux approximation methods"" "subject:"multipoints flux approximation methods""
1 |
Simulação Numérica de Escoamento Bifásico em reservatório de Petróleo Heterogêneos e Anisotrópicos utilizando um Método de Volumes Finitos “Verdadeiramente” Multidimensional com Aproximação de Alta OrdemSOUZA, Márcio Rodrigo de Araújo 22 September 2015 (has links)
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.
|
Page generated in 0.1283 seconds