Modelling and quantification of structural uncertainties in petroleum reservoirs assisted by a hybrid cartesian cut cell/enriched multipoint flux approximation approach

Efficient and profitable oil production is subject to make reliable predictions about reservoir performance. However, restricted knowledge about reservoir distributed properties and reservoir structure calls for History Matching in which the reservoir model is calibrated to emulate the field observed history. Such an inverse problem yields multiple history-matched models which might result in different predictions of reservoir performance. Uncertainty Quantification restricts the raised model uncertainties and boosts the model reliability for the forecasts of future reservoir behaviour. Conventional approaches of Uncertainty Quantification ignore large scale uncertainties related to reservoir structure, while structural uncertainties can influence the reservoir forecasts more intensely compared with petrophysical uncertainty. What makes the quantification of structural uncertainty impracticable is the need for global regridding at each step of History Matching process. To resolve this obstacle, we develop an efficient methodology based on Cartesian Cut Cell Method which decouples the model from its representation onto the grid and allows uncertain structures to be varied as a part of History Matching process. Reduced numerical accuracy due to cell degeneracies in the vicinity of geological structures is adequately compensated with an enhanced scheme of class Locally Conservative Flux Continuous Methods (Extended Enriched Multipoint Flux Approximation Method abbreviated to extended EMPFA). The robustness and consistency of proposed Hybrid Cartesian Cut Cell/extended EMPFA approach are demonstrated in terms of true representation of geological structures influence on flow behaviour. In this research, the general framework of Uncertainty Quantification is extended and well-equipped by proposed approach to tackle uncertainties of different structures such as reservoir horizons, bedding layers, faults and pinchouts. Significant improvements in the quality of reservoir recovery forecasts and reservoir volume estimation are presented for synthetic models of uncertain structures. Also this thesis provides a comparative study of structural uncertainty influence on reservoir forecasts among various geological structures

Numerical simulation of fracture pattern development and implications for fuid flow

Simulations are instrumental to understanding flow through discrete fracture geometric representations that capture the large-scale permeability structure of fractured porous media. The contribution of this thesis is threefold: an efficient finite-element finite-volume discretisation of the advection/diffusion flow equations, a geomechanical fracture propagation algorithm to create fractured rock analogues, and a study of the effect of growth on hydraulic conductivity. We describe an iterative geomechanics-based finite-element model to simulate quasi-static crack propagation in a linear elastic matrix from an initial set of random flaws. The cornerstones are a failure and propagation criterion as well as a geometric kernel for dynamic shape housekeeping and automatic remeshing. Two-dimensional patterns exhibit connectivity, spacing, and density distributions reproducing en echelon crack linkage, tip hooking, and polygonal shrinkage forms. Differential stresses at the boundaries yield fracture curving. A stress field study shows that curvature can be suppressed by layer interaction effects. Our method is appropriate to model layered media where interaction with neighbouring layers does not dominate deformation. Geomechanically generated fracture patterns are the input to single-phase flow simulations through fractures and matrix. Thus, results are applicable to fractured porous media in addition to crystalline rocks. Stress state and deformation history control emergent local fracture apertures. Results depend on the number of initial flaws, their initial random distribution, and the permeability of the matrix. Straightpath fracture pattern simplifications yield a lower effective permeability in comparison to their curved counterparts. Fixed apertures overestimate the conductivity of the rock by up to six orders of magnitude. Local sample percolation effects are representative of the entire model flow behaviour for geomechanical apertures. Effective permeability in fracture dataset subregions are higher than the overall conductivity of the system. The presented methodology captures emerging patterns due to evolving geometric and flow properties essential to the realistic simulation of subsurface processes

Automatic mesh representation of urban environments

A robust watertight mesh generation framework for urban cityscape and waterscape is proposed. The framework, consisting of a set of algorithms implemented in MATLAB, uses geospatial data available from OpenStreetMap and United States Geological Survey repositories, and incorporates Triangle - a popular two-dimensional Delaunay triangulation software - to develop the mesh. For the cityscape component, the facades of the buildings are meshed as structured triangular grids while the roofs and terrains are meshed as unstructured triangular grids using Triangle. For the waterscape component, quadrilateral cells are created based on the requirements of Environmental Fluids Dynamics Code (EFDC) model – a popular modeling platform for environmental fluid flow analysis. The resulting mesh generated is watertight with little human intervention and can serve as a significant preprocessing tool in environmental computational fluid dynamics. Although, there are a few existing methodologies in the literature, most are limited in capacity and are difficult to implement

