The representation of two phase fault-rock properties in flow simulation models

2002 ◽  
Vol 8 (2) ◽  
pp. 119-132 ◽  
Author(s):  
T. Manzocchi ◽  
A.E. Heath ◽  
J.J. Walsh ◽  
C. Childs
Author(s):  
John White

Modelling water vapour flow, heat transfer and porosity in porous adsorbent is somewhat challenging simulation problem. Primary macroscopic water vapour flow models, such as Darcy's law, fail to predict the pressure drop entirely correctly for the reason that many of flow parameters not considered because of the simplifications that remain made for the multi-scale structure of the porous adsorbents. For one to develop a good physical understanding of such water vapour flows and the accuracy of existing 3D simulation models, there is a need for some accurate 3D geometry to be studied. This present work describes two-phase water vapour flow and adsorption/ desorption performed on porous adsorbent by a Dynamic vapour sorption (DVS). The CFD simulation results are associated with experiments results. It is decided that for such complex porous adsorbent CFD simulation problems the use of COMSOL Multiphysics and SolidWorks flow simulation will be utilised.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Md Saiful Islam ◽  
Tom Manzocchi

AbstractFaults are generally represented in conventional upscaled models as 2D planar surfaces with transmissibility multipliers used to represent single-phase fault properties. However, faults are structurally complex 3D zones in which both single-phase and two-phase fault rock properties can be significant. Ignoring this structural and petrophysical complexity within faults may impart considerable inaccuracy on the predictive performance of upscaled models. This study has developed a two-phase flow-based geometrical upscaling method capable of representing simultaneously the complex geometry and saturation-dependent two-phase flow properties of realistic fault zones. In this approach, high-resolution sector models are built of small portions of the fault zones and assigned appropriate single-phase and two-phase fault rock properties. Steady state two-phase flow simulations at different fractional flows of oil and water are used to determine the saturation dependent upscaled pseudo relative permeability functions which are incorporated into upscaled models. The method is applied to an example model containing two 3D fault zone components and tested by comparing the flow results of upscaled model with those of a high-resolution truth model. Results show that two-phase flow-based geometrical upscaling is a promising method for representing the effects of two-phase fault rock properties and complex 3D fault zone structure simultaneously.


2008 ◽  
Vol 11 (06) ◽  
pp. 1071-1081 ◽  
Author(s):  
Amy Whitaker ◽  
C. Shah Kabir ◽  
Wayne Narr

