Submitted by Fabio Sobreira Campos da Costa (fabio.sobreira@ufpe.br) on 2017-07-20T12:25:34Z No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) documento.pdf: 6110083 bytes, checksum: e763b9e4b979081c4ada6fef0eb596a6 (MD5) / Made available in DSpace on 2017-07-20T12:25:34Z (GMT). No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) documento.pdf: 6110083 bytes, checksum: e763b9e4b979081c4ada6fef0eb596a6 (MD5) Previous issue date: 2015-08-03 / ANP-PRH26 / Aimportânciadageomecânicaedoestudodeesquemasdeacoplamentoentreageomecânica e fluxo multifásico têm sido cada vez mais importantes e utilizados pela indústria a medida que formações cada vez mais profundas vêem sendo descobertas e exploradas. O entendimento do comportamento do estado de tensão em um reservatório permite produzir um melhor entendimento das implicações geomecânicas que ocorrem durante a fase de explotação, isso porque durante esta fase, as alterações na poro-pressão conduzem perturbações no equilíbrio mecânico afetando o estado de tensão de formações profundas, de maneira a alterar as propriedades da rocha tais como permeabilidade e porosidade. No entanto, a simulação acoplada (hidromecânica) em um grande campo heterogêneo implica na solução de equações de fluxo e mecânica, associadas a um grande número de graus de liberdade que torna esse tipo de abordagem inviável e computacionalmente cara. Neste contexto, um simulador geomecânico-linhas de fluxoé apresentado dentro de um algoritmo sequencial iterativo. Neste trabalho, aplica-se o método de elementos finitos com volume de controle para o subproblema poro-mecânico que fornece um campo de velocidade de Darcy pós-processado e a porosidade como entradas para o subproblema de transporte. Este subproblema é resolvido através do método de decomposição de operador, no qual basea-se em um esquema preditor-corretor com os passos preditor e corretor discretizados pelos esquemas baseados em tempo de vôo e volumes finitos, respectivamente. Simulações numéricas de injeção de água foram comparadas com soluções encontradas na literatura, mostrando bons resultados. Em problemas dominados pela advecção, envolvendo um reservatório naturalmente fraturado, a abordagem implementada foi capaz de predizer a distribuição do campo de saturação ao longo de toda simulação. Além disso, para avaliar a resposta geomecânica, simulações numéricas foram realizadas em um grande sistema de reservatório-rocha capeadora em uma fase de recuperação primária de hidrocarboneto, mostrou que a formulação apresentada provou ser: uma alternativa promissora para simulação hidro-geomecânica tradicional; úteis para o modelo de fluxo de redução de ordem nos casos em que o comportamento geomecânico são mais importantes do que o comportamento de fluxo e de uma ferramenta complementar para simulação geomecânica convencional. / The importance of geomechanics and the study of coupling between geomechanics and multiphase flow have been increasingly recognized and used by the industry as deeper formations are discovered and exploited. The knowledge of the state of stress in a reservoir yields a better understanding of the geomechanical implications during exploitation stage, because during the primary recovery stage, changes in pore pressure leads to perturbations inthemechanicalequilibrium,affectingthestressstateintheformationsinawaythatalters the rock properties such as permeability and porosity. However, the coupled simulation (hydromechanical) in large field heterogeneous models involves stress and flow equations solving, associated with a large number of degrees-of-freedom which becomes infeasible and computationally costly. In this context, a geomechanical-streamline simulator is presented within a iteratively coupled framework algorithm. In the present work, we applied control volume finite element method for the poromechanics subproblem which provides a Darcy velocityfieldthroughapost-processingvelocityprocedureandporosityasinputfieldstothe transportsubproblem.Suchsubproblemissolvedbymeansofanoperatorsplittingmethod, which is based on a predictor-corrector scheme with the predictor and corrector steps discretized by a time-of-flight and a finite volume based schemes, respectively. Numerical simulations of water-flooding are compared to the numerical results available in literature, showing good results. In convection-dominated problems, involving a naturally fractured reservoir, the approach was able to predict the saturation distributions for the whole simulation correctly. Furthermore, to appraisal the geomechanical response, numerical simulation was performed in a large reservoir-caprock system in a primary hydrocarbon recovery stage, showing that the formulation presented proved be: an promising alternative to traditional hydro-geomechanical simulation; useful for flow model order reduction in cases where the geomechanical behavior are more important than the flow behavior and a complementary tool for conventional geomechanical simulations.

