A large percentage of oil and gas reservoirs in the most productive regions such as the Middle East, South America, and Southeast Asia are naturally fractured reservoirs (NFR). The major difference between conventional reservoirs and naturally fractured reservoirs is the discontinuity in media in fractured reservoir due to tectonic activities. These discontinuities cause remarkable difficulties in describing the petrophysical structures and the flow of fluids in the fractured reservoirs. Predicting fluid flow behavior in naturally fractured reservoirs is a challenging area in petroleum engineering. Two classes of models used to describe flow and transport phenomena in fracture reservoirs are discrete and continuum (i.e. dual porosity) models. The discrete model is appealing from a modeling point of view, but the huge computational demand and burden of porting the fractures into the computational grid are its shortcomings. The affect of natural fractures on the permeability anisotropy can be determined by considering distribution and orientation of fractures. Representative fracture permeability, which is a crucial step in the reservoir simulation study, must be calculated based on fracture characteristics. The diagonal representation of permeability, which is customarily used in a dual porosity model, is valid only for the cases where fractures are parallel to one of the principal axes. This assumption cannot adequately describe flow characteristics where there is variation in fracture spacing, length, and orientation. To overcome this shortcoming, the principle of the full permeability tensor in the discrete fracture network can be incorporated into the dual porosity model. Hence, the dual porosity model can retain the real fracture system characteristics. This study was designed to develop a novel approach to integrate dual porosity model and full permeability tensor representation in fractures. A fully implicit, parallel, compositional chemical dual porosity simulator for modeling naturally fractured reservoirs has been developed. The model is capable of simulating large-scale chemical flooding processes. Accurate representation of the fluid exchange between the matrix and fracture and precise representation of the fracture system as an equivalent porous media are the key parameters in utilizing of dual porosity models. The matrix blocks are discretized into both rectangular rings and vertical layers to offer a better resolution of transient flow. The developed model was successfully verified against a chemical flooding simulator called UTCHEM. Results show excellent agreements for a variety of flooding processes. The developed dual porosity model has further been improved by implementing a full permeability tensor representation of fractures. The full permeability feature in the fracture system of a dual porosity model adequately captures the system directionality and heterogeneity. At the same time, the powerful dual porosity concept is inherited. The implementation has been verified by studying water and chemical flooding in cylindrical and spherical reservoirs. It has also been verified against ECLIPSE and FracMan commercial simulators. This study leads to a conclusion that the full permeability tensor representation is essential to accurately simulate fluid flow in heterogeneous and anisotropic fracture systems. / text

The transport of components and phases plays a fundamental role in the success of an EOR process. Because many reservoirs have harsh conditions of salinity, temperature and rock heterogeneity, which limit process options, a robust system with flexibility is required. Systematic experimental study of formulations capable to transport surfactant as foam at 94°C, formulated in sea water, is presented. It includes methodology to conduct core floods in sand packs using foaming surfactants and to develop “surfactant blend ratio- salinity ratio maps” using equilibrium phase behavior to determine favorable conditions for oil recovery in such floods. Mathematical model able to reproduce the foam strength behavior observed in sand packs with the formulations studied is presented. Visualization of oil recovery mechanism from matrix is realized using a model system of micro-channels surrounded by glass beads to mimic matrix and fractures respectively. The observations illustrate how components may distribute within the matrix, thereby releasing oil into the fractures. The use of chemicals to minimize adsorption is required when surfactant adsorption is important. The presence of anhydrite may limit the use of sodium carbonate to reduce adsorption of carbonates. A methodology is presented to estimate the amount, if any, of anhydrite present in the reservoir. The method is based on brine software analysis of produced water compositions and inductively coupled plasma (ICP) analysis of core samples. X-ray powder diffraction (XRD) was used to verify the mineralogy of the rock. X-ray photoelectron spectroscopy (XPS) was used to obtain surface composition for comparison with bulk composition of the rock. Adsorption of surfactants was measured using dynamic and static adsorption experiments. Determining the flow properties of the rock samples via tracer analysis permitted the simulation of the dynamic adsorption process using a mathematical model that considers the distribution of adsorbed materials in the three different regions of pore space. Using this method allows one to predict adsorption in a reservoir via simulation.

