Differences in the Upscaling Procedure for Compositional Reservoir Simulations of Immiscible and Miscible Gas Flooding

2021 ◽  
Author(s):  
Victor de Souza Rios ◽  
Arne Skauge ◽  
Ken Sorbie ◽  
Gang Wang ◽  
Denis José Schiozer ◽  
...  

Abstract Compositional reservoir simulation is essential to represent the complex interactions associated with gas flooding processes. Generally, an improved description of such small-scale phenomena requires the use of very detailed reservoir models, which impact the computational cost. We provide a practical and general upscaling procedure to guide a robust selection of the upscaling approaches considering the nature and limitations of each reservoir model, exploring the differences between the upscaling of immiscible and miscible gas injection problems. We highlight the different challenges to achieve improved upscaled models for immiscible and miscible gas displacement conditions with a stepwise workflow. We first identify the need for a special permeability upscaling technique to improve the representation of the main reservoir heterogeneities and sub-grid features, smoothed during the upscaling process. Then, we verify if the use of pseudo-functions is necessary to correct the multiphase flow dynamic behavior. At this stage, different pseudoization approaches are recommended according to the miscibility conditions of the problem. This study evaluates highly heterogeneous reservoir models submitted to immiscible and miscible gas flooding. The fine models represent a small part of a reservoir with a highly refined set of grid-block cells, with 5 × 5 cm2 area. The upscaled coarse models present grid-block cells of 8 × 10 m2 area, which is compatible with a refined geological model in reservoir engineering studies. This process results in a challenging upscaling ratio of 32 000. We show a consistent procedure to achieve reliable results with the coarse-scale model under the different miscibility conditions. For immiscible displacement situations, accurate results can be obtained with the coarse models after a proper permeability upscaling procedure and the use of pseudo-relative permeability curves to improve the dynamic responses. Miscible displacements, however, requires a specific treatment of the fluid modeling process to overcome the limitations arising from the thermodynamic equilibrium assumption. For all the situations, the workflow can lead to a robust choice of techniques to satisfactorily improve the coarse-scale simulation results. Our approach works on two fronts. (1) We apply a dual-porosity/dual-permeability upscaling process, developed by Rios et al. (2020a), to enable the representation of sub-grid heterogeneities in the coarse-scale model, providing consistent improvements on the upscaling results. (2) We generate specific pseudo-functions according to the miscibility conditions of the gas flooding process. We developed a stepwise procedure to deal with the upscaling problems consistently and to enable a better understanding of the coarsening process.

Author(s):  
Guillaume Bonnaffoux ◽  
Cécile Melis

After having relied for decades totally on small scale model tests, the design of anchoring systems for moored floaters like FPSOs is now widely performed numerically. Estimation of design maxima during mooring analyses requires calculating system response for a large number of sea states in order to screen all possible scenarios between wind, waves and current parameters. In addition, the slow-drift response motions of the system constituted by the floater and the anchoring system are highly dependent of the wave elevation realization, which is not an input parameter of the simulations and can lead to extremely variable responses. This is generally addressed by designers by performing the analysis for the same sea state and varying the wave group spectrum (or wave components phases) a large number of times N (20–50 realizations is a typical range, see Refs [4], [5]). For these N realizations, N response maxima are extracted, and a distribution of response extremes is derived, from which the response level is extracted. In terms of computational cost, performing N 3-hour simulations to derive N values of extreme response is extremely expensive. The paper will focus on methods that can be employed to reduce the computational cost of analyses. In a first step, the rapidity of statistical convergence of response estimates depending on the system will be investigated. This will allow pre-determining the number of sea states realizations required to reach a satisfactory convergence of response. In a second step, a mean of improving the computational efficiency of calculations carried-out to reach the statistical convergence will be proposed.


Author(s):  
Jeong-Tae Kim ◽  
So-Young Lee ◽  
Dong-Soo Hong ◽  
Jin-Hak Yi ◽  
Yoon-Koo Kang

In this study, vibration responses of harbor caisson structures with various boundary conditions are experimentally examined as a fundamental study to develop a health assessment technique for harbor structures. To achieve the objective, four-step approach is implemented. Firstly, a target caisson structure is selected and a small-scale model of the caisson model is constructed in the laboratory. Secondly, a finite element model of the caisson model is generated to analyze dynamic responses of the structure. Thirdly, experimental tests are performed on the caisson model to obtain dynamic responses under various boundary conditions and impact locations. Four different boundary conditions, 1) ‘hanging by crane’, 2) ‘standing on styrofoam block’, 3) ‘standing on sand mat’, and 4) ‘standing on concrete floor’ are considered. Finally, variation of frequency response ratio assurance criterion and correlation coefficients of frequency response functions are analyzed.