In this research we consider developing a reservoir simulator capable of simulating complex coupled poromechanical processes on massively parallel computers. A variety of problems arising from petroleum and environmental engineering inherently necessitate the understanding of interactions between fluid flow and solid mechanics. Examples in petroleum engineering include reservoir compaction, wellbore collapse, sand production, and hydraulic fracturing. In environmental engineering, surface subsidence, carbon sequestration, and waste disposal are also coupled poromechanical processes. These economically and environmentally important problems motivate the active pursuit of robust, efficient, and accurate simulation tools for coupled poromechanical problems. Three coupling approaches are currently employed in the reservoir simulation community to solve the poromechanics system, namely, the fully implicit coupling (FIM), the explicit coupling, and the iterative coupling. The choice of the coupling scheme significantly affects the efficiency of the simulator and the accuracy of the solution. We adopt the fixed-stress iterative coupling scheme to solve the coupled system due to its advantages over the other two. Unlike the explicit coupling, the fixed-stress split has been theoretically proven to converge to the FIM for linear poroelasticity model. In addition, it is more efficient and easier to implement than the FIM. Our computational results indicate that this approach is also valid for multiphase flow. We discretize the quasi-static linear elasticity model for geomechanics in space using the continuous Galerkin (CG) finite element method (FEM) on general hexahedral grids. Fluid flow models are discretized by locally mass conservative schemes, specifically, the mixed finite element method (MFE) for the equation of state compositional flow on Cartesian grids and the multipoint flux mixed finite element method (MFMFE) for the single phase and two-phase flows on general hexahedral grids. While both the MFE and the MFMFE generate cell-centered stencils for pressure, the MFMFE has advantages in handling full tensor permeabilities and general geometry and boundary conditions. The MFMFE also obtains accurate fluxes at cell interfaces. These characteristics enable the simulation of more practical problems. For many reservoir simulation applications, for instance, the carbon sequestration simulation, we need to account for thermal effects on the compositional flow phase behavior and the solid structure stress evolution. We explicitly couple the poromechanics equations to a simplified energy conservation equation. A time-split scheme is used to solve heat convection and conduction successively. For the convection equation, a higher order Godunov method is employed to capture the sharp temperature front; for the conduction equation, the MFE is utilized. Simulations of coupled poromechanical or thermoporomechanical processes in field scales with high resolution usually require parallel computing capabilities. The flow models, the geomechanics model, and the thermodynamics model are modularized in the Integrated Parallel Accurate Reservoir Simulator (IPARS) which has been developed at the Center for Subsurface Modeling at the University of Texas at Austin. The IPARS framework handles structured (logically rectangular) grids and was originally designed for element-based data communication, such as the pressure data in the flow models. To parallelize the node-based geomechanics model, we enhance the capabilities of the IPARS framework for node-based data communication. Because the geomechanics linear system is more costly to solve than those of flow and thermodynamics models, the performance of linear solvers for the geomechanics model largely dictates the speed and scalability of the coupled simulator. We use the generalized minimal residual (GMRES) solver with the BoomerAMG preconditioner from the hypre library and the geometric multigrid (GMG) solver from the UG4 software toolbox to solve the geomechanics linear system. Additionally, the multilevel k-way mesh partitioning algorithm from METIS is used to generate high quality mesh partitioning to improve solver performance. Numerical examples of coupled poromechanics and thermoporomechanics simulations are presented to show the capabilities of the coupled simulator in solving practical problems accurately and efficiently. These examples include a real carbon sequestration field case with stress-dependent permeability, a synthetic thermoporoelastic reservoir simulation, poroelasticity simulations on highly distorted hexahedral grids, and parallel scalability tests on a massively parallel computer. / text