Störungszonen in den karbonatischen Wechselfolgen des Muschelkalks (Mittlere Trias) sind potenzielle Ziele für hydrothermale Projekte im Oberrheingraben (ORG). Es wurden verschiedene Methoden miteinander kombiniert, um die Permeabilitäts-Strukturen solcher Störungszonen, deren Spannungszustände und lokale Spannungsfelder abzuschätzen. Diese Arbeit kann damit zur Exploration von störungsgebundenen Muschelkalk-Reservoiren im ORG beitragen. Sie vergleicht diese mit einem bereits erfolgreich getesteten (störungsgebundenen) hydrothermalen Reservoir, den Sandsteinen des Buntsandsteins. Um für den Muschelkalk die Charakteristika von Störungszonen zu definieren und die zugehörige Permeabilitätsstruktur abzuschätzen, wurden verschiedene Störungszonentypen (z.B. Abschiebungen, invertierte Störungszonen, Schrägabschiebung) mit unterschiedlichen Maßstäben (Versatz: mittel-skalig 1-10 m, groß-skalig >10 m) detailliert untersucht. Eine Aufschlussanalogstudie zu einer groß-skaligen Störungszone (Schrägabschiebung) in dickbankigen Sandsteinen (Buntsandstein, Untere Trias) wurde zum Vergleich herangezogen. Der besondere Schwerpunkt lag jeweils auf der Charakterisierung von Bruchzone und Störungskern sowie des assoziierten Bruchsystems (Bruchdichte, Öffnungsweitenverteilung, Vernetzungsgrad, vertikale Ausdehnung). Die Daten zeigen, dass Bruchsysteme in eher homogenen Einheiten einen positiven Effekt auf die Reservoir-Permeabilität haben. Sie bieten, vor allem in der Nähe des Störungskerns, potenzielle Fließwege auch über mehrere Schichten und zeigen eine starke Vernetzung zwischen vergleichsweise kurzen Brüchen. Im Gegensatz dazu scheinen störungsgebundene Bruchsysteme in Einheiten mit einer starken mechanischen Schichtung einen geringeren Einfluss auf die Reservoirpermeabilität zu haben. Störungskerne zeigen stellenweise eine signifikante Komplexität, da sie vor allem abdichtende aber stellenweise auch durchlässige Strukturen aufweisen. Groß-skalige Störungszonen (im Reservoir-Maßstab) lassen sich für beide potenzielle Reservoirhorizonte am besten als kombinierte Barriere-Leiter-Systeme beschreiben. Diese Barriere-Leiter-Systeme zeichnen sich durch ein potenziell hydraulisch durchlässiges Bruchsystem in der Bruchzone (und im Buntsandstein zusätzlich in der Übergangszone) sowie einen schwach durchlässigen bis abdichtenden Störungskern aus. Um die Kenntnisse zum lokalen Spannungsfeld in Störungszonen im geschichteten Muschelkalk (Reservoirtiefe: 2.900 m) zu verbessern, wurden mit der Finite-Elemente-Software COMSOL Multiphysics® numerische 3D-Modelle erstellt. Es wurden deutliche Unterschiede des lokalen Spannungsfelds in Abhängigkeit von (1) der Orientierung, (2) dem Einfluss der maximalen Horizontalspannung SH im Spannungsregime, (3) Störungszonen-Maßstab und (4) den Kontrasten der mechanischen Eigenschaften festgestellt. In Spannungsregimen mit starker horizontaler Kompression wurde für Spannungsmagnituden und Versatz eine deutliche Abhängigkeit von der Orientierung der Störungszone festgestellt. Vor allem groß-skalige Störungszonen mit einem 30°-Winkel zu SH scheinen einen SH-induzierten Horizontalversatz innerhalb weicher Störungszoneneinheiten zu begünstigen; in dieser wurden die höchsten Versatzbeträge festgestellt. Die typische Abnahme von Spannungsmagnituden in weicheren Störungszoneneinheiten verringert sich in Richtung von Störungszonen, die senkrecht zu SH streichen. Der Einfluss der mechanischen Schichtung steigt mit zunehmender horizontaler Kompression, was zu einem vertikal heterogenen Spannungsfeld führt. Um Annahmen zu einer möglichen hydraulischen Aktivität einer Störungszone zu treffen, werden für die untersuchten Störungszonen analytische Abschätzungen zu Bewegungs- und Dehnungstendenzen unter Reservoir-Bedingungen präsentiert. Die Ergebnisse zeigen unterschiedliche Spannungszustände der analysierten Störungszonen, die auf das rezente Übergangsregime sowie unterschiedliche Orientierungen der maximalen Horizontalspannungen SH zurückzuführen sind. In dieser Arbeit präsentierte Ergebnisse von Aufschlussanalogstudien helfen dabei, fundierte Annahmen zur Permeabilitätsstruktur von Störungszonen zu treffen und damit vielversprechende Bohrziele im ORG zu definieren. In diesem Zusammenhang konnten sedimentäre Wechselfolgen identifiziert werden, die als potenzielle geothermische Reservoire ausgeschlossen werden können. Die Ergebnisse der numerischen 3D-Modellierungen können dazu beitragen, möglicherweise nötige Stimulationsmaßnahmen im Muschelkalk Reservoir optimal auszulegen. Darüber hinaus bieten die Ergebnisse die Möglichkeit, Einblick in potenzielle Probleme während der Bohrungsherstellung in Muschelkalk-Reservoiren zu bekommen, wie zum Beispiel die Wahrscheinlichkeit eines vertikal heterogenen Spannungsfeldes.