Modelling of steady-state fluid flow in 3D fractured isotropic porous media: Application to effective permeability calculation

International audienceIn this paper, 3D steady-state fluid flow in a porous medium with many intersecting fractures is derived numerically by using collocation method. Fluid flow, in the matrix and fractures, are described by Darcy's law and Poisseuille's law, respectively. The recent theoretical development presented a general potential solution to model the steady-state flow in fractured porous media under a far-field condition. This solution is a hypersingular integral equation with pressure field in the fracture surfaces as the main unknown. The numerical procedure can resolve the problem for any form of fractures and also takes into account the interactions and the intersection between fractures. Once the pressure field, and then the flux field in fractures have been determined, the pressure field in the porous matrix is computed completely. The basic problem of a single fracture is investigated and a semi-analytical solution is presented. Using the solution obtained for a single fracture, Mori-Tanaka and self-consistent schemes are employed for upscalling the effective permeability of 3D fractured porous media

Towards large-scale modelling of fluid flow in fractured porous media

To date, the complexity of fractured porous media still precludes the direct incorporation of small-scale features into field-scale modelling. These features, however, can be instrumental in shaping and triggering coarsening instabilities and other forms of emergent behaviour which need to be considered on the field-scale. Here we develop numerical simulation methods for this purpose and demonstrate their improved performance in single-and two-phase flow simulations with models of fractured porous media. Material discontinuities in fractured porous media strongly influence single-and multi-phase fluid flow. When continuum methods are used to model transport across such interfaces, they smear out jump discontinuities of concentration or saturation. To overcome this drawback, we “explode” hybrid finite-element node-centred finite-volume models along these introducing complementary finite-volumes along the material interfaces. With this embedded discontinuity discretization we develop a transport scheme that realistically represents the dependent variable discontinuities arising at these interfaces. The main advantage of this new scheme is its ability to honour the flow effects that we know that these discontinuities have in physical experiments. We have also developed a new time-stepping control scheme for the transport equation. It allows the user to specify the volume fraction of the model in which he/she is prepared to relax the CFL condition. This scheme is applied in a study of the impact of fracture pattern development on solute transport. These two-dimensional simulations quantify the effect of the fractures on macro-scale dispersion in geomechanically generated fracture geometries, as opposed to stochastically generated ones. Among other insights, the results indicate that fracture density, fracture spacing, and the fracture-matrix flux ratio control anomalous mass transport in such media. We also find that it is crucial to embed discontinuities into large-scale models of heterogeneous porous media

A Force-Balanced Control Volume Finite Element Method for Multi-Phase Porous Media Flow Modelling

Dr D. Pavlidis would like to acknowledge the support from the following research grants: Innovate UK ‘Octopus’, EPSRC ‘Reactor Core-Structure Re-location Modelling for Severe Nuclear Accidents’) and Horizon 2020 ‘In-Vessel Melt Retention’. Funding for Dr P. Salinas from ExxonMobil is gratefully acknowledged. Dr Z. Xie is supported by EPSRC ‘Multi-Scale Exploration of Multi-phase Physics in Flows’. Part funding for Prof Jackson under the TOTAL Chairs programme at Imperial College is also acknowledged. The authors would also like to acknowledge Mr Y. Debbabi for supplying analytic solutions.Peer reviewedPublisher PD

Development of the VHP-Female Full-Body Computational Model and Its Applications for Biomedical Electromagnetic Modeling