SPE Journal ◽  
2016 ◽  
Vol 21 (04) ◽  
pp. 1192-1203 ◽  
Author(s):  
A.. Alkhatib ◽  
M.. Babaei

Summary Reservoir heterogeneity can be detrimental to the success of surfactant/polymer enhanced-oil-recovery (EOR) processes. Therefore, it is important to evaluate the effect of uncertainty in reservoir heterogeneity on the performance of surfactant/polymer EOR. Usually, a Monte Carlo sampling approach is used, in which a number of stochastic reservoir-model realizations are generated and then numerical simulation is performed to obtain a certain objective function, such as the recovery factor. However, Monte Carlo simulation (MCS) has a slow convergence rate and requires a large number of samples to produce accurate results. This can be computationally expensive when using large complex reservoir models. This study applies a multiscale approach to improve the efficiency of uncertainty quantification. This method is known as the multilevel Monte Carlo (MLMC) method. This method comprises performing a small number of expensive simulations on the fine-scale model and a large number of less-expensive simulations on coarser upscaled models, and then combining the results to produce the quantities of interest. The purpose of this method is to reduce computational cost while maintaining the accuracy of the fine-scale model. The results of this approach are compared with a reference MCS, assuming a large number of simulations on the fine-scale model. Other advantages of the MLMC method are its nonintrusiveness and its scalability to incorporate an increasing number of uncertainties. This study uses the MLMC method to efficiently quantify the effect of uncertainty in heterogeneity on the recovery factor of a chemical EOR process, specifically surfactant/polymer flooding. The permeability field is assumed to be the random input. This method is first demonstrated by use of a Gaussian 3D reservoir model. Different coarsening algorithms are used and compared, such as the renormalization method and the pressure-solver method (PSM). The results are compared with running Monte Carlo for the fine-scale model while equating the computational cost for the MLMC method. Both of these results are then compared with the reference case, which uses a large number of runs of the fine-scale model. The method is then extended to a channelized non-Gaussian generated 3D reservoir model incorporating multiphase upscaling The results show that it is possible to robustly quantify spatial uncertainty for a surfactant/polymer EOR process while greatly reducing the computational requirement, up to two orders of magnitude compared with traditional Monte Carlo for both the Gaussian and non-Gaussian reservoir models. The method can be easily extended to other EOR processes to quantify spatial uncertainty, such as carbon dioxide (CO2) EOR. Other possible extensions of this method are also discussed.


2013 ◽  
Vol 14 (2) ◽  
Author(s):  
Noor Fachrizal

Biomass such as agriculture waste and urban waste are enormous potency as energy resources instead of enviromental problem. organic waste can be converted into energy in the form of liquid fuel, solid, and syngas by using of pyrolysis technique. Pyrolysis process can yield higher liquid form when the process can be drifted into fast and flash response. It can be solved by using microwave heating method. This research is started from developing an experimentation laboratory apparatus of microwave-assisted pyrolysis of biomass energy conversion system, and conducting preliminary experiments for gaining the proof that this method can be established for driving the process properly and safely. Modifying commercial oven into laboratory apparatus has been done, it works safely, and initial experiments have been carried out, process yields bio-oil and charcoal shortly, several parameters are achieved. Some further experiments are still needed for more detail parameters. Theresults may be used to design small-scale continuous model of productionsystem, which then can be developed into large-scale model that applicable for comercial use.


2021 ◽  
Vol 11 (6) ◽  
pp. 2652
Author(s):  
Jung Han Kim ◽  
Ick-Hyun Kim ◽  
Jin Ho Lee

When a seismic force acts on bridges, the pier can be damaged by the horizontal inertia force of the superstructure. To prevent this failure, criteria for seismic reinforcement details have been developed in many design codes. However, in moderate seismicity regions, many existing bridges were constructed without considering seismic detail because the detailed seismic design code was only applied recently. These existing structures should be retrofitted by evaluating their seismic performance. Even if the seismic design criteria are not applied, it cannot be concluded that the structure does not have adequate seismic performance. In particular, the performance of a lap-spliced reinforcement bar at a construction joint applied by past practices cannot be easily evaluated analytically. Therefore, experimental tests on the bridge piers considering a non-seismic detail of existing structures need to be performed to evaluate the seismic performance. For this reason, six small scale specimens according to existing bridge piers were constructed and seismic performances were evaluated experimentally. The three types of reinforcement detail were adjusted, including a lap-splice for construction joints. Quasi-static loading tests were performed for three types of scale model with two-column piers in both the longitudinal and transverse directions. From the test results, the effect on the failure mechanism of the lap-splice and transverse reinforcement ratio were investigated. The difference in failure characteristics according to the loading direction was investigated by the location of plastic hinges. Finally, the seismic capacity related to the displacement ductility factor and the absorbed energy by hysteresis behavior for each test were obtained and discussed.