Orientador: Osvaldo Luís Manzoli / Abstract: This work presents a numerical method based on Discrete Fracture Model (DFM) and the Finite Element Method (FEM), where the fractures are approximated by a reduced model. The flow along and across the fracture is described by a simplified set of equations considering both conductive fractures and barriers. The coupled hydromechanical model is composed of a linear poroelastic Biot medium and a nonlinear model based on damage mechanics for the fractures, which captures the nonlinear normal deformation and shear dilation according to the Barton-Bandis model. Both flow and geomechanical models are approximated using the finite element model. Fractures are explicitly represented by three-node standard finite elements with high aspect ratio (i.e. ratio between the largest and the smallest element dimensions) and appropriate constitutive laws. These interface high aspect ratio elements represent a regularization method which continuously approximate the discontinuous pressure and displacement fields on a narrow material band around the fracture. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. The efficiency of the proposed method is demonstrated through numerical examples, as well as the effects of fractures in the hydraulic properties of porous rocks and its dependency of the stress state. / Resumo: Este trabalho apresenta um método numérico baseado no Modelo de Fratura Discreta (MFD) e no Método dos Elementos Finitos (MEF), onde as fraturas são aproximadas por um modelo reduzido. O fluxo ao longo e através da fratura é descrito por um conjunto simplificado de equações, considerando tanto fraturas condutoras quanto barreiras. O modelo hidromecânico acoplado é composto por um meio poroelástico linear e um modelo não linear para fraturas, baseado na mecânica do dano e que captura a deformação normal não linear e a dilatância ao cisalhamento de acordo com o modelo de Barton-Bandis. Os modelos de fluxo e geomecânico são aproximados usando o método dos elementos finitos. As fraturas são explicitamente representadas por elementos finitos triangulares de três nós com elevada razão de aspecto (isto é, a razão entre a maior e a menor dimensão do elemento) e leis constitutivas apropriadas. Esses elementos de elevada razão de aspecto representam um método de regularização que aproxima de forma contínua os campos de pressão e deslocamento descontínuos em uma estreita faixa material ao redor da fratura. A formulação matemática completa é apresentada juntamente com o algoritmo sugerido para sua implementação numérica. A eficiência do método proposto é demonstrada através de exemplos numéricos, bem como os efeitos de fraturas nas propriedades hidráulicas de rochas porosas e sua dependência do estado de tensão. / Mestre

