Using an Integrated Multidimensional Scaling and Clustering Method to Reduce the Number of Scenarios Based on Flow-Unit Models Under Geological Uncertainties

2020 ◽  
Vol 142 (6) ◽  
Author(s):  
Seyed Kourosh Mahjour ◽  
Manuel Gomes Correia ◽  
Antonio Alberto de Souza dos Santos ◽  
Denis José Schiozer

Abstract Understanding the role of geological uncertainties on reservoir management decisions requires an ensemble of reservoir models that cover the uncertain space of parameters. However, in most cases, high computation time is needed for the flow simulation step, which can have a negative impact on a suitable assessment of flow behavior. Therefore, one important point is to choose a few scenarios from the ensemble of models while preserving the geological uncertainty range. In this study, we present a statistical solution to select the representative models (RMs) based on a novel scheme of measuring the similarity between 3D flow-unit models. The proposed method includes the integration of multidimensional scaling and cluster analysis (IMC). IMC can be applied to the models before the simulation process to save time and costs. To check the validity of the methodology, numerical simulation and then uncertainty analysis are carried out on the RMs and full set. We create an ensemble of 200 3D flow-unit models through the Latin Hypercube sampling method. The models indicate the geological uncertainty range for properties such as permeability, porosity, and net-to-gross. This method is applied to a synthetic benchmark model named UNISIM-II-D and proves to offer good performance in reducing the number of models so that only 9% of the models in the ensemble (18 selected models from 200 models) can be sufficient for the uncertainty quantification if appropriate similarity measures and clustering methods are used.

Author(s):  
Mitsugu Yamaguchi ◽  
Tatsuaki Furumoto ◽  
Shuuji Inagaki ◽  
Masao Tsuji ◽  
Yoshiki Ochiai ◽  
...  

AbstractIn die-casting and injection molding, a conformal cooling channel is applied inside the dies and molds to reduce the cycle time. When the internal face of the channel is rough, both cooling performance and tool life are negatively affected. Many methods for finishing the internal face of such channels have been proposed. However, the effects of the channel diameter on the flow of a low-viscosity finishing media and its finishing characteristics for H13 steel have not yet been reported in the literature. This study addresses these deficiencies through the following: the fluid flow in a channel was computationally simulated; the flow behavior of abrasive grains was observed using a high-speed camera; and the internal face of the channel was finished using the flow of a fluid containing abrasive grains. The flow velocity of the fluid with the abrasive grains increases as the channel diameter decreases, and the velocity gradient is low throughout the channel. This enables reduction in the surface roughness for a short period and ensures uniform finishing in the central region of the channel; however, over polishing occurs owing to the centrifugal force generated in the entrance region, which causes the form accuracy of the channel to partially deteriorate. The outcomes of this study demonstrate that the observational finding for the finishing process is consistent with the flow simulation results. The flow simulation can be instrumental in designing channel diameters and internal pressures to ensure efficient and uniform finishing for such channels.


Author(s):  
Manuel Gomes Correia ◽  
Célio Maschio ◽  
Denis José Schiozer

Super-giant carbonate fields, such as Ghawar, in Saudi Arabia, and Lula, at the Brazilian pre-salt, show highly heterogeneous behavior that is linked to high permeability intervals in thin layers. This article applies Local Grid Refinements (LGR) integrated with upscaling procedures to improve the representation of highly laminated reservoirs in flow simulation by preserving the static properties and dynamic trends from geological model. This work was developed in five main steps: (1) define a conventional coarse grid, (2) define LGR in the conventional coarse grid according to super-k and well locations, (3) apply an upscaling procedure for all scenarios, (4) define LGR directly in the simulation model, without integrate geological trends in LGR and (5) compare the dynamic response for all cases. To check results and compare upscaling matches, was used the benchmark model UNISIM-II-R, a refined model based on a combination of Brazilian Pre-salt and Ghawar field information. The main results show that the upscaling of geological models for coarse grid with LGR in highly permeable thin layers provides a close dynamic representation of geological characterization compared to conventional coarse grid and LGR only near-wells. Pseudo-relative permeability curves should be considered for (a) conventional coarse grid or (b) LGR scenarios under dual-medium flow simulations as the upscaling of discrete fracture networks and dual-medium flow models presents several limitations. The conventional approach of LGR directly in simulation model, presents worse results than LGR integrated with upscaling procedures as the extrapolation of dynamic properties to the coarse block mismatch the dynamic behavior from geological characterization. This work suggests further improvements for results for upscaling procedures that mask the flow behavior in highly laminated reservoirs.


2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Zhenglun Alan Wei ◽  
Zhongquan Charlie Zheng ◽  
Xiaofan Yang

A parallel implementation of an immersed-boundary (IB) method is presented for low Reynolds number flow simulations in a representative elementary volume (REV) of porous media that are composed of a periodic array of regularly arranged structures. The material of the structure in the REV can be solid (impermeable) or microporous (permeable). Flows both outside and inside the microporous media are computed simultaneously by using an IB method to solve a combination of the Navier–Stokes equation (outside the microporous medium) and the Zwikker–Kosten equation (inside the microporous medium). The numerical simulation is firstly validated using flow through the REVs of impermeable structures, including square rods, circular rods, cubes, and spheres. The resultant pressure gradient over the REVs is compared with analytical solutions of the Ergun equation or Darcy–Forchheimer law. The good agreements demonstrate the validity of the numerical method to simulate the macroscopic flow behavior in porous media. In addition, with the assistance of a scientific parallel computational library, PETSc, good parallel performances are achieved. Finally, the IB method is extended to simulate species transport by coupling with the REV flow simulation. The species sorption behaviors in an REV with impermeable/solid and permeable/microporous materials are then studied.


