Made available in DSpace on 2014-12-17T14:08:47Z (GMT). No. of bitstreams: 1 MartinhoQAF_DISSERT.pdf: 3517855 bytes, checksum: e2d79876500fdde55b949a6a832b6ff0 (MD5) Previous issue date: 2011-05-27 / Water injection is the most widely used method for supplementary recovery in many oil fields due to various reasons, like the fact that water is an effective displacing agent of low viscosity oils, the water injection projects are relatively simple to establish and the water availability at a relatively low cost. For design of water injection projects is necessary to do reservoir studies in order to define the various parameters needed to increase the effectiveness of the method. For this kind of study can be used several mathematical models classified into two general categories: analytical or numerical. The present work aims to do a comparative analysis between the results presented by flow lines simulator and conventional finite differences simulator; both types of simulators are based on numerical methods designed to model light oil reservoirs subjected to water injection. Therefore, it was defined two reservoir models: the first one was a heterogeneous model whose petrophysical properties vary along the reservoir and the other one was created using average petrophysical properties obtained from the first model. Comparisons were done considering that the results of these two models were always in the same operational conditions. Then some rock and fluid parameters have been changed in both models and again the results were compared. From the factorial design, that was done to study the sensitivity analysis of reservoir parameters, a few cases were chosen to study the role of water injection rate and the vertical position of wells perforations in production forecast. It was observed that the results from the two simulators are quite similar in most of the cases; differences were found only in those cases where there was an increase in gas solubility ratio of the model. Thus, it was concluded that in flow simulation of reservoirs analogous of those now studied, mainly when the gas solubility ratio is low, the conventional finite differences simulator may be replaced by flow lines simulator the production forecast is compatible but the computational processing time is lower. / A inje??o de ?gua ? o m?todo de recupera??o suplementar mais difundido na maioria dos campos de petr?leo, devido a v?rios motivos, como o fato de a ?gua ser um eficiente agente deslocante de ?leos com baixa viscosidade associado ao fato de ser relativamente f?cil de injetar e de existir grande disponibilidade de ?gua a um custo relativamente baixo de se obter. Para o dimensionamento de um projeto de inje??o de ?gua ? necess?rio efetuar um estudo de reservat?rio com o objetivo de se definir os v?rios par?metros necess?rios para aumentar a efic?cia do m?todo. Para este estudo podem ser utilizados v?rios modelos matem?ticos classificados como anal?ticos ou num?ricos. Este trabalho tem como objetivo efetuar uma an?lise comparativa entre os resultados apresentados por um simulador de linhas de fluxo e um simulador convencional por diferen?a finita, que s?o dois tipos de simuladores de fluxo baseados em m?todos num?ricos, para um modelo de um reservat?rio de petr?leo leve submetido ? inje??o de ?gua. Para tanto foram montados dois modelos de reservat?rio, sendo um modelo heterog?neo e outro homog?neo contendo valores m?dios das propriedades de rocha do primeiro modelo. Foram realizadas compara??es dos resultados destes modelos considerando sempre as mesmas condi??es de opera??o. Em seguida foram alterados alguns par?metros de fluido e de rocha em ambos os modelos e efetuadas novamente compara??es de resultados. A partir do planejamento fatorial, realizado para a an?lise de sensibilidade dos par?metros de reservat?rio, foram escolhidos alguns casos para efetuar an?lise da varia??o dos par?metros operacionais escolhidos que foram a vaz?o de inje??o de ?gua e a localiza??o dos canhoneados. Observou-se que os resultados entre os simuladores s?o bastante similares na maioria dos casos estudados, tendo sido encontrado diverg?ncias apenas nos casos em que h? aumento na raz?o de solubilidade de g?s do modelo. Desta maneira foi conclu?do que o simulador convencional por diferen?as finitas pode ser substitu?do por um simulador por linhas de fluxo em reservat?rios semelhantes ao estudado, quando a raz?o de solubilidade possui valores baixos, sem preju?zo aos resultados de produ??o e com menores tempos de processamento computacional. / 2020-01-01