In geothermal reservoirs and unconventional gas reservoirs with very low matrix permeability, fractures are the main routes of fluid flow and heat transport, so the fracture permeability change is important. In fact, reservoir development under this circumstance relies on generation and stimulation of a fracture network. This thesis presents numerical simulation of the response of a fractured rock to injection and extraction considering the role of poro-thermoelasticity and joint deformation. Fluid flow and heat transport in the fracture are treated using a finite difference method while the fracture and rock matrix deformation are determined using the displacement discontinuity method (DDM). The fractures response to fluid injection and extraction is affected both by the induced stresses as well as by the initial far-field stress. The latter is accounted for using the non-equilibrium condition, i.e., relaxing the assumption that the rock joints are in equilibrium with the in-situ stress state. The fully coupled DDM simulation has been used to carry out several case studies to model the fracture response under different injection/extractions, in-situ stresses, joint geometries and properties, for both equilibrium and non-equilibrium conditions. The following observations are made: i) Fluid injection increases the pressure causing the joint to open. For non-isothermal injection, cooling increases the fracture aperture drastically by inducing tensile stresses. Higher fracture aperture means higher conductivity. ii) In a single fracture under constant anisotropic in-situ stress (non-equilibrium condition), permanent shear slip is encountered on all fracture segments when the shear strength is overcome by shear stress in response to fluid injection. With cooling operation, the fracture segments in the vicinity of the injection point are opened due to cooling-induced tensile stress and injection pressure, and all the fracture segments experience slip. iii) Fluid pressure in fractures increases in response to compression. The fluid compressibility and joint stiffness play a role. iv) When there are injection and extraction in fractured reservoirs, the cooler fluid flows through the fracture channels from the injection point to extraction well extracting heat from the warmer reservoir matrix. As the matrix cools, the resulting thermal stress increases the fracture apertures and thus increases the fracture conductivity. v) Injection decreases the amount of effective stress due to pressure increase in fracture and matrix near a well. In contrast, extraction increases the amount of effective stress due to pressure drop in fracture and matrix.

Le Buntsandstein, localisé dans le graben du Rhin supérieur, apparait comme une cible intéressante pour la géothermie, associant une formation argilo-gréseuse à l’anomalie thermique régionale. Cette étude vise à caractériser les propriétés pétrophysiques de ces grès ainsi que la fracturation les affectant, dans le but de fournir un modèle conceptuel de la formation qui servira de guide pour son exploitation futur. Les faciès sédimentaires sont composés par cinq faciès pétrographiques (grès propres, grès à enrobage argileux, grès à matrice argileuse, grès silicifiés et grès à ciment carbonaté), qui se répartissent dans des proportions variables et contrôlent une partie des propriétés pétrophysiques mesurées à l’échelle matricielle. La comparaison des données pétrophysiques, des données macroscopiques issues d’une analyse de gradients de température, des données de modélisation et de la fracturation permet de construire un modèle de circulation dans le réservoir. Ces analyses mettent en avant le rôle de la zone endommagée des zones de faille pour le transfert de fluides à grande échelle, mais aussi celui de deux faciès sédimentaires : les grès déposés dans un environnement de Playa Lake et fluvio-éolien. L’analyse de différents affleurements montre que la fracturation évolue en fonction de la situation dans la pile sédimentaire et en fonction de la situation par rapport aux accidents tectoniques majeures. / The Buntsandstein, located in the Upper Rhine Graben, appears to be an easy target forgeothermal exploitation, linking sandstone and clay with the regional thermal anomaly. This study aims at characterizing petrophysical characteristics of these sandstones as well as the fracturation affecting them, with the intention of providing a conceptual model of the formation which will act as guide for future exploitation. The sedimentary facies are composed by five petrographical facies (clean sandstones, sandstones with clayed coating, clay matrix sandstones, silicified sandstones and carbonated matrix sandstones) which split with variable proportions and control a part of petrophysical properties measured at matrix scale. The comparison between petrophysical data, macroscopic data from temperature gradient analysis, modelling data and fracturing, allows the building of a Buntsandstein Sandstones fluids circulation conceptual model. This analysis points the role of the damage zone of fault zones for fluids transfer at large scale, but also that of two sedimentary facies: marginal erg and Playa Lake. The analysis of different outcrops shows that the fracturation evolves according to the situation in the sedimentary pile and according to the situation in comparison with major tectonic accidents.

