Analytical Rate-Transient Analysis and Production Performance of Waterflooded Fields with Delayed Injection Support

2021 ◽  
pp. 1-23
Author(s):  
Daniel O'Reilly ◽  
Manouchehr Haghighi ◽  
Mohammad Sayyafzadeh ◽  
Matthew Flett

Summary An approach to the analysis of production data from waterflooded oil fields is proposed in this paper. The method builds on the established techniques of rate-transient analysis (RTA) and extends the analysis period to include the transient- and steady-state effects caused by a water-injection well. This includes the initial rate transient during primary production, the depletion period of boundary-dominated flow (BDF), a transient period after injection starts and diffuses across the reservoir, and the steady-state production that follows. RTA will be applied to immiscible displacement using a graph that can be used to ascertain reservoir properties and evaluate performance aspects of the waterflood. The developed solutions can also be used for accurate and rapid forecasting of all production transience and boundary-dominated behavior at all stages of field life. Rigorous solutions are derived for the transient unit mobility displacement of a reservoir fluid, and for both constant-rate-injection and constant-pressure-injection after a period of reservoir depletion. A simple treatment of two-phase flow is given to extend this to the water/oil-displacement problem. The solutions are analytical and are validated using reservoir simulation and applied to field cases. Individual wells or total fields can be studied with this technique; several examples of both will be given. Practical cases are given for use of the new theory. The equations can be applied to production-data interpretation, production forecasting, injection-water allocation, and for the diagnosis of waterflood-performanceproblems. Correction Note: The y-axis of Fig. 8d was corrected to "Dimensionless Decline Rate Integral, qDdi". No other content was changed.

SPE Journal ◽  
2013 ◽  
Vol 18 (04) ◽  
pp. 795-812 ◽  
Author(s):  
C.R.. R. Clarkson ◽  
J.D.. D. Williams-Kovacs

Summary Early fluid production and flowing pressure data gathered immediately after fracture stimulation of multifractured horizontal wells may provide an early opportunity to generate long-term forecasts in shale-gas (and other hydraulically fractured) reservoirs. These early data, which often consist of hourly (if not more frequent) monitoring of fracture/formation fluid rates, volumes, and flowing pressures, are gathered on nearly every well that is completed. Additionally, fluid compositions may be monitored to determine the extent of load fluid recovery, and chemical tracers added during stage treatments to evaluate inflow from each of the stages. There is currently debate within the industry of the usefulness of these data for determining the long-term production performance of the wells. “Rules of thumb” derived from the percentage of load-fluid recovery are often used by the industry to provide a directional indication of well performance. More-quantitative analysis of the data is rarely performed; it is likely that the multiphase-flow nature of flowback and the possibility of early data being dominated by wellbore-storage effects have deterred many analysts. In this work, the use of short-term flowback data for quantitative analysis of induced-hydraulic-fracture properties is critically evaluated. For the first time, a method for analyzing water and gas production and flowing pressures associated with the flowback of shale-gas wells, to obtain hydraulic-fracture properties, is presented. Previous attempts have focused on single-phase analysis. Examples from the Marcellus shale are analyzed. The short (less than 48 hours) flowback periods were followed by long-term pressure buildups (approximately 1 month). Gas + water production data were analyzed with analytical simulation and rate-transient analysis methods designed for analyzing multiphase coalbed-methane (CBM) data. This analogy is used because two-phase flowback is assumed to be similar to simultaneous flow of gas and water during long-term production through the fracture system of coal. One interpretation is that the early flowback data correspond to wellbore + fracture volume depletion (storage). It is assumed that fracture-storage volume is much greater than wellbore storage. This flow regime appears consistent with what is interpreted from the long-term pressure-buildup data, and from the rate-transient analysis of flowback data. Assuming further that the complex fracture network created during stimulation is confined to a region around perforation clusters in each stage, one can see that fluid-production data can be analyzed with a two-phase tank-model simulator to determine fracture permeability and drainage area, the latter being interpreted to obtain an effective (producing) fracture half-length given geometrical considerations. Total fracture half-length, derived from rate-transient analysis of online (post-cleanup) data, verifies the flowback estimates. An analytical forecasting tool that accounts for multiple sequences of post-storage linear flow, followed by late-stage boundary flow, was developed to forecast production with flowback-derived parameters, volumetric inputs, matrix permeability, completion data, and operating constraints. The preliminary forecasts are in very good agreement with online production data after several months of production. The use of flowback data to generate early production forecasts is therefore encouraging, but needs to be tested for a greater data set for this shale play and for other plays, and should not be used for reserves forecasting.


Energies ◽  
2021 ◽  
Vol 14 (16) ◽  
pp. 4887
Author(s):  
Suyang Zhu ◽  
Alireza Salmachi