Made available in DSpace on 2014-12-17T15:01:45Z (GMT). No. of bitstreams: 1 JennysLMB.pdf: 3002293 bytes, checksum: 3b71773a14bee418fd80243a633d0cc8 (MD5) Previous issue date: 2008-02-22 / Steam assisted gravity drainage process (SAGD) involves two parallel horizontal wells located in a same vertical plane, where the top well is used as steam injector and the bottom well as producer. The dominant force in this process is gravitational. This improved oil recovery method has been demonstrated to be economically viable in commercial projects of oil recovery for heavy and extra heavy oil, but it is not yet implemented in Brazil. The study of this technology in reservoirs with characteristics of regional basins is necessary in order to analyze if this process can be used, minimizing the steam rate demand and improving the process profitability. In this study, a homogeneous reservoir was modeled with characteristics of Brazilian Northeast reservoirs. Simulations were accomplished with STARS , a commercial software from Computer Modelling Group, which is used to simulate improved oil recovery process in oil reservoirs. In this work, a steam optimization was accomplished in reservoirs with different physical characteristics and in different cases, through a technical-economic analysis. It was also studied a semi-continuous steam injection or with injection stops. Results showed that it is possible to use a simplified equation of the net present value, which incorporates earnings and expenses on oil production and expenses in steam requirement, in order to optimize steam rate and obtaining a higher net present value in the process. It was observed that SAGD process can be or not profitable depending on reservoirs characteristics. It was also obtained that steam demand can still be reduced injecting in a non continuous form, alternating steam injection with stops at several time intervals. The optimization of these intervals allowed to minimize heat losses and to improve oil recovery / O processo de drenagem gravitacional com inje??o cont?nua de vapor (SAGD) envolve dois po?os horizontais paralelos localizados em uma mesma vertical, onde o po?o superior ? usado como injetor de vapor e o inferior como produtor. A for?a dominante neste processo ? a gravitacional. Este m?todo de recupera??o avan?ada tem sido demonstrado ser economicamente vi?vel em projetos comerciais de recupera??o de petr?leo pesado e extra pesado, mas ainda n?o foi implementado no Brasil. O estudo desta tecnologia em reservat?rios com caracter?sticas das bacias regionais ? necess?rio para analisar como se ad?qua o processo para minimizar a demanda de vapor obtendo a maior rentabilidade do processo. Neste estudo foi usado um modelo homog?neo com caracter?sticas de reservat?rios do Nordeste Brasileiro. As simula??es foram realizadas em um programa comercial da Computer Modelling Group , o STARS , m?dulo usado para realizar estudos de m?todos de recupera??o avan?ada de reservat?rios de ?leo. Neste trabalho, foi realizada uma otimiza??o do vapor em reservat?rios com diferentes caracter?sticas f?sicas e em diferentes cen?rios, atrav?s de uma an?lise t?cnico-econ?mica. Tamb?m foi estudada a inje??o de vapor semi-cont?nua ou com paradas. Os resultados obtidos mostraram que ? poss?vel utilizar uma equa??o simplificada do valor presente l?quido, que incorpora os ganhos e gastos na produ??o de ?leo e os gastos na inje??o de vapor, para otimizar a demanda do vapor obtendo um maior valor presente l?quido no processo. Observou-se que o m?todo (SAGD) pode ser ou n?o rent?vel dependendo das caracter?sticas do reservat?rio. Encontrou-se tamb?m que a necessidade de vapor pode ainda ser diminu?da utilizando esquemas de inje??o de vapor com paradas em intervalos de tempo otimizados, e isto permitiu minimizar as perdas de calor e melhorar a recupera??o.