L’eau des circuits d’une centrale nucléaire est purifiée à l’aide de résines échangeuses d’ions. La prédiction de leurs performances constitue une aide importante pour l’exploitation de ces réacteurs. Les méthodes du transport réactif sont particulièrement adaptées pour cela et constituent la base du code OPTIPUR, dédié à la modélisation de ces résines. Le travail présenté comporte trois axes principaux. Le premier est l’intégration d’une limitation au transfert de masse dans une colonne de résines échangeuses d’ions, avec une mobilité spécifique à chaque espèce chimique, dans le cadre d’un découplage des calculs de chimie et de transport. Ce modèle permet, sans paramètre ajustable, de reproduire assez fidèlement une série d’expériences réalisées précédemment par le CEA. Le second axe concerne les aspects numériques du transport réactif, avec l’utilisation de la méthode d’Anderson pour accélérer la convergence du couplage chimie-transport dans un schéma itératif. En utilisant les informations issues des itérations précédentes et sans modification majeure du code, la robustesse et les temps de calcul ont pu être nettement améliorés. La troisième thématique abordée est celle de l’équilibre d’échange d’ions. Les bases d’un modèle prenant en considération l’évolution de l’humidité de la résine, ainsi que son élasticité sont proposées ; les interactions entre groupes fonctionnels, contre-ions et eau sont considérées comme des équilibres chimiques. Les constantes d’équilibre sont ajustées à partir de mesures de la teneur en eau de la résine à différentes pressions de vapeur d’eau. Finalement, des coefficients de sélectivité apparents peuvent être calculés et comparés aux mesures disponibles. / In nuclear power plants, the water contained in different circuits is purified by passing through ion exchange resins. Prediction of the performance of these resins is an important help to the plant operators. To this end, the method of reactive transport modeling are well suited and is the basis of the OPTIPUR code that was designed to model the resins. The work presented in this manuscript covers three main aspects. The first one is the integration of a limitation to mass transfer in a ion exchange deep bed, taking into account a specific mobility for each chemical species, in the context of separated calculations for chemistry and transport. This model was shown to reproduce experimental data, without adjustable parameters. The second part of this work deals with the numerical aspects of reactive transport modelling. A method developped by Anderson was used to accelerate the convergence of the chemistry transport coupling in an iterative scheme. Using the information from previous iterations, and without major changes in the code, calculation times were largely decreased, as well as the number of calculations failures. The third topic is ion exchange equilibrium. The basis of a model that takes into account the change in the water content of the resin and its elasticity are described. The interactions between the fonctional groups, the counterions and water are considered as chemical reactions. The corresponding equilibrium constants are fitted to measurements of the water content of the resin at different relative humidity. Finaly, the selectivity coefficients can be calculated and compared to litterature values.

Tight gas and shale oil play an important role in energy security and in meeting an increasing energy demand. Hydraulic fracturing is a widely used technology for recovering these resources. The design and evaluation of hydraulic fracture operation is critical for efficient production from tight gas and shale plays. The efficiency of fracturing jobs depends on the interaction between hydraulic (induced) and naturally occurring discrete fractures. In this work, a coupled reservoir-fracture flow model is described which accounts for varying reservoir geometries and complexities including non-planar fractures. Different flow models such as Darcy flow and Reynold's lubrication equation for fractures and reservoir, respectively are utilized to capture flow physics accurately. Furthermore, the geomechanics effects have been included by considering a multiphase Biot's model. An accurate modeling of solid deformations necessitates a better estimation of fluid pressure inside the fracture. The fractures and reservoir are modeled explicitly allowing accurate representation of contrasting physical descriptions associated with each of the two. The approach presented here is in contrast with existing averaging approaches such as dual and discrete-dual porosity models where the effects of fractures are averaged out. A fracture connected to an injection well shows significant width variations as compared to natural fractures where these changes are negligible. The capillary pressure contrast between the fracture and the reservoir is accounted for by utilizing different capillary pressure curves for the two features. Additionally, a quantitative assessment of hydraulic fracturing jobs relies upon accurate predictions of fracture growth during slick water injection for single and multistage fracturing scenarios. It is also important to consistently model the underlying physical processes from hydraulic fracturing to long-term production. A recently introduced thermodynamically consistent phase-field approach for pressurized fractures in porous medium is utilized which captures several characteristic features of crack propagation such as joining, branching and non-planar propagation in heterogeneous porous media. The phase-field approach captures both the fracture-width evolution and the fracture-length propagation. In this work, the phase-field fracture propagation model is briefly discussed followed by a technique for coupling this to a fractured poroelastic reservoir simulator. We also present a general compositional formulation using multipoint flux mixed finite element (MFMFE) method on general hexahedral grids with a future prospect of treating energized fractures. The mixed finite element framework allows for local mass conservation, accurate flux approximation and a more general treatment of boundary conditions. The multipoint flux inherent in MFMFE scheme allows the usage of a full permeability tensor. An accurate treatment of diffusive/dispersive fluxes owing to additional velocity degrees of freedom is also presented. The applications areas of interest include gas flooding, CO₂ sequestration, contaminant removal and groundwater remediation. / text