2017 ◽  
Vol 2017 ◽  
pp. 1-9 ◽  
Author(s):  
Yiwei Lu ◽  
Hanlong Liu ◽  
Changjie Zheng ◽  
Xuanming Ding

X-section cast-in-place concrete pile is a new type of foundation reinforcement technique featured by the X-shaped cross-section. Compared with a traditional circular pile, an X-section pile with the same cross-sectional area has larger side resistance due to its larger cross-sectional perimeter. The behavior of static loaded X-section pile has been extensively reported, while little attention has been paid to the dynamic characteristics of X-section pile. This paper introduced a large-scale model test for an X-section pile and a circular pile with the same cross-sectional area subjected to cyclic axial load in sand. The experimental results demonstrated that cyclic axial load contributed to the degradation of shaft friction and pile head stiffness. The dynamic responses of X-section pile were determined by loading frequency and loading amplitude. Furthermore, comparative analysis between the X-section pile and the circular pile revealed that the X-section pile can improve the shaft friction and reduce the cumulative settlement under cyclic loading. Static load test was carried out prior to the vibration tests to investigate the ultimate bearing capacity of test piles. This study was expected to provide a reasonable reference for further studies on the dynamic responses of X-section piles in practical engineering.


Author(s):  
Segen F. Estefen ◽  
Paulo Roberto da Costa ◽  
Eliab Ricarte ◽  
Marcelo M. Pinheiro

Wave energy is a renewable and non-polluting source and its use is being studied in different countries. The paper presents an overview on the harnessing of energy from waves and the activities associated with setting up a plant for extracting energy from waves in Port of Pecem, on the coast of Ceara State, Brazil. The technology employed is based on storing water under pressure in a hyperbaric chamber, from which a controlled jet of water drives a standard turbine. The wave resource at the proposed location is presented in terms of statistics data obtained from previous monitoring. The device components are described and small scale model tested under regular waves representatives of the installation region. Based on the experimental results values of prescribed pressures are identified in order to optimize the power generation.


2007 ◽  
Vol 46 (7) ◽  
pp. 1067-1079 ◽  
Author(s):  
M. Kanda ◽  
M. Kanega ◽  
T. Kawai ◽  
R. Moriwaki ◽  
H. Sugawara

Abstract Urban climate experimental results from the Comprehensive Outdoor Scale Model (COSMO) were used to estimate roughness lengths for momentum and heat. Two different physical scale models were used to investigate the scale dependence of the roughness lengths; the large scale model included an aligned array of 1.5-m concrete cubes, and the small scale model had a geometrically similar array of 0.15-m concrete cubes. Only turbulent data from the unstable boundary layers were considered. The roughness length for momentum relative to the obstacle height was dependent on wind direction, but the scale dependence was not evident. Estimated values agreed well with a conventional morphometric relationship. The logarithm of the roughness length for heat relative to the obstacle height depended on the scale but was insensitive to wind direction. COSMO data were used successfully to regress a theoretical relationship between κB−1, the logarithmic ratio of roughness length for momentum to heat, and Re*, the roughness Reynolds number. Values of κB−1 associated with Re* for three different urban sites from previous field experiments were intercompared. A surprising finding was that, even though surface geometry differed from site to site, the regressed function agreed with data from the three urban sites as well as with the COSMO data. Field data showed that κB−1 values decreased as the areal fraction of vegetation increased. The observed dependency of the bulk transfer coefficient on atmospheric stability in the COSMO data could be reproduced using the regressed function of Re* and κB−1, together with a Monin–Obukhov similarity framework.


2011 ◽  
Vol 25 (4) ◽  
pp. 1485-1496 ◽  
Author(s):  
Min Kuang ◽  
Zhengqi Li ◽  
Pengfei Yang ◽  
Jinzhao Jia ◽  
Qunyi Zhu

Sign in / Sign up

Export Citation Format

Share Document