Orientadores: Denis Jose Schiozer, Celio Maschio / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica e Instituto de Geociencias / Made available in DSpace on 2018-08-13T03:31:41Z (GMT). No. of bitstreams: 1 Sousa_SergioHenriqueGuerrade_M.pdf: 1592526 bytes, checksum: a308d06cf11fb891b71f7b3396952356 (MD5) Previous issue date: 2007 / Resumo: O problema do ajuste de histórico é uma das tarefas que mais demandam tempo em um estudo de reservatório baseado em simulações de fluxo, porque é um problema inverso onde os resultados (dados de produção) são conhecidos, porém os valores de entrada (a caracterização do reservatório) não são integralmente conhecidos. Adicionalmente, as funções objetivo que medem a qualidade do ajuste costumam ser expressões compostas por uma série de componentes que tornam a topologia do espaço de soluções complexa e repleta de não linearidades. A metodologia adotada neste trabalho foi a modelagem do problema de ajuste de histórico como um problema de otimização combinatória de modo que ele pudesse ser abordado através de processos metaheurísticos. Em particular, a metaheurística Busca Dispersa (Scatter Search) foi acoplada a um algoritmo de Busca Direta baseado no método de Hooke e Jeeves para resolver o problema do ajuste de histórico. Reservatórios sintéticos de solução conhecida foram utilizados para fazer a validação da metodologia e, em seguida, ela foi aplicada a outro reservatório, também sintético, mas com características de reservatórios reais onde a solução do ajuste é desconhecida. São discutidos ao longo do texto o uso da metodologia de forma automática e assistida e também os benefícios do uso da computação distribuída na execução do método. As maiores contribuições deste trabalho em relação à questão do ajuste de histórico são: a introdução de uma nova metodologia versátil para uso automático ou assistido, a discussão de algumas características que dificultam o processo de ajuste e de que forma eles podem ser contornados e também a abordagem do tema do ajuste automático vs. o ajuste assistido ilustrado com exemplos. / Abstract: The history matching problem is one of the most demanding tasks in a reservoir simulation study; because it's an inverse problem where the results (production data) are known but the input data (the reservoir characterization data) are not entirely known. Moreover, the objective function that guide the match is usually made out of a series of components that make the topology of the objective function both complex and full of non-linearities. The methodology adopted in this work was to model the history matching problem as a combinatorial optimization problem in order for it to be solved by metaheuristic processes. In particular, the Scatter Search metaheuristic was coupled with a direct search method based on Hooke and Jeeve's method to solve the history matching problem. Synthetic reservoirs of known solutions where used to validate the methodology and then the methodology was applied to another reservoir, also synthetic, but with characteristics of real reservoirs where the solution is not known in advance. Throughout the text, the mixed use of the methodology on both an assisted and automatic fashion is discussed along with the benefits attained by the use of distributed computing resources. The greatest contributions of this work related to the history matching problem are: the introduction of a new versatile methodology for both automatic and assisted matches, the discussion of some characteristics that burden the entire process and some ways to overcome the difficulties, and also the discussion of some tradeoffs between automatic versus assisted history matching with examples to illustrate the matter. / Mestrado / Reservatórios e Gestão / Mestre em Ciências e Engenharia de Petróleo