Submitted by Renata Lopes (renatasil82@gmail.com) on 2018-06-29T15:51:42Z No. of bitstreams: 0 / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2018-07-03T15:49:10Z (GMT) No. of bitstreams: 0 / Made available in DSpace on 2018-07-03T15:49:10Z (GMT). No. of bitstreams: 0 Previous issue date: 2018-03-27 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho será proposto um método de acoplamento iterativo de interface com um esquema de subcycling no tempo eficiente e preciso. Este será aplicado a pro-blemas do tipo sólido-fluido discretizados, respectivamente, pelos métodos dos elementos finitos clássico (MEF) e espectral (MEFE). Adicionalmente, será proposta uma melhoria no esquema de subcycling, de modo que para convergir não sejam necessários métodos de relaxação. Aplicando o MEFE em subdomínios com geometrias pouco distorcidas, pode-se usufruir da alta precisão numérica com baixo custo de armazenamento oferecidos pelo método ao mesmo tempo em que é possível aplicar o MEF aos subdomínios com geometrias complexas, acrescentando versatilidade ao método. Diferentes exemplos nu-méricos são apresentados e analisados para demonstrar a precisão e a potencialidade das formulações numéricas propostas. / In this work, an iterative interface coupling method with an efficient and precise time subcycling scheme is proposed. It is applied to solid-fluid type problems discretized respectively by classical finite element method (FEM) and spectral finite element method (SFEM), additionally, an improvement in the subcycling scheme is proposed so as not to require relaxation methods to converge. Applying the SFEM in subdomains with low distorted geometries one can take advantage of the high numerical precision with low cost of storage offered by the method, while it is possible to apply the FEM in subdomains with complex geometries, adding versatility to the method. Many numerical examples are presented and analyzed here to show the accuracy and potentiality of the proposed numerical formulations.

Les activités humaines dans la subsurface se développent rapidement (stockage de déchets,nouvelles techniques minières, stockage à haute fréquence de l’énergie), alors que dans le même temps les attentes du public et des autorités s’intensifient. L’évaluation de chaque étape de ces opérations souterraines repose sur des études détaillées de la sûreté et des impacts environnementaux.Elles reposent sur des simulateurs élaborés et sur de la modélisation multiphysique. Avec leur approche orientée processus, les simulations en transport réactifs proposent une méthode efficace pour comprendre et prévoir le comportement de ces systèmes complexes, à différentes échelles de temps et d’espace.Le but de ce travail est d’intégrer la résolution de l’écoulement diphasique compressible dans le cadre de codes de transport réactifs à l’aide d’une méthode de séparation d’opérateurs. Un module multiphasique a été créé dans le code de transport réactif HYTEC. Une nouvelle approche a ensuite été développée pour coupler écoulement multicomposant multiphasique compressible, description de propriétés thermo-dynamiques complexes pour les fluides, avec des codes de transport réactif. La méthode a été intégrée dans HYTEC. Des cas de validation sont proposés, puis des exemples d’application pour la simulation du stockage souterrain de CO2 et des impuretés associées. / Human activity in the subsurface has rapidly been expanding and diversifying (waste disposal, new mining technologies, high-frequency storage of energy), while the public and regulatory expectations keep growing. The assessment of each step of underground operations requires careful safety and environmental impact evaluations. They rely on elaborate simulators and multiphysics modeling. With its process-based approach, reactive transport simulation provides an effective way to understand and predict the behavior of such complex systems at different time and spatial scale.This work aims at incorporating a compressible multiphase flow into conventional reactive transport framework by an operator splitting approach. A multiphase flow module is developed in the HYTEC reactive transport software. A new approach is then developed to fully couple multiphase multicomponent compressible flow, the complex thermodynamic description of the fluid properties, with existing reactive transport codes. The method is implemented in HYTEC. Some validation is provided, before application to the simulation of underground storage of CO2 and associated impurities.