Computational modeling offers better insight into a wide range of bioelectrical and biomechanical problems with improved tools for the design of medical devices and the diagnosis of pathologies. Electromagnetic modeling at low and high frequencies is particularly necessary. Modeling electromagnetic, structural, thermal, and acoustic response of the human body to different internal and external stimuli is limited by the availability of numerically efficient computational human models. This study describes the development to date of a computational full-body human model - Visible Human Project (VHP) - Female Model. Its unique feature is full compatibility both with MATLAB and specialized FEM computational software packages such as ANSYS HFSS/Maxwell 3D. This study also describes progress made to date in using the newly developed tools for segmentation. A visualization tool is implemented within MATLAB and is based on customized version of the constrained 2D Delaunay triangulation method for intersecting objects. This thesis applies a VHP - Female Model to a specific application, transcranial Direct Current Stimulation (tDCS). Transcranial Direct Current Stimulation has been beneficial in the stimulation of cortical activity and treatment of neurological disorders in humans. The placement of electrodes, which is cephalic versus extracephalic montages, is studied for optimal targeting of currents for a given functional area. Given the difficulty of obtaining in vivo measurements of current density, modeling of conventional and alternative electrode montages via the FEM has been utilized to provide insight into the tDCS montage performance. An insight into future work and potential areas of research, such as study of bone quality have been presented too

Kinetic Intersections as a Source of Information on Rate Coefficients

C

GIS platform for management of shallow geothermal resources