2021 ◽  
Vol 345 ◽  
pp. 00015
Author(s):  
Matěj Jeřábek ◽  
Michal Volf ◽  
Daniel Duda

The article describes a numerical simulation of flow in the cooling system of an electromagnetic calorimeter by analysing the temperature and pressure fields. Two fundamentally different approaches were used to analyse the pressure field - analytical 1D calculation and numerical 3D flow simulation. The article contains a detailed evaluation and description of individual analyses using the commercial software ANSYS 2020 R1.


Energies ◽  
2019 ◽  
Vol 12 (19) ◽  
pp. 3663
Author(s):  
Lindsey Rasmussen ◽  
Tianguang Fan ◽  
Alex Rinehart ◽  
Andrew Luhmann ◽  
William Ampomah ◽  
...  

The efficiency of carbon utilization and storage within the Pennsylvanian Morrow B sandstone, Farnsworth Unit, Texas, is dependent on three-phase oil, brine, and CO2 flow behavior, as well as spatial distributions of reservoir properties and wettability. We show that end member two-phase flow properties, with binary pairs of oil–brine and oil–CO2, are directly dependent on heterogeneity derived from diagenetic processes, and evolve progressively with exposure to CO2 and changing wettability. Morrow B sandstone lithofacies exhibit a range of diagenetic processes, which produce variations in pore types and structures, quantified at the core plug scale using X-ray micro computed tomography imaging and optical petrography. Permeability and porosity relationships in the reservoir permit the classification of sedimentologic and diagenetic heterogeneity into five distinct hydraulic flow units, with characteristic pore types including: macroporosity with little to no clay filling intergranular pores; microporous authigenic clay-dominated regions in which intergranular porosity is filled with clay; and carbonate–cement dominated regions with little intergranular porosity. Steady-state oil–brine and oil–CO2 co-injection experiments using reservoir-extracted oil and brine show that differences in relative permeability persist between flow unit core plugs with near-constant porosity, attributable to contrasts in and the spatial arrangement of diagenetic pore types. Core plugs “aged” by exposure to reservoir oil over time exhibit wettability closer to suspected in situ reservoir conditions, compared to “cleaned” core plugs. Together with contact angle measurements, these results suggest that reservoir wettability is transient and modified quickly by oil recovery and carbon storage operations. Reservoir simulation results for enhanced oil recovery, using a five-spot pattern and water-alternating-with-gas injection history at Farnsworth, compare models for cumulative oil and water production using both a single relative permeability determined from history matching, and flow unit-dependent relative permeability determined from experiments herein. Both match cumulative oil production of the field to a satisfactory degree but underestimate historical cumulative water production. Differences in modeled versus observed water production are interpreted in terms of evolving wettability, which we argue is due to the increasing presence of fast paths (flow pathways with connected higher permeability) as the reservoir becomes increasingly water-wet. The control of such fast-paths is thus critical for efficient carbon storage and sweep efficiency for CO2-enhanced oil recovery in heterogeneous reservoirs.


1974 ◽  
Vol 39 (1) ◽  
pp. 51-74 ◽  
Author(s):  
R. G. Matson ◽  
D. L. True

AbstractThis study is a comparison of the results of a variety of clustering methods and 2 multidimensional scaling techniques on data from sites in northern Chile. While differences do occur, the similarities among the results are strong in spite of differing inputs. In general, results of relative frequency analysis appear to be superior to those of presence/absence, and the techniques used seem to be viable additions to existing archaeological tools.


2005 ◽  
Author(s):  
Nobuyuki Yomoda ◽  
Masahiko Kubo ◽  
Norihiko Watanabe

2015 ◽  
Vol 32 (2) ◽  
pp. 289-307 ◽  
Author(s):  
Yexiang Xiao ◽  
Zhengwei Wang ◽  
Jidi Zeng ◽  
jintai Zheng ◽  
Jiayang Lin ◽  
...  

Purpose – The purpose of this paper is to experimentally and numerically investigate the interference characteristics between two ski-jump jets on the flip bucket in a large dam spillway when two floodgates are running. Design/methodology/approach – The volume of fluid (VOF) method together with the Realizable k-ε turbulence model were used to predict the flow in two ski-jump jets and the free surface motion in a large dam spillway. The movements of the two gates were simulated using a dynamic mesh controlled by a User Defined Function (UDF). The simulations were run using the prototype dam as the field test to minimize errors due to scale effects. The simulation results are compared with field test observations. Findings – The transient flow calculations, accurately predict the two gate discharges compared to field data with the predicted ski-jump jet interference flow pattern similar to the observed shapes. The transient simulations indicate that the main reason for the deflected nappe is the larger opening difference between the two gates as the buttress side gate closes. When both gates are running, the two ski-jump jets interfere in the flip bucket and raise the jet nappe to near the buttress to form a secondary flow on this jet nappe surface. As the gate continues to close, the nappe surface continues to rise and the surface secondary flow become stronger, which deflects the nappe over the side buttress. Originality/value – A dynamic mesh is used to simulate the transient flow behavior of two prototype running gates. The transient flow simulation clarifies the hydraulics mechanism for how the two ski-jump jets interfere and deflect the nappe.


2011 ◽  
Vol 44 ◽  
pp. 16
Author(s):  
Erke Arıbaş ◽  
Şenol Pişkin ◽  
M. Serdar Çelebi
Keyword(s):  

Sign in / Sign up

Export Citation Format

Share Document