Orientador: Osvair Vidal Trevisan / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica, Instituto de Geociencias / Made available in DSpace on 2018-08-12T20:49:17Z (GMT). No. of bitstreams: 1 MateoHernandez_JuanAlberto_M.pdf: 12401820 bytes, checksum: 705c41d03544d222731ce816469f47c2 (MD5) Previous issue date: 2006 / Resumo: Neste trabalho são investigados os impactos individual e coletivo dos mecanismos de gás em solução, geração de CO2, destilação, embebição capilar e drenagem gravitacional, sobre a recuperação de óleo e gás, durante a injeção continua de vapor num reservatório naturalmente fraturado contendo óleo pesado. A investigação é feita através de simulação numérica dos fenômenos em modelos padrões de reservatórios. Dois modelos numéricos semelhantes são usados para representar o processo de aquecimento da matriz. O primeiro descreve o aquecimento de uma seção horizontal bidimensional de um bloco da matriz circundado por uma fratura na qual circula vapor. O segundo modelo descreve o aquecimento de um bloco de matriz semelhantemente circundado por uma fratura em que circula vapor, porém na direção vertical, visando agregar o efeito da ação da gravidade. Os estudos foram conduzidos para rochas saturadas com óleo vivo. As propriedades da rocha são as de um reservatório carbonático fraturado real e as propriedades dos fluidos se referem também ao mesmo caso real. Alem disso, as condições operacionais adotadas de pressão e temperatura são as observadas no campo, tornando o estudo e suas conclusões como próprias de um estudo de caso. Os resultados mostram que os principais mecanismos de recuperação de óleo da matriz durante o intervalo de aquecimento de 10 anos, foram os mecanismos de gás em solução e de destilação por arraste de vapor. Este último é o mecanismo de maior importância e é responsável pelo melhoramento da qualidade do óleo produzido / Abstract: In this work, the individual and collective impacts of the mechanisms solution gas drive, CO2 generation, steam distillation, capillary imbibition and drainage gravitational, on the oil and gas recovery, were investigated during the steamflooding of a naturally fractured reservoir containing heavy oil. The investigation was performed for standard reservoir models through numeric simulation. Two similar numerical models represent the matrix heating process. The first describes the heating of a horizontal cross-section of a matrix block surrounded by a fracture, in which the steam is flooding. The second model describes the same method of matrix heating, which was represented in the first model, but in the vertical direction, investigating the action of gravity. The studies were performed for a rock saturated with live oil. The rock properties are the same of a real fractured carbonate reservoir and the fluid properties also refer to the same real case. In addition, the adopted field operational parameters (pressure and temperature) refer to field conditions, turning the study and its conclusions as proper of a case study. The results show that the main mechanisms of oil recovery for the matrix block during the heating interval of 10 years were the integrated action of solution gas and steam distillation. The latter is the dominant mechanism and it is responsible for the improvement in the quality of the produced oil / Mestrado / Reservatórios e Gestão / Mestre em Ciências e Engenharia de Petróleo

