Modeling transient flow behavior with the high velocity non-Darcy effect in composite naturally fractured-homogeneous gas reservoirs

Author(s):  
Ren-Shi Nie ◽  
Xiaohui Fan ◽  
Min Li ◽  
Zhangxin Chen ◽  
Qi Deng ◽  
...  
2019 ◽  
Vol 2019 ◽  
pp. 1-13 ◽  
Author(s):  
Kongjie Wang ◽  
Lian Wang ◽  
Caspar Daniel Adenutsi ◽  
Zhiping Li ◽  
Sen Yang ◽  
...  

To improve the carbonate gas reservoir development and production, highly deviated wells (HDW) are widely used in the field. Production decline analysis of HDW is crucial for long-term gas reservoir development. However, it is a new challenge to incorporate the complex pore structure of naturally fractured-vuggy carbonate gas reservoirs and evaluate the production performance of HDW. This paper presents a semianalytical model to analyze the pressure and production behavior of HDW in naturally fractured-vuggy carbonate gas reservoirs, which consist of fractures, vugs, and matrix. The primary flow occurs only through the fracture and the outer boundary is closed. Introducing pseudopressure and pseudotime, the Laplace transformation, Fourier transformation, and its inverse and Stehfest numerical inversion were employed to establish a point source and line source solutions. Furthermore, the validity of the proposed model was verified by comparing a field data from the Arum River Basin in Turkmenistan. Finally, the effects of major parameters on the production decline curves were analyzed by using the proposed model and it was found that they had influences at different stages of gas production history and the sensitivity intensity of each parameter was different. With its high efficiency and simplicity, this semianalytical model will serve as a useful tool to evaluate the well production behavior for the naturally fractured-vuggy carbonate gas reservoirs.


Author(s):  
Atheer Dheyauldeen ◽  
Omar Al-Fatlawi ◽  
Md Mofazzal Hossain

AbstractThe main role of infill drilling is either adding incremental reserves to the already existing one by intersecting newly undrained (virgin) regions or accelerating the production from currently depleted areas. Accelerating reserves from increasing drainage in tight formations can be beneficial considering the time value of money and the cost of additional wells. However, the maximum benefit can be realized when infill wells produce mostly incremental recoveries (recoveries from virgin formations). Therefore, the prediction of incremental and accelerated recovery is crucial in field development planning as it helps in the optimization of infill wells with the assurance of long-term economic sustainability of the project. Several approaches are presented in literatures to determine incremental and acceleration recovery and areas for infill drilling. However, the majority of these methods require huge and expensive data; and very time-consuming simulation studies. In this study, two qualitative techniques are proposed for the estimation of incremental and accelerated recovery based upon readily available production data. In the first technique, acceleration and incremental recovery, and thus infill drilling, are predicted from the trend of the cumulative production (Gp) versus square root time function. This approach is more applicable for tight formations considering the long period of transient linear flow. The second technique is based on multi-well Blasingame type curves analysis. This technique appears to best be applied when the production of parent wells reaches the boundary dominated flow (BDF) region before the production start of the successive infill wells. These techniques are important in field development planning as the flow regimes in tight formations change gradually from transient flow (early times) to BDF (late times) as the production continues. Despite different approaches/methods, the field case studies demonstrate that the accurate framework for strategic well planning including prediction of optimum well location is very critical, especially for the realization of the commercial benefit (i.e., increasing and accelerating of reserve or assets) from infilled drilling campaign. Also, the proposed framework and findings of this study provide new insight into infilled drilling campaigns including the importance of better evaluation of infill drilling performance in tight formations, which eventually assist on informed decisions process regarding future development plans.


2020 ◽  
Vol 10 (8) ◽  
pp. 3333-3345
Author(s):  
Ali Al-Rubaie ◽  
Hisham Khaled Ben Mahmud

Abstract All reservoirs are fractured to some degree. Depending on the density, dimension, orientation and the cementation of natural fractures and the location where the hydraulic fracturing is done, preexisting natural fractures can impact hydraulic fracture propagation and the associated flow capacity. Understanding the interactions between hydraulic fracture and natural fractures is crucial in estimating fracture complexity, stimulated reservoir volume, drained reservoir volume and completion efficiency. However, because of the presence of natural fractures with diffuse penetration and different orientations, the operation is complicated in naturally fractured gas reservoirs. For this purpose, two numerical methods are proposed for simulating the hydraulic fracture in a naturally fractured gas reservoir. However, what hydraulic fracture looks like in the subsurface, especially in unconventional reservoirs, remain elusive, and many times, field observations contradict our common beliefs. In this study, the hydraulic fracture model is considered in terms of the state of tensions, on the interaction between the hydraulic fracture and the natural fracture (45°), and the effect of length and height of hydraulic fracture developed and how to distribute induced stress around the well. In order to determine the direction in which the hydraulic fracture is formed strikethrough, the finite difference method and the individual element for numerical solution are used and simulated. The results indicate that the optimum hydraulic fracture time was when the hydraulic fracture is able to connect natural fractures with large streams and connected to the well, and there is a fundamental difference between the tensile and shear opening. The analysis indicates that the growing hydraulic fracture, the tensile and shear stresses applied to the natural fracture.


