Numerical Modeling of Coupled Fluid Flow and Geomechanical Stresses in a Petroleum Reservoir

2020 ◽  
Vol 142 (6) ◽  
Author(s):  
Manojkumar Gudala ◽  
Suresh Kumar Govindarajan

Abstract A fully coupled hydro and geomechanical model has been used to predict the transient pressure disturbance, reservoir deformation, and effective stress distribution in both homogeneous and heterogeneous reservoirs. The heterogeneous reservoir is conceptualized by explicitly considering the spatial distributions of porosity and permeability as against assuming it as constant values. The finite element method was used in the coupled model in conjunction with the poroelasticity. Transient pressure disturbance is significantly influenced by the overburden during the production in both homogeneous and heterogeneous reservoirs for all the perforation schemes. Perforation scheme 2 provides the optimum reservoir performance when compared with other three schemes in terms of transient pressure distribution and reservoir subsidence. It also has the ability to overcome both the water and gas coning problems when the reservoir fluid flow is driven by both gas cap and water drive mechanisms. A Biot–Willis coefficient is found to significantly influence both the pressure and stress distribution right from the wellbore to the reservoir boundary. Maximum effective stresses have been generated in the vicinity of the wellbore in the reservoir at a high Biot–Willis coefficient of 0.9. Thus, the present work clearly projects that a Biot–Willis coefficient of 0 cannot be treated to be a homogeneous reservoir by default, while the coupled effect of hydro and geomechanical stresses plays a very critical role. Therefore, the implementation of the coupled hydro and geomechanical numerical models can improve the prediction of transient reservoir behavior efficiently for the simple and complex geological systems effectively.

2021 ◽  
Author(s):  
Murat Zeybek ◽  
Lei Jiang ◽  
Hadrien Dumont

Abstract The radius of investigation (ROI) of pressure transient analyses has been traditionally assessed using analytical formulations with basic reservoir parameters for homogenous systems. Numerous studies aimed to improve ROI formulations to incorporate all reservoir and testing parameters such as gauge resolution and rate for more accurate ROI assessments. However, new generation wireline formation testers aim to improve deep transient tests with significant developments in gauge resolution and increasing rate. Challenges still remain in heterogeneous formations such as shaly sands and carbonate reservoirs. In this study, detailed conceptual high-resolution numerical models are set up, including comprehensive reservoir and measurement parameters, to investigate more realistic ROI assessments in layered heterogeneous systems without and with hydraulic communication. Several conceptual examples are presented in layered systems with permeability contrasts. In addition, deviation from infinite-acting radial flow (IAFR) and pressure propagation in highly heterogeneous layered systems are investigated to detect the presence of geological features, including closed boundary systems and the presence of a fault in the proximity of wellbore.


2021 ◽  
Vol 2021 ◽  
pp. 1-24
Author(s):  
Jiaming Lei ◽  
Jianmin Zhang ◽  
Lifang Zhang

The aerator can reduce erosion by mixing a large amount of air into the water in the solid wall area. The effectiveness of erosion reduction is mainly based on air concentration and its bubble size distribution. However, simultaneous simulation of the air concentration and its bubble size distribution in numerical simulations is still a hot and difficult area of research. Aiming at the downstream aerated flow of hydraulic aeration facilities, several numerical models, such as VOF, mixture, Euler, and Population Balance Model (PBM), are compared and verified by experiments. The results show that the CFD-PBM coupled model performs well compared to other conventional multiphase models. It can not only obtain the evolution law of the bubble distribution downstream of the aerator but also accurately simulate the recombination and evolution process of bubble aggregation and breakage. The Sauter mean diameter of the air bubbles in the aerated flow decreases along the way and eventually reaches a stable value. The bubble breakage is the main process in the development of the bubbles. It reveals the aeration law that the small air bubbles are closer to the bottom plate, while the large bubbles float up along the aerated flow, which provides a powerful support for the basic research on the mechanism of aeration and erosion reduction.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yousif M. Makeen ◽  
Xuanlong Shan ◽  
Mutari Lawal ◽  
Habeeb A. Ayinla ◽  
Siyuan Su ◽  
...  

AbstractThe Abu Gabra and Bentiu formations are widely distributed within the interior Muglad Basin. Recently, much attention has been paid to study, evaluate and characterize the Abu Gabra Formation as a proven reservoir in Muglad Basin. However, few studies have been documented on the Bentiu Formation which is the main oil/gas reservoir within the basin. Therefore, 33 core samples of the Great Moga and Keyi oilfields (NE Muglad Basin) were selected to characterize the Bentiu Formation reservoir using sedimentological and petrophysical analyses. The aim of the study is to de-risk exploration activities and improve success rate. Compositional and textural analyses revealed two main facies groups: coarse to-medium grained sandstone (braided channel deposits) and fine grained sandstone (floodplain and crevasse splay channel deposits). The coarse to-medium grained sandstone has porosity and permeability values within the range of 19.6% to 32.0% and 1825.6 mD to 8358.0 mD respectively. On the other hand, the fine grained clay-rich facies displays poor reservoir quality as indicated by porosity and permeability ranging from 1.0 to 6.0% and 2.5 to 10.0 mD respectively. A number of varied processes were identified controlling the reservoir quality of the studies samples. Porosity and permeability were enhanced by the dissolution of feldspars and micas, while presence of detrital clays, kaolinite precipitation, iron oxides precipitation, siderite, quartz overgrowths and pyrite cement played negative role on the reservoir quality. Intensity of the observed quartz overgrowth increases with burial depth. At great depths, a variability in grain contact types are recorded suggesting conditions of moderate to-high compactions. Furthermore, scanning electron microscopy revealed presence of micropores which have the tendency of affecting the fluid flow properties in the Bentiu Formation sandstone. These evidences indicate that the Bentiu Formation petroleum reservoir quality is primarily inhibited by grain size, total clay content, compaction and cementation. Thus, special attention should be paid to these inhibiting factors to reduce risk in petroleum exploration within the area.


