(Formation Evaluation and the Analysis of Reservoir Performance) Module for: Analysis of Reservoir Performance Introduction T.A. Blasingame, Texas A&M U. Department of Petroleum Engineering Texas A&M University College Station, TX 77843-3116 (979) 845-2292 t-blasingame@tamu.edu Introduction Analysis of Reservoir Performance Slide 1
Executive Summary Module Summary The overall objective of this module is to familiarize the student with the modern methods for the analysis, interpretation, and modelling (rate/pressure prediction) of hydrocarbon reservoir systems. Module Structure: Module 1: Introductory Concepts History of data analysis, phase behavior, flow concepts. Module 2: Material Balance/Late-Time Flow Behavior Material balance methods and gas flow models. Module 3: Wellbore Phenomena/Near-Well Behavior Module 4: Well Test Analysis Deliverability and well test analysis, well test design. Module 5: Analysis of Production Data Data acquisition, decline curves, EUR, decline type curves. Introduction Analysis of Reservoir Performance Slide 2
Module 1: Introductory Concepts History of Data Analysis Correlation plots cumulative-ip, rate-time, etc. Rate-time correlations (log-log, semilog, etc.). Arps rate-time relations (exponential/hyperbolic). Phase Behavior Gas z-factor (Standing-Katz correlation plot(s)). Gas compressibility. Gas viscosity. Fundamentals of Fluid Flow in Porous Media General form of the gas diffusivity equation. Pseudopressure-pseudotime form. Pressure-squared form. Pseudotime and pressure-squared criteria plots. Introduction Analysis of Reservoir Performance Slide 3
Module 1: History of Data Analysis Data Plots a. From: Manual for the Oil and Gas Industry Arnold (1919). Production decline analysis: Over 80 years old! Objective was economic, not technical production extrapolations were even referenced to the tax year! Very humble origins "whatever worked" plots seemed to be popular (e.g., Cartesian, log-log, and semilog). b.from: Estimation of Underground Oil Reserves by Oil-Well Production Curves Cutler (1924). Introduction Analysis of Reservoir Performance Slide 4
Module 1: History of Data Analysis Rate Plots a. The "engineer's solution" (i.e., the log-log plot this plot did not stand the test of time). b.the "gee it works" plot "I wonder if there is some theory?"... (yes). From: Estimation of Underground Oil Reserves by Oil-Well Production Curves Cutler (1924). c. The "scratch your head" plot... interesting, but... how does it work? Introduction Analysis of Reservoir Performance Slide 5
Module 1: History of Data Analysis Arps Relations a. From: SPE-Transactions Arps (1944). "Arps" rate decline analysis: Introduction of exponential and hyperbolic families of "decline curves" (Arps, 1944) Introduction of log-log "type curve" for the "Arps" family of "decline curves" (Fetkovich, 1973). Empirical... but seems to work as a general tool. Is this more coincidence or theory? (... theory) b.from: SPE 04629 Fetkovich (1973). Introduction Analysis of Reservoir Performance Slide 6
Module 1: Phase Behavior z-factor (SK-data) a. SK base plot (z vs. p pr ) Poettmann-Carpenter Data (5960 data points). b. SK "revised" plot (z vs. p pr /T pr ) Poettmann- Carpenter Data (5960 data points). c. SK "revised" plot (z vs. ρ pr ) Poettmann- Carpenter Data (5960 data points). z-factor (SK-data): z vs. p pr original Standing-Katz (SK) formulation (data only). z vs. p pr /T pr "modified" Standing- Katz formulation (data only). z vs. ρ pr "modified" Standing-Katz formulation (data only). This serves as the basis for later reduced density correlations (development of an equation of state (EOS)). Introduction Analysis of Reservoir Performance Slide 7
Module 1: Phase Behavior z-factor (DAK-EOS) a. SK base plot (z vs. p pr ) Original Dranchuk- Abou-Kassem data fit (DAK-EOS). b. SK "revised" plot (z vs. p pr /T pr ) Original Dranchuk-Abou-Kassem data fit (DAK-EOS). c. SK "revised" plot (z vs. ρ pr ) Original Dranchuk- Abou-Kassem data fit (DAK-EOS). z-factor (DAK-EOS): z vs. p pr Original Dranchuk-Abou- Kassem data fit (DAK-EOS). z vs. p pr /T pr "modified" Standing- Katz formulation, original Dranchuk- Abou-Kassem data fit (DAK-EOS). z vs. ρ pr "modified" Standing-Katz formulation, original Dranchuk- Abou-Kassem data fit (DAK-EOS). Introduction Analysis of Reservoir Performance Slide 8
Module 1: Phase Behavior Gas Compressibility a. Definition of gas compressibility. b. Definition of reduced gas compressibility. d. "Reduced compressibility" plot (c r vs. p pr, 1.4< T pr <3.0) Mattar, Brar, and Aziz (1975). c. "Reduced compressibility" plot (c r vs. p pr, 1.05< T pr <1.4) Mattar, Brar, and Aziz (1975). Gas Compressibility "Reduced compressibility" concept is used to correlate data. c r function is computed using the Dranchuk-Abou-Kassem EOS convenient analytical form. Introduction Analysis of Reservoir Performance Slide 9
Module 1: Phase Behavior Gas Viscosity a. Jossi, et al plot pure component data. Note the correlation with ρ r. b.jossi, et al plot various data. Note the corelation with ρ r. c. Jossi, et al plot database. Note the temperature (T r ) dependance. Correlations must be created using ρ r and T r. From: The Viscosity of Pure Substances in the Dense Gaseous and Liquid Phases Jossi, Stiel, and Thodos (1962). Introduction Analysis of Reservoir Performance Slide 10
Module 1: Phase Behavior Gas Viscosity a. Jossi, et al revised correlation for gas viscosity. c. New correlation for gas viscosity. New Gas Viscosity Correlation: b. Lee, et al revised correlation for gas viscosity. Introduction Analysis of Reservoir Performance Slide 11
Module 1: Fundamentals of Fluid Flow in Porous Media a. General definition of the gas diffusivity equation. d. Illustration of the pressure-squared criteria. b. Pseudopressure-pseudotime form of the gas diffusivity equation. c. Pressure-squared form of the gas diffusivity equation. e. Illustration of the pseudotime criteria. Introduction Analysis of Reservoir Performance Slide 12
Module 2: Material Balance Gas Material Balance Dry gas reservoir systems. Water drive gas reservoir systems. "Abnormally-pressured" gas reservoir systems. Generalized gas material balance equation. Late-Time Flow Behavior Simplified flow behavior (empirical relations). Exponential decline (liquid flow solution). Simplified (approximate) gas flow relation. Introduction Analysis of Reservoir Performance Slide 13
Module 2: Late-Time Flow Behavior a. Simplified flow behavior (empirical relations (Arps)). b. Exponential decline (liquid flow solution, constant p wf case). d. Log-log "type curve" Arps rate relations. c. Simplified (approximate) gas flow relation (valid for p i <6000 psia, constant p wf case). e. Log-log "type curve" simplified (approximate) gas flow relation. Introduction Analysis of Reservoir Performance Slide 14
Module 2: Gas Material Balance Concepts a. Material balance concept for a dry gas reservoir with water influx and pore/water compressibility. c. Illustration of the "pore collapse concept for p/z data (Fetkovich, et al (SPE 22921)). b. Illustration of the "pore collapse concept (Fetkovich, et al (SPE 22921)). d. Field analysis using generalized material balance relation (Fetkovich, et al (SPE 22921)). Introduction Analysis of Reservoir Performance Slide 15
Module 2: Gas Material Balance Relations a. Material balance relation for a dry gas reservoir. b. Generalized material balance relation (Dake). c. Generalized material balance relation (Fetkovich, et al (SPE 22921)). (Warning) Do not use statistical methods (e.g., regression analysis) as the sole mechanism to solve material balance problems physically inconsistent parameter estimates are likely. Introduction Analysis of Reservoir Performance Slide 16
Module 3: Wellbore Phenomena/Near-Well Behavior Calculation of Bottomhole Pressures General relation (energy balance) Static (non-flowing) bottomhole pressure (dry gas). Flowing bottomhole pressure (dry gas). Near-Well Reservoir Flow Behavior Steady-state "skin factor" concept used to represent damage or stimulation in the near-well region. "Variable" skin effects: non-darcy flow, well cleanup, and gas condensate banking. Introduction Analysis of Reservoir Performance Slide 17
Module 3: Calculation of Bottomhole Pressure a. Basic energy balance for flow in inclined pipes. b. Solution assuming average T and z-values. d. Schematic illustration of wellbore configuration. c. Cullender-Smith solution of the energy balance. e. Wellbore diagram surface pressure measurement. Introduction Analysis of Reservoir Performance Slide 18
Module 3: Near-Well Behavior Skin Concept a. Simple skin concept steady-state flow of a liquid in the "altered" zone. b. Governing relations for steady-state flow of a liquid in the "altered" zone (k s =permeability in the "altered" zone). c. Schematic pressure behavior in a "steady-state" skin zone. Introduction Analysis of Reservoir Performance Slide 19
Module 3: Extensions of the Skin Concept a. Concepts of "infinitesimal" skin as well as use of the effective wellbore radius these schematics seek to represent the concept of "skin" as a physical phenomena, as well as a (simple) mathematical model. b. 2-zone, radial "composite" model (condensate bank). Skin Concept Generally speaking, use of the skin concept (or skin "factor") tends to "isolate" pressure behavior that cannot be directly attributed to the reservoir. This is an oversimplification, but it is convenient, and is widely used. Introduction Analysis of Reservoir Performance Slide 20
Module 4: Well Test Analysis Orientation This module will focus specifically on the analysis and interpretation of deliverability test data and pressure transient test data. The issues must be clear: test design, data acquisition/data quality control, and test execution are critical activities. Deliverability Testing: Keep it simple a "4-point" test is appropriate. Isochronal testing is very difficult to implement. Pressure Transient Test Analysis/Interpretation: Conventional analysis is consistent/appropriate. Model identification (log-log (or type curve) analysis). Test design keep it simple. Introduction Analysis of Reservoir Performance Slide 21
Module 4: Deliverability Testing Basics a. "Standard" 4-point test deliverability test note that the rates increase (to protect the reservoir). c. Modified "Isochronal" test sequence note that each "buildup" is not required to achieve p i. b. "Isochronal" test sequence note that each "buildup" is required to achieve p i. d. Governing equations for "deliverability" test analysis/interpretation. Introduction Analysis of Reservoir Performance Slide 22
Module 4: Deliverability Testing Orientation a. Basic "pressure-squared" relation that is presumed to describe gas flow analogous form can be derived from steady-state flow theory (Darcy's law). c.traditional "deliverability" plot probably derived from empirical plotting of data. b."rate-squared" (or velocity-squared) formulation analogous form can be derived from steady-state flow theory (Forchheimer Eq.). d. Modified "deliverability" plot note that bq sc 2 must be known (... need alternative approach). Introduction Analysis of Reservoir Performance Slide 23
Module 4: Multirate Testing Example Case a. Multirate (4-point) rate sequence (note pressure match (solid trend through the data). b. Log-log "summary plot" note good agreement in comparison of data and model. c. Results summary note that non-darcy flow, changing wellbore storage, and an infiniteacting reservoir system were considered in this analysis. Introduction Analysis of Reservoir Performance Slide 24
Module 4: "Well Interference" Example Case a. "Well Interference" plot note the linear trend through the data functions (confirms interference). b. Log-log "summary plot" note the corrected and uncorrected data (well interference). c. Horner semilog plot note the two semilog trends confirm the radial composite model. Discussion: "Well interference" is much more common than previously thought and we must recognize the characteristic behavior on each plot: Log-log plot (b) Semilog plot (c) Specialized plot (a) Introduction Analysis of Reservoir Performance Slide 25
Module 4: Pressure Transient Testing Basics a. Log-log "preliminary analysis" plot wellbore storage and radial flow (C s, k). c. Semilog "middle-time" plot used to analyze radial flow behavior (k, s). e. Cartesian "Arps" plot used to estimate average reservoir pressure. b. Cartesian "early-time" plot used to analyze wellbore storage (p 0, C s ). d. Horner "middle-time" plot used to analyze radial flow behavior (k, s, p*). f. Log-log "summary" plot summary of all analysis (C s, k, s, A, etc). Introduction Analysis of Reservoir Performance Slide 26
Well Test Analysis WBS Type Curves a. Type Curve: Radial flow with wellbore storage and skin effects (p D, p Dd ). b. Type Curve: Radial flow with wellbore storage and skin effects (p D, p Dd, p Dr1 ). c. Type Curve: Radial flow with wellbore storage and skin effects (p D, p Ddd ). d. Type Curve: Radial flow with wellbore storage and skin effects (p Di, p Did ). e. Type Curve: Radial flow with wellbore storage and skin effects (p Di, p Did, p Dir1 ). f. Type Curve: Radial flow with wellbore storage and skin effects (p Di, p Dir2 ). Introduction Analysis of Reservoir Performance Slide 27
Well Test Analysis Bounded Reservoir a. Type Curve for sealing faults (p Dd ). b. Type Curve for conductive (leaky) faults (p Dd ). c. Type Curve for pressure buildup test in a closed rectangular reservoir (p Dd ). d. Type Curve for pressure buildup test in a closed rectangular reservoir (p Did ). Introduction Analysis of Reservoir Performance Slide 28
Well Test Analysis Composite Systems a. Composite Reservoir (η r =1x10-3 ). b. Composite Reservoir (η r =1x10-2 ). c. Composite Reservoir (η r =1x10-1 ). d. Composite Reservoir (η r =1x10 0 ). e. Composite Reservoir (all η r cases). Introduction Analysis of Reservoir Performance Slide 29
Well Test Analysis Fractured Wells a. Type Curve: C fd = various, no C Df cases. b. Type Curve: C fd =1, C Df = various. c. Type Curve: C fd =2, C Df = various. d. Type Curve: C fd =5, C Df = various. e. Type Curve: C fd =10, C Df = various. f. Type Curve: C fd =1x10 3, C Df = various. Introduction Analysis of Reservoir Performance Slide 30
Module 5: Analysis of Production Data Production data "low frequency" (taken at large (or even random) intervals) and "low resolution" (data quality (i.e., accuracy) is minimal). Simple Analysis: Rate-time decline curve analysis. EUR analysis (usually rate-cumulative or variation). New stuff advanced analysis based on a simple, yet robust model (e.g., Knowles q g -G p analysis) Decline type curve analysis: Systematic, model-based analysis approach. Model identification transient data analysis. Volume estimates pseudosteady-state behavior is dictated by material balance (very consistent). Introduction Analysis of Reservoir Performance Slide 31
Module 5: Production Analysis Example Case Rate and pressure profile for a mid-continent (U.S.) gas well, note the daily and seasonal fluctuations in the data. Introduction Analysis of Reservoir Performance Slide 32
Module 5: Production Analysis EUR Analysis Estimated ultimate recovery (EUR) profile for a mid-continent (U.S.) gas well there is considerable variation in the data. Introduction Analysis of Reservoir Performance Slide 33
Module 5: Production Analysis WPA Approach Note that the WPA approach provides a unique analysis/ interpretation of the well performance history. Introduction Analysis of Reservoir Performance Slide 34
Module 5: Production Analysis the hard way... "Van Everdingen-Meyer Method: "Analysis by simulation" (use analytical solution to define x- axis plotting function). Considers all of the data, needs a complete model to generate an appropriate analysis/interpretation. Theoretically simple, practical. Pro: Theoretically simple and practical (can use field data). Con: Limited by solution model as well as data quality. From: SPE 15482 Whitson and Sognesand (1988). Introduction Analysis of Reservoir Performance Slide 35
Module 5: Production Analysis History Lessons From: SPE 04629 Fetkovich (1973). Composite Transient Type Curve: Collapses the transient flow trends into "stems" related to reservoir size and skin factor (Fetkovich, 1973). Composite Total Type Curve: Addition of the "Arps" empirical trends for "boundary-dominated flow behavior (Fetkovich, 1973)." Assumptions: Constant bottomhole pressure. "Liquid" flow (not gas). From: SPE 04629 Fetkovich (1973). Introduction Analysis of Reservoir Performance Slide 36
Module 5: Production Analysis History Lessons From: SPE 28688 Doublet, et al (1994). Fetkovich Derivative Type Curve: Good concept, but just try to take the derivative of production data... Fetkovich-McCray Type Curve: Concept is to generate "integral" functions for data analysis, much better performance than simply using rate. Still Need: Variable pressure/rate methods. Other models fractured wells, horizontal wells, etc... From: SPE 25909 Palacio, et al (1993). Introduction Analysis of Reservoir Performance Slide 37
Module 5: Production Analysis History Lessons From: SPE 25909 Palacio, et al (1993). UNFRACTURED Well Case Variable Rate/Pressure Approach: Use "material balance time" (xaxis) and "pressure drop normalized rate" (y-axis) functions. Good news: New concept provides unique behavior during boundarydominated flow regime. Not-So-Good-News: Wellbore pressure data are critical. From: SPE 28688 Doublet, et al (1994). Introduction Analysis of Reservoir Performance Slide 38
Module 5: Production Analysis History Lessons a. Doublet, et al (1996): "Fetkovich-McCray" format INFINITE conductivity vertical fracture (F cd = ). c. Pratikno (2002): "Fetkovich-McCray" format FINITE conductivity vertical fracture (F cd =0.5). b. Pratikno (2002): "Fetkovich-McCray" format FINITE conductivity vertical fracture (F cd =10). Decline Type Curves: Fractured Wells Infinite fracture conductivity: Less complex solution, but somewhat ideal for use in practice. Finite fracture conductivity: F cd =10: Moderate to high fracture conductivity case. F cd =0.5: Low fracture conductivity case. Introduction Analysis of Reservoir Performance Slide 39
Module 5: Production Analysis History Lessons a. INFINITE conductivity vertical fracture (F cd = ) (1996). b. FINITE conductivity vertical fracture (F cd =0.1). c. FINITE conductivity vertical fracture (F cd =1). d. FINITE conductivity vertical fracture (F cd =10). e. FINITE conductivity vertical fracture (F cd =100). f. FINITE conductivity vertical fracture (F cd =1000). Introduction Analysis of Reservoir Performance Slide 40
Module 5: Production Analysis History Lessons Agarwal, et al Methodology: Basically the same as Blasingame, et al work. More like pressure transient test analysis/interpretation. From: SPE 57916 Agarwal, et al (1998). From: SPE 57916 Agarwal, et al (1998). From: SPE 57916 Agarwal, et al (1998). Introduction Analysis of Reservoir Performance Slide 41
Module 5: Production Analysis History Lessons MULTIWELL Analysis Multiwell case can be "recast" into single well case using cumulative production for entire field. Homogeneous reservoir example shows that all cases (9 wells) align same behavior observed for heterogeneous reservoir cases. From: SPE 71517 Marhaendrajana (2001). From: SPE 71517 Marhaendrajana (2001). From: SPE 71517 Marhaendrajana (2001). Introduction Analysis of Reservoir Performance Slide 42
Module 5: Production Analysis New Stuff (Gas) This work presents an analysis and interpretation sequence for the estimation of reserves in a volumetric dry-gas reservoir. This is based on the "Knowles" ratecumulative production relation for pseudosteady-state gas flow given as: qg = qgi 1 pwf pi 2qgi / / zwf zi G 2 p G + 1 qgi pwf / zwf pi / zi 2 2 G 2 Gp "Knowles" relations for gas flow: q g G p follows quadratic "rate-cumulative" relation. Approximation valid for p i <6000 psia. Assumes p wf = constant. Introduction Analysis of Reservoir Performance Slide 43
Module 5: Production Analysis New Stuff (Gas) a. West Virginia Well A: log(q g ) versus log(t) (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). b. West Virginia Well A: log(q g ) versus t (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). c. West Virginia Well A: log(g p ) versus log(t) (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). q g vs. t and G p vs. t: Good data matches (a. and c.) data quality provides clear trends. "High" and "low" q i cases are +/- 10 percent assist in orienting analysis in the spreadsheet. Note that production performance (i.e., rate data) is very consistent (semilog rate-time plot). Introduction Analysis of Reservoir Performance Slide 44
Module 5: Production Analysis New Stuff (Gas) a. q g vs. G p : West Virginia Well A (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). b. q gi,gp vs. G p : West Virginia Well A (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). q g vs. G p : Good data trend (as noted earlier) model fits data quite well. The location of the minimum of the q g vs. G p trend is the gas-inplace (G). This analysis should not be performed using regression regression will favor statistics, rather than the physical problem. q gi,gp vs. G p : Similar to the q q vs. G p plot smoother than q q data. This function serves to validate/confirm the q g vs. G p behavior. The comparison is very clear in this perspective (useful in the distinction of the model trends). Introduction Analysis of Reservoir Performance Slide 45
Module 5: Production Analysis New Stuff (Gas) a. (q gi -q g )/G p vs. G p : West Virginia Well A (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). b. (q gi,gp -q gi )/G p vs. G p : West Virginia Well A (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). c. (q gi,gp -q)/g p vs. G p : West Virginia Well A (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). Comments: (q gi -q g )/G p vs. G p : Extrapolates to 2G good straight-line trend. (q gi,gp -q gi )/G p vs. G p : "Smoother" than the rate function. Extrapolates to 3G note the clarity of "late-time" portion of the trend. (q gi,gp -q g )/G p vs. G p : Reasonably good linear trend extrapolates to 3/2G. Introduction Analysis of Reservoir Performance Slide 46
Module 5: Production Analysis New Stuff (Gas) a. "Quadratic" rate-cumulative production decline type curve. West Virginia Well A (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). b. "Hyperbolic" rate-cumulative production decline type curve. West Virginia Well A (p i = 4175 psia, p wf =710 psia, G Arps =2.14 BSCF). Rate-Cumulative Decline Type Curve Analysis: Note excellent agreement of data/quadratic model trends for this example. "Hyperbolic" analysis appears plausible (as an empirical approach) upon close inspection, the data appear to cross two (2) decline curve stems. Introduction Analysis of Reservoir Performance Slide 47
Module 5: Production Analysis New Stuff (Gas) Rate-Time Decline Type Curve: (Rigorous Analysis) West Virginia Well A (p i = 4175 psia, p wf =710 psia, G quad =3.29 BSCF). Very clear trends all data functions. The estimate of gas reserves using rate-time decline type curve analysis is conservative, possibly/probably a result of the lack of pressure data. Introduction Analysis of Reservoir Performance Slide 48
Module 5: Production Analysis Scaling Production data analyses and pressure transient analyses "see" the reservoir as a volumeaveraged set of properties. New solutions/models will also have this view of the reservoir, but quantifying heterogeneity may be possible by the analysis of data at the "local" level. Scaling will remain a major issue regardless of the mechanism used to analyze reservoir performance. From: Simulator Parameter Assignment and the Problem of Scaling in Reservoir Engineering Halderson (1986). Introduction Analysis of Reservoir Performance Slide 49
Module 5: Production Analysis Future Future of Production Data Analysis: Evolutionary changes: Better data acquisition (major issue). Multiwell analysis (major issue). Better software (major issue). More reservoir/well solutions (minor issue). Revolutionary changes: Direct dataflow into integrated packages for analysis/simulation (5-10 years). Real-time rate-pressure optimization, simultaneous monitoring and control (5-10 years). Introduction Analysis of Reservoir Performance Slide 50
(Formation Evaluation and the Analysis of Reservoir Performance) Module for: Analysis of Reservoir Performance Introduction End of Presentation T.A. Blasingame, Texas A&M U. Department of Petroleum Engineering Texas A&M University College Station, TX 77843-3116 (979) 845-2292 t-blasingame@tamu.edu Introduction Analysis of Reservoir Performance Slide 51