2016 ◽  
Vol 26 (3/4) ◽  
pp. 1187-1225 ◽  
Author(s):  
Nicola Massarotti ◽  
Michela Ciccolella ◽  
Gino Cortellessa ◽  
Alessandro Mauro

Purpose – The purpose of this paper is to focus on the numerical analysis of transient free convection heat transfer in partially porous cylindrical domains. The authors analyze the dependence of velocity and temperature fields on the geometry, by analyzing transient flow behavior for different values of cavity aspect ratio and radii ratio; both inner and outer radius are assumed variable in order to not change the difference ro-ri. Moreover, several Darcy numbers have been considered. Design/methodology/approach – A dual time-stepping procedure based on the transient artificial compressibility version of the characteristic-based split algorithm has been adopted in order to solve the transient equations of the generalized model for heat and fluid flow through porous media. The present model has been validated against experimental data available in the scientific literature for two different problems, steady-state free convection in a porous annulus and transient natural convection in a porous cylinder, showing an excellent agreement. Findings – For vertically divided half porous cavities, with Rayleigh numbers equal to 3.4×106 for the 4:1 cavity and 3.4×105 for the 8:1 cavity, the numerical results show that transient oscillations tend to disappear in presence of cylindrical geometry, differently from what happens for rectangular one. The magnitude of this phenomenon increases with radii ratio; the porous layer also affects the stability of velocity and temperature fields, as oscillations tend to decrease in presence of a porous matrix with lower value of the Darcy number. Research limitations/implications – A proper analysis of partially porous annular cavities is fundamental for the correct estimation of Nusselt numbers, as the formulas provided for rectangular domains are not able to describe these problems. Practical implications – The proposed model represents a useful tool for the study of transient natural convection problems in porous and partially porous cylindrical and annular cavities, typical of many engineering applications. Moreover, a fully explicit scheme reduces the computational costs and ensures flexibility. Originality/value – This is the first time that a fully explicit finite element scheme is employed for the solution of transient natural convection in partially porous tall annular cavities.


2018 ◽  
Vol 81 (1) ◽  
Author(s):  
Shiferaw Regassa Jufar ◽  
Tareq M Al-Shami ◽  
Ulugbek Djuraev ◽  
Berihun Mamo Negash ◽  
Mohammed Mahbubur Rahman

A numerical simulation of flow of oil blob through a capillary tube constriction is presented. The simulation was run in a 2D axisymmetric model. Water is injected at the inlet to mobilize oil blob placed near the capillary tube constriction. Transient flow images were used to understand the flow evolution process. Results from the study show that pulsed injection effectively assisted to squeeze out the oil blob through the capillary tube constriction with shorter time compared to continuous injection.  Pulsed injection reduced the time required for the first droplet to cross the capillary tube constriction by about 3 folds compared to continuous injection. In addition, the droplet that crossed the constriction is larger when the flow was pulsed. In both cases, there is a reverse flow in the opposite direction of the injection. However, the severity of the reverse flow is stronger in the case of continuous injection. Immediately downstream the constriction, there is an adverse pressure gradient zone during continuous injection which limits the mobility of droplet that crossed the constriction. However, in the case of pulsed injection, there is a favorable pressure gradient zone immediately downstream the constriction. This zone expedites mobility of droplets that cross the constriction by transporting them further downstream through suction effect. Apparently, pulsed injection eases off the adverse pressure gradient and allowed more volume of oil to pass through the constriction. Within about two periods of pulsation, 84% of original oil placed at the beginning crossed the constriction compared to only 35% in the case of continuous injection. Even though the same amount of water was injected in both cases, pulsed injection clearly altered the flow behavior. The observation from this study may be extended to more complex flows in order to tailor the method for certain specific applications, such as flow of residual oil through a reservoir.


2014 ◽  
Vol 105 (3) ◽  
pp. 539-558 ◽  
Author(s):  
Lei Wang ◽  
Xiaodong Wang ◽  
Erhui Luo ◽  
Junlei Wang

2018 ◽  
Vol 2018 ◽  
pp. 1-20 ◽  
Author(s):  
Mingxian Wang ◽  
Zifei Fan ◽  
Xuyang Dong ◽  
Heng Song ◽  
Wenqi Zhao ◽  
...  

This study develops a mathematical model for transient flow analysis of acid fracturing wells in fractured-vuggy carbonate reservoirs. This model considers a composite system with the inner region containing finite number of artificial fractures and wormholes and the outer region showing a triple-porosity medium. Both analytical and numerical solutions are derived in this work, and the comparison between two solutions verifies the model accurately. Flow behavior is analyzed thoroughly by examining the standard log-log type curves. Flow in this composite system can be divided into six or eight main flow regimes comprehensively. Three or two characteristic V-shaped segments can be observed on pressure derivative curves. Each V-shaped segment corresponds to a specific flow regime. One or two of the V-shaped segments may be absent in particular cases. Effects of interregional diffusivity ratio and interregional conductivity ratio on transient responses are strong in the early-flow period. The shape and position of type curves are also influenced by interporosity coefficients, storativity ratios, and reservoir radius significantly. Finally, we show the differences between our model and the similar model with single fracture or without acid fracturing and further investigate the pseudo-skin factor caused by acid fracturing.


Sign in / Sign up

Export Citation Format

Share Document