Orientadores: Denis Jose Schiozer, Celio Maschio / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica e Instituto de Geociencias / Made available in DSpace on 2018-08-12T07:35:35Z (GMT). No. of bitstreams: 1 Becerra_GustavoGabriel_M.pdf: 16760750 bytes, checksum: 0609c24d13d46b9121f71356ce9d42a1 (MD5) Previous issue date: 2007 / Resumo: A escassez de informações de qualidade introduz risco ao processo de previsão da produção de petróleo tornando imprescindível o ajuste de histórico de produção, que é a calibração do modelo a partir da resposta produtiva registrada. O ajuste é um problema inverso, em que diferentes combinações dos valores dos parâmetros do reservatório podem conduzir a respostas aceitáveis, especialmente quando o grau de incerteza desses parâmetros é elevado. A integração do ajuste de histórico com a análise probabilística dos cenários representativos conduz à obtenção de uma metodologia para detecção dos modelos calibrados dentro de uma faixa de aceitaçãodefinida. O tratamento de atributos interdependentes de influência global e local e o avanço por etapas são necessários. Desta forma, o objetivo deste trabalho é apresentar uma metodologia que integra a análise de incertezas com o ajuste de histórico em modelos de reservatórios complexos. Este procedimento auxilia a detectar os atributos incertos críticos e sua possível variação com o intuito de estimar a faixa representativa das reservas a desenvolver. Não é alvo obter o melhor ajuste determinístico, mas refletir como o histórico possibilita uma mitigação das incertezas. Assim, a meta é usar modelos mais complexos e aprimorar a metodologia iniciada por Moura Filho (2006), desenvolvida para um modelo teórico simples. São utilizados dois casos de estudo de complexidade similar. Um deles referente ao reservatório do Campo de Namorado, utilizado para verificar e validar, em nível global, a aplicação da metodologia. Na etapa de aplicação, é usado um modelo sintético construído a partir de dados de afloramentos reais no Brasil e compreendendo informações de campos análogos com sistemas turbidíticos depositados em águas profundas. Os métodos aplicados, mediante a redefinição das probabilidades associadas e níveis dos atributos incertos, permitem: (1) reduzir a faixa de ajustes possíveis e obter modelos mais confiáveis; (2) identificar e condicionar à incerteza presente em função dos dados registrados; (3) diminuir os intervalos de incerteza dos parâmetros críticos identificados; (4) demarcar os limites seguros do desempenho futuro do reservatório. A conseqüência é um aumento da confiança no uso da simulação como ferramenta auxiliar do processo decisório. Além disso, procura-se fornecer à equipe multidisciplinar uma metodologia para reduzir o tempo empregado no gerenciamento de múltiplos atributos incertos na etapa de ajuste do modelo. / Abstract: The lack of reliable data or with high degree of uncertainty yields risk to the process of production prediction making the history matching, the model calibration from the registered field production indispensable. History matching is an inverse problem and, in general, different combinations of reservoir attributes can lead acceptable solutions, especially whit high degree of uncertainty of these attributes. The integration of history matching with a probabilistic analysis of representative models yields a way to detect matched models inside an acceptance interval, providing more efficient framework for predictions. It is necessary to consider dependences between global and local attributes. The scope of this work is to present a methodology that integrates the uncertainty analysis with the history matching process in complex models. This procedure helps to detect critical subsurface attributes and their possible variation, in order to estimate a representative range of the additional reserves to be developed. . It is not an objective to obtain the best deterministic model, but to mitigate uncertainties by using observed data. The objective is to improve the methodology initiated by Moura Filho (2006), applied to a simple model. The methodology presented in this work is applied in two study cases with similar complexity. Firstly, the methodology is verified and validated, on global scale, in Namorado Field. Then, at the application stage, it is chosen a synthetic reservoir model made from real outcrop data of Brazil and involving information from analog fields with turbiditic systems deposited in deep waters. The methodology allows the redefinition of the probability and levels of the dynamic and static attributes in order: (1) to reduce the group of possible history matching obtaining more realistic models; (2) to identify the existent uncertainty as a function of observed data; (3) to decrease the uncertainty range of critical reservoir parameters; (4) to increase the confidence in production forecast. One contribution of this work is to present a quantitative approach to increase the reliability on the use of reservoir simulation as an auxiliary tool in decision processes. Another purpose of this work is to provide a procedure to reduce the consumed time to handle multiples uncertainty attributes during the history matching. / Mestrado / Reservatórios e Gestão / Mestre em Ciências e Engenharia de Petróleo

Orientadores: Denis José Schiozer, Célio Maschio / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica e Instituto de Geociências / Made available in DSpace on 2018-08-23T20:38:54Z (GMT). No. of bitstreams: 1 Davolio_Alessandra_D.pdf: 20299340 bytes, checksum: 39bd0e5f79603d95c1cb4c74b2616103 (MD5) Previous issue date: 2013 / Resumo: O foco deste trabalho é usar dados de engenharia de reservatórios gerados através da simulação de fluxos para melhorar a interpretação quantitativa da sísmica 4D. A ideia é usar os conhecimentos de engenharia para minimizar possíveis informações incorretas geradas pela sísmica 4D, antes de usá-la para atualizar modelos de simulação de reservatórios em um procedimento de ajuste de histórico. Neste trabalho a integração entre estes dois conjuntos de dados é feita no domínio da pressão e saturação. Dessa forma, a primeira parte do trabalho apresenta uma metodologia para estimar variações de pressão e saturação a partir da sísmica 4D através de um procedimento de inversão petro-elástica. Este procedimento pode ser visto como uma ferramenta para aplicar a metodologia de integração da segunda parte, que é uma das principais contribuições deste trabalho. A metodologia de integração usa múltiplos modelos de simulação para restringir a estimativa de pressão e saturação a partir da sísmica 4D. Como resultado, observaram-se mapas menos ruidosos que permitem, portanto, uma melhor interpretação das variações ocorridas no reservatório. Seguindo a sequencia tradicional de trabalho, o mapa de saturação estimado foi então usado como dado de entrada em um processo de ajuste de histórico apresentado na terceira parte da tese. Foi mostrado que o procedimento de ajuste gera melhores resultados quando o dado de entrada, neste caso mapa de saturação, respeita o balanço de massa, o que não acontecia para os dados considerados. Assim, uma metodologia para calibrar o volume de água injetada associado ao mapa de saturação gerado pela sísmica 4D é apresentada na última parte do trabalho, que é outra contribuição importante a se destacar. Para melhor controlar os resultados de todas as metodologias aqui apresentadas foram usados dados sintéticos em todo o trabalho. Embora resultados satisfatórios tenham sido observados para este conjunto de dados, é importante destacar que as principais contribuições deste trabalho não são apenas os resultados observados, mas as metodologias propostas, que apresentam uma perspectiva inovadora para integração entre dados de sísmica 4D e engenharia de reservatórios / Abstract: The focus of this work is to use reservoir engineering data from numerical flow simulation to improve the quantitative interpretation of 4D seismic signals. The idea is to use engineering knowledge to minimize possible incorrect information provided by 4D seismic before using it to update reservoir simulation models in a history matching procedure. In this work the integration between the two dataset is done in the pressure and saturation domain. So, the first part of the work presents a methodology to estimate pressure and saturation changes from 4D seismic through a petro-elastic inversion procedure. This procedure can be seen as a tool to apply the integration methodology of the second part which is one of the main contributions of this work. The integration methodology uses multiple simulation models to constrain the estimation of pressure and saturation from 4D seismic. As a result, less noisy maps were obtained, allowing a better interpretation of the reservoir changes. Following the traditional sequence, the estimated saturation map is then used as input to the history matching process presented in the third part of the work. It was shown that the history matching procedure provides better results if the input data, in this case the saturation map, respects the expected mass balance, which was not the case for the dataset considered. Thus, a methodology to calibrate the volume of injected water associated to the saturation map provided by 4D seismic is presented in the last part of the work, which is another important contribution to be highlighted. To better control the results of all the methodologies here presented, a synthetic dataset was used in the entire work. Although satisfactory results were observed for this dataset, it is important to highlight that the main contributions of this work are not only the results, but the methodologies proposed, that present an innovative perspective for 4D seismic and reservoir engineering data integration / Doutorado / Reservatórios e Gestão / Doutora em Ciências e Engenharia de Petróleo

