A Multiple-Continuum Approach For Modeling Multiphase Flow in Naturally Fractured Vuggy Petroleum Reservoirs
By Yu-Shu Wu, Guan Qin, R. Ewing, Y. Efendiev, Zhijiang Kang, Yulin Ren. All Days, 2006
87 citations recorded by the source corpus.
The existence of vugs or cavities in naturally fractured reservoirs has long been observed. Even though these vugs can be largely attributed to reserves of oil, natural gas, and groundwater, few investigations of vuggy fractured reservoirs have been conducted. In this paper, a new multiple-continuum conceptual model is developed, based on geological data and observations of core samples from carbonate formations in China, to investigate multiphase flow behavior in such vuggy fractured reservoirs. The conceptual model has been implemented into a three-dimensional, three-phase reservoir simulator with a generalized multiple-continuum modeling approach. The conceptual model considers vuggy fractured rock as a triple- or multiple-continuum medium, consisting ofhighly permeable fractures,low-permeability rock matrix, andvarious-sized vugs. The matrix system may contain a large number of small or isolated cavities (of centimeters or millimeters in diameter), whereas vugs are larger cavities, with sizes from centimeters to meters in diameter, indirectly connected to fractures through small fractures or microfractures. Similar to the conventional double-porosity model, the fracture continuum is responsible for the occurrence of global flow, while vuggy and matrix continua, providing storage space, are locally connected to each other (and interacting with globally connecting fractures). For practical application of the multi-continuum concept in reservoir simulation, we propose a novel upscaling method for computing equivalent gridblock permeabilities of coarse blocks containing large isolated vugs, in which the local problems consisting of Darcy and Stokes flows are solved. In addition, we describe an efficient boundary condition for accurate computation of upscaled permeabilities. In the numerical implementation, a control-volume, integral finite-difference method is used for spatial discretization, and a first-order finite-difference scheme is adapted for temporal discretization of governing flow equations in each continuum. The resulting discrete nonlinear equations are solved fully implicitly by Newton iteration. The numerical scheme is verified and applied to simulate water-oil flow through the fractured vuggy reservoirs of the Tahe Oil Field in China. Since the 1960s, investigation of flow and transport processes in fractured reservoirs has received much attention with significant progress being made. Driven by the increasing need to develop petroleum and geothermal reservoirs (as well as to resolve subsurface contamination problems), many numerical modeling approaches and techniques have been developed [6, 25, 15, 21]. The petroleum industry is currently facing a growing demand for oil and natural gas, while at the same time fewer new oil reserves exist worldwide. The efficient development of naturally fractured reservoirs, possessing a large portion of current world oil and gas reserves, has become a top priority. Thanks to the known low oil recovery rates from these reservoirs, interest in enhancing oil and gas recovery from such reservoirs has grown, with more investigations conducted for multiphase flow and transport phenomena in fractured reservoirs [2, 19, 4, 12]. Also, environmental concerns over subsurface contamination have motivated many related, similar studies [31, 26, 7]. Even though significant progress has been made towards understanding and modeling of flow and transport processes in fractured rock since the 1960s [6, 25, 15, 16, 21], most studies have focused primarily on naturally fractured reservoirs without taking into consideration of cavities. Recently, characterizing vuggy fractured rock has received attention, because a number of fractured vuggy reservoirs have been found worldwide that can significantly contribute to reserves and the production of oil and gas [17, 23, 18, 13, 9].
Subfields: Computational Theory and Mathematics, Environmental Engineering, Mechanical Engineering
View Paper: https://www.semanticscholar.org/paper/e476406107bea2ddccb8d476ec4d9e2c5a024b0b
Citations in the corpus
- Numerical Calculation of Equivalent Permeability Tensor for Fractured Vuggy Porous Media Based on Homogenization Theory (2011, unknown): “Email addresses: huangzhqin@gmail.com (Z. Huang), yaojunhdpu@126.com (J. Yao), yajunok@126.com (Y. Li), wcc1220@163.com (C. Wang), xinruiabc@gmail.com (X. Lv) http://www.global-sci.com/ 180 c©2011 Global-Science Press Z. Huang et al. / Commun.”