Mathematical modeling of thermal stabilization of vertical wells on high performance computing systems

Temperature stabilization of oil and gas wells is used to ensure stability and prevent deformation of a subgrade estuary zone. In this work, we consider the numerical simulation of thermal stabilization using vertical seasonal freezing columns. A mathematical model of such problems is described by a time-dependent temperature equation with phase transitions from water to ice. The resulting equation is a standard nonlinear parabolic equation. Numerical implementation is based on the finite element method using the package Fenics. After standard purely implicit approximation in time and simple linearization, we obtain a system of linear algebraic equations. Because the size of freezing columns are substantially less than the size of the modeled area, we obtain mesh refinement near columns. Due to this, we get a large system of equations which are solved using high performance computing systems.Comment: 9 pages, 5 figure

Parallel Splitting and Decomposition Method for Computations of Heat Distribution in Permafrost

A mathematical model, numerical algorithm and program code for simulation and long-term forecasting of changes in permafrost as a result of operation of a multiple well pad of northern oil and gas field are presented. In the model the most significant climatic and physical factors are taken into account such as solar radiation, determined by specific geographical location, heterogeneous structure of frozen soil, thermal stabilization of soil, possible insulation of the objects, seasonal fluctuations in air temperature, and freezing and thawing of the upper soil layer. A parallel algorithm of decomposition with splitting by spatial variables is presented

Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers

Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and investigate numerically the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geologocal configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems

Modeling Thermal Turbulence Using Implicit Large Eddy Simulation

A general description of a thermally coupled fluid flow is given by the incompressible Navier-Stokes equations coupled with the heat equation using Boussinesq approximation, whose mathematical structure is much well understood. A variational multiscale finite element approximation has been considered for the formulation of incompressible Navier-Stokes equation and heat equation. The complexity of these problems makes their numerical solution very difficult as the standard finite element method is unstable. In the incompressible Navier Stokes equations, two well known sources of numerical instabilities are the incompressibility constraint and the presence of the convective term. Many stabilization techniques used nowadays are based on scale separation, splitting the unknown into a coarse part induced by the discretization of the domain and a fine subgrid part. The modeling of the subgrid scale and its influence leads to a modified coarse scale problem providing stability. In convection-diffusion problem once global instabilities have been overcome by a stabilization method, there are still local oscillations near layers due to the lack of monotonicity of the method. Shock capturing techniques are often employed to deal with them. Proper choice of stabilization and shock capturing techniques can eliminate the local instabilities near layers of convection-diffusion equation. A very important issue of the formulation presented in this thermally coupled incompressible flow is the possibility to model turbulent flows. Some terms involving the velocity subgrid scale arise from the convective term in the Navier-Stokes equations which can be understood as the contribution from the Reynolds tensor of a LES approach and the contribution from the cross stress tensor. This opens the door of modeling thermal turbulence using LES automatically inherited by the formulation used in this work. Different classical benchmark problems are numerically solved in this thesis work for the convection-diffusion equation to show the capabilities of different combination of stabilization and shock capturing methods. In the case of thermally coupled incompressible flows some numerical and industrial examples are exhibited to check the performance of the different combination of stabilization and shock capturing methods and to compare them. The objective is to conclude which method works better to approximate the exact solution and eliminate instabilities and local oscillations

Finite element simulations: computations and applications to aerodynamics and biomedicine.

171 p.Las ecuaciones en derivadas parciales describen muchos fenómenos de interés práctico y sus solucionessuelen necesitar correr simulaciones muy costosas en clústers de cálculo.En el ámbito de los flujos turbulentos, en particular, el coste de las simulaciones es demasiado grande sise utilizan métodos básicos, por eso es necesario modelizar el sistema.Esta tesis doctoral trata principalmente de dos temas en Cálculo Científico.Por un lado, estudiamos nuevos desarrollos en la modelización y simulación de flujos turbulentos;utilizamos un Método de Elementos Finitos adaptativo y un modelo de ¿número de Reynolds infinito¿para reducir el coste computacional de simulaciones que, sin estas modificaciones, serían demasiadocostosas.De esta manera conseguimos lograr simulaciones evolutivas de flujos turbulentos con número deReynolds muy grande, lo cual se considera uno de los mayores retos en aerodinámica.El otro pilar de esta tesis es una aplicación biomédica.Desarrollamos un modelo computacional de Ablación (Cardiaca) por Radiofrecuencia, una terapiacomún para tratar varias enfermedades, por ejemplo algunas arritmias.Nuestro modelo mejora los modelos existentes en varias maneras, y en particular en tratar de obteneruna aproximación fiel de la geometría del sistema, lo cual se descubre ser crítico para simularcorrectamente la física del fenómeno

Analytical design method for cold production of heavy oil with bottom water using Bilateral Sink Wells

Few heavy oil reservoirs with strong bottom water drives have been developed successfully because of severe water coning. Water coning tends to cause low ultimate recovery, low well productivity, and high water production. Although thermal and gravity-assisted methods might improve recovery in oil reservoirs, such methods are widely perceived as either economically unfavorable or technologically infeasible. This study proposes a new, cold production technique, called Bilateral Water Sink (BWS), to meet those challenges. The BWS method suppresses water cresting by producing oil and water simultaneously from separate, horizontal wells completed in the oil and water zones; the oil and water completions are parallel, with the oil well directly above the water well. In conventional horizontal well production, water cresting causes water to bypass oil, making the water drive mechanism ineffective. BWS controls water invasion by altering the pressure distribution in the near-well area. With cresting suppressed, the oil completion remains water-free, allowing water to displace oil from the edges of the well drainage area to the oil completion, increasing ultimate recovery. Unlike existing heavy oil recovery methods, BWS exploits the natural reservoir energy in the bottom water drive. This makes BWS economically, technically, and environmentally appealing – especially for offshore applications, where cold production is currently the only option and oil-water separation is a problem. In this study, BWS oil recovery is investigated analytically and numerically. A new mathematical model identifies controlling variables and project design parameters, and describes the relationships among them. The design model is used to select rates of water and oil in BWS wells for best performance. The analytical model is verified by a comparison to numerical simulations. These two approaches together provide the quantitative account of the BWS’s effect on avoiding water cresting and improving oil recovery. The results show that BWS can increase oil recovery from 10 percent to over 40 percent in a conventional case, while avoiding the problem of oil-contaminated water production. As a result, the mathematical model of BWS well behavior is shown to be a practical reservoir management tool to guide development of heavy oil reservoirs with bottom water drives

Aeronautical Engineering. A continuing bibliography, supplement 115

This bibliography lists 273 reports, articles, and other documents introduced into the NASA scientific and technical information system in October 1979