Two phase flow and horizontal well completion pose additional challenges for rate-transient analysis (RTA) techniques in under-saturated coal seam gas (CSG) reservoirs. To better obtain reservoir parameters, a practical workflow for the two phase RTA technique is presented to extract reservoir information by the analysis of production data of a horizontal well in an under-saturated CSG reservoir. This workflow includes a flowing material balance (FMB) technique and an improved form of two phase (water + gas) RTA. At production stage of a horizontal well in under-saturated CSG reservoirs, a FMB technique was developed to extract original water in-place (OWIP) and horizontal permeability. This FMB technique involves the application of an appropriate productivity equation representing the relative position of the horizontal well in the drainage area. Then, two phase (water + gas) RTA of a horizontal well was also investigated by introducing the concept of the area of influence (AI), which enables the calculation of the water saturation during the transient formation linear flow. Finally, simulation and field examples are presented to validate and demonstrate the application of the proposed techniques. Simulation results indicate that the proposed FMB technique accurately predicts OWIP and coal permeability when an appropriate productivity equation is selected. The field application of the proposed methods is demonstrated by analysis of production data of a horizontal CSG well in the Qinshui Basin, China.


SPE Journal ◽  
2020 ◽  
Vol 25 (04) ◽  
pp. 1636-1656
Author(s):  
Qian Sun ◽  
Luis F. Ayala

Summary Considerable research has been focused on the development of rate-transient-analysis (RTA) models to estimate the reserves of gas/condensate reservoirs. Currently, broadly deployed RTA tools rely on multiphase pseudopressure concepts to enable multiphase production-data analysis. In any multiphase pseudopressure calculation, the determination of the saturation/pressure (So/p) relationship plays a vital role because it directly influences the ability of multiphase RTA methods to reliably forecast original gas in place (OGIP). In this work, we present a thermodynamics-based So/p model that provides a better understanding of the phase behavior for the boundary-dominated gas/condensate reservoirs. The proposed So/p model is derived from the thermodynamic nature of extended black-oil formulations. A noniterative flash-calculation protocol is used to establish the So/p path in the condensate-buildup region. The developed method can be coupled with RTA tools and services for the calculation of multiphase pseudopressure. In this work, we present case studies of three gas/condensate reservoirs with different types of fluids. Two RTA multiphase analysis models are used to scrutinize the production data using the newly proposed So/p relationship, and results are compared with the use of a traditional steady-state method coupled with constant-volume-depletion (CVD) data. Results of the case studies show that RTA models that use the proposed So/p consistently yield more accurate OGIP estimation. Thus, this work presents a practical approach to remove commonly used yet potentially faulty assumptions in multiphase RTA applications for liquid-rich gas/condensate reservoirs.


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):  
Ruslan Rubikovich Urazov ◽  
Alfred Yadgarovich Davletbaev ◽  
Alexey Igorevich Sinitskiy ◽  
Ilnur Anifovich Zarafutdinov ◽  
Artur Khamitovich Nuriev ◽  
...  

Abstract This research presents a modified approach to the data interpretation of Rate Transient Analysis (RTA) in hydraulically fractured horizontal well. The results of testing of data interpretation technique taking account of the flow allocation in the borehole according to the well logging and to the injection tests outcomes while carrying out hydraulic fracturing are given. In the course of the interpretation of the field data the parameters of each fracture of hydraulic fracturing were selected with control for results of well logging (WL) by defining the fluid influx in the borehole.


SPE Journal ◽  
2020 ◽  
Vol 25 (06) ◽  
pp. 3280-3299
Author(s):  
Hongyang Chu ◽  
Xinwei Liao ◽  
Zhiming Chen ◽  
W. John John Lee

Summary Because of readily available production data, rate-transient analysis (RTA) is an important method to predict productivity and reserves, and for reservoir and completion characterization in unconventional reservoirs. In addition, multihorizontal well pads are a common development method for unconventional reservoirs. Close well spacing between multifractured horizontal wells (MFHWs) in the multiwell pads makes interference from adjacent MFHWs especially significant. For RTA of production data from multihorizontal well pads, the influence of adjacent MFHWs cannot be ignored. In this work, we propose a semianalytic RTA model for the multihorizontal well pad with arbitrary multiple MFHW properties and starting-production times. Combining Laplace transformation and finite-difference analysis, we obtained a general solution of a multiwell mathematical model to use in RTA. Our model is applicable to cases of multiple MFHWs with different bottomhole pressures (BHPs), varying hydraulic-fracture properties, and different starting-production times. In the solutions, we observe bilinear flow, linear flow, transition flow, and multi-MFHW flow. Rate-normalized pressure (RNP) and its derivative are also affected by multi-MFHW flow. Two case studies revealed that the negative effect of interwell interference on the parent-well productivity is closely related to the pressure distribution caused by the production of child wells.