Orientador: Osvair Vidal Trevisa / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica e Instituto de Geociências / Made available in DSpace on 2018-08-28T00:16:54Z (GMT). No. of bitstreams: 1 MercadoSierra_DianaPatricia_D.pdf: 5821830 bytes, checksum: 86d0611dd821cb67544d7463dfed3a39 (MD5) Previous issue date: 2015 / Resumo: A combustão in-situ é um processo multiescala, multifísico que envolve simultaneamente o escoamento de fluidos no meio poroso, o equilíbrio de fases e a cinética das reações químicas. A simulação desse processo tem alcançado um elevado grau de desenvolvimento, no entanto, mecanismos básicos ainda são representados de maneira incompleta, impondo inúmeros desafios na modelagem. A dificuldade de modelar fenômenos relacionados com a combustão tem a ver com a representação do efeito da frente de combustão e a modelagem do consumo de combustível. Na combustão in-situ as reações químicas acontecem em uma zona delgada de menos de um metro de espessura, que é pequena quando comparada com a escala do reservatório de centenas ou milhares de metros. Na simulação na escala de campo, o uso de células de tamanho maior do que a zona de reação leva a erros na distribuição da temperatura. Consequentemente, a velocidade das reações não pode ser bem representada. De outro lado, os simuladores não permitem controlar a ocorrência das reações a partir da energia de ativação. Como resultado, o início das reações se torna independente da temperatura. O objetivo desta tese é desenvolver um modelo pseudocinético para a simulação numérica da combustão in-situ na escala de campo. Com o modelo pseudocinético pretende-se representar os fenômenos na zona de combustão, reduzindo o efeito do tamanho de célula. O trabalho foi desenvolvido em etapas. Primeiro foram estabelecidas as condições que o simulador deveria atender e definida a estratégia de abordagem, que foi a de desenvolver um modelo pseudocinético. Depois foi definida a metodologia de obtenção do modelo pseudocinético. Após o modelo pseudocinético concluído, este foi utilizado para a simulação de um campo de óleo pesado brasileiro submetido à combustão in-situ. O modelo pseudocinético proposto consiste em expressar a energia de ativação das reações em função da temperatura. Através do modelo, é possível restringir a ocorrência da reação de craqueamento, de modo que o início da formação do coque aconteça somente para temperaturas acima dos valores observados na zona de craqueamento. Note-se que neste cenário a quantidade de coque depositado pode ser modelada usando a reação de craqueamento, o que se constitui numa das principais contribuições do trabalho. O modelo permite manter a dependência da taxa de reação com a temperatura mediante o uso de valores de energia de ativação apropriados. Além disso, consegue-se reduzir o efeito da distribuição de temperatura mediante o controle da taxa de reação em função dos valores médios de temperatura observados nas células do modelo de simulação na escala de campo. Na simulação do piloto de combustão in-situ, o modelo pseudocinético foi obtido do ajuste progressivo dos parâmetros cinéticos das reações químicas, partindo da simulação do processo na escala de laboratório até a escala de campo. Os dados experimentais utilizados na simulação na escala de laboratório foram obtidos de um ensaio em tubo de combustão seca realizado no Laboratório de Métodos Térmicos de Recuperação do Departamento de Energia da UNICAMP. O fluido utilizado foi um óleo pesado de 15,3 °API proveniente da Bacia do Espírito Santo / Abstract: The in-situ combustion is a multi-scale, multi-physics process, involving fluid flow in porous media, thermodynamic equilibrium of the phases involved and chemical kinetics of reactions. The simulation of this process has achieved a high degree of development, however basic mechanisms are still represented incompletely, imposing numerous challenges in modeling. The issues in the combustion modeling are related with the representation of the combustion front effect and the fuel consumption modeling. Chemical reactions of the in-situ combustion process take place in a thin zone of less than a meter thick, which is small compared to the field scale of hundreds or thousands of meters. Numerical simulations at the field scale typically use grid blocks that are at least two orders of magnitude greater than that. Such divergence leads to improper representations of key aspects of the process, as the temperature distribution and the reaction kinetics. In accordance with that the reaction occurrence is not controlled by the activation energy in the simulation models. The major shortcome is on fuel deposition, a key issue in in-situ combustion, which will happen from the start, since the cracking reaction may proceed even at reservoir temperature. The objective of this thesis is to develop a new pseudokinetic model for field-scale simulation of in-situ combustion. With the pseudokinetic model meant to improve the representation of the combustion zone effects reducing the gridblock size effect. The work was carried out in stages. First establishes the conditions that the simulator should meet and defined the strategy to develop a pseudokinetic model. Then a methodology was defined for obtaining the pseudokinetic model. After the pseudokinetic model is completed, it is applied to the in-situ combustion modeling of a Brazilian heavy oil field. The models pursue the idea of making the activation energy a function of the grid block temperature. The model allows restricting the cracking reaction occurrence by the temperature, so that the beginning of the coke deposition occurs at temperatures greater than the temperature observed in the cracking zone. Note that in this scenario the cracking reaction can be used to represent the coke deposition, which constitutes one of the main contributions of this work. The model allows maintaining the dependence of reaction rate with temperature through the use of appropriate activation energy values. Furthermore, the model reduces the temperature distribution effect by controlling the reaction rate based on average temperature values observed in the field simulation model. In the simulation of the in-situ combustion pilot, the pseudokinetic model was obtained from the progressive tuning of the kinetic parameters of chemical reactions, based on the simulation of the process from the laboratory to field scale. The experimental data used in the laboratory scale simulation were obtained from a dry combustion tube test carried out at the Thermal Recovery Methods Laboratory of the Energy Department at UNICAMP. The fluid used was a 15.3 ° API heavy oil from the Espírito Santo Basin / Doutorado / Reservatórios e Gestão / Doutora em Ciências e Engenharia de Petróleo