Quantifying the influence of pressure on the effective elastic rock properties is important for applications in rock physics and reservoir characterization. Here I investigate the relationship between effective pressure and seismic velocities by performing inversion on the laboratory-measured data from a suite of clastic, carbonate and igneous rocks, using different analytic and discrete inversion schemes. I explore the utility of a physical model that models a natural fracture as supported by asperities of varying heights, when an effective pressure deforms the tallest asperities, bringing the shorter ones into contact while increasing the overall fracture stiffness. Thus, the model is known as the ?asperity-deformation? (ADM) or ?bed-of-nails? (BNM) model. Existing analytic solutions include one that assumes the host rock is infinitely more rigid than the fractures, and one that takes the host-rock compliance into account. Inversion results indicate that although both solutions can fit the data to within first-order approximation, some systematic misfits exist as a result of using the rigid-host solution, whereas compliant-host inversion returns smaller and random misfits, yet out-of-range parameter estimates. These problems indicate the effects of nonlinear elastic deformation whose degree varies from rock to rock. Consequently, I extend the model to allow for the pressure dependence of the host rock, thereby physically interpreting the nonlinear behaviors of deformation. Furthermore, I apply a discrete grid-search inversion scheme that generalizes the distribution of asperity heights, thus accurately reproduces velocity profiles, significantly improves the fit and helps to visualize the distribution of asperities. I compare the analytic and numerical asperity-deformation models with the existing physical model of elliptical ?pennyshape? cracks with a pore-aspect-ratio (PAR) spectrum in terms of physical meaning and data-fitting ability. The comparison results provide a link and demonstrate the consistency between the use of the two physical models, making a better understanding of the microstructure as well as the contact mechanism and physical behaviors of rocks under pressure. ADM-based solutions, therefore, have the potential to facilitate modeling and interpretation of applications such as time-lapse seismic investigations of fractured reservoirs.

[pt] Meios que apresentam descontinuidades como fraturas e falhas em um maciço rochoso ou reservatórios de petróleo impõem algumas dificuldades na simulação numérica pelo Método dos Elementos Finitos. Uma dessas dificuldades é a necessidade de geração de malhas muito refinadas, principalmente na região próxima à descontinuidade, até a obtenção de uma resposta confiável do modelo, o que pode consumir um tempo significativo. Ao mesmo tempo, a discretização da descontinuidade com elementos muito pequenos, quando comparados ao restante do modelo, pode conduzir a um aumento no tempo de simulação. Neste trabalho é apresentada a formulação de um elemento finito cortado por uma descontinuidade para aplicação em problemas com acoplamento fluido-mecânico em meios saturados com um único fluido. A inserção da descontinuidade no elemento é obtida pela adição de novos termos à função de interpolação, dispensando a sua discretização. Esses termos adicionais conseguem reproduzir a mudança no campo de deslocamento e poro-pressão no elemento devido à presença da descontinuidade. A resposta do elemento é verificada através da comparação com uma solução analítica unidimensional e com exemplos simples simulados em um programa comercial. / [en] Media that present discontinuities as fractures and faults in a rock mass or oil reservoirs impose some difficulties in numerical simulation using standard Finite Element Method. One of these difficulties is the need for very refined mesh generation, especially in the region near the discontinuities, to obtain a reliable answer of the model, which can consume significant time. At the same time, the discretization of the discontinuity with very small elements compared to the rest of the model may lead to an increase in simulation time. This work presents the formulation of an element that is crossed by a discontinuity for use in coupling fluid-mechanical problems in single fluid saturated mediums. The insertion of the discontinuity in the element is obtained by adding new terms in the interpolation function, which eliminates its discretization. These additional terms can reproduce the change in the displacement and pore pressure field in the element due to the presence of discontinuity. The response of the element is validated by comparing it with one dimensional analytical solution and simple examples simulated in a commercial program.

No description available.