Premi extraordinari doctorat UPC curs 2015-2016, àmbit d’Enginyeria CivilThis thesis promotes an efficient use of shallow geothermal energy by means of an integrated management system to organize its exploitation. Shallow geothermal energy is a renewable resource based on thermal energy exchange with the ground. Due to the growth in demand for this energy, the development of management techniques to organize the exploitation of this resource is mandatory to protect both groundwater and the users' rights. Shallow geothermal performance of underground is closely related to groundwater behavior, so it is necessary to understand and improve the knowledge about it. Thus, an integrated methodology is proposed for the 3D visualization of underground resources related to groundwater. A set of tools named HEROS3D was developed in a GIS environment to support the generation of 3D entities representing geological, hydrogeological, hydrochemical and geothermal features. The GIS technology also gives a wide-ranging support to environmental modeling, either conceptual or numerical, especially to groundwater modeling. However, there is a scarcity of tools to implement the conceptual model in numerical modeling platforms. This transition needs of specific methodologies to adapt the geometries and alpha-numerical data from the conceptual model to the numerical model to get optimal numerical results. Although most necessities can be satisfied with inherent GIS tools, there are particular steps in the implementation of hydrogeological conceptual model into the numerical modeling software that have not been solved yet. To overcome this gap, a set of tools is presented, named ArcArAz. It focuses on the configuration of geometry and parameterization for groundwater numerical models. Once both the hydrogeological conceptual model and the numerical model are defined, a solid basis for management of Shallow Geothermal energy is available. This thesis proposes two methodologies for the management of this energy resource at two different scales: for a regional scale and for a metropolitan scale. The first GIS methodology provides a response to the need for a regional quantification of the geothermal potential that can be extracted by Boreholes Heat Exchangers and its associated environmental impacts. For the first time, advection and dispersion heat transport mechanisms and the temporal evolution from the start of operation of the BHE are considered in the regional estimation of the variables of interest. A sensitivity analysis leads to the conclusion that the consideration of dispersion effects and temporal evolution of the exploitation prevent significant differences up to a factor 2.5 in the heat exchange rate accuracy and up to several orders of magnitude in the impacts generated. To deepen the management of Shallow Geothermal Energy, this thesis proposes to establish a market of shallow geothermal energy use rights which would allow managing this resource at a metropolitan scale. This methodology is based on a GIS framework and is composed of a geospatial database to store the main information required to manage the installations and a set of GIS tools used to define, implant and control this use rights market. Thermal impacts derived from the exploitation of this resource can also be registered geographically, by taking into account the groundwater flow direction and adjusting the thermal impact to the available plot.Esta tesis promueve el uso eficiente de la geotermia somera a través de un sistema integrado de gestión de este recurso. La geotermia somera es un recurso renovable que se basa en el intercambio de energía con el suelo. Los Intercambiadores de calor, o Borehole Heat Exchangers (BHEs) se están popularizando como sistema para explotarla. Debido al crecimiento en la demanda de geotermia somera, es imprescindible establecer una gestión integrada de este recurso para organizar su explotación y proteger tanto a las aguas subterráneas como a los beneficiarios de esta energía renovable. Debido a que la geotermia somera está íntimamente relacionada con el comportamiento de las aguas subterráneas, es imprescindible ahondar y mejorar su conocimiento. Para ello, se propone una metodología para la visualización tridimensional de los recursos subterráneos relacionados con la hidrogeología. Se ha desarrollado un conjunto de herramientas, llamado HEROS3D, en un entorno SIG. Estas herramientas facilitan la creación de entidades tridimensionales que representan datos geológicos, hidrogeológicos, hidrogeoquímicos y geotermales. Están relacionadas con una base de datos donde tanto la información bruta como la interpretada se encuentran almacenadas. La tecnología SIG también da soporte, no sólo a la modelación conceptual, sino también a la numérica, especialmente en el caso de la hidrogeología. Para facilitar la implementación de los modelos conceptuales en las plataformas de modelación numérica, esta tesis presenta un segundo conjunto de herramientas llamado ArcArAz. Estas herramientas ofrecen soluciones a los problemas más comunes relacionados con la configuración de la geometría de entrada al modelo numérico, así como su parametrización. Las bases para una gestión eficiente de la geotermia somera se establecen llamado ArcArAz. Estas herramientas ofrecen soluciones a los problemas más comunes relacionados con la configuración de la geometría de entrada al modelo numérico, así como su parametrización. Las bases para una gestión eficiente de la geotermia somera se establecen una vez que hemos definido y están disponibles tanto el modelo hidrogeológico conceptual como el modelo numérico. En relación a este aspecto, en esta tesis se proponen dos metodologías de gestión enfocadas a escalas diferentes: escala regional y escala metropolitana o local. La primera metodología SIG ofrece una respuesta a la necesidad de una cuantificación regional del potencial geotérmico somero que puede extraerse con intercambiadores de calor o Borehole Heat Exchangers, así como sus impactos térmicos asociados. Por primera vez pueden tenerse en cuenta en la estimación regional de las variables de interés la advección y dispersión de calor, como mecanismos de transporte de calor, así como la evolución temporal desde el inicio de la explotación. Un análisis de sensibilidad demuestra que la consideración de los efectos de dispersión así como el régimen temporal de la explotación supone diferencias de hasta 2.5 veces el potencial extraído y hasta de varios ordenes de magnitud en los impactos térmicos generados. Para profundizar en la gestión de la geotermia somera a escala local, esta tesis propone establecer un mercado de derechos de uso de este recurso. Esta metodología se ha implementado en un ambiente SIG y está compuesta de una base de datos donde se almacena la información principal necesaria para gestionar las instalaciones y de un conjunto de herramientas para definir, implantar y controlar este mercado de derechos de uso de geotermia somera. Los impactos térmicos derivados de la explotación de este recurso pueden quedar registrados geográficamente, teniendo en cuenta la dirección de flujo de las aguas subterráneas y ajustando estos impactos a la superficie de la parcela disponible una vez que hemos definido y están disponibles tanto el modelo hidrogeológico conceptual como el modelo numérico. En relación a este aspecto, en esta tesis se proponen dos metodologías de gestión enfocadas a escalas diferentes: escala regional y escala metropolitana o local. La primera metodología SIG ofrece una respuesta a la necesidad de una cuantificación regional del potencial geotérmico somero que puede extraerse con intercambiadores de calor o Borehole Heat Exchangers, así como sus impactos térmicos asociados. Por primera vez pueden tenerse en cuenta en la estimación regional de las variables de interés la advección y dispersión de calor, como mecanismos de transporte de calor, así como la evolución temporal desde el inicio de la explotación. Un análisis de sensibilidad demuestra que la consideración de los efectos de dispersión así como el régimen temporal de la explotación supone diferencias de hasta 2.5 veces el potencial extraído y hasta de varios ordenes de magnitud en los impactos térmicos generados. Para profundizar en la gestión de la geotermia somera a escala local, esta tesis propone establecer un mercado de derechos de uso de este recurso. Esta metodología se ha implementado en un ambiente SIG y está compuesta de una base de datos donde se almacena la información principal necesaria para gestionar las instalaciones y de un conjunto de herramientas para definir, implantar y controlar este mercado de derechos de uso de geotermia somera. Los impactos térmicos derivados de la explotación de este recurso pueden quedar registrados geográficamente, teniendo en cuenta la dirección de flujo de las aguas subterráneas y ajustando estos impactos a la superficie de la parcela disponibleAward-winningPostprint (published version