Unlike the traditional two-stage methods, a conditional and inverse-conditional simulation approach may directly generate independent, identically distributed realizations to honor both static data and state data in one step. The Markov chain Monte Carlo (McMC) method was proved a powerful tool to perform such type of stochastic simulation. One of the main advantages of the McMC over the traditional sensitivity-based optimization methods to inverse problems is its power, flexibility and well-posedness in incorporating observation data from different sources. In this work, an improved version of the McMC method is presented to perform the stochastic simulation of reservoirs and aquifers in the framework of multi-Gaussian geostatistics. First, a blocking scheme is proposed to overcome the limitations of the classic single-component Metropolis-Hastings-type McMC. One of the main characteristics of the blocking McMC (BMcMC) scheme is that, depending on the inconsistence between the prior model and the reality, it can preserve the prior spatial structure and statistics as users specified. At the same time, it improves the mixing of the Markov chain and hence enhances the computational efficiency of the McMC. Furthermore, the exploration ability and the mixing speed of McMC are efficiently improved by coupling the multiscale proposals, i.e., the coupled multiscale McMC method. In order to make the BMcMC method capable of dealing with the high-dimensional cases, a multi-scale scheme is introduced to accelerate the computation of the likelihood which greatly improves the computational efficiency of the McMC due to the fact that most of the computational efforts are spent on the forward simulations. To this end, a flexible-grid full-tensor finite-difference simulator, which is widely compatible with the outputs from various upscaling subroutines, is developed to solve the flow equations and a constant-displacement random-walk particle-tracking method, which enhances the com / Fu, J. (2008). A markov chain monte carlo method for inverse stochastic modeling and uncertainty assessment [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/1969 / Palancia

En simulation dynamique des réservoirs, un des artéfacts les plus gênants pour la prédiction de production est l’effet de l’orientation du maillage. Bien que celui-ci soit « normal » pour tout schéma numérique, il se trouve amplifié par l’instabilité du modèle physique, ce qui a lieu lorsque le contraste de mobilités entre l’eau (fluide poussant, utilisé dans les procédés de récupération secondaires) et l’huile (fluide poussé, contenant les hydrocarbures) dépasse un certain seuil critique. On parle alors d’écoulements à rapport de mobilités défavorable. Connu depuis longtemps, ce problème a fait l'objet de nombreux travaux dans les années 1980 ayant abouti au schéma dit à neuf points. Actuellement implanté dans PumaFlow, logiciel développé et commercialisé par IFPEN, ce schéma fonctionne relativement bien en maillages carrés et dépend d’un paramètre scalaire dont le réglage varie selon les auteurs sur la base de considérations heuristiques. Dans cette thèse, nous proposons une nouvelle démarche méthodologique afin non seulement d’ajuster ce paramètre libre de manière optimale mais aussi de généraliser le schéma aux maillages rectangulaires. La stratégie que nous préconisons repose sur une analyse d’erreur du problème, à partir de laquelle il est possible de définir une notion d’erreur angulaire et de garantir que le comportement du schéma obtenu soit le « moins anisotrope » possible via une minimisation de son écart par rapport à un comportement idéal. Cette procédure de minimisation est ensuite appliquée à deux autres familles de schémas numériques~ : (1) un schéma multidimensionnel proposé par Kozdon, dans lequel le paramètre libre est une fonction~ ; (2) un autre schéma à neuf points faisant intervenir deux paramètres scalaires. C’est ce dernier qui réduit le mieux l’effet de l’orientation lorsque le rapport des pas de maillage s’éloigne de 1. Enfin, une extension de la méthode à des modèles physiques plus complets est envisagée. / In dynamic reservoir simulation, one of the most troublesome artifacts for the prediction of production is the grid orientation effect. Although this normally arises from any numerical scheme, it happens to be amplified by the instability of the physical model, which occurs when the mobility contrast between the water (pushing fluid, used in the processes of secondary recovery) and the oil (pushed fluid, containing the hydrocarbons) exceeds a some critical threshold. We then speak of flows with adverse mobility ratio. This GOE issue has received a lot of attention from the engineers. Numerous works dating back to the 1980s have resulted in the so-called nine-point scheme. Currently implemented in the IFPEN software PumaFlow, this scheme performs relatively well in square meshes and depends on a scalar parameter whose value varies from one author to another, on the grounds of heuristic considerations. In this thesis, we propose a new methodological approach in order not only to optimally adjust this free parameter, but also to extend the scheme to rectangular meshes. The strategy that we advocate is based on an error analysis of the problem, from which it is possible to define a notion of angular error and to guarantee that the behavior of the obtained scheme is the "least anisotropic" possible through a minimization of its deviation from some ideal behavior. This minimization procedure is then applied to two other families of numerical schemes: (1) a multidimensional scheme proposed by Kozdon, in which the free parameter is a function; (2) another nine-point scheme involving two scalar parameters. The latter provides the best results regarding GOE reduction when the ratio of the mesh steps is far away from 1. Finally, an extension of the method to more sophisticated physical models is envisaged.

[pt] Este trabalho tem por objetivo investigar a redução do custo computacional associado ao cálculo das principais estatísticas das saídas dos modelos de propagação de incertezas. Para tal, apresentamos uma implementação alternativa ao método tradicional de Monte Carlo, chamado Caos Polinomial; que é adequado a problemas onde o número de variáveis de incerteza não é muito alto. No método Caos Polinomial, o valor esperado e a variância das saídas do simulador são diretamente estimados, como funções de distribuições de probabilidade de variáveis de incerteza na entrada do simulador. A principal vantagem do método de Caos Polinomial é que o número de pontos necessários para uma boa estimativa das estatísticas da saída de um simulador, comparado com Monte Carlo, é menor. Aplicações de Caos Polinomial em reservatórios de petróleo serão apresentadas para a propagação de até quatro variáveis, apesar do método poder ser aplicado a problemas de dimensões maiores. Nossos principais resultados são aplicados a dois modelos de reservatórios de petróleo sintéticos. / [en] In this work we investigate the reduction of the computational cost of the calculus of statistical moments of simulator s output in uncertainties propagation s models. For do that, we present an alternative s implementation to the traditional Monte Carlo s Method, called Polynomial Chaos; that is adequate to problems where the number of uncertain variables is not so high. In the Polynomial Chaos method, the expectation and the variance of the simulator s output are directly estimated, as functions of the probability distribuition of the uncertain variables in simulator input. The great advantage of Polynomial Chaos is that number of points necessary for a good estimation of the output statistics have smaller magnitude, compared to the Monte Carlo Method. Applications of Polynomial Chaos on oil reservoir simulations will be presented. As it is just a preliminar implementation, we just treat propagation s problems with at most four uncertainties variables, despite of the method being applicable to problems with more dimensions. Our main results are applied to two models of synthetic oil reservoirs.

Les réservoirs de type "tight oil" et "shale gaz" ont une partie importante de leur volume poreux occupée par des micropores (< 2nm) et des mesopores (entre 2 et 50 nm). Ce type d'environnement crée de fortes forces d’interaction dans le fluide confiné avec les parois du pores et entre ses propres molécules, ce qui change fortement la thermodynamique du fluide. Un travail important doit donc être effectué sur le développement de méthodes de mise à l'échelle de la distribution de pore pour effectuer des simulations réservoir à grande échelle. Premièrement, des simulations moléculaires sont effectuées sur des fluides confinés afin d'obtenir des propriétés thermodynamiques de référence à l'équilibre liquide/vapeur pour différentes tailles de pore. Ensuite, une comparaison des données de simulation moléculaire avec les résultats issus des équation d'état utilisées dans la littérature a permis de mettre en valeur la méthode de flash avec pression capillaire et changement du point critique comme la meilleure méthode existante pour décrire la physique du fluide confiné. Des simulations fines d'écoulement matrice/fracture ont donc été effectuées pour différentes tailles de pore. Des modèles de mise à l'échelle en maillage grossier ont été ensuite construits à partir du même cas synthétique et les résultats ont été comparés avec ceux des simulations de référence en maillage fin. Un nouveau modèle de triple porosité considérant fracture, petit pores et grand pores avec une approche MINC a donné des résultats très proches du maillage fin. Finalement un réservoir stimulé hydrauliquement à grande échelle a été simulé pour différentes distributions de pores avec le modèle développé. / Tight oil and shale gas reservoirs have a significant part of their pore volume occupied by micro (below 2nm) and mesopores (between 2 and 50nm). This kind of environment creates strong interaction forces in the confined fluid with pore walls as well as between its own molecules and then changes dramatically the fluid phase behavior. An important work has therefore to be done on developing upscaling methodology of the pore size distribution for large scale reservoir simulations. Firstly, molecular simulations are performed on different confined fluids in order to get reference thermodynamic properties at liquid/vapor equilibrium for different pore sizes. Then, the comparison with commonly used modified equation of state (EOS) in the literature highlighted the model of flash with capillary pressure and critical temperature and pressure shift as the best one to match reference molecular simulation results. Afterwards fine grid matrix/fracture simulations have been built and performed for different pore size distributions. Then, coarse grid upscaling models have then been performed on the same synthetic case and compared to the reference fine grid results. A new triple porosity model considering fracture, small pores and large pores with MINC (Multiple Interacting Continua) approach, has shown very good match with the reference fine grid results. Finally a large scale stimulated reservoir volume with different pore size distribution inside the matrix has been built using the upscaling method developed here.