2014 ◽  
Vol 17 (02) ◽  
pp. 209-219 ◽  
Author(s):  
H.. Luo ◽  
G.F.. F. Mahiya ◽  
S.. Pannett ◽  
P.. Benham

Summary The evaluation of expected ultimate recovery (EUR) for tight gas wells has generally relied upon the Arps equation for decline-curve analysis (DCA) as a popular approach. However, it is typical in tight gas reservoirs to have limited production history that has yet to reach boundary-dominated flow because of the low permeability of such systems. Commingled production makes the situation even more complicated with multiboundary behavior. When suitable analogs are not available, rate-transient analysis (RTA) can play an important role to justify DCA assumptions for production forecasting. The Deep-basin East field has been developed with hydraulically fractured vertical wells through commingled production from multiple formations since 2002. To evaluate potential of this field, DCA type curves for various areas were established according to well performance and geological trending. Multiple-segment DCA methodology demonstrated reasonable forecasts, in which one Arps equation is used to describe the rapidly decreasing transient period in early time and another equation is used for boundary-dominated flow. However, a limitation of this approach is the uncertainty of the forecast in the absence of extended production data because the EUR can be sensitive to adjustments in some assumed DCA parameters of the second segment. In this paper, we used RTA to assess reservoir and fracture properties in multiple layers and built RTA-type well models around which uncertainty analyses were performed. The distributions of the model properties were then used in Monte Carlo analysis to forecast production and define uncertainty ranges for EUR and DCA parameters. The resulting forecasts and EUR distribution from RTA modeling generally support the DCA assumptions used for the type curves for corresponding areas of the field. The study also showed how the contribution from the various commingled layers changes with time. The proposed workflow provides a fit-for-purpose way to quantify uncertainties in tight gas production forecasting, especially for cases when production history is limited and field-level numerical simulation is not practicable.


2021 ◽  
Author(s):  
Michael B. Vasquez ◽  
Pedro M. Adrian

Abstract Analysis of modern production data also known as Rate Transient Analysis (RTA) is a technique to perform reservoir characterization using the combination of bottomhole flowing pressure and flow rate data without the need to close wells. These methods allow the estimation of the Hydrocarbon Initially In-Place (HIIP), production forecast and main reservoir parameters. Several RTA methods have already been developed to analyze different reservoir models such as homogeneous, naturally fractured, geopressurized, hydraulically fractured, however, in the case of layered reservoirs the studies are almost null although there are several studies conducted in the area of pressure transient analysis. This paper presents the analytical derivation of the Palacio-Blasingame type curves to analyze production data of a two-layered reservoir model without crossflow or hydraulic communication between them. A new set of type curves were generated by applying the Gaver Stehfest algorithm with Matlab to achieve the solution of the inverse of the Laplace space considering a constant flow of production flow and a flow regime in the radial pseudosteady-state, then applying the definitions dimensionless the proposed method was derived. Synthetic data were generated with a commercial simulator to validate the method. Furthermore this paper presents a field case study application. The results were compared to the type curve for homogenous reservoirs, volumetric method as well as well testing results. Results confirmed the applicability of rate transient analysis technique in a two-layered reservoir without crossflow with a single drainage area and the same initial pressure for all layers (same pressure gradient of formation), and different values of thickness of the layers, permeability and porosity.


SPE Journal ◽  
2019 ◽  
Vol 25 (01) ◽  
pp. 369-390 ◽  
Author(s):  
Yang Wang ◽  
Luis F. Ayala

Summary Current rate-transient-analysis tools for gas wells producing under boundary-dominated-flow (BDF) conditions largely rely on the deployment of the Arps empirical decline models (Arps 1945), or liquid-based analytical models rewritten in terms of pseudofunctions. Recently, Stumpf and Ayala (2016) demonstrated that, contrary to common practice, decline exponents (b) used in Arps’ hyperbolic equations when applied to gas-well analysis can be rigorously estimated before any field-production data are collected. This determination is solely dependent on gas pressure/volume/temperature (PVT) properties and prevailing constant-bottomhole-pressure (BHP) specification for volumetric, single-phase gas-flow conditions. In the study, we extend that work to a more-realistic variable-BHP condition, which is the most common production-specification condition, in terms of the ratio of changing BHP to average reservoir pressure. The decline exponent (b) is thus rederived, and it is shown that under such conditions, variable BHP hyperbolic decline coefficients become solely dependent on fluid PVT properties and take their largest possible magnitude compared with constant-BHP production. Step-by-step analysis procedures are presented that enable explicit and straightforward estimation of original gas in place (OGIP) and other reservoir properties by universal-type-curve and straight-line analysis. Finally, several cases using simulated and field data are discussed in detail to validate the capabilities of the proposed approach.


Sign in / Sign up

Export Citation Format

Share Document