GeoArabia ◽  
1996 ◽  
Vol 1 (2) ◽  
pp. 267-284
Author(s):  
John L. Douglas ◽  

ABSTRACT The North ‘Ain Dar 3-D geocellular model consists of geostatistical models for electrofacies, porosity and permeability for a portion of the Jurassic Arab-D reservoir of Ghawar field, Saudi Arabia. The reservoir consists of a series of shallow water carbonate shelf sediments and is subdivided into 10 time-stratigraphic slices on the basis of core descriptions and gamma/porosity log correlations. The North ‘Ain Dar model includes an electrofacies model and electrofacies-dependent porosity and permeability models. Sequential Indicator Simulations were used to create the electrofacies and porosity models. Cloud Transform Simulations were used to generate permeability models. Advantages of the geostatistical modeling approach used here include: (1) porosity and permeability models are constrained by the electrofacies model, i.e. by the distribution of reservoir rock types; (2) patterns of spatial correlation and variability present in well log and core data are built into the models; (3) data extremes are preserved and are incorporated into the model. These are critical when it comes to determining fluid flow patterns in the reservoir. Comparison of model Kh with production data Kh indicates that the stratigraphic boundaries used in the model generally coincide with shifts in fluid flow as indicated by flowmeter data, and therefore represent reasonable flow unit boundaries. Further, model permeability and production estimated permeability are correlated on a Kh basis, in terms of vertical patterns of distribution and cumulative Kh values at well locations. This agreement between model and well test Kh improves on previous, deterministic models of the Arab-D reservoir and indicates that the modeling approach used in North ‘Ain Dar should be applicable to other portions of the Ghawar reservoir.


Hydrology ◽  
2021 ◽  
Vol 8 (4) ◽  
pp. 146
Author(s):  
Javier Fernández-Pato ◽  
Pilar García-Navarro

Numerical simulation of flows that consider interaction between overland and drainage networks has become a practical tool to prevent and mitigate flood situations in urban environments, especially when dealing with intense storm events, where the limited capacity of the sewer systems can be a trigger for flooding. Additionally, in order to prevent any kind of pollutant dispersion through the drainage network, it is very interesting to have a certain monitorization or control over the quality of the water that flows in both domains. In this sense, the addition of a pollutant transport component to both surface and sewer hydraulic models would benefit the global analysis of the combined water flow. On the other hand, when considering a realistic large domain with complex topography or streets structure, a fine spatial discretization is mandatory. Hence the number of grid cells is usually very large and, therefore, it is necessary to use parallelization techniques for the calculation, the use of Graphic Processing Units (GPU) being one of the most efficient due to the leveraging of thousands of processors within a single device. In this work, an efficient GPU-based 2D shallow water flow solver (RiverFlow2D-GPU) is fully coupled with EPA’s Storm Water Management Model (SWMM). Both models are able to develop a transient water quality analysis taking into account several pollutants. The coupled model, referred to as RiverFlow2D-GPU UD (Urban Drainge) is applied to three real-world cases, covering the most common hydraulic situations in urban hydrology/hydraulics. A UK Environmental Agency test case is used as model validation, showing a good agreement between RiverFlow2D-GPU UD and the rest of the numerical models considered. The efficiency of the model is proven in two more complex domains, leading to a >100x faster simulations compared with the traditional CPU computation.


2020 ◽  
Author(s):  
Philipp Eichheimer ◽  
Marcel Thielmann ◽  
Wakana Fujita ◽  
Gregor J. Golabek ◽  
Michihiko Nakamura ◽  
...  

Abstract. Fluid flow on different scales is of interest for several Earth science disciplines like petrophysics, hydrogeology and volcanology. To parameterize fluid flow in large-scale numerical simulations (e.g. groundwater and volcanic systems), flow properties on the microscale need to be considered. For this purpose experimental and numerical investigations of flow through porous media over a wide range of porosities are necessary. In the present study we sinter glass bead media with various porosities. The microstructure, namely effective porosity and effective specific surface, is investigated using image processing. We determine flow properties like hydraulic tortuosity and permeability using both experimental measurements and numerical simulations. By fitting microstructural and flow properties to porosity, we obtain a modified Kozeny-Carman equation for isotropic low-porosity media, that can be used to simulate permeability in large-scale numerical models. To verify the modified Kozeny-Carman equation we compare it to the computed and measured permeability values.


2011 ◽  
Vol 223 ◽  
pp. 431-438 ◽  
Author(s):  
Aldo Attanasio ◽  
Elisabetta Ceretti ◽  
Cristian Cappellini ◽  
Claudio Giardini

In cutting field, residual stress distribution analysis on the workpiece is a very interesting topic. Indeed, the residual stress distribution affects fatigue life, corrosion resistance and other functional aspects of the workpiece. Recent studies showed that the development of residual stresses is influenced by the cutting parameters, tool geometry and workpiece material. For reducing the costs of experimental tests and residual stress measurement, analytical and numerical models have been developed. The aim of these models is the possibility of forecasting the residual stress distribution into the workpiece as a function of the selected process parameters. In this work the residual stress distributions obtained simulating cutting operations using a 3D FEM software and the corresponding simulation procedure are reported. In particular, orthogonal cutting operations of AISI 1045 and AISI 316L steels were performed. The FEM results were compared with the experimental residual stress distribution in order to validate the model effectiveness.


Sign in / Sign up

Export Citation Format

Share Document