Summary The extent to which fractures affect fluid pathways is a vital component of understanding and modeling fluid flow in any reservoir. We examined the Wafra Ratawi grainstone for which production extending for 50 years, including recent horizontal drilling, has provided some clues about fractures, but their exact locations, intensity, and overall effect have been elusive. In this study, we find that a limited number of total fractures affect production characteristics of the Ratawi reservoir. Although fractures occur throughout the Wafra field, fracture-influenced reservoir behavior is confined to the periphery of the field where the matrix permeability is low. This work suggests that for the largest part of the field, explicit fractures are not necessary in the next-generation Earth and flow-simulation models. The geologic fracture assessment included seismic fault mapping and fracture interpretation of image logs and cores. Fracture trends are in the northeast and southwest quadrants, and fractures are mineralized toward the south and west of the field. Pressure-falloff tests on some peripheral injectors indicate partial barriers, and most of these wells lie on seismic-scale faults in the reservoir, suggesting partial sealing. A few wells show fractured-reservoir production characteristics, and rate-transient analysis on a few producers indicates localized dual-porosity behavior. Producers proximal to dual-porosity wells display single-porosity behavior, however, to attest to the notion of localized fracture response. The spatially restricted fracture-flow characteristics appear to correlate with fracture or vug zones in a low-permeability reservoir. Presence of fracture-flow behavior was tested by constructing the so-called flow-capacity index (FCI), the ratio of khwell (well test-derived value) to khmatrix (core-derived property). Data from 80 wells showed khmatrix to be consistently higher than khwell, a relationship that suggests insignificant fracture production in these wells. Introduction The Wafra field is in the Partitioned Neutral Zone (PNZ) between Kuwait and Saudi Arabia, as shown in Fig. 1. The field has been producing since the 1950s and has seen renewed drilling activity since the late 1990s, including horizontal drilling and implementation of peripheral water injection (Davis and Habib 1999). The Lower Cretaceous Ratawi formation contains the most reserves of the producing intervals at Wafra. The Ratawi oolite (a misnomer--it is a grainstone) reservoir has variable porosity (5 to 35%) and permeability that ranges from tens to hundreds of md (Longacre and Ginger 1988). The main Wafra structure is a gentle (i.e., interlimb angle >170°), doubly plunging anticline trending north-northwest to south-southeast, which culminates near its northern end. The East Wafra spur is a north-trending branch that extends from the center of the main Wafra structure. As seen in Fig. 1, relief on the Main Wafra structure exceeds that on East Wafra. The Ratawi oolite in the Wafra field has been studied at length, and various authors have reported geologic and engineering elements, leading to reservoir characterization and understanding of reservoir performance. Geologic studies are those of Waite et al. (2000) and Sibley et al. (1997). In contrast, Davis and Habib (1999) presented implementation of peripheral water injection, whereas Chawathé et al. (2006) discussed realignment of injection pattern owing to lack of pressure support in the reservoir interior. Previous studies considered the reservoir to behave like a single-porosity system. But recent image-log fracture interpretations indicate high fracture densities, suggesting that the implementation of a dual-porosity model may be necessary because the high impact of fractures during field development has been recognized in some Middle East reservoirs for more than 50 years (Daniel 1954). Static and dynamic data are required to characterize fracture reservoir behavior accurately (Narr et al. 2006). Geologic description of the fracture system, by use of cores, borehole images, seismic data, and well logs, does not in itself determine whether fractures affect reservoir behavior. While seismic and some image logs were available to locate fractures in the Wafra Ratawi reservoir, no dynamic testing with the specific objective of understanding fracture impact has occurred. So, to determine whether fractures influence oil productivity significantly, we used diagnostic analyses of production data and well tests of available injectors. The assessment of fracture effects in the Ratawi reservoir will be used to guide the next generation of geologic and flow-simulation models. Dynamic data involving pressure and rate have the potential to reveal the influence of open fractures in production performance. Unfortunately, pressure-transient testing on single wells does not always provide conclusive evidence about the presence of fractures with the characteristic dual-porosity dip on the pressure-derivative signature (Bourdet et al. 1989). That is because a correct mixture of matrix/fracture storativity must be present for the characteristic signature to appear (Serra et al. 1983). In practice, interference testing (Beliveau 1989) between wells appears to provide more-definitive clues about interwell connectivity, leading to inference about fractures. In contrast to pressure-transient testing, rate-transient analysis offers the potential to provide the same information without dedicated testing. In this field, all wells are currently on submersible pumps. Consequently, the pump-intake pressure and measured rate provided the necessary data for pressure/rate convolution or rate-transient analysis. We provide the Ratawi-reservoir case study primarily as an example of the integration of diverse geologic and engineering data to develop an assessment of fracture influence on reservoir behavior. It illustrates the use of production-data diagnostic tests to determine fracture influence in the absence of targeted fracture-analysis testing. The workflow can be applied to similar static/dynamic problems, such as fault-transmissivity determination. Secondly, this analysis illustrates the process of deciding that fractures, although present throughout the reservoir, may not lead to widespread fractured-reservoir characteristics (e.g., Allan and Sun 2003).


2021 ◽  
Author(s):  
Mokhles Mezghani ◽  
Mustafa AlIbrahim ◽  
Majdi Baddourah

Abstract Reservoir simulation is a key tool for predicting the dynamic behavior of the reservoir and optimizing its development. Fine scale CPU demanding simulation grids are necessary to improve the accuracy of the simulation results. We propose a hybrid modeling approach to minimize the weight of the full physics model by dynamically building and updating an artificial intelligence (AI) based model. The AI model can be used to quickly mimic the full physics (FP) model. The methodology that we propose consists of starting with running the FP model, an associated AI model is systematically updated using the newly performed FP runs. Once the mismatch between the two models is below a predefined cutoff the FP model is switch off and only the AI model is used. The FP model is switched on at the end of the exercise either to confirm the AI model decision and stop the study or to reject this decision (high mismatch between FP and AI model) and upgrade the AI model. The proposed workflow was applied to a synthetic reservoir model, where the objective is to match the average reservoir pressure. For this study, to better account for reservoir heterogeneity, fine scale simulation grid (approximately 50 million cells) is necessary to improve the accuracy of the reservoir simulation results. Reservoir simulation using FP model and 1024 CPUs requires approximately 14 hours. During this history matching exercise, six parameters have been selected to be part of the optimization loop. Therefore, a Latin Hypercube Sampling (LHS) using seven FP runs is used to initiate the hybrid approach and build the first AI model. During history matching, only the AI model is used. At the convergence of the optimization loop, a final FP model run is performed either to confirm the convergence for the FP model or to re iterate the same approach starting from the LHS around the converged solution. The following AI model will be updated using all the FP simulations done in the study. This approach allows the achievement of the history matching with very acceptable quality match, however with much less computational resources and CPU time. CPU intensive, multimillion-cell simulation models are commonly utilized in reservoir development. Completing a reservoir study in acceptable timeframe is a real challenge for such a situation. The development of new concepts/techniques is a real need to successfully complete a reservoir study. The hybrid approach that we are proposing is showing very promising results to handle such a challenge.


Sign in / Sign up

Export Citation Format

Share Document