Welltest Help PDF
Welltest Help PDF
Welltest Help PDF
Trademarks
IHS Markit and the IHS Markit logo are trademarks of IHS Markit. Other trademarks appearing in this
publication are the property of IHS Markit or their respective owners.
About WellTest 39
Licensing 40
Best practice 41
Standalone licensing 43
Network licensing 43
Test Types 50
Workflows 61
Conventional Workflow 62
IPR Workflow 64
PITA Workflow 69
Procedures 72
Exporting 75
Exporting to CMG 76
Multiphase Options 81
Starting a Project 82
Performing an Analysis 83
Creating an Analysis 86
Identifying Closure 90
Identifying ISIP 92
Modeling 131
Running Numerical Models Created in WellTest 2013 v2 (v 7.7.0) and Earlier 157
Reporting 233
Pseudo-Pressure 262
References 266
References 268
Reference 271
Summary of Equations for Bilinear Fracture Flow Derivative Analysis (Field Units) 281
Introduction 283
References 284
References 287
Derivative 292
Summary of Equations for Linear Channel Flow Derivative Analysis (Field Units) 294
References 294
Summary of Equations for Linear Fracture Flow Derivative Analysis (Field Units) 298
Derivative 302
Summary of Equations for Linear Channel Flow Derivative Analysis (Field Units) 304
Derivative 307
Summary of Equations for Pseudo-steady State Flow Derivative Analysis (Field Units) 309
Permeability 311
Skin 311
Derivative 313
Permeability 314
Skin 315
Summary of Equations for Radial Flow Derivative Analysis (Field Units) 315
References 318
Introduction 322
Nolte 333
Introduction 349
Theory 350
Analysis 352
References 355
Theory 357
Analysis 359
References 361
References 363
Models 364
Prerequisites 369
Parameters 379
References 381
Parameters 384
Prerequisites 392
References 397
References 398
References 401
References 403
References 407
References 409
References 414
References 417
References 420
References 421
References 424
References 426
References 428
References 430
Tabs 437
Toolbar 446
Export 448
Filter 449
Hotkeys 475
Toolbar 484
Filtering 487
Toolbar 490
Forecasting 498
Toolbar 503
Forecasting 509
Saturations 522
References 551
References 583
Skin 586
References 586
References 592
References 592
References 594
References 596
Undersaturated 612
Reference 612
Undersaturated 614
Undersaturated 615
Heavy-Oils 615
Heavy-Oils 615
Undersaturated 616
Reference 616
Undersaturated 618
Reference 618
Undersaturated 619
Undersaturated 620
Undersaturated 620
Reference 621
Reference 622
Gas-Saturated 622
Undersaturated 623
Reference 623
Undersaturated 624
Undersaturated 624
Reference 625
Undersaturated 625
Reference 626
Undersaturated 627
Undersaturated 628
Reference 628
Undersaturated 631
Reference 631
Field 639
Marquardt-Levenberg 644
Gas 647
References 648
References 658
Minifrac 671
Introduction 671
PAS 682
Recombination 688
Wettability 690
Corey 694
Honarpour 695
Sandstone 695
Limestone 696
Stone I 698
Stone II 698
Baker 699
References 704
Salinity 731
Quality 734
Conclusions 739
References 768
Background 772
Differences between Black Oil, Gas Condensate, and Volatile Oil Models 781
Pseudo-Time 811
References 815
Nomenclature 823
References 835
AOF 839
Deconvolution 841
IPR 843
Pseudo-Pressure 845
Pseudo-Time 845
Contact us 849
Resources 850
WellTest Wizard
Starting a Project
WellTest Wizard: Startup
Test Types
Drawdown / Buildup
Injection / Falloff
Minifrac
PITA & Closed Chamber Tests
Workflows
Conventional
IPR
Minfrac
PITA
Step Down Test
Procedures
Performing an Analysis
Filtering
Generating & Submitting PAS AER Files
Importing Data
Modeling
Preparing Data for Analysis
Reporting
Models
Overview
Advanced Models: Analytical
Advanced Models: Numerical
Legacy Models
User Interface
Tabs
Reference Materials
Calculations & Correlations
General Concepts
Modeling Theory
Time Functions
Deliverability forecasts, at varying line pressures, can be generated for any model, and conversion
of wellhead pressures to sandface conditions is a standard feature. WellTest generates Pressure
ASCII Standard (PAS) files for electronic pressure data submission to the Alberta Energy Regulator
(AER).
Video Tutorials
See the Training Videos web page (opens in a new window) for addi-
tional information.
A standalone license is installed on your local machine. For information on standalone licensing, see
activating / deactivating your license.
A network license is installed on a network server and WellTest must be configured to the network
server. If you want to work offline, a network license can be checked out as a commuter, and is
accounted for in the network license seat count. After you are finished working offline, you can check
your commuter in. For information on setting up network licensing, see the IHS Markit Installation
Guide and IHS Markit License Manager help (both links open in a new window / tab).
To set up a network licensing system, you need to install the License Service on your server and the
License Manager on one or more machines with connectivity to the server. For more information
(including activation / deactivation information), see the IHS Markit Installation Guide (opens in a
new window / tab).
Best practice
For network entitlements, activate / deactivate licenses through the License Manager. For stan-
dalone entitlements, activate / deactivate licenses through WellTest. Do not activate licenses using
the EMS Customer Portal unless instructed to do so by a customer care representative.
1. Open your EMS - Entitlement Certificate email and copy your entitlement ID (EID).
2. Start WellTest.
3. Open the Licensing dialog box by clicking the Help menu and selecting Licensing.
5. Paste your EID (from the email) in the Entitlement Id field and click Connect.
6. To activate your license, in the Activate tab, select the product you want to use and click the
Activate button.
Your license is now active.
Note: You can also queue up your renewed license as soon as you receive your new enti-
tlement email. There is no need to remove your current license, or to wait until your
current license expires.
7. To deactivate your license, click the Deactivate tab, select the product you want to deac-
tivate and click the Deactivate button.
The product is now displayed in the Activate tab and can be activated on a new machine.
Standalone licensing
Active License -- Standalone — indicates that your license is ready for use. Note that if
your license is not active, errors are displayed for troubleshooting purposes.
Configure — opens the Licensing Settings dialog box where you can configure your licens-
ing (for example, you can install or uninstall a license that has already been activated).
Activate / Deactivate — opens the IHS Markit Online Activation dialog box where you can
activate and deactivate your licenses. For more information, see activating / deactivating
licenses.
Network licensing
Check Out Commuter button — enables you to work offline by commuting a network seat.
(This button is only visible if you have not checked out a commuter.)
Days field — use this field to set the number of days your commuter can exist offline. By
default, you can check out the commuter for a maximum of 90 days. Note that this field is
only visible when the Check Out Commuter button is displayed.
Check In Commuter button — click this button when you are finished working offline to
check your commuter back in. (This button is only visible if you have checked out a com-
muter.)
Current Users icon — opens the Current Users dialog box, which lists all current users,
machine names, start dates, and if users' licenses are commuters. (Commuter licenses
enable you to work offline.) This list of current users can be useful if you are trying to track
down a seat that could be made available. If the Commuter column lists "Yes", the user is
likely away from the office.
Configure — opens the Licensing Settings dialog box where you can configure your licens-
ing (for example, you can set the name of your network license server).
Activate / Deactivate — opens the IHS Markit Online Activation dialog box where you can
activate and deactivate your licenses.
At the top, the name of your licensed module is displayed. Beneath this, the following information is
displayed:
Application Version — the build number for WellTest.
Feature Version — indicates whether or not a license is version locked. A version of ALL
indicates that a license is not version locked and can use any version of WellTest com-
patible with the Gemalto licensing system.
License Start Date — the first day the license can be used.
License Expiry Date — the last day the license can be used.
License Library Version — the build number for the license library.
License Service Version — the build number for the license server (only visible for net-
work systems).
Network Licensing
Use this section to specify the location of your network license server. If you are hosting all
of your IHS Markit network licenses (for applications using the Gemalto licensing system) on
the same server, use the Common License Server field. If you have separate servers for
WellTest and other IHS Markit applications, use the WellTest License Server field.
Note: If you fill in both fields, the WellTest License Server field takes precedence.
Installed Licenses — lists all of the standalone licenses currently activated. This section
shows active licenses, as well as inactive licenses (expired or future licenses that are
queued up).
Install button — click this button to install standalone licensing that you have activated using
the customer portal.
Uninstall button — if you want to uninstall standalone licensing on your local computer,
click this button. Note that standalone licenses that are uninstalled with this button are not
available for activation on other computers. If you want to move a standalone license to
another computer, it must first be deactivated, and then activated on the new computer. For
more information, see activating / deactivating licenses.
In the unlikely event of licensing issues, Customer Care may ask you to select one or more
of these logging options. We recommend that these options remain deselected, unless you
are actively investigating an issue.
Entitlement Id — enter the entitlement ID (sent to you via email) and click Connect. This
enables you to activate licenses or deactivate / move licenses from one machine to another.
Remember Me — if you click this checkbox, your entitlement ID is saved, so that you can
quickly reconnect to the EMS
Save Troubleshooting Data button — if you do not see the products you want in the list,
Customer Care may ask you to click this button to generate a text file for troubleshooting pur-
poses. If you need assistance, please contact us. Note that this button only applies to the
Activate tab.
Drawdown / Buildup
Injection / Falloff
Minifrac
Pressure buildup tests are conducted by producing the well, and then shutting it in while measuring
changes in pressure with time. Ideally, the flow preceding the shut-in would at least be long enough
for the rates to stabilize. The quality of data obtained from buildup tests is far superior to that of draw-
down tests because the rates are stable (zero). When analyzing data from buildup tests, it is valu-
able to analyze the drawdown and buildup data.
When the test type is Injection / Falloff, imported rates are converted to negative values, and the
average reservoir pressure on the Properties tab defaults to the lowest measured pressure. Wells
with injection and production, such as gas storage wells, can be analyzed using either Drawdown /
Buildup or Injection / Falloff test types provided the rate history is properly specified.
It is extremely difficult and impractical to obtain estimates of initial pressure in horizontal multifrac
wells following the hydraulic fracture treatments. Minifrac tests are a convenient solution. In com-
parison to Perforation Inflow tests, Minifrac tests are also attractive because they can give estimates
of initial pressure and reservoir flow characteristics within a shorter time frame (provided injection
volume is small), and because the mini-fracture can cut through near wellbore damage and contact
the native reservoir.
The Minifrac test type supports pre-closure analysis, after-closure analysis, modeling, and step-
down test analysis. Minifrac after-closure data may also be analyzed in the same way as falloff data
by changing the test type to Injection / Falloff. However, changing the test type will delete any exist-
ing minifrac analyses. If you would like to have a combination of minifrac analyses and conventional
analyses, we recommend saving the file with a new name prior to changing the test type.
Perforation Inflow
Perforation Injection
Slug
Slug Injection
Perforation Injection tests are analogous to this and occur when the well is perforated in an over-bal-
anced condition. These tests have been referenced as: Slug tests, Surge tests, Perforation Inflow
Diagnostics (PID), or Closed Chamber tests.
Perforation Inflow tests are attractive in unconventional (low-permeability) reservoirs because estim-
ates of initial reservoir pressure (pi), permeability (k), and skin (s') can be obtained in a shorter dur-
ation than buildup tests, and because they don’t require any venting / flaring into the atmosphere.
However, the radius of investigation achieved from these tests is less than that of a buildup test.
In this case, wellbore fluid is injected into the reservoir during the test. The objective is to extract key
reservoir parameters (initial reservoir pressure (pi), permeability (k), and skin (s')) from the pressure
data.
When applying this testing method to formations with moderate to high permeability (1 - 100 md),
estimates of key reservoir properties (initial reservoir pressure (pi), permeability (k), and skin (s'))
can be obtained without shutting in downhole, and resembles a perforation inflow test analysis
(PITA). However, to obtain estimates of key reservoir properties in low permeability formations (< 1
md), downhole shut-in is required to reduce the wellbore volume (Vw) and decrease the time
required to see reservoir-dominated (radial) flow. The downhole shut-in usually occurs after some
formation fluid has been produced into the wellbore. Hence, this type of closed chamber test is com-
prised of two periods: a flow period when the wellbore volume (Vw) is large (PITA), and a buildup
period when the wellbore volume (Vw) is small. To account for this reduction in wellbore storage, we
have developed a different analysis technique, which is referred to as Closed Chamber Test (CCT)
analysis.
Closed Chamber Injection tests refer to the same test conducted when perforating over-balanced,
where the initial cushion pressure (pwo) is greater than the initial reservoir pressure (pi).
Note: In WellTest, the PITA analysis has been developed further than slug.
PITA facilitates the estimation of p*, reservoir permeability, and skin whereas Slug's Impulse Radial
Analysis only provides an estimation of p*. Therefore, PITA is recommended for analyzing data from
slug tests. Both PITA and Slug tests have identical modeling capabilities in WellTest.
For the same reasons described in the Slug Test Type, we recommend using the Perforation Injec-
tion Test Type to analyze slug injection data.
We recommend that you use the wizard. Beside many of the steps in the workflows, you'll see the
Wizard icon ( ), which you can click to jump to the appropriate step in the wizard.
Conventional
IPR
Minifrac
PITA
Step Down
Test Design
1. Depending on the type of data you have, start a Drawdown / Buildup, or Injection / Falloff pro-
ject. See Starting a Project for details.
3. Import your rate data. If you have multiple files containing rate data, such as cleanup flow
and test production, import each file and then merge the rates together.
5. Click the Production Editor tab. The Select Data to Analyze dialog box opens. Ensure that
the imported pressure dataset is displayed in the Active Pressure drop-down list, and that
the imported rate dataset is displayed in the Gas, Oil, or Water Rate drop-down lists; click
OK. See Loading Data into the Production Editor for additional information.
6. Remove unnecessary data points (see Deleting Data at the Start or End of the Production
Editor for details), such as pressure data before and after the test.
7. Specify shut-in points by entering a zero rate in the row of the data table with the first buildup
(or falloff) pressure. Check the data chart to ensure that the rates go to zero immediately
before the pressure buildup (or falloff) begins.
9. Filter the data in the Production Editor. We recommend filtering the data down to fewer than
5,000 points to speed up model calculations.
10. Click the Properties tab and enter the relevant fluid properties.
11. Convert pressures to sandface. If you have downhole gauge data, adjust the pressures
from run depth to sandface conditions. If you have wellhead pressures, calculate sandface
pressures.
13. Create a model. See Creating an Advanced Model or Creating a Legacy Model for addi-
tional information.
1. Launch WellTest and start a new Drawdown / Buildup project.See Starting a Project for
details.
The WellTest Wizard: Startup dialog box opens.
Note: You can also create an IPR by clicking the IPR menu and selecting Create Liquid IPR.
In the Liquid IPR 1 tab, on the left side, there is a section to input values for pressure and
rate parameters. On the right side, there is an empty plot that is generated based on your
inputs.
Note: If you have detailed injection rates, you should import them and synchronize them to
the pressures.
3. Click the Production Editor tab and load the imported data into the Production Editor.
4. If you did not import rates in Step 2, enter the injection rates manually. Most operators report
rates in units of bbl/min. You can change units by right-clicking the units of the Water Rate
column.
5. Enter a zero rate in the row with the first falloff pressure. Check the data chart to ensure the
8. Filter the data. We recommend filtering the data down to fewer than 5,000 points to speed up
model calculations.
a. If KCL water was used, specify the water gravity. This will be used to calculate sandface
pressures.
15. Collate the pre-closure, Nolte after-closure, Soliman after-closure, and model reports.
Then, save the collated report as a pdf file.
1. If the pressures build up (i.e., the pressure in the wellbore prior to perforating was less than
the reservoir pressure), start a new PITA inflow project. If the pressures fall off, start a new
PITA injection project. See Starting a Project for details.
3. Click the Production Editor tab. The Select Data to Analyze dialog box will open. Ensure that
the imported pressure dataset is displayed in the Active Pressure drop-down list and click
OK. See Loading data into the Production Editor for more information.
4. Remove all data prior to perforating and after the build up (or fall off). See Deleting Data at
the Start or End of the Production Editor for additional information.
6. Filter the data. We recommend filtering the data down to fewer than 5,000 points to speed up
model calculations.
7. Click the Properties tab. Depending on the type of reservoir fluid specified earlier, click
either the Gas Properties, Oil Properties, or Water Propertiessub-tabs. Change the
Reservoir Temperature value and enter the fluid properties.
12. Collate the default Late Time, Early Time, and model reports. Then, save the collated
report as a pdf file.
Note: Wellhead pressures must be converted to sandface pressures before conducting an ana-
lysis.
2. Click the Production Editor tab and load the imported data into the Production Editor by
selecting OK on the Data Selection dialog box.
3. Convert wellhead pressures to sandface pressures by clicking the Calc SFP icon on the Pro-
duction Editor toolbar and specifying wellbore parameters.
4. Create a Step Down Analysis. In the main menu, select Step Down Test Analysis from the
Analysis menu.
5. Select one data point per step rate. Click the white arrow to the left of each row in the grid.
This puts the cursor in the point selection mode. Click a desired data point in the History or
Zoom plot. This places an annotation arrow on the History plot and populates the appro-
priate row on the grid with data. You can use the Zoom plot to fine-tune the point selection.
6. On the bottom-right section, click the Step Down Plot tab to see a pressure vs. rate plot of
the selected data points.
7. In the Model Parameters section (top-right), enter initial estimates of model parameters with
which to history match the analysis data. A model is generated and displayed in the Step
Down Plot when all input parameters have been entered.
8. Modify the model parameters, especially those with low confidence levels, to improve the
match between the selected data points and the model's calculated values. Alternatively,
select parameters for automatic estimation by checking the checkbox beside each para-
meter, and then click the Automatic Parameter Estimation ( ) icon.
10. To generate a report, click the Reports tab and Preview Report sub-tab. A step down test
report is generated.
Changing Global Unit Settings Generating & Submitting PAS AER Files
Starting a Project
1. Click the annotation arrow icon located on the toolbar on any of the plots ( ).
b. To add synchronized arrow annotations to the same point on each plot, select the syn-
chronized arrow ( ).
When a synchronized arrow is moved, the corresponding arrows on other curves and plots
will also move.
2. Click the plot near the data point you want the annotation arrow to point to.
3. Click-and-drag the arrow to a different data point, or click the arrow and use the arrow keys
to fine-tune your selection.
1. Click the Edit menu; then click Defaults; then File Defaults… The New File Defaults dialog
box opens.
Note: This feature is available for the black oil and modified black oil models, but it is not available
for coalbed methane (CBM) numerical models.
To export to CMG:
2. Save your .dat file with the default filename, or rename it if you want.
After this is done, you can run the .dat file in CMG IMEX. If you wish to modify the .dat file,
you can load it into the CMG Builder and edit it there. Or, you can use any text editor.
l Rock properties
l Fluid properties
l Well control information (e.g., if the Use Oil Rate calculation method is selected for
the model, oil rates are written to the appropriate section of the .dat file)
Note: The current implementation of the Export to CMG feature does not support some of the
advanced model parameters. If you try to export a model that has such parameters, an error
message is displayed in the Status Window, and the CMG .dat file is not created. To open
the Status Window, click the View menu and select Toolbars and Docking Windows,
Status Window.
2. In the Entity Viewer - Hierarchy pane, select one or more wells or a group for exporting; then
click the Export icon ( ) in the main toolbar and select To WellTest. Note that scenarios
in Harmony Enterprise are not supported in WellTest.
You can export several wells to the same file; however, you are only able to open one well at
a time in WellTest.
3. In the Save As dialog box, browse to your preferred location and enter a filename; then click
Save.
4. Open WellTest.
5. Click Open Existing File and select Open File. Navigate to your previously saved .hmexp
file.
l click the File menu and select Open after opening a WellTest project.
l
click the Harmony Enterprise icon ( ) in the main toolbar after opening a WellTest
project.
6. If your .hmexp file contains more than one well, the Well Selection dialog box opens. Select
the well you want to open in WellTest, and then click Open.
1. Find the Start of Injection line. It is a red vertical line on the Data Chart on the Production
Editor tab. By default, it is placed at time zero. If you can't find it, verify that the test type is
set to Minifrac. .
2. Scale the time axis of the Data Chart, so that you can see the Start of Injection line and injec-
tion portion of the minifrac test data.
3. Click and drag the Start of Injection line to the breakdown point.
For conventional tests, fluid types are only enabled if flow rates of that fluid are entered in the Pro-
duction Editor. For Minifrac or PITA / Closed Chamber Tests, Gas, Oil, and Water are always
enabled.
3. Click OK.
Note: For drawdown / buildup tests, fluid types will be grayed-out unless a rate for the fluid type
has been specified in the Production Editor.
Note: Detailed rock and fluid properties are required to use the multiphase fluid types. These prop-
erties may be entered in the Properties / Multiphase Properties tab.
When you click New File, and select your type of test, the Startup Wizard opens.
The wizard guides you through common procedures, from entering data to completing an analysis.
Each wizard dialog box contains instructions, or requires your input to complete that step. If you
close a dialog box, you can continue with the wizard steps using the Wizard Menu.
l a Harmony Enterprise project. For more information, see interoperability with Har-
mony Enterprise.
When you click Help, you can open this help system, or open the Licensing dialog box to configure
or activate / deactivate licenses. For more information, see licensing.
Creating an Analysis
Identifying Closure
Identifying ISIP
2. Click the Create Deconvolution Case button located on the Deconvolution Manager sub-
tab.
Note: If you don't have an analysis set up, the Create Deconvolution Case button will be
grayed-out.
A new tab is created, using the current data from the Production Editor tab.
To create an analysis:
2. Click the Analysis menu, and select your analysis type. Your analysis will open in a new tab
named "Diagnostic #".
3. Rename this tab by right-clicking it and selecting Rename. Use a name that you can easily
remember.
For your analysis, a Parameters dialog box is displayed in the window.
4. Enter values for all of the fields marked Required. After you have done this, the Lines dialog
box is displayed.
5. Use these dialog boxes to customize your analysis and determine your desired reservoir
parameters.
6. After you are satisfied with your analysis, create a model by clicking the Create Model but-
ton in the Lines dialog box, or by clicking the Models menu.
2. Right-click inside the arrow's parameter display box, and select Radius of Investigation.
Note: The calculated Radius of Investigation updates automatically when the arrow is
moved to another point on the plot.
4. If you want to modify the permeability value after you have already calculated the radius of
investigation, right-click inside the arrow's parameter display box and select Radius of
Investigation Options.
1. On the G-Function plot, rotate the Fracture Closure line through the straight-line portion of
the Semilog Derivative data by clicking the line and dragging it.
2. Click the closure arrow icon on the G-Function plot toolbar and select the synchronized
arrows.
3. Click the G-Function plot where the Semilog Derivative data starts to deviate downward from
the straight line to add the closure arrow annotation.
4. Click the closure arrow and use the left and right arrow keys to fine-tune your selection.
1. Position the ISIP line through the final injection point and the early falloff pressure data. The
line may be rotated around its center point by clicking and dragging the line, or by moving
the line by clicking and dragging the center point.
4. Click the ISIP arrow to select it, and then use the left and right arrow keys to position the
arrow where the pressure data starts to deviate from the ISIP line.
1. Create a diagnostic analysis. If the test was performed on a vertical well, click the Analysis
menu and select Vertical Diagnostic. Or, if it was performed on a horizontal well, click the
Analysis menu and select Horizontal Diagnostic.
2. Enter pay (h), porosity (Φ), and saturations in the Parameters dialog box. Enter effective hori-
zontal well length (Le) for a Horizontal Diagnostic.
4. Identify flow regimes. On the Derivative plot (i.e., bottom-right plot), look at the slope of the
semi-log derivative data. A well-established zero slope (horizontal trend) indicates radial
flow has developed; a slope of 1/2 indicates linear flow; and a slope of 1/4 indicates bilinear
flow. See Flow Regimes for more information.
If radial flow was not identified, you may want to place a radial flow line at the end of the data
to get estimates of permeability and p* to initialize the reservoir model.
6. Position the lines through the appropriate data. Click and drag the center point of the line on
the Derivative plot. On the specialized plots (e.g., Radial plots), you may also click and drag
the center point of the line, or click the line to rotate it.
Note: The radial line yields estimates for flow capacity (kh), skin, and p*.
2. Enter Datum Depth (TVD) in the Parameters dialog box. This is used to calculate the ISIP
(Instantaneous Shut-In Pressure) gradient, which is displayed on the parameter annotations
on the pre-closure analysis plots after ISIP has been identified.
3. Identify ISIP.
4. Identify closure.
2. Enter net pay (h) and porosity (phi) in the Parameters dialog box.
a. If radial flow was identified, click the Impulse Radial button in the Lines dialog box.
b. If linear flow was identified, click the Impulse Linear button in the Lines dialog box.
5. Position the lines through the appropriate after-closure data. Click and drag the center point
of the line on the Derivative plot. On the specialized plots (Minifrac Radial or Minifrac Linear
plots), you may also click and drag the center point of the line, or click the line to rotate it.
Note: If radial flow was not identified, place a radial flow line at the end of the after-closure data to
get upper limits of permeability and p*.
2. Ensure net pay (h) and porosity (phi) have been specified in the Parameters dialog box.
3. Identify flow regimes. On the Derivative plot (i.e., bottom-left plot), look at the slope of the
semi-log derivative data. A well established zero slope (horizontal trend) indicates radial
flow has developed, and a slope of +1/2 indicates linear flow.
b. If linear flow was identified, click the Linear button in the Lines dialog box.
5. Position the lines through the appropriate after-closure data. click and drag the center point
of the line on the Derivative plot. On the specialized plots (Minifrac Radial or Minifrac Linear
plots), you may also click and drag the center point of the line, or click the line to rotate it.
Radial flow gives an estimate of p* and permeability. Linear flow yields Xf*sqrt(k).
Note: If radial flow was not identified, place a radial flow line at the end of the after-closure data to
get upper limits of permeability and p*.
a. If radial flow was identified, select Vertical from the Model menu.
b. If linear flow was identified, select Fracture with Boundaries from the Model menu.
b. Fracture with Boundaries model: initialize Xf and k if estimates were not obtained from the
Soliman After-Closure Analysis.
1. If the PITA test was performed on a vertical well, click the Analysis menu and select Ver-
tical Diagnostic. Or, if it was performed on a horizontal well, click the Analysis menu and
select Horizontal Diagnostic.
2. Enter pay (h) and porosity (Φ) in the Parameters dialog box. Enter effective horizontal well
length (Le) for a Horizontal Diagnostic.
Note: The pay thickness (h) for a PITA test in a vertical well is usually entered as the height
of the perforated interval.
3. Decide whether there is a changing liquid level in the wellbore. If there is a changing liquid
level, it dominates the changes in pressure in the wellbore, and you should click the
Changing Liquid Level checkbox within the Parameters dialog box. If a single-phase
exists in the wellbore, changes to pressure in the wellbore are controlled by compressibility.
Note: The Changing Liquid Level option is not available for PITA inflow tests in gas reser-
voirs.
4. If the Changing Liquid Level checkbox is not selected (e.g., for a PITA inflow test in a gas
reservoir):
a. Click the button beside the Vw field.
b. Select the casing size from the list of common casing IDs (inner diameters). Enter the meas-
ured depth of the well and click Select. A calculated value is displayed in the Vw field in the
Parameters dialog box.
b. Select the casing size from the list of common casing IDs and click Select. A calculated
value is displayed in the Vu field in the Parameters dialog box.
6. Identify flow regimes. On the Derivative plot (i.e., bottom-right plot), look at the slope of the
semi-log derivative data. A well established zero slope (horizontal trend) indicates radial
flow has developed, a slope of 1/2 indicates linear flow, and a slope of 1/4 indicates bilinear
flow.
a. If radial flow was identified, click the Late Time-Impulse Radial button in the Lines dialog
box.
b. After an impulse radial line has been added, the Early Time +ve and -ve Skin buttons are
enabled. If you do not know if the skin is positive or negative, you can add both and see
which skin value makes the most sense.
8. Position the lines through the appropriate data. Click and drag the center point of the line on
the Derivative plot. On the specialized plots (e.g., Late Time), you may also click and drag
the center point of the line, or click the line to rotate it.
Note: The impulse radial line yields estimates of permeability and p*. The early time lines
yield estimates of skin.
9. Select the Early Time +ve Skin line from the Analysis Lines box.
This line is displayed on the Early Time plot, and you can use it to estimate skin.
1. Click the Model menu and select a model. In most PITA tests performed on vertical wells, a
vertical model is used for PITA tests.
A new tab opens (e.g., Vertical 1).
2. History match the PITA data. Focus on the late-time data to try to get estimates of per-
meability and pi.
You may decide to do a history match to get more accurate results. For additional inform-
ation, see History Matching Overview. For a PITA test, we recommend doing the history
match on the Late Time plot.
1. Memory management: After importing files with millions of data points, data reduction is
required to prevent slow, sluggish behaviour and "out of memory" related errors. Although,
Data Management is partly a placeholder of the raw data, there is no need to keep millions
of data points, and it is more important to free up computer memory for other modules within
WellTest. We recommend filtering gauges with millions of points down to 100,000 points
immediately after importing. Data reduction for memory management is primarily accom-
plished by Filtering in Data Management. Although this filter has advanced options, in most
cases, a simple arithmetic filter can be used for the entire dataset, and will keep enough res-
olution for analysis. For additional information, see Handling Large Datasets.
2. Calculation efficiency: When preparing data for analysis, another level of filtering is per-
formed to reduce calculation times. Filtering in the Production Editor drastically reduces con-
vergence times when performing automatic parameter estimation and running numerical
simulations. Although 20 points could capture the trend of any flow period, it is impossible to
know which 20 points to pick before scrutinizing the data. Therefore, we find a balance
between the amount of data we keep to analyze and computational time. For analytical ana-
lysis / modeling, it is highly recommended to filter the Production Editor data down to less
than 5,000 data points. For numerical modeling, 1,000 points is more than sufficient.
To filter data:
1. Click the Data Management tab; then click the Filter sub-tab.
If you have two or more gauges, the Select Gauge / Dataset dialog box opens. Select a
gauge; then click OK.
The default filter settings use an Arithmetic filter for flow periods, and a Logarithmic filter for
shut-in periods. The default filter settings may be adjusted, or a unique filter may be applied
to each section of data.
2. To apply different levels of filtering to different sections of the data, create partitions by click-
ing in the partition area at the top of the plot. A partition marker is displayed at each point
where you click. Drag a partition marker to a new location if needed, or drag it off the slider to
remove it.
There is a table underneath the plot that displays different options for filtering your data.
1. Identify the flow period (row) of interest and deselect the Use Default checkbox.
3. If the filter type chosen is Arithmetic, you can specify the Time Unit and Arithmetic Mode.
Note that the Time Unit changes the units of the Arithmetic Time Step column.
1. Click the Filter icon ( ) on the Production Editor toolbar and select Primary Data Fil-
ter. The Filter Data dialog box opens.
2. Use the checkboxes to select the number of data points to be used for each flow / shut-in
period. We recommend that you use Logarithmic checkboxes for the Buildup data points,
and Arithmetic checkboxes for the Drawdown data points.
3. You can edit the numbers in the Logarithmic Points Per Cycle and Arithmetic Time Step
columns to control the amount of filtering.
1. Populate the mandatory fields in each of the sub-tabs within the AER PAS tab.
Note: While the certification feature provides a sanity check, there are some known issues
where our Certify feature misses rules checked by the AER (false positives), or
detects problems the AER does not (false negatives). The most reliable way to
ensure a PAS file meets AER requirements is to validate through the AER DDS web-
site. Continue to the next step after resolving the obvious errors detected by this fea-
ture.
l Specify a PAS file name. Click the drop-down arrow in the PAS file name field to
browse to the directory where you'd like to save your PAS file.
l Specify a sub-zip file name. Click the drop-down arrow in the Sub zip file name field
to lto browse to the directory where you'd like to save your zip file.
l Specify the location of the PDF copy of the report. Click the drop-down arrow in the
Image (pdf) file name field to navigate to the PDF report that you’ve already gen-
erated.
Note: The report is not generated from the PAS tab. If you have not yet created a PDF copy
of the report, you must generate a report and then create a PDF of this report1.
l Click the Save PAS and ZIP file(s) button. This will generate the PAS file and a .zip
file that contains the PAS file and the PDF report.
l Click the AER PAS Help link, and sign in to their system. Navigate to AER / Sub-
missions / Well Test Data in the left-side table of contents. Validate the .zip file to see
if there are errors, and submit the .zip file if it's error-free.
1 use the Print functionality, except select Adobe PDF as your printer
2. Click the Enable Alberta AER PAS Options box to remove its checkmark. This will hide
the PAS tab.
3. Generate and submit a PAS TRG file paying special attention to the fields in the AER PAS /
AOF - IPR Results Summary tab.
*If a single-point deliverability test, such as a flow and buildup, is performed and the AOF of the well
is greater than 300,103 m3/d, an application for an exception to the rule can be made to the AER. If
the application is approved, and the exception granted, the AER provides a Single Point Author-
ization Number. This number is then entered in the Single-point Authorization field in the AOF / IPR
Results Summary tab, and the PAS file is accepted, and the initial deliverability requirement for that
well is met.
Note: If the well is producing inline, a minimum of 14 days (336 hrs) of initial production with meas-
ured wellhead pressures must be entered in the Production Editor tab. There must be at
least 28 lines of data spanning the (minimum of) 14 days of production. The resultant PAS
TRG file also fulfills the initial production requirement.
1. Conduct a well test, such as a buildup or falloff test, and analyze the data. Analysis line extra-
polations, or preferably, modeling can provide estimates of initial pressure. See Creating an
Analysis for additional information.
3. Generate and submit a PAS TRG file paying special attention to the following:
l Verify that the Test Purpose field on the Test Data tab is set to Initial Test.
l Verify that the Test Interpretation Present field on the Test Data tab is set to Yes.*
* – this field can be set to "No" in the rare case where the shut-in pressure changes
less than 2 kPa/hr over a 6 hour period.
l In the Analysis Input Parameters tab, select the most appropriate analysis or model
from the drop down box. This action populates fields within the Analysis Input Para-
meters and Pressure Results Summary tabs.
l Ensure all gauge information is specified within the Gauges / Gauge Header tabs.
Note: Acoustic well sounder (AWS) data cannot be used to determine the initial pool pressure (for
the first well in a pool).
FLD
PAS
Note: For Excel files, WellTest only imports from the tab that was viewed when the file was saved.
If the last tab viewed and saved does not contain a data table, an error message is dis-
played.
If you have date, clock time, and cumulative time columns within the same file, either
import date and clock time, OR cumulative time. If "date" is selected for one column,
cumulative time is not available for selection in another column and vice versa.
If date and clock time are in the same column, specify "date" as the data type, and the
units of the date for the units. After selecting your columns and clicking Next, a dialog
box opens where you can specify the clock time units (e.g., hh:mm:ss).
If the date or clock time columns are in an unsupported unit format, use cumulative
time if it's available. Otherwise, the format may need to be changed in MS Excel.
b. In the second row of the header, select the units of the data in that column.
d. Click Next.
The data is displayed in a new tab within the Data Management tab, Data sub-tab.
Data may be imported using the wizard, or by opening the Import data dialog box, which may be
accessed by clicking File, Import Data. The wizard takes you step-by-step through the import pro-
cess and this topic describes each step.
3. Click the Skip Blank Data Line checkbox (if it is not checked).
4. Select the data that you'd like to bring in to WellTest, and click OK.
Your data will open in a new tab in the Data Management, Data sub-tab area.
Data may be imported using the wizard, or by opening the Import data dialog box, which may be
accessed by clicking File, Import Data. The wizard takes you step-by-step through the import pro-
cess and this topic describes each step.
PRD PAS
To import a PRD PAS file:
4. Select the data you'd like to bring in to WellTest, and click OK.
Your data opens in a new tab in the Data Management, Data sub-tab area.
PAS AER
To import a PAS AER file:
History Matching
Advanced Models
Legacy Models
The choice of model can drastically change the outcome of a forecast. For example, pressure tran-
sient analysis of radial flow data can be used to find the reservoir flow capacity (kh). If heterogeneity
is observed in the data, one might use the composite model with different values for permeability in
each zone. Another analyst might decide that using a two-layered multilayer model, having different
net pays, is more probable. The forecast of the composite and multilayered models for this situation
can differ drastically, even though both models can match the test data. The analyst must decide
which model is better.
For assistance in selecting a model, refer to basic model typecurves to help identify which model
best represents your pressure response.
1. Adjust the model input parameters until a close match is obtained between the simulated
and measured pressures.
2. Perform automatic parameter estimation (APE) to fine-tune the results. See Performing APE
in Advanced Models and Performing APE in Legacy Models.
3. If an acceptable match is not obtained, you should manually adjust the parameters and
repeat the APE process, or consider a different reservoir model.
In the advanced suite of analytical models within WellTest, an anchor point is added by clicking the
Anchor icon on the History Plot toolbar.
Clicking the Anchor icon automatically places an anchor arrow on the measured pressure preceding
the first selected flow period. At the same time, the mouse pointer changes to an anchor, so that you
can relocate the anchor from the default position by clicking any measured pressure on the history
plot. The anchor point is annotated with an arrow, and if you change the position of the anchor arrow,
the model will automatically recalculate.
Note: The anchor option is currently only available for analytical models in Calculate Pressure
mode.
1. Click the Advanced Models tab. The Model Manager sub-tab will be displayed by default.
3. Customize the model by specifying the model name and fluid type(s).
The forecasting and test design module within advanced analytical models can also be used without
production history, to investigate the potential deliverability of a new reservoir given the uncertainties
of reservoir flow capacity and volumetric parameters. It can also be used to investigate which flow
regimes would be investigated during a test for a specified set of reservoir properties, flow, and shut-
in periods.
To create a forecast:
1. After creating an Advanced Model, click the Forecast and Test Design sub-tab.
3. In the Duration column of the forecast table, enter a value for the length of the forecast.
Changing the number of steps may change the calculated expected ultimate recovery as
cumulative production is dependent on the length of each timestep. Flowing pressure for the
forecast is automatically populated to the last historical value; however, this value can also
be edited.
4. (Optional) After a value is entered in the forecast table, an additional row is displayed where
a second forecast period can be entered. This row does not have to be populated for a fore-
cast to run. Rows can also be added or deleted by right-clicking a row and selecting Insert /
Delete row.
Forecasted data is displayed on the Results plot. This data consists of a fluid rate, fore-
casted pressure, and the average reservoir pressure. If more than one forecast case is
present, the results of each are displayed on the Results plot.
5. View the results of your forecast, displayed as a plot, in the Forecast Results section.
Forecasting Injection
Running forecasts or test designs is similar for injection and production wells. However, for injection
wells:
l Calculate the recombined gas rate factor (RGRF), which is dependent on the instant-
aneous value of the condensate gas ratio (CGR)
Note: RGRF is calculated from separator conditions, condensate gravity, gas gravity, and CGR.
There are two options that control how CGR is calculated over the forecast period:
l To enter a CGR that varies linearly over the forecasted period, select Ramp for the
Ratio Interpolation method; then enter a Final CGR value.
l in the Properties tab / Gas Properties Input pane, select Gas Liquid Content from
the Gas Type drop-down list, and Ovalle et al from the Vaporized Oil Ration Cor-
relation drop-down list.
The resulting forecast graph and table show recombined gas, gas (dry), and condensate
rates.
In splitting the forecast of analytical models, the assumption is that these effects / properties in the
reservoir can be ignored:
Note: The effect of liquid blockage on well productivity can be modeled by increasing skin with
time.
For gas analytical models, pseudo-pressure and pseudo-time are used to account for the change in
gas properties.
For oil analytical models, fluid properties are assumed to be constant: properties are estimated at
the initial pressure using correlations or tables set in the Properties tab, and these properties are
used for calculations throughout the production history and the forecast.
This assumption is acceptable if the pressure in the reservoir remains above the bubble point, since
the variations for under-saturated oil properties with pressure are negligible. However, if the pres-
sure drops below the bubble point, oil properties (in particular – oil compressibility) change sig-
nificantly.
Example: Consider an under-saturated oil reservoir with: pi= 5000 psi and pbp = 4500 psi.
The well was producing at a constant rate (drawdown) and then was shut in (buildup). During the
drawdown, reservoir pressure dropped to 4400 psi; therefore during the buildup, the reservoir was
filled with saturated oil. As a result, saturated oil properties should be used to history match the
buildup.
2. Expand the PVT Pressure group and set ppropto 4400 psi.
Note that with this approach, the drawdown portion of the data will not be matched. However, match-
ing just the buildup portion of the data is useful, as it helps to estimate current reservoir and com-
pletion properties.
3. Run APE.
Note: Automatic parameter estimation (APE) is not available for numerical models.
Deselecting Points
If your data has some outliers, or some periods of unreliable data, you can deselect these points. To
specify data points to be ignored by APE:
1. On the toolbar of any plot within the Plots sub-tab of an advanced analytical model, click the
Display a toolbar with options to select data points icon, and select the Select Points
option from the sub-toolbar.
2. Right-click and drag your mouse over the data on plot. (Clicking-and-dragging reselects any
deselected data.)
Weighting Points
You may want to weight some data points more heavily than others when performing APE. For
example, when analyzing buildup data to get an estimate of average reservoir pressure, it is more
important to match the data at the end of the buildup, than the start of the buildup. Similarly, if a well-
bore effect is observed, you may want to ignore those points when history matching.
To weight points:
1. Click the Weight Data Points icon on the history plot toolbar.
2. Move the mouse over the pressure data on the history plot. The cursor changes to a spray
can.
3. Click-and-drag the cursor to spray the points you would like to add weighting to. Right-click
and drag to remove weighting from weighted points.
1. In the Analytical Model pane, select the parameters you'd like to vary by clicking the Auto-
matically calculate variable checkbox to the right of the parameter's field.
WellTest adjusts the parameters in order to improve the history match between the sim-
ulated and measured pressures on the data plot. The process stops automatically, or if you
want to stop it manually, click the APE icon again.
After you have achieved your desired result in CalcP mode (calculates sandface flow-
ing pressure based on measured rate), we recommend that you switch to CalcR mode
(calculates the rate based on given sandface flowing pressures), and verify if the calculated
rate matches the actual rate. If these rates are not close, additional history matching should
be done.
Note: If an acceptable match is not obtained, you can manually adjust the parameters and repeat
the APE process, or consider a different analytical model.
1. Some of the model parameters (i.e., initial pressure, fluid saturations, net pay etc.) are
already specified in the Properties tab. To copy parameters from the Properties tab to the
model, click the Defaults icon ( ) on the toolbar, and select Reservoir from the drop-
down list.
2. If you have an analytical model, copy parameters from the analytical to numerical model by
clicking the Defaults icon ( ) and select the analytical model from the drop-down list.
3. In the Numerical Model pane, select your calculation method . You can use the rates of any
of the fluids presented in the model to calculate rates for the remaining fluids, and the sand-
face flowing pressure. Alternatively, you can use the sandface flowing pressure to calculate
rates for all the fluids presented in the model.
Note: Before the model starts a simulation, it performs multiple consistency checks on the input
properties.
While the model is calculating, a rotating gear is displayed on the header of the tab. To stop
the simulation, click the Stop Simulation ( ) icon.
2. You can replace one or more plots on the dashboard by clicking the Add an available view
icons.
3. After the calculation is finished, you can see how pressures or saturations change through
time by clicking the following icons on the toolbar of any Color Shading plot.
6. To history match the model, change the parameters manually. Most of the history matching
parameters are located on the History Model pane. However, you may also want to change
some parameters located in the Properties tab (e.g., relative permeability curves, Rs curve,
Rv curve, etc.).
After the model has been populated and calibrated to the production and pressure data, pro-
duction forecasts can be created and compared under a variety of different constraints. The
procedure for forecasting using a numerical model is similar to Forecasting Using Analytical
Models.
Note: Wellhead forecasting is not currently available for the numerical models.
In WellTest 2013 v2 (v 7.7.0) and earlier: it is possible to create a (black oil) numerical model for a
liquid-rich gas well. However, liquid drop-out is not modeled directly in such models. Instead, gas
and condensate rates measured at surface conditions are combined into a richer single-phase gas
using recombination; then rates and properties of the recombined gas are used in numerical mod-
eling.
This approach is appropriate if condensation takes place mostly at the separator and in the wellbore.
In cases where there is a significant liquid drop-out in the reservoir, the model may be inaccurate,
because the liquid drop-out in the reservoir and its influence on the flow hindrance are not taken into
account. Another disadvantage of numerical modeling for liquid-rich gas wells in previous versions
of WellTest, is that only the recombined gas rate is reported, rather than reporting a separate dry-
gas rate and condensate rate.
In WellTest 2014 v1 (v 7.8.0) and later: liquid-rich gas wells are modeled using Gas Condensate
numerical models. In Gas Condensate numerical models, liquid drop-out in the reservoir is modeled
directly. Both phases (i.e., gas and liquid) are taken into account in the model; therefore flow
hindrance due to interaction between phases is also honoured.
If you have a (black oil) numerical model created in WellTest 2013 v2 (v 7.7.0) or earlier for a liquid-
rich gas well, this model will not run in WellTest 2014 v1 (v 7.8.0) or later. This happens because
starting from WellTest 2014 v1 (v 7.8.0), liquid-rich gas wells are modeled using Gas Condensate
numerical models. You have two options for converting such older models:
Option 1: If you want to consider liquid drop-out in the reservoir, you have to create a new
Gas Condensate numerical model. After the model is created, copy all parameters from the
existing black oil numerical model, using the Defaults Manager. The gas condensate model
then runs in WellTest 2014 v1 (v 7.8.0).
The advantage of this option is that it takes into account liquid drop-out in the reservoir and
honours flow hindrance due to interaction between phases. However, because the nature of
calculation has changed compared to the earlier version, you have to re-match the model.
1. On the Production Editor tab, copy the Recombined Gas Rate column to the
Gas Rate column.
2. In Properties tab, navigate to the Gas Properties group and copy the value for
Recombined Gas Gravity (Gr).
4. In the Properties tab, paste the copied value for Gr into the Gas Gravity (G)
field.
The advantage of Option 2 is that the calculation essentially follows the one performed in
earlier versions of WellTest: instead of modeling liquid drop-out directly, we use recom-
bination, and then use recombined gas rates and properties in numerical modeling. Thus,
the calculation results do not change, so you do not have to re-match the model. However, it
is important to remember that this model is only applicable if liquid drop-out in the reservoir
is negligible.
l Forecasting
Note: Some of the model parameters (i.e., initial pressure, fluid saturations, net pay, etc.) are
already specified in the Properties tab.
1. Copy parameters from the Properties tab to the model by clicking the Defaults icon ( )
on the toolbar; then select Properties Tab from the drop-down list.
2. Copy all available parameters from other models by clicking the Defaults icon ( ) again.
If you have an analytical gas horizontal multifrac model, we recommend that you use it to
populate the hybrid model's parameters.
1. To calculate the sandface flowing pressure based on measured rate, click the CalcP icon (
).
2. To calculate the rate based on a given sandface flowing pressure, click the CalcR icon ( ).
3. To have the model calculate automatically, and recalculate automatically every time you
If the AutoCalc icon is deselected, you will need to click the Synthesize icon ( ) every
time you want to recalculate the model.
Note: If the CaIcP icon ( ) is selected, WellTest uses rates to calculate pressures.
The purpose of history matching is to get the calculated pressures to reproduce actual pressures as
closely as possible by changing the model's parameters. You can vary model parameters manually,
or you can use automatic parameter estimation.
1. Click the Display a toolbar with options to select data points icon on the Plots tab, and
select the appropriate sub-option to remove unwanted data.
2. If it is important to match some specific portion of the data (e.g., the last part of the pro-
duction), you can weight data points by clicking the Weight Data Points icon ( ).
1. To select which parameters to vary, click the Automatically calculate variable checkbox to
the right of the parameter's field.
3. Enter the appropriate values, and then click the Accept button.
While APE is running, selected parameters are automatically varied to minimize the error between
the simulated pressure and the actual pressure. The updated simulated pressure is shown on the
plots while the calculation is running.
While APE is running, you can wait until it finishes, or if the history match looks good, you can stop
After you have achieved your desired result in CalcP mode, we recommend that you switch to CalcR
mode and verify if the calculated rate matches the actual rate. If these rates are not close, additional
history matching should be done.
To speed up your calculation times, you can enable the sparse data option, which is available for the
hybrid model. By doing this, you will run simulations with a smaller number of timesteps.
2. Click the Use Sparse Data icon ( ) on the hybrid model toolbar.
The Use Sparse Production Data dialog box opens. When using sparse data, the hybrid
model uses coarse timesteps for calculations, and each coarse timestep contains n points.
l If you select CalcP mode, the sparse gas rate used for calculations is displayed as a
thin red line
l If you select CalcR mode, the sparse flowing pressure used for calculations is dis-
played as a thin brown line
4. To ensure that the sparse dataset represents the actual production history, adjust the value
for n, as required.
Note: Your weighting and point selection are honored when running APE using sparse data.
While calculation results are slightly less accurate, the calculations run faster. (After APE
finds a match using sparse data, you can run one calculation using the entire production
dataset (n = 1) to ensure that the match is still good.)
For CalcR mode: Drawdowns are divided in coarse timesteps. Coarse timesteps contain n points,
and the pressure for each timestep is calculated as "last measured pressure of the timestep". For
buildups, all of the production data is used.
1. Enter the pressure data obtained at a shut-in well (observation well) in the Production Editor
tab.
2. Enter the production / injection rates measured at the producer / injector (active well) in the
Production Editor tab.
Note: Any combination / number of production rates and injection rates can be accom-
modated, but injection rates should be entered as production rates with a negative
value.
4. Enter model parameters based on your available information. For example, initial estimates
of permeability are usually obtained from pressure buildup (PBU) tests and/or from pressure
falloff (PFO) tests performed with active well data.
5. Select the Observation Well checkbox, and set the location of the observation well relative
to the active well.
The legacy suite of models within WellTest, which are accessible via the Models menu, provide
three possible anchor options: initial pressure, final flowing pressure, or any other specified pres-
sure.
1. Initial Pressure Anchor: This selection is used when analyzing drawdown data with a
known (measured) initial pressure. The synthetic initial pressure is set to the reservoir pres-
sure specified in pressure, volume, temperature (PVT) properties.
2. Final Flowing Pressure Anchor: This anchor point is convenient when analyzing buildup
or falloff data. The initial pressure is usually not measured, and with anchoring to the final
flowing pressure, you can focus on history matching the good quality, shut-in data, and then
checking to ensure the synthetic initial pressure is reasonable.
3. Specified Pressure Anchor: In situations where neither the initial pressure nor the final
flowing pressure is known, it is possible to specify any other measured pressure.
An anchor point is selected by clicking one of three icons on the model’s toolbar.
Note: Legacy models may also be created by clicking the Wizard menu and selecting Model
selection. This opens the WellTest Wizard: Model Selection dialog box. Click the radio but-
ton beside the type of model you wish to create and click Next.
To create a forecast:
1. Select the model input parameters you want to vary by clicking the Auto checkbox beside
each parameter in the Parameters dialog box.
2. Click the start Auto button on the model toolbar to start the auto match. WellTest adjusts the
parameters to improve the history match between the simulated and measured pressures.
Weighting Points
You may want to weight some data points more heavily than others when performing APE. For
example, when analyzing buildup data to get an estimate of average reservoir pressure, it is more
important to match the data at the end of the buildup, than the start of the buildup. Similarly, if a well-
bore effect is observed, you may want to ignore those points when history matching.
APE tries to minimize the error between the synthetic (model) pressures and the measured pres-
sures. The error terms for weighted points are amplified in the convergence algorithm. In other
words, APE tries really hard to converge on weighted points, less hard converging on points without
weighting, and ignores deselected points.
The model toolbar contains a number of buttons above the Auto button that are used for weighting
points.
l
The yellow icon ( ) is for weighting points
l
The white icon ( ) is for excluding points
To weight points:
2. Move the mouse over any plot with measured data, such as the History plot, the Derivative
plot, or the Typecurve plot.
The cursor changes to a spray can.
3. Spray the data points you would like to add weighting to by clicking over the measured pres-
sures. Right-clicking changes the weighting of data points back to a normal weight.
l
This icon ( ) selects all points.
l
This icon ( ) excludes all points outside the selected flow period, and selects all
points within the selected flow period.
l
This icon ( ) icon selects 10 equally spaced points within the selected flow
period, and excludes all other points.
l
This icon ( ) weights or unweights the point immediately preceding the selected
flow point (e.g., the final flow pressure).
Entering Custom Data for Dry Gas / Wet Gas Renaming Gauges
2. Ensure that you have values in either the Tubing Pressure or Casing Pressure columns of
the table.
3. If appropriate, specify liquid levels in the Liquid Level column of the table. Static oil and water
gradients may be entered in the Oil and Water Gradient columns. Alternatively, you can opt
to calculate these gradients by clicking the Calc checkbox in the header of the Oil and Water
Gradient columns.
If you do not see these columns on the table, see Grid Options.
l casing specifications
6. If you are dealing with more than one phase, click the Pressure Loss Correlation tab and
select the most appropriate correlation.
7. If the well is deviated, click the Deviation Survey tab. This table can be populated manually,
or by copying and pasting from the deviation survey.
1. Find the first or last data point you want to keep by clicking that point in the Data Chart.
You may fine-tune your selection by clicking the annotation arrow on the Zoomed Data
Chart, and then use the right and left arrow keys to select neighboring points. This highlights
the corresponding row on the data table.
2. Right-click the number of the highlighted row in the data table and select Delete All Rows
Before or Delete All Rows After.
2. Spray the data you would like to remove by holding down the left-mouse button and drag-
ging your mouse over the points.
3. Click the Delete selected points toolbar icon in the Production Editor Data Chart.
2. On the top left plot, position the line through the general trend of the drawdown (or injection)
pressure data.
3. On the top right plot, position the line through the first few shut-in points.
4. Fine-tune the line positions on the Data Chart plot on the bottom left.
6. After you are satisfied with the shut-in point, click the Click here to select shut-in point
box in the Shut-in Point plot, or click the Select Shut-in Point button in the top left of the
module window. Then click Yes in the popup that opens.
The shut-in point will be inserted into the Production Editor tab.
7. Check the shut-in point on the Production Editor tab. Sometimes zero rates that were in the
Production Editor tab prior to inserting the shut-in point need to be deleted manually.
Note: When manually entering rates, make sure the rates displayed in the Data Chart match your
expectations.
3. Select correlations from the PVT Correlation and Viscosity Correlation drop-down lists.
5. Enter data for pcal and Tcal in the Calibration Properties section.
6. Enter Rs cal, Bo cal, co cal, and μo cal in the Calibration Properties section.
7. Check the plots to ensure that the calibration is reasonable. If the calculated oil com-
pressibility goes negative, check the gas and oil properties, as well as the calibration to
ensure that everything has been entered correctly.
4. Click the PVT Correlation drop-down list and select Custom Table.
5. Enter data in the p column; then enter data into either the Z or Bg columns of the Custom
Table. The other data will be calculated.
4. Select the properties you want to display (i.e., numerical or non-numerical) and click OK.
5. Click the PVT Correlation drop-down list and select Custom Table.
6. Enter data into the p column; then enter data into either the Z or Bg columns of the Custom
Table. The other data will be calculated.
9. Click the Gas Type drop-down list and select Liquid-Rich Gas.
10. Click the Vaporized Oil Ratio Correlation drop-down list and select Custom Table: Liq.
Drop-Out.
3. Enter your Oil / Condensate properties. If applicable, enter your Calibration properties.
4. Click the PVT Correlation drop-down list and select Custom Table.
5. Enter data into the p, Rso, and Bo columns of the Saturated Oil Custom Table.
2. Insert rows into the Production Editor tab and enter the time and rate for each new row.
l Double-clicking the Date / Clock Time column launches a Calendar dialog box. Altern-
atively, you can enter the date and time manually provided the format is consistent
with the units of the column.
l In many cases, the Time column is the same as the Cumulative Time column.
However, delta times may be entered as negative times. To enter daily production
data, enter -24 hrs in the Time column, and the Cumulative Time and Date / Clock
Time columns automatically update.
l When calculating sandface pressures, all fluid rates must be entered to establish the
most representative sandface flowing pressures.
l When using measured sandface pressures, or when there are no pressures available
for the previous production period, the following rates are used:
l Dry gas reservoirs
l Gas-condensate reservoirs:
-gas rate only if condensate gas ratio (CGR) is constant
-gas and condensate rates if CGR is changing
l Oil reservoirs
-oil and water rates if above bubble point
-gas, oil, and water rates if below bubble point
Click OK. Blank tables are displayed in place of the data plots.
2. Manually enter values in the table. Or, click the Import Data icon ( ) in the Inputs pane
toolbar to populate the table.
The maximum amount of data that can be imported into WellTest depends on the amount of random
access memory (RAM) the computer has, as well as how many other programs are using up the
RAM.
l When extracting data from a database, use appropriate queries to extract enough
data to obtain a good interpretation, but not so much that it affects performance.
Extracting hundreds of thousands of points from millions of points yields a large
enough dataset to perform further interactive filtering in WellTest, without losing sig-
natures in the data.
l When dealing with flat files, such as ASCII files, most wireline companies can provide
lower frequency data than the gauge sampling frequency.
l Use IHS ValiData™, a separate program, to reduce data before bringing it into
WellTest. ValiData can handle approximately five times more data than WellTest.
l Filter gauges. The filter in WellTest's Data Management tab is meant to be used for
data reduction and should be used when dealing with super large datasets. See Fil-
tering in Data Management for additional information.
l In extreme cases, where re-drawing issues or sluggish behaviour prevents the inter-
pretation of data, deleting data management gauges reduces the memory footprint of
the file. This can be done by right-clicking the gauge tab above the grid and selecting
"Delete Gauge".
1. Click the Data Selection icon ( ). The Select data to analyze dialog box opens.
2. Use the drop-down menus to select the datasets that you wish to use for analysis and mod-
eling.
1. Right-click a heading in the data table and select Manipulate. The Manipulate dialog box
opens.
Specify how you would like to change the data under the Manipulate Options section, and to
which row(s) these changes will apply under the Rows section.
On the left-side pane, one gauge is set as the 'Primary' gauge. If merging two gauges with like data-
sets, the value from the primary gauge will be kept when dates from the two gauges are identical.
Note: Prior to using this function, you must import data into at least two gauges.
To merge gauges:
1. Open the Data Management tab; then click the Merge tab. The Select Gauges to Merge dia-
log box opens.
2. Select the primary gauge and secondary gauge(s); then click OK.
Note: Multiple gauges may be selected in the Secondary Gauge section of the dialog box
by holding down the Ctrl or Shift key while clicking the gauges to be merged.
3. Click Apply.
To smooth data:
1. Click the Data Management tab; then click the Smooth sub-tab.
2. If you have two or more gauges, the Select Gauge / Dataset dialog box opens. Select a
gauge; then click OK.
Gauges with rate and pressure data are automatically partitioned into flow and shut-in peri-
ods. You can merge pressures and rates prior to smoothing. Alternatively, you can manually
partition your data by clicking in the partition area at the top of the Smoothing Plot. Drag a
partition marker to move it, or drag it off the plot to remove it.
4. Click-and-drag the middle slider to the right to smooth more, or to the left to smooth less. The
slider bars on the top and bottom control the degree to which smoothing is tapered off on the
left and right ends of the selected region.
5. Repeat steps 3 and 4 for all datasets and regions you would like to smooth.
6. Click Apply.
To synchronize data:
1. Click the Data Management tab; then click the Synchronize sub-tab.
2. Select the primary and secondary datasets from the drop-down menus.
4. Click Apply.
If you want to add multiple sync points, click the Options button in the top toolbar. The data
between two sync points can be compressed / expanded, or the first / last point from the sec-
ondary data point can be shifted.
collating a report
designing a report
previewing a report
To collate a report:
1. Click the Reports tab. By default, the Preview Report sub-tab is selected.
A pane opens to the right of the Reports pane with an folder named Collated Reports.
3. From the Reports tree, click the name of a Report Template and drag it to the Collated
Reports pane. To select multiple templates, drag a selection box around the template
names.
4. To print your collated report, click the Print icon ( ) on the toolbar.
Additional Options
Additional options include the following:
1. To create a new report category, right-click the empty space within the Collated Reports
pane and select Create New Report.
To create a report:
1. Click the Reports tab; then click the Design Report sub-tab.
3. In the Template Name field, type in your template's name, or use the default name of Tem-
plate1.
5. In the Design Tools section, add different components to the report. Additional information is
provided in the next section of this topic.
6. To save the template for later use, click the Save Report icon ( ).
2. Add or remove plots, grids, annotations, lines, frames, images, text, or schematics.
3. To save the template for later use, click the Save Report icon ( ).
Note: You can modify reports or templates. If you modify and save a report, a new file-specific
report is generated. If you modify and save a template, a new custom template is generated.
The custom template is used to generate custom reports for other files, which contain the
same model and fluid type.
1. From the Existing Reports and Templates section, select Default Templates.
2. From the drop-down list, select the type of report template you would like to modify.
1. Click the Add plot to current report page icon ( ) in the Design Tools section of the
Design Report sub-tab.
2. Double-click the plot rectangle to open the Select Plot dialog box.
Or, right-click the rectangle and select Select Plot….
Note: Clicking the Show Plot Headers checkbox adds header information to the plot. For
additional information, see WellTest Wizard: Headers.
After a plot has been added, right-clicking the plot displays additional options:
l Copy to Clipboard — you can paste the plot into another application.
l Scale Text — changes the size of the text on the plot relative to the curves that
have been drawn.
l Send to Back / Front — with multiple plots you can layer them with some over-
laying portions of each other with the most important one being on top.
l Specify Size... — specify the height and width of your plot in inches.
l Fit to Page — scales the plot to automatically fit your page size.
l Repeat on Every Page — places the plot of every page of your report.
Adding a Grid
A grid often contains parameter information or key results, which are displayed in a table format. To
add a grid:
1. Click the Add data grid to current report page icon ( ) in the Design Tools section of
the Design Report sub-tab.
2. Double-click the plot rectangle to open the Select Grid dialog box.
Or, right-click the rectangle and select Change Grid….
l
Single Page — you can view the report one page at a time. This option is
meant for cases where there are multiple page reports, such as a collated report.
This option is grayed-out if the report template being viewed is a single-page
report.
l
Multiple Pages — you can view the thumbnail view of a multi-page report.
This option is grayed-out if the report template being viewed is a single-page
report.
After a grid has been added, right-clicking the grid displays additional options:
l Copy to Clipboard — you can paste the plot into another application.
l Send to Back / Front — with multiple grids you can layer them with some over-
laying portions of each other with the most important one being on top.
l Specify Size... — specify the height and width of your grid in inches.
l Fit to Page — scales the grid to automatically fit your page size.
l Specific Fit — you can select columns, rows, or both to be fitted in the grid's rect-
angle.
l Dynamic Row Sizing — you can have your grid / table dynamically grow if the num-
ber of rows changes. When this option is selected, columns are always fitted in the
rectangle provided.
l Make Grid Same Size On All Pages — applies to a grid that spreads across sev-
eral pages and uses the same rectangle as the first- page rectangle to do the sizing
of the grid.
Adding an Annotation
To add an annotation:
1. Click the Add annotation to current report page icon ( ) in the Design Tools section of
the Design Report sub-tab.
Note: To change the default font settings, click the Options button.
l Save As a block — opens the Provide Block Information dialog box where you
can enter your block name, and associate the block with your application data.
l Scale Annotation — sets the size of the annotation, which can be enlarged or
reduced.
l Copy To Clipboard — you can paste the annotation into another application.
l Send to Back / Front — with multiple annotations you can layer them with some
overlaying portions of each other with the most important one being on top.
l Specify Size — specify the height and width of your annotation in inches.
l Fit to Page — scales the annotation to automatically fit your page size.
l Repeat on every page — places the annotation on every page of your report.
1. Click the Add line to current report page icon ( ) in the Design Tools section of the
Design Report sub-tab.
l Format… — opens a dialog box where you can set the thickness, style, and
color of the line to be selected.
l Make Vertical / Make Horizontal — toggles the line from horizontal to vertical,
and vertical to horizontal.
l Send to Back / Front — with multiple lines you can layer them with some over-
laying portions of each other with the most important one being on top.
l Specify Size — specify the height and width of your line in inches.
l Fit to Page — scales the line to automatically fit your page size.
l Repeat on every page — places the line on every page of your report.
Adding a Frame
This option creates a blank frame that can be re-sized, placed around objects, and moved, so that
you can decide which portions of a report / template to group together.
To add a frame:
3. Specify the frame's thickness and style from the drop-down lists; then click OK.
You can click-and-drag a plot or image into the frame you just created.
l Format… — opens a dialog box where you can set the thickness and style of the
frame to be selected.
l Send to Back / Front — with multiple lines you can layer them with some overlaying
portions of each other with the most important one being on top.
l Specify Size — specify the height and width of your line in inches.
l Fit to Page — scales the line to automatically fit on the available space on your
page.
l Repeat on every page — places the line on every page of your report.
Adding an Image
To add an image:
l Restore Original Aspect Ratio... — restores the object's original aspect ratio.
l Scale Image
l Copy to Clipboard — you can paste the image into another application.
l Send to Back / Front — with multiple lines you can layer them with some overlaying
portions of each other with the most important one being on top.
l Specify Size — specify the height and width of your line in inches.
l Fit to Page — scales the line to automatically fit your page size.
l Repeat on every page — places the image on every page of your report.
Adding Text
With this option you can create labels or static text in the report template.
To add text:
l Edit / Format Text... — used to change the text in the selected object.
l Copy to Clipboard — you can paste the image into another application.
l Select text from annotation — copies the text from an existing annotation.
l Vertical Justification — select top, center, or bottom for your text's vertical jus-
tification.
l Border Width / Border Style / Border Color — sets the width, style, and color of
the border.
l Background Color — displays a selection of colors for you to change the color of
the text box.
l Send to Back / Front — with multiple text boxes, you can layer them with some over-
laying portions of each other with the most important one being on top.
l Specify Size — specify the height and width of your text box in inches.
l Fit to Page — scales the text box to automatically fit your page size.
l Repeat on every page — places the text box on every page of your report.
Adding a Schematic
To add a schematic:
l Send to Back / Front — with multiple text boxes, you can layer them with some over-
laying portions of each other with the most important one being on top.
l Specify Size — specify the height and width of your schematic in inches.
l Fit to Page — scales the schematic to automatically fit your page size.
l Repeat on every page — places the schematic on every page of your report.
To preview a report:
1. Click the Reports tab. By default, the Preview Report sub-tab will be selected.
2. Click a report template in the Reports pane. The first page of your report is displayed in the
middle of the screen.
3. If your report has more than one page, click the Multiple Pages icon ( ) to view a thumb-
nail view of the report.
4. If you are in the single page view of a multi-page report, click the Next Page ( ) and Pre-
Rearranging Pages
To rearrange the pages of your report:
Deconvolution
Note: The AOF and deliverability plots can be generated at both the wellhead and sandface.
Isochronal Test
A fundamental reason that the conventional test is theoretically sound is that the radius of invest-
igation is constant for each flow period. In order to uphold this principle, the isochronal test takes
An isochronal test is performed by flowing a well at several different flow rates for periods of equal
duration, normally much less than the time required for stabilization. A shut-in, long enough for the
pressure to reach essentially static conditions, is performed between each flow period. In addition,
an extended flow rate, long enough to reach pressure stabilization, is required.
Note: In tight reservoirs, the length of time required to reach pressure stabilization between flow
periods could make the isochronal test impractical.
Extended Flow
Normally an isochronal test includes one flow rate that is extended to stabilization, and a stabilized
pressure and flow rate point is determined. This point is the extended flow pressure and flow rate for
the test. Single point tests do not include the multi-rate portion of a test and consist of only an exten-
ded rate and pressure.
Stabilized Shut-in
Stabilized generally refers to a test in which the pressure no longer changes significantly with time.
For AOF tests, the stabilized shut-in pressure is a pressure that reflects the average reservoir pres-
sure at the time. It is either measured during the test, or it is determined from the interpretation of the
data.
Many tests however, are not flowed to stabilization because of time constraints (especially in tight
reservoirs). An extended flow and stabilized shut-in are still performed at the end of these deliv-
erability tests, so that the buildup data can be analyzed, and from that the stabilized rate calculated.
Stabilized flow can be determined by calculation, or by creating a model of the reservoir, doing a fore-
cast at a specified pressure, and finding the point when the rate has stabilized (usually at 3 months,
6 months, or 1 year).
Types of Analyses
Two types of analysis are available: the simplified analysis and the Laminar-Inertial-Turbulent (LIT)
analysis.
LIT analysis is more rigorous than the simplified analysis, and is usually only used in tests where tur-
bulence is dominant, and the extrapolation to the AOF is large. However, in most cases the sim-
plified analysis is sufficient to determine the AOF and deliverability.
Pressure Method
For both the simplified and LIT analysis, two pressure options are available, the pressure squared or
the pseudo-pressure approach.
Pressure Squared
The pressure-squared approach is the more traditional method, and is often used because it is
easier to understand and calculate. However, it is only valid for medium-to-low pressure ranges, but
is just as accurate as the pseudo-pressure approach in this range.
Pseudo-Pressure
Using pseudo-pressure is more accurate than the pressure-squared approach, especially when deal-
ing with a high-pressure system, where gas viscosity (µg) and compressibility (cg) cannot be
assumed to be constant. Thus, pseudo-pressure works for all pressure ranges, although it is more
difficult to calculate and requires more computational time.
The analysis of a modified isochronal test using the simplified method is shown below. For the mod-
ified isochronal test, pws must be used instead of pR because the duration of each shut-in period is
too short to reach static conditions.
The data is plotted on a log-log plot of ∆p2 versus qst where ∆p2 is defined as:
The flow and shut-in periods of equal duration provide the information required to plot four points. A
straight line, called the transient deliverability line, is drawn through these four points.
The duration of the last flow rate is extended until the pressure response has stabilized. This inform-
ation is used to plot another point called the stabilized point. A line parallel to the transient deliv-
erability line is drawn through the stabilized point. This is called the stabilized deliverability line.
If the extended flow period does not reach pressure stabilization, a stabilized point can be found by
calculation from a buildup test.
The other parameter, C, can be determined using n and the coordinates (qst and pR) of any point on
the stabilized deliverability line (e.g., the stabilized point) as follows:
Note that C and n are considered to be constant for a limited range of flow rates. In theory, it is expec-
ted that this form of the deliverability relationship is used only for the range of flow rates used during
the test. However, in practice it is used indiscriminately for a wide range of rates and pressures.
Note that the pseudo-pressure squared terms (a qst and b qst2) are equivalent to skin due to dam-
age (sd) and skin due to turbulence (sturb). The coefficients a and b are defined in the example
below.
The analysis of an isochronal test using the LIT method is shown below.
As in the simplified analysis, the transient deliverability line is drawn through the four isochronal
points and a parallel stabilized deliverability line is drawn through the stabilized point.
From this equation, the slope of the line is equal to b. The parameter a is determined by rearranging
the above equation to solve for a and then substituting b and the coordinates (qst and ψR) of any
point on the line.
References
"Theory and Practice of the Testing of Gas Wells", L. Mattar, G. Brar, and M. Mumby, Energy
Resources Conservation Board (1978) Third Edition, Chapter 3.
"Gas Reservoir Engineering", J. Lee and R. Wattenbarger, Society of Petroleum Engineers Inc.
(1996) Volume 5, 73 - 74, 173 - 181.
l Superimpose this typecurve with historical rate data to calculate synthetic pres-
sures
l Calculate the error between the calculated pressures and the measured pres-
sures
l Generate a new typecurve and repeat the process until the error between cal-
culated and measured pressures is minimized
The deconvolved typecurve can then be analyzed using conventional diagnostic analysis tech-
niques to determine various reservoir characteristics such as permeability (k), skin (s'), reservoir
size, well location, and more. Deconvolution provides an alternative to conventional diagnostic ana-
lysis, and can show additional flow regime information that would not normally be seen within the
specified time-frame of the buildup test.
Note: Deconvolution is a purely mathematical process and should be used with caution.
The following guidelines show the ideal conditions for deconvolving data to obtain the best possible
results:
l Data should be free of a lot of noise and outliers in both rate and pressure
l Buildup typecurves used for deconvolution need to be consistent with each other
l Use as few flow periods as possible and assess whether or not the typecurve
selected as the initial typecurve (Init TC) is a good initial guess.
l When selecting more flow periods it is important that they are as consistent as
possible. To determine this, select each flow period one at a time and compare
their Deconvolved Typecurves.
l If you have good quality measured point(s) in a deselected flow period, try select-
ing these points individually.
l The overall results can be drastically affected by pi, especially with high rates. If
you have a more reliable value for pi we recommend using it instead of the
defaults provided.
l Avoid using "auto" on either pi or Pu1, especially at the same time, as this can sig-
nificantly reduce performance and produce poor results.
l For noisy data, increase the curvature weighting (wcurv) to make it easier for the
algorithm to match the data.
Common Questions
Why doesn’t the typecurve for a certain flow period show up on the Data Typecurve plot?
l Check to see that the Show TC checkbox on the task bar has been selected for
the flow period.
l The range is too short, or has very few measured pressure points in it, and thus
the typecurve cannot be calculated.
How do I select a different flow period to be the initial guess for the Deconvolved Type-
curve?
Go to the task bar and select the appropriate flow period using the radio buttons in the Init
TC column, or go to the Total Test plot, right-click the desired flow period; then select Initial
Typecurve from the pop-up menu.
References
This deconvolution method is based on the work of Thomas Von Schroeter et al. and Michael M.
Levitan et al. For more information on the concepts of deconvolution, refer to the following papers:
SPE 71574, SPE 77688, SPE 84290, SPE 90680.
Note: This relationship was developed from Darcy’s law for the steady state radial flow of a single,
incompressible liquid.
Vogel IPR
Vogel showed that as depletion proceeds in a solution-gas drive reservoir, the productivity of a typ-
ical well decreases, primarily because the reservoir pressure is reduced, and because the increas-
ing gas saturation causes greater resistance to oil flow.
The result is a progressive deterioration of the inflow performance relationship. This is shown in the
following diagram, which compares the straight line IPR with the Vogel IPR.
The purpose of analyzing afterflow data is to determine the wellbore storage constant (C), and to
determine when afterflow ends and reservoir-dominated data begins. Thus, the wellbore storage
constant (C) is a measure of the storage capacity associated with the wellbore volume (Vw).
In the fluid-filled case, the wellbore is presumed to be full of either all liquid or gas, which gets com-
pressed as afterflow continues into the wellbore. The wellbore storage constant (C) is a function of
the wellbore volume (VW) and the wellbore fluid compressibility (cwb) as shown below:
When both gas and liquids co-exist in the wellbore and the liquid level changes, the wellbore storage
constant (C) becomes a function of the wellbore area (Awb) and liquid gradient (Gliq) as shown
below.
or
The wellbore storage constant (C) is often converted to a dimensionless wellbore storage constant
(CD) for use in analytical reservoir models.
For gas:
These equations are linear with respect to time and, as a result, afterflow (wellbore storage) data will
form a straight line when plotted on a Cartesian plot of pressure versus delta time (p vs. ∆t) as shown
below.
Using the slope of this line (m), the wellbore storage constant (C) can be determined. The table
below summarizes the equations used to calculate C for different fluid types and flow periods.
Buildup
Oil
Buildup
Gas
Starting with the constant rate solution (see above), the derivative with respect to the logarithm of
time is:
For gas:
For gas:
Using any derivative point on this line, the wellbore storage constant (C) can be determined. The
table below summarizes the equations used to calculate C for different fluid types and flow periods.
Oil
Buildup
Gas Drawdown
For gas:
Bilinear fracture flow data forms a straight line when placed on a semi-log plot of pressure versus a
quad root time function (p vs. (time function)0.25) for a buildup as shown below.
The slope (m) of this line is used to calculate root fracture flow capacity (sqrt(kfwf)) as follows:
For gas:
The table below summarizes the equations used to calculate root fracture flow capacity (sqrt(kfwf))
for different fluid types and flow periods.
Oil
Buildup
Drawdown
Ga-
s
Buildup
Derivative Analysis
The signature of Bilinear fracture flow data on a derivative plot is a straight line with a slope of 1/4.
The position of this line is used to calculate root fracture flow capacity (sqrt(kfwf)). Note that root frac-
ture flow capacity increases as the position of this line moves to the right.
For gas:
For gas:
This result is linear with time, and as a result, the derivative of bilinear fracture flow data falls on a straight
line with a slope equal to 1/4 on a log ∆p versus log ∆t plot as shown below.
For gas:
Oil
Drawdown
Gas
Buildup
Due to the fact that the derivative is often noisy and that information can be lost from over smooth-
ing, it is recommended that the standard analyses are also used to fine-tune the parameter values
obtained from the derivative analysis.
Introduction
By definition, the derivative is the slope (m) of the data when plotted on semi-log coordinates (see
radial flow analysis) as follows:
There are several numerical techniques available to calculate a derivative. Two methods, standard
and Bourdet, are available within the software and are described below. Both methods incorporate a
smoothing algorithm to reduce noise in the derivative. Be aware that smoothing should be used with
caution, since over smoothing the derivative can completely change the shape, and mask or show
flow regimes that may not exist. Therefore, it is recommended that a minimum amount of smoothing
be applied, so that results are in a recognizable shape of the derivative curve.
Standard Derivative
The standard definition of the semi-log derivative is as follows:
This definition assumes two points along the semi-log data curve are used. By applying a least
squares fitting technique, this definition can be expanded to allow for more than two points (up to the
total number of data points) to be used to calculate the derivative.
Using this definition, ∆X represents the log-cycle fraction used to control the amount of smoothing.
Note this value is typically quite small (0.01 to 0.2) and a small increase represents a significant
increase in smoothing.
References
1. "Derivative Analysis Without Type Curves", Louis Mattar, JCPT Special Edition (1999) vol.
38.
2. "Use of Pressure Derivative in Well-Test Interpretation", D. Bourdet, J.A. Ayoub, and Y.M.
Pirard, SPE Formation Evaluation, (June 1989).
The permeability (k) and skin (s') are obtained from a radial analysis and reservoir shape and drain-
age area are specified. Using this information, a vertical model is used along with a material balance
equation to rigorously solve for the final average reservoir pressure.
The elliptical flow regime exhibits a slope of 1/3 on a log-log plot as shown below.
The position of this line is used to calculate the square root of the horizontal permeability ratio (sqrt
(ky / kx)) as follows:
The semi-log pressure derivative during the elliptical flow period in dimensionless form is:
In dimensional parameters, for an oil reservoir, the semi-log pressure derivative becomes:
Where:
Upon re-arranging the second equation shown above, the square root of the horizontal permeability
ratio can be expressed as follows:
For gas:
Where:
To calculate the square root of the horizontal permeability ratio, the above equations are used along
with a ∆pw' or ∆ψw' value read from the 1/3 slope line on the log-log plot at a corresponding ∆t
value. kxy is obtained from the horizontal radial analysis line.
References
"Determination of Horizontal Permeability Anisotropy from Horizontal Well Tests", M.B. Issaka, K.
Zaoral, A.K. Ambastha, L. Mattar, Paper presented at the 2000 SPE Saudi Arabia Section Technical
Symposium, Dahran, Saudi Arabia, October 21 - 23.
The equations and theory for analyzing horizontal radial flow data are identical to that for radial flow.
The only difference is the permeability obtained and shown in the equations is replaced by horizontal
permeability (kxy).
It is important to note that in a horizontal well, two radial flow regimes occur. One at early time which
occurs in the vertical plane of the formation, and the second at late time in the horizontal plane.
Thus, it is important to analyze the proper region of data to obtain correct results. The plots below
show the distinction between vertical and horizontal radial flow on the semi-log and derivative plots,
as well as the correct placement of the analysis lines.
For gas:
Linear channel flow data will form a straight line when placed on a semi-log plot of pressure versus a
square root time function (p vs. sqrt (time function)) for a buildup as shown below.
For gas:
The permeability can be obtained from the radial flow analysis, or estimated from core data, or other
tests. After permeability is determined, channel width is found using the following equation:
Drawdown
Oil
Buildup
Drawdown
Gas
Buildup
Derivative
The signature of linear channel flow data on a derivative plot is a straight line with a slope of 1/2,
between the middle and late time flow regions. The position of this line is used to calculate channel
width and the square root of permeability (w sqrt(k)).
Note: The channel width and the square root of permeability increases as the position of this line
moves to the right.
Starting with the constant rate solution (see above) the derivative with respect to the logarithm of
square-root time is:
For gas:
At long flow times, this result is linear with time, and as a result, the derivative of linear channel flow
data falls on a straight line with a slope equal to 1/2 on a log ∆p versus log ∆t plot as shown below.
Note that for shorter flow times linear equivalent time should be used to maintain a 1/2 slope.
Using any derivative point on this line, the product of channel width and root permeability (w sqrt(k)).
can be determined as follows:
For gas:
Oil
Buildup
Drawdown
Gas
Buildup
References
"Pressure Transient Analysis of Wells in Very Long Narrow Reservoirs", R. Nutakki, L. Mattar; Soci-
ety of Petroleum Engineers Paper 11221, 1982.
For gas:
Linear fracture flow data will form a straight line when placed on a semi-log plot of pressure versus a
square root time function (p vs. sqrt (time function)) for a buildup as shown below.
For gas:
The permeability can be obtained from the radial flow analysis, or estimated from core data, or other
tests. After permeability is determined, fracture half-length is found using the following equation:
Oil
Buildup
Drawdown
Gas
Buildup
Starting with the constant rate solution (see above), the derivative with respect to the logarithm of
square-root time is:
For gas:
For gas:
Using any derivative point on this line, the product of fracture half-length and root permeability
(Xfsqrt(k)). can be determined as follows:
For gas:
Drawdown
Oil
Buildup
Drawdown
Gas
Buildup
For gas:
Linear horizontal flow data forms a straight line when placed on a semi-log plot of pressure versus a
square root time function (p vs. sqrt (time function)) for a buildup as shown below.
For gas:
Permeability in the y-direction (ky) can then be directly calculated by taking the square of the result
shown above.
Drawdown
Oil
Buildup
Drawdown
Gas
Buildup
Derivative
The signature of linear horizontal flow data on a derivative plot is a straight line with a slope of 1/2,
between the vertical radial and elliptical, or horizontal radial flow regions. The position of this line is
used to calculate the square root of permeability in the y-direction (sqrt(ky)).
Note: The square root of permeability in the y-direction increases as the position of this line moves
to the right.
Starting with the constant rate solution (see above) the derivative with respect to the logarithm of
square-root time is:
For gas:
This result is linear with time and, as a result, the derivative of linear horizontal flow data falls on a
straight line with a slope equal to 1/2 on a log ∆p versus log ∆t plot as shown below.
Using any derivative point on this line, the square root of permeability in the y-direction (sqrt(ky)) can
be determined as follows:
Oil
Buildup
Drawdown
Gas
Buildup
For gas:
These equations are linear with respect to time, and as a result, pseudo-steady state flow data forms
a straight line when plotted on a Cartesian plot of delta pressure over normalized rate versus delta
time (∆p / q vs. ∆t) as shown below.
For gas:
The reservoir pore volume is then used to calculate the original hydrocarbons in place as follows:
Draw-
down
Gas
Derivative
The signature of pseudo-steady state flow data on a derivative plot is a straight line with a unit slope
at late time. The position of this line is used to calculate reservoir pore volume (Vp) and the original
hydrocarbons in place (OOIP or OGIP). Note that reservoir pore volume (Vp) and thus original hydro-
carbons in place (OOIP or OGIP) increases as the position of this line moves to the right. The
pseudo-steady state flow analysis cannot be applied to buildup or falloff test data.
Starting with the constant rate solution (see above) the derivative with respect to the logarithm of
time is:
For gas:
This result is linear with time, and as a result, the derivative of pseudo-steady state flow data on a
log-log plot is a straight line with unit slope on a log ∆p/q versus log ∆t plot as shown below.
Notice that both wellbore storage (early time) and pseudo-steady state (late time) exhibit tank beha-
vior and have a signature of a unit slope on the derivative plot.
Using any derivative point on this line, the pore volume (Vp in ft3) and original hydrocarbons in place
(OOIP (stbbls) or OGIP (MMscf)) can be determined as follows:
Oil
Drawdown
Ga-
s
For gas:
Radial flow data will form a straight line when placed on a semi-log plot of pressure versus a log time
function (p vs. log(time function)) for a buildup as shown below.
For gas:
Skin
Drawdown Tests
Permeability along with the slope (m) and position of the line is used to calculate apparent or total
skin (s'). The equation for total skin (s') is found by rearranging the radial flow equation as follows:
The values of pwf and t can be taken from any point on the semi-log straight line (not a measured
data point).
Buildup Tests
For a buildup test, total skin (s') cannot be determined from shut-in pressure alone. It can still be
found however, by making use of the information provided by any flowing pressure point. By con-
vention, the final flowing pressure (pwfo) is used.
Any pressure off the semi-log line (not the measured data) can be used, corresponding to any spe-
cified shut-in time. In most texts the pressure at 1 hour is used. In the software, the extrapolated
buildup pressure (p*) at infinite shut-in time (t = 1.0) is used in the total skin (s') calculation. Using this
convention, the equation used for calculating total skin (s') from a buildup test is as follows:
For gas:
O-
il
Buildup
Draw-
down
G-
a-
s Buildup
Derivative
The signature of radial flow data on a derivative plot is a horizontal straight line. The position of this
line is used to calculate permeability (k), and apparent or total skin (s'). Note that an increase in the
vertical position of the horizontal straight line corresponds to a decrease in permeability.
Starting with the constant rate solution (see above) the derivative with respect to the logarithm of
time is:
At long flow times, this result is linear with time, and as a result, the derivative of radial flow data falls
on a horizontal straight line with a slope equal to 0 on a log ∆p versus log ∆t plot as shown below.
Note that for shorter flow times, radial equivalent time should be used to maintain a horizontal slope.
Permeability
The value of the derivative of radial flow data corresponds to the vertical position of the horizontal
straight line and is used to calculate permeability (k) as follows:
For gas:
O-
il
Buildup
Drawdown
G-
a-
s Buildup
Spherical flow analysis can be used prior to the start of the infinite-acting radial flow regime at very
early times. To apply the spherical flow model in the case of production through a cylindrical well-
bore, an equivalent spherical wellbore radius (rs) needs to be applied, which acts as the wellbore
radius (rw) for all calculations involving spherical flow. This is calculated as follows:
During the spherical flow regime, the pressure drawdown occurring with time can be expressed as:
The above equation can be manipulated to show the expression for spherical permeability as:
Where:
Spherical flow data forms a straight line when placed on a semi-log plot of pressure versus an
inverse square root time function (p vs. 1 / sqrt (time function)) for a buildup as shown below. m is
defined as the slope of this line. Note that to perform a spherical flow analysis, the value of hp must
be specified to calculate the values of ks and rs.
The equations and theory for analyzing vertical radial flow data are identical to that for radial flow.
The only difference is the permeability obtained and shown in the equations is replaced by vertical
permeability (kxy), and net pay (h) is replaced with effective wellbore length (Le).
It is important to note that in a horizontal well, two radial flow regimes occur. One at early time, which
occurs in the vertical plane of the formation, and the second at late time in the horizontal plane.
Thus, it is important to analyze the proper region of data to obtain correct results. The plots below
show the distinction between vertical and horizontal radial flow on the semi-log and derivative plots,
as well as the correct placement of the analysis lines.
Pre-closure
After-closure
Step down
Pre-closure analysis consists of identifying closure and analyzing the early pressure falloff period
while the induced fracture is closing. One of the most critical parameters in fracture treatment design
is the closure pressure. One specialized analysis technique used to identify closure is, in pressure
transient parlance, the G-Function Analysis.
l Instantaneous Shut-In Pressure (ISIP) = Final injection pressure - Pressure drop due
to friction
l Net Fracture Pressure (Δpnet) – Net fracture pressure is the additional pressure
within the frac above the pressure required to keep the fracture open. It is an indic-
ation of the energy available to propagate the fracture.
l Fluid efficiency – Fluid efficiency is the ratio of the stored volume within the fracture to
the total fluid injected. A high fluid efficiency means low leakoff and indicates the
energy used to inject the fluid was efficiently utilized in creating and growing the frac-
ture. Unfortunately, low leakoff is also an indication of low permeability. For minifrac
after-closure analysis, high fluid efficiency is coupled with long closure durations and
even longer identifiable flow regime trends
A straight-line trend of the G-function derivative (Gdp/dG) is expected where the slope of the deriv-
ative is still increasing. Position the Fracture Closure Identification line, which is anchored to the ori-
gin by default, through the straight-line portion of the G-Function derivative. Fracture closure is
identified as the point where the G-Function derivative starts to deviate downward from the straight
line as shown in the following figure.
Note: The shape of the G-function derivative prior to closure qualitatively describes how fluid
moves from the fracture into the formation (see Fluid Leakoff Types).
l The equation for α = 1.0 is for low leakoff, or high efficiency where the fracture area
open after shut-in varies approximately linearly with time.
l The equation for α = 0.5 is for high leakoff, or low efficiency fluids where the fracture
surface area varies with the square-root of time after shut-in.
One of the key variables identified by Nolte is the difference between a high efficiency (upper limit)
and a low efficiency (lower limit) leak-off condition. These two conditions have no significant effect
on the overall shape of the curves.
Normal Leakoff
Normal leakoff occurs when the fracture area is constant during shut-in and the leakoff occurs
through a homogeneous rock matrix.
l The G-Function derivative (G dP/dG) lies on a straight line that passes through
the origin.
The fracture closure point can be identified when the G-Function derivative (G dP/dG) starts to devi-
ate downward from the straight line. The time and pressure corresponding to this point are identified
from a marker set at this point.
l The portion of the normal leakoff lies on a straight line passing through the origin.
As shown below, the fracture closure point is identified from the G-Function derivative G dP/dG
when it starts to deviate downward from the straight line.
1. Transverse storage occurs when the main fracture intercepts a secondary fracture network,
which could be natural or induced. This differs from PDL in that the dominant effect of the
secondary fractures is to provide pressure support to the main fracture, rather than addi-
tional surface area for leakoff. There can be cases where transverse storage (pressure sup-
port) dominates, followed by a period of PDL before closure of the main fracture occurs.
2. Fracture height recession occurs if the fracture propagates through adjoining impermeable
layers (above or below the pay zone) during injection. In the normal leakoff scenario, fluid
can leak off from the entire surface area of the fracture. For fracture height recession, leakoff
can only occur in the portion of the fracture which is in communication with the permeable
zone. As a result, the leakoff rate is slower than the normal case. Eventually, the fracture
area in the impermeable layer(s) starts closing (height recession), and during this period the
rate of pressure decline increases. After the fracture height recedes to the edge of the per-
meable zone, the entire area of the frac contributes to leak off, and a period of normal leakoff
ensues.
The three characteristic signatures for height recession during shut-in are:
l The G-Function derivative G dP/dG initially exhibits a large positive slope that
continues to decrease with shut-in time, yielding a concave-down curvature.
l Any straight line fit through the G-Function derivative G dP/dG intersects the y -
axis above the origin.
Until the main fracture closes, the G-Function derivative behaves similarly to PDL, and it is difficult to
distinguish between PDL and fracture tip extension. The following plot shows a typical response for
fracture tip extension after shut-in.
l ISIP = Final injection pressure - Pressure drop due to friction in the wellbore and per-
forations or slotted liner
ISIP can be a difficult value to quantify and, as a result, there are numerous ways to obtain an estim-
ate of ISIP in WellTest.
l If sandface pressures have been calculated, a default value for ISIP is calculated by
subtracting the friction component of the sandface calculation from the final injection
pressure. This tends to overestimate the value of ISIP because it doesn't account for
friction through the perforations or near the wellbore.
l It is also common practice to estimate ISIP by placing a straight line on the early fal-
loff portion of the history plot. After shut-in, the friction decreases rapidly and this can
appear to cause a step drop in pressure, or a brief linear trend in the data (prior to the
expected concave-up-trend of the falloff).
After an ISIP arrow has been added to the plot, the pressure value from the arrow will be
added to the list of ISIP defaults.
If a closure arrow is placed on the G-function plot, ISIP can be calculated from the following
equation, which assumes normal leakoff.
Closure pressure is equal to the minimum in-situ stress because the pressure required to open a
fracture is the same as the pressure required to counteract the stress in the rock perpendicular to the
fracture. Closure pressure is determined from the G-Function or the Sqrt(t) plot.
It is possible to obtain reasonable estimates of reservoir pressure and formation permeability from
diagnostic analyses if the pressure falloff is recorded long enough to achieve radial flow. Minifrac
analytical models provide the significant advantage of verifying and improving on results obtained
from diagnostic analyses, similar to the analysis workflow used for many years in traditional PTA.
With modeling, you can match transitions between flow regimes and you do not need to depend on
the identification of fully developed flow regimes.
Nolte
This after-closure analysis method is based on the work of K.G. Nolte, and expanded on by R.D. Bar-
ree. It is based on the solution of a constant pressure injection followed by a falloff. The impulse
equations are obtained by approximating the injection duration as very small. Nolte's techniques use
injected volume as the impulse volume and consider the closure point as the beginning of the falloff.
If impulse radial flow is identified, the pressure data overlays the derivative data. This diagnostic may
be used to help fine-tune the estimate of initial pressure.
There are two MiniFrac ACA (Nolte) Linear Time functions. FL1 is an approximation of FL2. The
above equation holds regardless of the linear time function used.
Impulse linear flow shows up as a straight line on a plot of pressure vs. FL (FL1 or FL2).
After closure pressure and permeability have been specified, the position of the impulse linear line is
used to calculate p* and fluid loss coefficients, CR and CT. Both MiniFrac ACA (Nolte) Linear Time
functions are bounded by 1.0 (at early shut-in time) signifying the time of closure and 0.0 (at late
time) representing an infinite shut-in time. P* is the extrapolation of the linear line to FL= 0, which cor-
responds to an infinite shut-in time.
p@tc is the extrapolation of the linear flow line at closure time, tc.
Note: FL = 1.0 at tc
The MiniFrac Nolte ACA Radial Time function, FR1, has a minimum value of 0.0 (at late time) rep-
resenting an infinite shut-in time. p* is therefore obtained from the extrapolation of the linear line to
FR1= 0. As described in Reference 1 (Barree, 2007), the equation governing impulse radial flow is
defined as:
The radial flow portion of the after-closure falloff data has a straight-line trend on a plot of pressure
vs. FR1.
The slope of the line on the MiniFrac Radial plot is mR1 and permeability may be calculated as fol-
lows:
With the position of the impulse radial line on the Derivative plot, you can form a similar calculation:
D. Craig developed an analytical model that accounts for fracture growth, leakoff, closure, and after-
closure. The late-time approximation of his model produced impulse radial and bilinear flow equa-
tions that are the same as Soliman's solutions. However, the constant in Craig's impulse linear flow
equation differs from Soliman's. The linear flow equation implemented in the Soliman / Craig ana-
TM
lysis method in WellTest uses the constant obtained by D.Craig.
By default, the delta pressure curve is not displayed on this plot. It can be added by right-clicking the
legend and selecting to view the delta pressure curve from the customization menu. The delta pres-
sure curve on the Derivative plot is defined as the measured falloff pressure minus the initial pres-
sure. The initial reservoir pressure defaults to the p* (extrapolated pressure) of the first late-time
analysis line added. This default may be overridden manually by typing a value into the Initial Pres-
sure field of the Parameters dialog box, or by using the default value from another analysis line by
clicking the Default button beside the Initial Pressure field in the Parameters dialog box. If the initial
pressure is linked to an analysis line, and if delta pressure is displayed on the typecurve, the delta
pressure updates as the line is moved.
A bilinear analysis line can be added from the Analysis Options dialog box. The bilinear line should
then be positioned through the bilinear flow portion of the falloff data (identified by a straight line on
the MiniFrac bilinear plot, or a 1/4 slope on the Derivative plot). When the bilinear line is selected, the
Analysis Options dialog box displays a field for entering permeability, as shown below. If radial flow
developed, permeability could be estimated from the Radial analysis. Otherwise, a best estimate of
permeability must be entered. After permeability is entered, fracture conductivity is calculated.
A linear analysis line can be added from the Analysis Options dialog box. The linear line should then
be positioned through the linear flow portion of the falloff data (identified by a straight line on the Min-
iFrac linear plot, or a 1/2 slope on the Derivative plot). When the linear line is selected, the Analysis
Options dialog box displays a field for entering permeability, as shown below. If radial flow
developed, permeability could be estimated from the Impulse Radial analysis. Otherwise, a best
estimate of permeability must be entered. After permeability is entered, fracture half-length is cal-
culated.
The MiniFrac Soliman / Craig ACA Radial Time function has a minimum value of 0.0 (at late time)
representing an infinite shut-in time. p* is the extrapolation of the radial line to the time of 0, which
corresponds to an infinite shut-in time.
The radial flow portion of the after-closure falloff data has a straight-line trend on a plot of
pressure vs.
A radial analysis line can be added from the Analysis Options dialog box. The radial line should then
be positioned through the radial flow portion of the falloff data (identified by a straight line on the Min-
iFrac radial plot, or a zero slope on the Derivative plot).
To perform this test, a fluid of known properties (for example, water) is injected into the formation at a
rate high enough to initiate a small frac. The injection rate is then reduced in a stair-step fashion,
each rate lasting an equal time interval, before the well is finally shut-in. The resulting pressure
response caused by the rate changes is influenced by perforation and near-wellbore friction. Tor-
tuosity and perforation friction pressure losses vary differently with rate. By analyzing the pressure
losses experienced at different rates, we can differentiate between pressure losses due to tortuosity
and due to perforation friction.
Pressure drops across perforations and due to tortuosity are given mathematically by the following
equations:
where
For step-down tests, it is essential to keep as many variables controlled as possible, so that the pres-
sure response during the rate changes is due largely to perforations and tortuosity, and not some
Step-down test analysis is done by plotting the pressure / rate data points with the same time since
the last rate change on a pressure-rate plot, and matching the pressure loss model (given by the
equations above) to these points. On the basis of the model, the perforation and tortuosity com-
ponents of the pressure loss are calculated, and the defining parameters are also estimated.
The step-down test analysis module in WellTest is as shown in the figure below. The module com-
prises four sections. The bottom-left section houses the history plot. The top-left section contains the
data grid where the user inputs data points to use in the analysis. Data points can be entered either
manually in the grid, or by using the “Point Selection” arrow to pick desired points in the history plot.
The top-right section is where you enter initial estimates of the parameters used in the pressure loss
model. The pressure fracture (pfrac) parameter initial estimate should be set to the instantaneous
shut-in pressure (ISIP). The bottom-right section serves two purposes, depending on the active tab;
the Zoom tab is used to fine-tune point selection, while the Step Down Plot tab displays the analysis
results graphically.
PITA
Slug Analysis
Introduction
Conventional well testing has served the industry faithfully for decades as the primary and most reli-
able means of:
l Quantifying deliverability
l Characterizing a reservoir
Conventional well tests are performed by producing reservoir fluid either into a gathering system, or
by venting / flaring to the atmosphere. The production period, which can be a few minutes to several
days or even years, is usually followed by a buildup period. Typically, the analysis is focused on the
pressure buildup data, which requires the input of the production history prior to the buildup test.
An alternative to the conventional test is one that simply allows the well to flow into the closed well-
bore after perforating (closed chamber test), which results in a pressure buildup. The pressure
buildup data is collected for a period of hours or days, depending on the reservoir’s flow
potential. These tests have been variously referenced to as: slug test, surge test, perforation inflow
diagnostic (PID), or closed chamber test (CCT). Following a critical review of the literature published
on these unconventional tests, we have developed a complete and systematic analysis that yields
an estimate of initial reservoir pressure (pi), permeability (k), and skin (s'). We have called this pro-
cedure perforation inflow test analysis (PITA).
The following subsections illustrate the basis and the process of how PITA works, for a gas well per-
forated in an under-balanced condition (initial cushion pressure (pwo) < initial reservoir pressure
(pi)). In this case, the formation fluid enters the wellbore (inflow). However, the same methodology
applies to a perforation test done in an over-balanced condition, where initial cushion pressure (pwo)
> initial reservoir pressure (pi). In such a case, wellbore fluid is actually injected into the reservoir dur-
ing the test. The objective of PITA is to extract key reservoir parameters (initial reservoir pressure
(pi), permeability (k), and skin (s')) from the pressure data.
In gas wells, the pressure is usually measured at the wellhead and converted to sandface. This con-
version is primarily due to hydrostatic head because the influx rate into the wellbore diminishes rap-
idly (and friction is negligible). The influx flow rate is not measured but can be estimated using closed
chamber calculations, provided the assumption of single-phase flow can be justified. The figure
below shows the typical profiles of measured pressure and the calculated gas influx rate for a per-
foration test of a gas well. As shown here, the influx rate declines rapidly.
If a reservoir is expected to produce liquid, the pressure must be measured downhole (near the per-
forations). For oil reservoirs, wet coals or aquifers that have little or no associated gas production at
the time of perforating, single-phase flow can be assumed, and the influx flow rate can be estimated
using closed chamber calculations.
A typical derivative using the above equation is shown in Figure 2. It illustrates that the early time
data (wellbore storage) has a slope of 2 for a positive skin (well test derivative has a slope of 1), and
the late time data (reservoir flow) has a slope of 0 (flat line, same as the conventional well test deriv-
ative).
Note: For a zero or negative skin, the derivative will yield an initial slope of 1.75.
Using this convention for the derivative, it is easy to recognize if reservoir flow exists. If it does, initial
reservoir pressure (pi) and permeability (k) can be determined. At least some of the late time data
should fall on the flat part of the derivative, as shown in Figure 2, to obtain a reliable analysis. If this
data exists, then skin (s') can also be calculated from the early time data. If reservoir flow is not
apparent, then a unique interpretation for these reservoir characteristics is not possible.
If a formation is partially penetrated, and has low permeability (k), only the pay thickness exposed to
the wellbore should be considered in the analysis. However, in some cases where perforation plug-
ging restricts communication with the formation, using the total pay thickness can yield a low estim-
ate of permeability (k).
Note: The fundamental equations were originally derived for liquid flow in terms of time and pres-
sure. To apply these equations to gas flow, time and pressure are replaced with pseudo-
time (ta) and pseudo-pressure (ψ) respectively.
Analysis
In traditional well test interpretation, we start analyzing the data from early time to late time. In PITA,
we start with the late time analysis first to obtain initial reservoir pressure (pi) and permeability (k).
Then the early time data (where the derivative slope is 2 or 1.75) is analyzed and skin (s') is cal-
culated using this analysis along with the permeability (k) obtained from the late time analysis.
A summary of the procedure for analyzing a perforation inflow test for a gas well is presented below:
4. Calculate the impulse derivative and plot versus pseudo-time on a log-log scale (see Figure
2).
5. Determine the start of late time reservoir dominated flow (Figure 2).
6. Analyze the late time data by plotting pseudo-pressure (ψ) versus inverse pseudo-time (1 /
∆ta) on reverse Cartesian coordinates (see Figure 3). A straight line through the valid late
time data gives permeability (k) (from slope (m)) and initial reservoir pseudo-pressure (ψi)
(from intercept) using the following equation:
8. Analyze the early time data by plotting pseudo-pressure versus pseudo-time on Cartesian
coordinates (see Figure 4). Based on the early time slope on the derivative (Figure 2), use
the appropriate equation below along with a ψw and ∆ta read from a point on the early time
analysis line to calculate skin (s') (Figure 4).
In terms of the slope, of the "early time negative skin" analysis line on the Early Time plot.
After rwa is calculated, the skin (s') is found from the following equation:
After the key reservoir parameters (initial reservoir pressure (pi), permeability (k), and skin (s')) are
estimated from PITA, a PITA model is used to confirm the estimated parameters. This is especially
critical when reservoir dominated (radial) flow is not achieved within a test period, or when data scat-
ter aggravates the analysis.
Experience shows that it can be difficult to differentiate between a small positive skin (s') and a zero
or negative skin (s'), based on the early time analysis alone. Under such circumstances, the PITA
model provides a more reliable indication of the skin (s') value, if reservoir dominated flow is
achieved, or at least, if a transition to reservoir dominated flow is developed.
References
1. "Use of PITA for Estimating Key Reservoir Parameters", N.M. Anisur Rahman, Mehran
Pooladi-Darvish, Martin S. Santo and Louis Mattar, Paper CIPC 2006 - 172 presented at
2006 Canadian International Petroleum Conference, Calgary, Alberta, Canada, June 13 -
15.
2. "Development of Equations and Procedure for Perforation Inflow Test Analysis (PITA)", N.M.
Anisur Rahman, Mehran Pooladi-Darvish and Louis Mattar, Paper SPE 95510 presented at
2005 ATCE, Dallas, TX, October 9 - 12.
Note: For additional information, see Closed Chamber Inflow Test Type.
Theory
Just like in PITA, closed chamber test analysis is based on the principles of slug tests, and incor-
porates procedures for both liquid and gas flow, using the concept of a stepped change in wellbore
storage. Early time data (during the flow period) is used to estimate skin (s'), and late time data (dur-
ing the buildup period) is used to estimate initial reservoir pressure (pi) and permeability (k), similar
to PITA. Pressures are measured downhole, below the shut-in device (see Figure 1).
Flow Regimes
Just like in PITA, the critical part of the analysis is distinguishing between the data that is dominated
by wellbore storage (afterflow), and the data that is dominated by the reservoir response (radial
flow). This is accomplished by constructing a derivative plot similar to that used for PITA. If downhole
shut-in is achieved early in the test (during the wellbore storage period), the CCT derivative will
exhibit an immediate decrease in wellbore storage, and a quick transition to reservoir (radial) flow fol-
lows, as shown in below.
Figure 4 shows the late time analysis plot used to calculate initial reservoir pressure (pi) and per-
meability (k) for a gas zone. After estimates of initial reservoir pressure (pi) and permeability (k) are
obtained, an early time analysis is performed to calculate skin (s'), as depicted in Figure 5.
Note: The skin (s') calculations are performed using the same equations developed for PITA.
After the key reservoir parameters (initial reservoir pressure (pi), permeability (k), and skin (s')) are
estimated from the CCT analysis, a CCT model is performed to confirm the estimated
parameters. This is especially critical when reservoir-dominated (radial) flow is not achieved within a
test period, or when data-scatter aggravates the analysis.
In WellTest, the PITA analysis has been developed further than slug. PITA facilitates the estimation
of p*, reservoir permeability and skin whereas Slug's Impulse Radial Analysi only provides an estim-
ation of p*. Therefore, PITA is recommended for analyzing data from slug tests. Both PITA and Slug
tests have (identical) modeling capabilities in WellTest.
Transition period flow data will form a straight line when placed on a Cartesian plot of pressure
versus reciprocal time (p vs. 1 / ∆t) as shown below.
The basis of this approach is illustrated with the following mathematical development.
The wellbore flowing pressure (pwf) (during fluid influx) can be derived as:
Thus, a plot of p versus 1 / ∆t will result in a straight line, and the intercept of which on the pressure
axis provides the initial pressure (pi).
In order to select or confirm the range of data that can be analyzed by the impulse radial analysis, a
log-log plot of the semi-log derivative versus delta time (∆t) can be used as shown below.
Note: Any data that falls within the range of the -1 slope line can be used for analysis.
References
"A New Method for Determination of Average Reservoir Pressure", F.J. Kuchuk, SPE Paper 56418
Presented at 1999 ATCE, Houston, TX (October 3 - 6).
Overview
Advanced Models
Legacy Models
Models are used to simulate your data using the process of history matching pressure transient data
based on mathematical principles. There are many different models available to match data, depend-
ing on the situation. Thus, it is important to analyze the pressure transient data before modeling
because it forces you to think about the probable reservoir configurations, and provide good estim-
ates of reservoir parameters. Models are not unique (different model types can match the same set
of data), and as a result, we recommend that the choice of model type occur after the analysis step.
l Modeling makes use of all the information within a dataset. While analyzing data,
the analyst might try and determine permeability and skin by analyzing the data
points which make up the zero slope on the derivative plot and the semi-log
straight line on the radial plot, but ignore the data points in the transition period
between wellbore storage and radial flow. Models make use of the information
contained in transition periods.
l Modeling takes all flow regimes into account. In multi-rate situations, analyses
depend on the superposition of the equation for a single-flow regime. For
example, the derivation of Horner time includes the assumption that all flow,
including the entire drawdown, is radial. Modeling doesn’t assume that only one
flow regime has occurred.
l With modeling, you can simultaneously analyze multiple flow periods, so that a
single set of parameters can be found.
Parameter values obtained during the analysis step provide a good starting point for an appro-
priately chosen model type. Parameters can then be optimized by automatic parameter estimation
(APE). Before using the APE method, corrupted data should be removed from the dataset to prevent
the attempted match of invalid points.
Model Assumptions
There are many assumptions that go into the model itself, which can lead or mislead the analyst.
Most models assume that a reservoir is homogeneous (dual porosity excluded). In nature, there isn’t
a single reservoir that is actually homogeneous, but many reservoirs behave as homogeneous reser-
voirs.
Analytical models assume single-phase flow. They support the simultaneous analysis of multiple
flow periods. Data from different flow periods may be viewed simultaneously on the derivative and
specialized plots. The analytical models support drawdown / buildup and injection falloff test types.
The numerical models have automatic grids, designed to provide an ease-of-use similar to the ana-
lytical models, while accounting for the complexities of multiphase flow. Currently, the numerical
models are only functional for drawdown / buildup test types.
The hybrid model is essentially a numerical model, but with certain modifications to significantly
reduce computation time, so that it is almost as fast as an analytical model.
Performing an analysis using a Hybrid Model is very similar to performing an analysis using an ana-
lytical model. However, the calculations used by the hybrid model are more rigorous, as they accur-
ately account for changes in gas properties. To view the calculations used in the hybrid model, see
Hybrid Model Theory.
Prerequisites
Before you start modeling reservoirs, ensure that you have all the required information in your
WellTest project:
l Well production data (rates and sandface pressures) must be populated in the Pro-
duction Editor tab
With the composite model, you can add any number of different cylindrical zones. With an unlimited
number of zones, virtually any pressure or rate transient response can be matched; therefore, it is
important to exercise good judgment to determine when this model is appropriate. No reservoir is
perfectly cylindrically concentric composite in nature; however, many reservoirs do behave the same
way as composite reservoirs do. Some common situations where the composite model is useful
include injection cases (which cause changes in viscosity and compressibility), reservoir het-
erogeneities (such as changes in flow capacity (kh)), and cases where the well was drilled into a nat-
urally fractured reservoir with varying fracture distribution.
Note that the following assumptions apply when using this model:
l Both the inner and outer regions can be assumed to be homogeneous or dual poros-
ity
Note that the effective wellbore length (Le) defines the wellbore area open to fluid flow. The cyl-
indrical source solution is used at very early times, which is followed by Green’s function solutions
for horizontal wells, as developed by Thompson et al. (1991). No-flow boundaries are modeled using
the method of images. The result is superposed in time based on the rate history provided. The fol-
lowing flow regimes can be handled by this model:
l Wellbore storage
l Elliptical flow
The water drive model is an analytical radial composite reservoir model. The reservoir is rep-
resented by the inner zone, and is given a radius of "re". The outer zone represents the aquifer, with
a radius of "raq". The aquifer radius may be set to any value greater than "re".
The model accounts for the mobility difference between the reservoir and aquifer by giving you the
option to enter a mobility ratio (M). This is defined as the ratio of the aquifer mobility to that of the
reservoir.
It is a rectangular model that contains an inactive horizontal well fed by multiple identical and
equally-spaced transverse fractures. The portion of the reservoir between the fracture tips and the
entire reservoir length is defined as the inner reservoir, and the rest is the outer reservoir, as illus-
trated in the figure below. The permeabilities of the inner and outer regions can differ, making this
model useful for modeling a stimulated reservoir volume (SRV), created by hydraulic fracturing,
which is fed by an unstimulated outer region.
This model calculates the reservoir's response from early-time storage and fracture flow, through the
transition into boundary-dominated flow. The fundamental building block of this model is the tri-linear
fracture model for a vertical well. The outer reservoir feeds the inner reservoir via linear flow, the
inner reservoir feeds the fractures via linear flow, and the fluid within the fractures travels linearly
towards the wellbore. However, for a horizontal well, the fluid within transverse vertical fractures
The reservoir dimensions, number of fractures, and fracture half-length may be specified by the
user, provided the entire wellbore and all fractures fit within the reservoir boundaries. The dimen-
sionless fracture conductivity must also be specified. Dual-porosity behavior can be modeled within
the inner-reservoir, but the outer-reservoir remains strictly homogeneous. A turbulence factor, D,
can also be specified.
Parameters
Pi = initial pressure
h = net pay
Φt = total porosity
Sg = gas saturation
So = oil saturation
Sw = water saturation
cf = formation compressibility
PSS = Transfer from secondary porosity to primary porosity is modeled as pseudo-steady state
(PSS).
Transient Slabs = Transfer from secondary porosity to primary porosity is modeled as transient,
using a slab shape factor.
Transient Cubes = Transfer from secondary porosity to primary porosity is modeled as transient,
using a cube shape factor.
Transient Sticks = Transfer from secondary porosity to primary porosity is modeled as transient,
using a stick shape factor.
ω = storativity ratio
λ = interporosity coefficient
Xe = reservoir length
Ye = reservoir width
rw = wellbore radius
References
1. “Comparison of Fractured Horizontal-Well Performance in Conventional and Uncon-
ventional Reservoirs", E. Ozkan, SPE, Colorado School of Mines, M. Brown, SPE, Colorado
School of Mines, R. Raghavan, SPE, Phillips Petroleum Co. (Retd.), and H. Kazemi, SPE,
Colorado School of Mines. Paper SPE 121290, Presented at the 2009 SPE Western
Regional Meeting held in San Jose, California, USA, 24-26 March 2009.
It is a rectangular model that contains an inactive horizontal well fed by multiple identical and
equally-spaced transverse fractures. X1 describes the extent of the improved permeability region.
The width of the reservoir is defined by the distance between the tips of the fractures, and the length
of the reservoir is defined by the horizontal well length.
Note: A detailed description of the model is given by Stalgorova and Mattar (2012).
Damage skin is applied along the length of the horizontal wellbore and a turbulence factor may be
specified.
h = net pay
Φt = total porosity
Sg = gas saturation
So = oil saturation
Sw = water saturation
cf = formation compressibility
PSS = Transfer from secondary porosity to primary porosity is modeled as pseudo-steady state
(PSS).
Transient Slabs = Transfer from secondary porosity to primary porosity is modeled as transient,
using a slab shape factor.
Transient Cubes = Transfer from secondary porosity to primary porosity is modeled as transient,
using a cube shape factor.
Transient Sticks = Transfer from secondary porosity to primary porosity is modeled as transient,
using a stick shape factor.
ω = storativity ratio
λ = interporosity coefficient
Xe = reservoir length
rw = wellbore radius
Because all fractures follow a repeating pattern, after you specify fracture locations, fracture half-
lengths, and fracture conductivities for one stage, the parameters of the rest of the stages are pop-
ulated automatically.
Prior to calculating the whole system’s solution, the model generates the solution for one single
stage through superposition of individual infinite conductivity fracture solutions in space. Then the
whole system’s solution is calculated by multiplying the single stage’s solution by the number of
stages. Consequently, the speed of calculation is improved and depends on the number of fractures
in each stage, rather than the total number of fractures.
In this model, (xf)y is always half of Ye. Changing the value of one dynamically changes the value of
the other.
Each of these configurations can be launched for single-phase (gas or water), two-phase (gas-oil or
gas-water), or three-phase (gas-oil-water) systems. Additionally, you can model gas condensate
and volatile oil systems with or without accounting for water.
1. Validation of Analytical models — Numerical models are used to ensure that complex
analytical models exhibit expected behaviours. It should be noted that analytical and numer-
ical solutions never yield exactly the same answer. However, for a successful validation
case, the difference ought to be less than 10% when comparing cumulative production over
a fixed period, or when comparing simulated pressure responses to a defined rate history.
(In some very complex cases, such as multi-frac’d shale gas completions with pressure-
dependent permeability and/or adsorbed gas, we have noted greater than a 10% difference
between the analytical and numerical solutions. Thus, the numerical models provide a valu-
able backup, in the event that the analytical models have been stretched beyond their cap-
abilities.) In general, it is a good idea to approach well performance analysis from several
independent perspectives. Agreement between analytical and numerical models usually
leads to the most robust interpretation possible.
Prerequisites
Before you start modeling, ensure that all required information is added to the WellTest project. Well
production data (e.g., rates and sandface pressures) should be populated in the Production Editor
Tab.
We recommend that you create and match an analytical model prior to numerical modeling. This
gives you an idea as to which reservoir and completion parameters to use, and makes it easier to his-
tory match.
To create a numerical model that accurately captures the reservoir behaviour, it is essential that the
fluid properties are populated correctly. See the section below for details.
Important: If the Gas Type is Liquid-Rich Gas, properties used for numerical modeling
are slightly different from the properties used for other types of modeling and analysis.
To be able to display and edit properties used for numerical modeling, click the Plot
Options icon ( ) on the toolbar and select the Display numerical analysis prop-
erties option in the Plot Options dialog box.
If you select Liquid-Rich Gas as the Gas Type, you have to set the Vaporized Oil Ratio (Rv).
If you have lab data available, you can set Rv as a table. Alternatively, you can use the cor-
relation suggested by Ovalle et al.
1. Estimate the gas-oil ratio (GOR) at the beginning of the production history (when the GOR is
fairly flat).
2. Enter your estimate as the Solution GOR (Rsol). (Pbp is calculated based on Rsol.)
By following this workflow, you ensure that the modeled GOR matches the production GOR
(at least at the beginning of the production).
Relative Permeability
For accurate numerical modeling of the multiphase system, relative permeabilities should be set.
By default, relative permeabilities use the Generalized-Corey correlation with preset default para-
meters. If you have relative permeability data (from lab studies, or based on an analog well), we
recommend that you adjust correlation parameters, or set permeabilities as a custom table.
Horizontal Slug
References
"A New Approach for Interpreting Pressure Data to Estimate Key Reservoir Parameters from
Closed-Chamber Tests", N. M. Anisur Rahman, SPE, Schlumberger, M.S. Santo, and L. Mattar,
SPE, Fekete Assoc. Inc., Paper SPE 109860, presented at 82nd Annual Technical Conference and
Exhibition of the SPE, Anaheim, CA, November 11 - 14, 2007.
With the composite model, you can add any number of different cylindrical zones. With an unlimited
number of zones, virtually any pressure transient response could be matched; therefore, it is import-
ant to exercise good judgement to determine when this model is appropriate. No reservoir is per-
fectly cylindrically concentric composite in nature; however, many reservoirs do behave the same
way as composite reservoirs do. Some common situations where the composite model is useful
include injection cases (which cause changes in viscosity and compressibility), reservoir het-
erogeneities (such as changes in flow capacity (kh)), and cases where the well was drilled into a nat-
urally fractured reservoir with varying fracture distribution.
Note that the following assumptions apply when using this model:
l Both the inner and outer regions can be assumed to be homogeneous or dual
porosity
References
1. "An Analytical Model for Composite Reservoirs Produced at Either Constant Bottomhole
Pressure or Constant Rate", Olarewaju, J.S. and Lee, W.J., SPE 16763 (1987).
At early times, this model uses the concept of Lee and Brockenbrough (1986) of tri-linear flow to rep-
resent a finite conductivity fracture (see figure below). Three linear-flow zones that dominate the
pressure behavior are:
Fracture diffusivity has been assumed constant at 1.0X106, as suggested by Cinco-Ley et al. (1978).
The tri-linear fracture flow results merge into the solution for infinite-acting radial flow in the middle
times. Thus, the tri-linear flow solution is truncated as soon as the flow becomes pseudo-radial.
Occasionally, the merging of these two solutions is not smooth, and the derivative exhibits dis-
continuities (spikes). These are localized aberrations and can be ignored as they do not affect the
rest of the results. Ultimately at late times, the model uses the solution for pseudo-steady state for a
no-flow outer boundary, or continues to use the solution for infinite-acting radial flow.
At very early times, the cylindrical source solution is used, which is followed by Green’s function solu-
tions. The Green’s function solution, as developed by Thompson et al. (1991) with slight modi-
fications, is used to simulate an infinite-conductivity vertical fracture. No-flow boundaries are
modeled using the method of images. The result is superposed in time based on the rate history
provided. Note that the dual porosity model provided is based on pseudo-steady state interporosity
flow.
This well may be at any location within the reservoir and the model supports infinite, no-flow, and
constant pressure boundaries. Thus, classical configurations like a well near a sealing fault, or a con-
stant pressure boundary near intersecting faults can be easily modeled.
At very early times, the cylindrical source solution is used, which is followed by Green’s function solu-
tions. No-flow boundaries are modeled using the method of images. The result is superposed in time
based on the rate history provided. Note that the dual porosity model provided is based on pseudo-
steady state interporosity flow.
The horizontal well is oriented in the x-direction and may be at any location within the reservoir (see
figure below), and the model supports infinite, no-flow, and constant pressure boundaries. Note that
the effective wellbore length (Le) defines the wellbore area open to fluid flow. Thus, classical con-
figurations like a well near a sealing fault, or a constant pressure boundary near intersecting faults
can be easily modeled.
At very early times, the cylindrical source solution is used, which is followed by Green’s function solu-
tions for horizontal wells, as developed by Thompson et al. (1991). No-flow boundaries are modeled
using the method of images. The result is superposed in time based on the rate history
provided. Note that the dual porosity model provided is based on pseudo-steady state interporosity
flow. The following flow regimes can be handled by this model:
l Wellbore storage
l Elliptical flow
l Boundary effects
The model accounts for transient flow within the fault as well. The formulation considers a line sink
well located in a two-zone reservoir, separated by a finite-conductivity fault. The width of the fault is
neglected for simplicity when developing the solution, except in the definition of the fault
conductivity. The fault is characterized by two parameters:
Based on the definitions mentioned above, the pressure drop along the fault length is due to the fault
conductivity (FCD), and the pressure drop across the fault is due to skin (sfault). Thus, various fluid
flow behaviour between the reservoir zones can be modeled by adjusting these parameters accord-
ingly. For example, fluid flow from a reservoir region far away from the wellbore to the other reservoir
region can be modeled using a highly conductive fault (high FCD), or a no-flow boundary could be
modeled by setting FCD to zero and sfault to a high value.
Note: Due to the nature of the mathematics involved in this model, the computational time is sub-
stantially more than in other models.
References
"Analytical Solution to the Transient-Flow Problems for a Well Located near a Finite-Conductivity
Fault in Composite Reservoirs", N.M. Anisur Rahman, M.D. Miller and L. Mattar, SPE Paper 84295
for Presentation in SPE Annual Technical Meeting and Exhibition, Denver, CO (October 2003).
Minifrac models yield results consistent with the Soliman / Craig after-closure analyses when clearly
defined flow regimes are developed. However, as with conventional pressure transient analysis, the
models provide more realistic estimates when clearly developed flow regimes have not developed,
which highlight misinterpretation of data in straight-line or derivative analysis.
Minifrac models should not be used to match the early-time data because the models do not account
for fracture closure, and are only designed for analyzing after-closure data. Changing wellbore stor-
age will be incorporated at a later date to assist in matching the early time data.
The injected volume in minifrac theory is intended to represent the total volume of fluid injected into
the reservoir. The value of Vinj, calculated in WellTest, is the volume injected from the Start of Injec-
tion line (located on the Data Chart on the Production Editor tab) to shut-in. In order for Vinj to equal
the total injected volume, the Start of Injection line would need to be placed at the point-in-time when
the wellbore had been pressured up to the initial reservoir pressure (see figure below). Any fluid
injected beyond this time would pressure up the wellbore, generate a fracture, and eventually leak-
off into the formation. If the Start of Injection line is placed at the point where the wellbore had been
pressured up to the initial reservoir pressure, a value for Vw is not required.
When the Start of Injection line is placed at breakdown, Vinj does not represent the total injected
volume. To account for the added energy of pressuring up the wellbore, a fluid expansion term is
added to Vinj, which requires an estimate of Vw. (VT is the value used in the model calculations.)
2. "A Study of the Behaviour of Bounded Reservoir Composed of Stratified Layers", H.C.
Lefkovits, P. Hazebroek, E.E. Allen and C.S. Matthews, SPEJ (March 1961) 43 - 58.
Each layer is considered to have a cylindrical geometry with an unequal initial pressure (pi) to other
layers, as well as its own skin factor, reservoir properties, and outer boundary condition. No
When the layers have unequal initial pressures, crossflow in the wellbore can occur from one layer to
another upon perforating the layers, or after shut-in as the pressures equalize. The following figure
illustrates the crossflow in the wellbore when the well is shut-in (pre-production period).
This shows that fluid is entering Layer 2 and Layer 3 from Layer 1 (the top layer) during crossflow
through the wellbore. The model calculates the sandface rate from each layer with time, where a pos-
itive rate means fluid is leaving the layer, and a negative rate means fluid is entering the layer. The
amount of crossflow at a given time is a function of many factors including the initial pressure (pi),
skin damage (sd), flow capacity (kh), storativity, boundary condition, and pore volume (Vp) of the
individual layers, as well as wellbore storage.
Hydrostatic pressure differences between layers can also be accounted for in this model by entering
a datum depth (Ddatum) and individual layer depths (Dlayer). All depths are referenced from the sur-
face.
References
1. "The New Analytical Solution to Pressure Transient Problems in Commingled, Layered
Zones with Unequal Initial Pressures Subject to Stepped Changes in Production Rates",
2. "A Study of the Behaviour of Bounded Reservoir Composed of Stratified Layers", H.C.
Lefkovits, P. Hazebroek, E.E. Allen and C.S. Matthews, SPEJ (March 1961) 43 - 58.
3. "Wells Producing Commingled Zones with Unequal Initial Pressures and Reservoir Prop-
erties", L. Larsen, Paper SPE 10325 for Presentation in AFTCE, San Antonio, TX, October 5
- 7, 1981.
4. "An Efficient Algorithm for Computation of Well Responses in Commingled Reservoirs", J.B.
Spath, E. Ozkan and R. Raghavan, SPEFE (June 1994) 115 - 121.
Each layer is considered to have a rectangular geometry with an identical initial pressure (pi) to other
layers, as well as its own skin factor, reservoir properties, and outer boundary conditions. This well
At very early times, the cylindrical source solution is used, which is followed by Green’s function solu-
tions. No-flow boundaries are modeled using the method of images. The result is superposed in time
based on the rate history provided. No crossflow between the layers can occur except at the well-
bore. The dual porosity flow within each layer can be modeled either as pseudo-steady state or tran-
sient interporosity flow.
References
1. "The Application of the Laplace Transformation to Flow Problems in Reservoirs", A.F. van
Everdingen and W. Hurst, Trans. AIME (1949) vol. 186, 305 - 324.
2. "The Use of Source and Green’s Functions in Solving Unsteady-Flow Problems in Reser-
voirs", A.C. Gringarten and H.J. Ramey, Jr., SPEJ (October 1973) 285 - 296.
3. "Transient Pressure Behaviour of Commingled Reservoirs", F.J. Kuchuk and D.J. Wilkinson,
SPEFE (1991).
At very early times, the cylindrical source solution is used, which is followed by Green’s function solu-
tions. No-flow boundaries are modeled using the method of images. The result is superposed in time
based on the rate history provided. Note that the dual porosity model provided is based on pseudo-
steady state interporosity flow.
References
"The Use of Source and Green’s Functions in Solving Unsteady-Flow Problems in Reservoirs", A.C.
Gringarten and H.J. Ramey, Jr., SPEJ (October 1973) 285 - 296.
At very early times, the cylindrical source solution is used, which is followed by Green’s function solu-
tions. No-flow boundaries are modeled using the method of images. The result is superposed in time
based on the rate history provided. Note that the dual porosity model provided is based on pseudo-
steady state interporosity flow.
A perforation inflow test is performed by shutting in the well at the surface and perforating wellbore
casing in an under-balanced condition (i.e., the cushion pressure (pwo) is significantly less than the
formation pressure). This creates a continuous pressure increase in the wellbore as fluid flows into
the wellbore at the sandface. As the wellbore fills with fluid, this increase in pressure slows down and
stabilizes near the initial pressure (pi) as shown in the plot below. This increase in pressure is meas-
ured over time, and the data obtained can by modeled to determine permeability (k), skin (s'), and ini-
tial pressure (pi) of the reservoir.
If a formation is perforated over-balanced (i.e., the cushion pressure (pwo) is greater than the form-
ation pressure), the wellbore pressure decreases as the wellbore fluid flows into the formation.
When the cushion pressure (pwo) exceeds the formation fracture pressure (during a fracture
Pressure data during the test can be obtained at the wellhead when inflow is expected to be single-
phase gas, and a fairly rapid pressure response is anticipated. The wellhead pressures can then be
converted to sandface conditions for analysis. When single-phase liquid inflow is expected, sand-
face pressures must be measured directly.
The late time analysis is used to estimate the initial pressure (pi), permeability (k), and skin (s') for
the model. These initial estimates can then be verified by the model. In cases where reservoir
(radial) flow behaviour is not fully developed within the test period, better estimates of these key
reservoir parameters are obtained directly from modeling.
The wellbore storage effect is accounted for by specifying the volume of the compressible fluid in the
wellbore (Vw for liquid or gas) or by specifying a changing fluid (liquid) level condition with a wellbore
volumetric capacity (Vu for liquid only). These values directly affect the results and must be estim-
ated as accurately as possible.
Note: Measured rates are not required for interpretation since the model calculates the fluid influx
rates using the pressure data only.
Because only the allowable fluid inflow into the wellbore is needed, the PITA test can be run over a
much shorter period of time than a conventional test. Thus, it is important to note that the reservoir
information obtained may only represent a small portion of the reservoir near the wellbore, espe-
cially when the permeability (k) is low, or the skin (s') is high. Under these conditions, if a transition
from wellbore storage to reservoir- dominated (radial) flow is sufficiently developed, meaningful
reservoir parameters can still be obtained. If reservoir-dominated (radial) flow is not sufficiently
developed, then a downhole shut-in is required to minimize wellbore storage and reduce the time
required to achieve reservoir-dominated flow. In this case, a step-change in wellbore storage results
when downhole shut-in occurs, and analysis / modeling is recommended to be undertaken using the
Closed Chamber Test (CCT) analysis or CCT model.
A test design PITA model can be used to determine the appropriate test duration, and the need for
down-hole shut-in. For a gas reservoir, down-hole shut-in can be simulated by simply reducing the
wellbore volume (Vw). Similarly, for a liquid-filled system, the wellbore volumetric capacity (Vu) can
be reduced to represent a small wellbore volume (Vw).
References
1. "Use of PITA for Estimating Key Reservoir Parameters", N. M. Anisur Rahman, Mehran
Pooladi-Darvish, Martin S. Santo and Louis Mattar, Paper CIPC 2006 - 172, presented at
2. "Development of Equations and Procedure for Perforation Inflow Test Analysis (PITA)", N.
M. Anisur Rahman, Mehran Pooladi-Darvish and Louis Mattar, Paper SPE 95510, presen-
ted at 80th Annual Technical Conference and Exhibition of the SPE, Dallas, TX, October 9 -
12, 2005.
This model supports infinite, no-flow, and constant pressure boundaries. Note that the effective well-
bore length (Le) defines the wellbore area open to fluid flow. Thus, classical configurations like a well
near a sealing fault, or a constant pressure boundary near intersecting faults can be easily modeled.
At very early times, the cylindrical source solution is used, which is followed by Green’s function solu-
tions for horizontal wells as developed by Thompson et al. (1991). No-flow boundaries are modeled
using the method of images. The result is superposed in time based on the rate history provided.
Note that the dual porosity model provided is based on pseudo-steady state interporosity flow. The
following flow regimes can be handled by this model:
l Wellbore storage
l Elliptical flow
l Boundary effects
Due to the additional complexity of a slanted wellbore, the solution method requires that the wellbore
be split up into specified number of segments (Nw). Increasing the number of segments will result in
higher precision in the solution, but with slower computational performance. A value of 100 seg-
ments is recommended as it handles most situations adequately.
References
"Efficient Algorithms for Computing the Bounded Reservoir Horizontal Well Pressure Responses",
L.G. Thompson, J.L. Manrique and T.A. Jelmert, Paper SPE 21827 presented at 1991 Rocky Moun-
tain Regional Meeting and Low-Permeability Reservoirs Symposium, Denver, CO, April 15 - 17.
The slug test is performed by shutting in the well at surface and introducing an instantaneous pres-
sure drop across the sandface either by perforating the wellbore casing, or opening a downhole
valve. This creates a continuous pressure increase in the wellbore as fluid flows into the wellbore at
sandface. As the wellbore fills with fluid, this increase in pressure slows down and stabilizes near the
initial pressure (pi) as shown in the figure below. This increase in pressure is measured over time,
and the data obtained can be modeled to determine permeability (k), skin (s'), and initial pressure
(pi) of the reservoir.
Because only the allowable fluid inflow into the wellbore is needed, the slug test can be run over a
much shorter period of time than a conventional test. Thus, it is important to note that the reservoir
The impulse radial analysis is used to estimate the initial pressure (pi) for the model. Permeability (k)
and skin (s') cannot be determined uniquely from the impulse radial analysis; therefore an initial
guess for these values is required.
The wellbore storage effect is accounted for by specifying the volume of the compressible fluid in the
wellbore (Vw for liquid or gas), or by specifying a changing fluid (liquid) level condition with a well-
bore volumetric capacity (Vu for liquid only).
Note: Measured rates are not required for interpretation since the model calculates the fluid influx
rates using the pressure data only.
References
1. "Analysis of Slug Test or DST Flow Period Data", H.J. Ramey, Jr., R.G. Agarwal and I.
Martin, JCPT (July - September 1975) 37 - 47.
2. "Annulus Unloading Rates as Influenced by Wellbore Storage and Skin Effect", H.J. Ramey,
Jr. and R.G. Agarwal, SPEJ (October 1972) 453 - 462.
3. "Analysis of Slug and Drillstem Tests", A.M.M. Peres, Ph.D. Thesis, The University of Tulsa
(1989).
With this model, you can observe the pressure response at an observation point (Zo) along the well-
bore to determine the horizontal permeability (kh or kxy) and vertical permeability (kv or kyz).
The anisotropy is handled using a conformal mapping procedure that adjusts the boundary sizes
accordingly to mimic the effect of increased or decreased permeability in each direction. This well
may be at any location within the reservoir and the model supports infinite, no-flow, and constant
pressure boundaries. Thus, classical configurations like a well near a sealing fault, or a constant
pressure boundary near intersecting faults can be easily modeled.
At very early times, the cylindrical source solution is used, which is followed by Green’s function solu-
tions. No-flow boundaries are modeled using the method of images. The result is superposed in time
based on the rate history provided. Note that the dual porosity model provided is based on pseudo-
steady state interporosity flow.
At early times, the wedge boundaries are modeled using the method proposed by Chen and
Raghavan (1997), while late time is modeled using the infinite-acting radial solution. The result is
superposed in time based on the rate history provided. Note that the dual porosity model provided is
based on pseudo-steady state interporosity flow.
References
1. "The Application of the Laplace Transformation to Flow Problems in Reservoirs", A.F. van
Everdingen and W. Hurst, Trans. AIME (1949) vol. 186, 305 - 324.
2. "Computing Pressure Distribution in Wedges," C.-C. Chen and R. Raghavan, SPEJ (March
1997) 24 - 32.
Data Tables
Tabs
Wizard Menu
This menu includes standard grid options, such as copy-and-paste, as well as options where you
can set precision, find a value within the column, or manipulate data. Data manipulation is applied to
any selected cells, which gives you the flexibility to select the entire column, or a portion of the data.
l Print and Preview data table — You can print the table. There is currently no
print filter, and this print operation can easily print hundreds of pages for sand-
face pressure recorders.
l Create / Delete Calculated Dataset — This option is only available in the tool-
bar for the All Data tab. Datasets can be calculated based on the difference or
average of similar data types from two different gauges. For example, a pressure
comparison dataset can be generated by calculating the difference between the
pressure values from tandem sandface gauges. The primary pressure derivative
of pressure datasets can also be calculated.
l Export gauge data — This option is not available for the All Data tab. Data can
be exported to a csv format.
l Find next / previous shut-in — This option is not available for the All Data tab
and is only active in the toolbar when gauges have rate information. It is often
necessary to navigate to shut-in points to ensure data is valid and has been syn-
chronized properly. With this option, you can efficiently navigate to the next shut-
in point.
l Zero gauge start time — This option is not available for the All Data tab and is
only active in the toolbar when the first row of a gauge has a non-zero cumulative
time. This operation changes the gauge start time to the date in the first row of
the data table for that gauge. All the cumulative times update to reflect the
change in gauge start time.
Data Management
Production Editor
Properties
AER PAS
Reports
Comparison Plot
Data Management
Deconvolution
Advanced Models
There are five sub-tabs within the Data Management tab, which are described below.
Data Sub-tab
The Data sub-tab stores a copy of the imported gauges, as well as any datasets created from mer-
ging multiple gauges, or smoothing data within a gauge. All imported data is automatically displayed
in a table with a corresponding plot. The data tables and plots work together to display information
efficiently. When you click a point on a plot, the plot's data marker moves to the nearest data point,
and the corresponding row in the data table is highlighted in blue. Likewise, when you select a row in
the data table, the marker on the plot moves to the data point, which corresponds to that row. Both
the data tables and plots have tools to help navigate to specific points, or to manipulate data.
The plots within the Data sub-tab enable you view and manipulate data.
l Switch Active Dataset — This option is active on plots displaying more than one
curve. It controls which dataset the zoomed plot displays and the plots marker it is
attached to.
l Select / Delete range — Range markers may also be used to select data. The
Delete option enables you to remove all points within or outside a selected range.
l Delete selected points — The Delete option has a sub-toolbar with options to
remove data selected on the active dataset, on all datasets within the gauge, or
inside / outside a selected range.
The top left pane contains tabs of data tables for each imported gauge, and an all data table of all
imported data.
The bottom pane contains tabs of plots for each imported gauge and an all data plot.
l
To add new plots, click the New Plots icon ( ) on the toolbar.
l
To calculate datasets, click the Create calculated dataset ( ) icon. Note that it is
common to create a gauge comparison plot if tandem recorders are run downhill.
The top right pane displays the active dataset in the Zoomed Data Chart.
Synchronize Sub-tab
In the Synchronize sub-tab, you can synchronize data from different gauges, or datasets contained
within a single gauge.
Note: To open this sub-tab, at least two valid gauges, or two datasets within one gauge must exist.
The top left pane displays the primary plot, and the data in this plot stays where it is.
The bottom left pane displays the secondary dataset. The data in this plot will shift when the Apply
The bottom right pane displays the results plot, which illustrates what the data would look like when
the Apply icon on the primary dataset is clicked.
Merge Sub-tab
Data from different gauges (i.e., imported from different files) can be merged into one continuous file
by clicking the Merge sub-tab.
Note: To open this sub-tab, at least two valid gauges must exist.
In the Merge tab, you can select the datasets to be merged and view the result of the merge on the
plot.
The Data Table / Grid, Data Chart, and Zoomed Data Chart are all linked. If you click a row in the
table, the annotation arrow moves to the corresponding point on the data chart. Similarly, if you click
a point in the data chart, the corresponding row on the table is highlighted.
Note: The data table is also known as a grid, so the Grid Options icon corresponds to the data in
the table.
Note: The notation WellTest uses for rate data is "end of period". For additional information, see
End of Period - Rate Data.
The Total Gas Rate column is not editable and is used for analysis of gas wells. This column is for
the recombined gas rate for condensate wells, and is equal to the gas rate for wells with no con-
densate. When the Use RGRF checkbox in the header of the Cond. Rate column is checked, the
rates in the Gas Rate column are multiplied by the recombined gas rate factor (RGRF) in the Recom-
bination section of the Properties / Gas Properties tab. This option should be used when the con-
densate rates aren’t available, or when the produced condensate gas ratio (CGR) is sporadic due to
challenges with measurement and reporting. When this option is unchecked, a different recombined
gas rate factor is calculated for each row. The inputs to this calculation are the gas and condensate
rates, used to obtain a CGR, the separator pressure and temperature, and condensate gravity in the
Recombination section of the Gas Properties tab.
The oil and water gradient columns are used for calculating sandface pressures for the hydrostatic
head of the fluid below a specified fluid level. These gradients may be calculated using the PVT prop-
erties when the Calc checkbox in the header of these columns is selected, or entered manually,
when the Calc checkbox is unchecked.
Data Chart
The pressure and rate data is plotted on the Data Chart. To focus on the test data, the Data Chart
autoscales to display all the pressures. However, you may wish to customize the scale to view the
previous production.
Note: To find a shut-in point quickly, right-click the Data Chart within the first few buildup points,
and then press the left arrow key to move the arrow annotation back a few points. This is
much more efficient than trying to click on the final flow point. This is because the data dens-
ity is much lower at the start of the buildup because the pressures are changing quickly with
time.
Grid Options
To hide or show different columns on the data table, click the Grid Options icon ( )
on the toolbar and select your desired options from the drop-down list:
l Hide Show Columns – can also be accessed by right-clicking anywhere on the table
and selecting Show Columns. When you select this option, the Hide/Show Columns
dialog box opens, where you can click the checkboxes for the columns you want to
show (or hide) on the table.
l Show Sandface Pressure Columns – shows the time, rate, and sandface pressure
columns, and hides the wellhead pressure, liquid level, and gradient columns.
l Set Drag Scroll Speed – sets the speed at which the mouse wheel scrolls through
the data table / grid.
l Show Errors / Warnings – shows Production Editor errors / warnings, such as time
out of sequence errors, if there are any.
Data Selection
This icon launches a dialog box where you can select data from the Data Management tab, and
bring it into the Production Editor for analysis. Pressures and rates may be brought it from different
gauges, and are automatically merged into the Production Editor.
Previous Production
This icon launches the WellTest Wizard: Previous Production dialog box, which prompts you to
either enter previous production data manually, or specify an effective producing time. Note that this
icon is only available for test types that require rates.
Undo / Redo
These are standard Undo / Redo operations for any changes made to the Production Editor dataset.
Clicking the Filter icon ( ) opens the Filter data dialog box, where you have the option of
conducting a primary data filter. The primary data filter is a standard step in preparing data for ana-
lysis. For additional information, see Filtering.
The data is segmented into flow and shut-in periods and different filtering parameters may be
applied to each period. It is very important that the flow and shut-in periods have been correctly spe-
cified before filtering. With the checkboxes you can select the filter mode (e.g., logarithmic, arith-
metic, or no filter) for each period. The values in the Logarithmic Points Per Cycle and Arithmetic
Time Step columns are editable. The Current and Remaining Data Points columns keep track of how
many points there originally were, and how many points remain after the filter is applied. The final
row shows the total number of points before and after filtering.
Logarithmic filtering keeps more data at the start of the flow period, when the pressure is changing a
lot with time, and less data at the end, when the change in pressure is more gradual. (This is the
default for all shut-in periods.) Arithmetic filtering takes evenly spaced points (e.g., 1 point per hour)
and is the default for all flowing (or injecting) periods.
Determine Shut-in
Clicking this icon opens the Determine Shut-in window where you can insert a time, pressure, and
zero-rate point into the middle of your dataset. This is useful when the final flow point is not clearly
defined.
In the Properties tab, you can enter and view reservoir and fluid parameters. Within the tab there is a
Property Type View, a Reservoir Properties View, an Inputs Pane, and a Main View.
An exclamation mark beside a property type indicates that the properties have not been fully spe-
cified. In the screenshot below, Oil has an exclamation mark beside it. This is only a problem if your
reservoir has oil. If this is the case, click Oil and enter the required parameters.
l
View Data — view your data in the Main view in a table format. This icon is only
displayed when the Main view is displaying plots.
l
View Plot — view your data in the Main view represented by plots. This icon is
only displayed when the Main view is displaying tables.
l
Import Data — bring in your custom property data.
l
Help for Correlations — after clicking this icon, select either PVT Correlation
Help, or Viscosity Correlation Help. A dialog box opens, which displays inform-
ation about each of the available correlations. This icon is disabled for property types
that don't have correlations.
l
Change Bottom Left Plot — while viewing your plots, there may be an addi-
tional plot that can be viewed by clicking this icon. This icon is only displayed when
more than four plots are available. For example, when liquid-rich gas is selected,
there is an additional plot to view the vaporized oil ratio.
Gas Properties
This view displays drop-down lists where you can select the PVT correlation (for Z, Bg, and cg) and
viscosity correlation.
Gas Properties, Gas Composition, and Critical Properties can be entered as input parameters in
each of their sections. The calculated properties at pi and TR are also displayed.
Gas properties vs pressure are plotted, including Z, ψ, Bg or Bgd, ρg, cg, Rv, VL / Vsat, and μg.
l pdew — dew point pressure. It should be set to the same value as the bubble point
pressure defined under Oil / Condensate properties.
l γcond — condensate gravity. It should be set to the same value as the oil gravity
defined under Oil / Condensate properties.
l RGRF — recombined gas rate factor. We recommend keeping the default value.
l Gr — recombined gas gravity. All gas property functions are generated based on this
gas gravity (as well as critical temperature and pressure).
After you select Oil / Condensate in the selection tree, you can select a PVT correlation, constant
properties, or custom table from the PVT Correlation drop-down list. In the Viscosity Correlation
drop-down list, you can select a correlation, constant properties, or custom table.
Oil properties can be entered as input parameters. The calculated properties at pi and TR are also
displayed. Calibration properties can be entered to calibrate your oil / condensate properties.
Oil properties vs pressure are plotted, including Rso, Bo, ρo, co, and μo.
Water Properties
This view displays a drop-down list and fields where you can enter your water properties.
After you select Water in the selection tree, you can select a correlation, constant properties, or cus-
tom table from the General Correlation drop-down list.
You can enter Water Composition as input parameters, and you can select the Gas Saturated check-
box to take into account gas saturation in water. The calculated properties at pi and TR are dis-
played.
Water properties vs pressure are plotted, including Rsw, Bw, ρw, cw, and μw.
CBM Properties
WellTest supports coalbed methane (CBM) analytical models (appropriate for dry CBM wells) and
numerical models (for wet or dry CBM wells). The pseudo-steady state (PSS) models, designed to
Enter the Langmuir Isotherm and Coal properties to account for the gas adsorbed on the coal.
CBM properties vs pressure are plotted, including Gas Content, k / ki, and Φ / Φi.
Adsorption Properties
Adsorption properties should be entered to account for adsorbed gas in shales. For adsorption in
CBM reservoirs, the CBM properties should be populated.
After you select Adsorption in the selection tree, Langmuir Isotherm, and Shale Properties can be
entered as input parameters.
After you select Geomechanical in the selection tree, you can select a correlation or custom table
from the Permeability Ratio Correlation drop-down list. In the Compressibility Ratio Correlation drop-
down list, you can select a correlation or custom table.
After you select Capillary Pressure in the selection tree, you can select a correlation or custom
table from the Capillary Pressure Correlation drop-down list.
Capillary pressures vs pressure are plotted, including pcgo, pcgw and pcow.
Relative permeabilities vs pressure are plotted for gas-water, gas-oil, and oil-water. The 3-phase rel-
ative permeability is also presented in a triangle graph.
Main View
This view consists of plots or tables displaying how properties vary with pressure or saturation, and
depends on the property type selected in the Property Type view. To switch between viewing plots
or tables, click the View Data or View Plot icons on the toolbar of the Inputs pane.
Note: The tables are read-only unless a Custom correlation is selected in the Inputs pane.
Note: If you do not deal with wells in Alberta, Canada you can hide the PAS tab.
This tab houses a number of sub-tabs that facilitate the entry of all well, test, and gauge data
required to generate the PAS file.
1. The left section displays a tree structure of available reports generated for the current file.
The reports displayed on the left are organized into Default, Custom and File Specific folders.
l Custom — contains reports generated from templates that were created or cus-
tomized in the Design Report tab.
l File Specific — contains reports not generated from templates, but rather cre-
ated or modified from existing default or custom reports. These reports are saved
to the file, and are not available to other files.
Ico-
Name Description
n
Open You can navigate to and open an existing report template (.PRT exten-
Exist- sion).
ing If the report template being opened already exists in the application, it is
Report automatically selected in the preview tree. If the report template being
opened is not associated with the current application, then the template is
displayed in the temporary folder.
Chang- Launches the Report Paths window, where you can change the location
e the where custom report templates are saved. By default, the report templates
report are located in:
paths
In a network installation, this path may not be found on the C drive. In that
case, the IT department may be required to help locate this folder. After
the location of the report templates are determined, the path can be set by
clicking the Browse icon. When this setting has been changed, WellTest
must be restarted for the changes to take effect.
Note that the Reset button can be used to clear all saved paths, and reset
the save path to the default setting.
Chang- Opens the Report Settings window, where you can change whether or not
e the filename, file path, or date is displayed in the report footer.
Report
Set-
tings
Reload Refreshes the view of the tree to reflect changes to the preview tree that
report- may not have been displayed yet. In most cases, this is automatically
s from done.
hard-
drive
Restor- You can recover deleted reports from the current session of the application
e by opening the Retrieve Deleted Templates dialog box.
delete- Note: If a deleted report template is not recovered before the application is
d closed, it is permanently deleted. A warning is displayed when you exit the
report application to inform you if there are templates available to be recovered.
tem- Also, report templates in the default folder cannot be deleted.
plate
(s)
Collate You can create groups of templates, which gives you the option of printing
Report- multiple templates at once.
s Clicking this option displays another tree next to the preview window.
Collated Reports
With a collated report, you can group several reports together, so you can print several reports at
once. There are several options available when working with collated reports.
To display the current list of collated reports, click the Collate Reports icon ( ).
To create a new report folder, right-click the empty space within the Collated Reports window, and
then select Create New Report.
You can add reports from the preview tree to the Collated Reports folder. To do this, click the report
name, then drag it to the Collated Reports folder. To select multiple templates, drag a selection box
around the template names.
1. The left section contains all of the tools required to modify and customize reports and tem-
plates.
l Reports are generated from templates, or by modifying existing reports in the Design
Report tab. For example, there is a pre-defined default template for a fracture model
for a gas reservoir. If a specific file has the reservoir fluid type set to gas, and has
three fracture models, three different reports are generated from that single report
template. The values in each report correspond to the values in the model associated
to it. Reports are file-specific. When reports are modified in the Design Report tab,
the saved report file (with a .prt extension) is only available for that specific WellTest
file.
l Templates refer to the report layout structure for a specific fluid type, and a specific
analysis or model. Templates are used to generate reports for any file with a match-
ing fluid type, and analysis or model type. In fact, if a file has multiple instances of the
same analysis or model, the corresponding template is used to generate multiple
reports. Templates are not file-specific. They are saved (with an .lof extension) to a
specified location on the hard drive, or network. Custom templates may be created
from scratch, or by modifying existing templates.
Note: You are likely to modify existing reports, and customize existing templates, rather than cre-
ating them from scratch.
In the left section of the Design Reports tab, is the Existing Reports and Templates list where you
may select the report or template to start from. After you have selected a report or template, a drop-
down list displays all the reports or templates available for modification. If you start from a report (a
file in the Default Reports, Custom Reports, File Specific Reports, or All Reports section), the mod-
ified file is saved as a report, and is displayed in the File-Specific section of the Preview Report tab,
and is not available to other project files. If you start from a template (a file in the Default Reports,
Custom Reports, or All Templates sections), the modified file is saved as a custom template and is
used to create custom reports for similar WellTest files.
Toolbar Icons
The Design Report toolbar contains the following icons:
Rearrange Displays all the pages in one view, so they can be rearranged by
Pages dragging- and-dropping pages.
Preview Tem- Switches the view from the Design Report tab to the Preview
plate Report tab.
Save Template Saves the currently active report template.
Add/Edit Tem- You can change the default location for saved files.
plate Locations
Hotkeys
When you are working in the report, the following hotkeys are available:
l Control D — Displays the order of objects on a report page. The order of reports
can be changed by clicking the objects in a specific order.
l Control R — Saves the relative paths for all image files. (By default, absolute
paths are used.)
In a network installation, this path may not be found on the C drive. In that case, the IT department
may be required to help locate this folder. After the location of the report templates are determined,
If the software is installed on a location on a network with restricted user rights, the first time you cre-
ate a report and try to save it, you are notified that you do not have access to the default location,
and you are asked to provide an alternate location. However, you still have access to the original
default location.
To change the location where files are saved, click the Browse... button.
By default, Fekete reports are included with the software. These reports are displayed under the
Default folder in the preview tree. Dynamic reports (reports that are automatically created) reside in
a temporary folder, and are not affected by changing the default folder's location.
Custom Templates
The following section describes the Existing Reports and Templates section of the Design Report
tab.
l Modified Reports — displays reports that have been modified from Default
Reports, Custom Reports, and Modified Reports. The modified reports can be
accessed by clicking Modified Reports in the Existing Reports and Templates
section. It is displayed in the File Specific folder of the Preview Report tab. These
reports are saved as part of the current file and do not get displayed in any other
files.
l All Reports — displays the default, custom and modified reports available to the
current file. All reports can be accessed by clicking All Reports in the Existing
Reports and Templates section.
l All Templates — displays the default and custom templates. They can be
accessed by clicking All Templates in the Existing Reports and Templates sec-
tion.
The easiest way to get data into the comparison plot is to copy-and- paste a plot into the comparison
plot.
2. Click the Comparison Plot tab; right-click below the tab and select Paste Plot.
Now the datasets from both plots are displayed on the same comparison plot. If two data-
sets with the same color and symbol are pasted onto the same plot, the symbols for the
second dataset will automatically be changed, so that the curves may be easily dis-
tinguished from one another.
3. (Optional) Modifications to dataset names, symbols, and line types can be made by right-
clicking the comparison plot and selecting Customize.
Deconvolution Plots
Deconvolution plots are described below.
Deconvolved Typecurve
This plot displays the deconvolved typecurve (unit rate pressure response) for the reservoir. Note
that until the deconvolution is run, this plot only shows the initial typecurve specified in the task bar,
and does not show the final result.
Data Typecurve
This plot displays the match between the deconvolved and selected typecurve data. Additional flow
periods can be displayed (overlayed) on this plot in one of two ways:
Total Test
This plot displays the pressure and rate history. Pressure points can be selected (light blue) or
deselected (white) for fitting on this plot.
Task Bar
The task bar provides an easy way of displaying flow periods, selecting points, and selecting the ini-
tial TypeCurve (TC). Each row corresponds to either a DrawDown (DD) or a BuildUp (BU) period.
The task bar can be docked, undocked, or hidden as needed.
l Show TC — Toggle these checkboxes to show / hide any flow period on the
Data Typecurve plot. (TC = Type Curve)
l Select Periods — Toggle these checkboxes to select / deselect any flow period
pressure data for fitting.
l Initial TC — Click a radio button to select which flow period typecurve to use as
the initial Deconvolved Typecurve. By default, the final buildup is selected. The
initial typecurve can also be manually adjusted by clicking-and-dragging any
derivative point (red) on the Deconvolved Typecurve plot. Alternatively, you can
use the Linear Point Alignment option from the toolbar on the Deconvolved Type-
curve plot by clicking a derivative point (red) and drawing a trend line that the
derivative points should follow. When the left mouse button is released the deriv-
ative points is aligned along this trend line.
Note: After any change to the derivative, the Pu curve is updated accordingly.
Note: For more information on deconvolution parameters and how they affect the fit, see Decon-
volution Theory.
Toolbar
1. pi — Use the specified pi (no anchor); note that the specified pi can be entered in the field, or
set by clicking-and-dragging the calculated pressure curve on the Total Test plot up or
down. Selecting the checkbox for this parameter includes it during the fitting process. We
recommend that when pi is not known, that the other anchoring options be used to minimize
instability during the fitting process.
2. pwfo — Anchor to the final flowing pressure before the last buildup is selected; pi is back-cal-
culated so the calculated pressure curve on the Total Test plot matches this point.
3. Any Point — Anchor on any pressure point. After selecting this option, click the pressure
point on the Total Test plot to anchor to. The A # symbol is displayed over the pressure point
you are anchored to when this option is on. To change the anchor point at a later time, click
the anchor icon on the Total Test plot toolbar and click another pressure point.
Typecurve Options
These parameters directly affect the Deconvolved Typecurve, as well as the calculated pressure:
l Pu1 — Sets the value of the first point of the pu curve on the Deconvolved Type-
curve. Selecting the checkbox for this parameter includes it in the fitting process.
Note that the fitting process may take longer to run when this parameter is set to
auto. We do not recommend fitting on this parameter and pi at the same time, so
as to avoid instability in the fitting process.
l Wpres — Sets the weighting factor that the pressure error (Epres) is multiplied
by during the fit.
l Wcurv — Sets the weighting factor that the curvature error (Ecurv) is multiplied
by during the fit.
l Wrate — Sets the weighting factor that the rate error (Erate) is multiplied by dur-
ing the fit. Note that this option is only visible when the Adjust Rates option is
selected.
l Adjust Rates — Select this checkbox to include rate corrections in the fit. If this
checkbox is not selected, the Wrate and Erate parameters are hidden.
Deconvolution Results
These fields display the following calculated errors during the fit.
l Epres — The overall error between the measured and calculated pressures.
l Ecurv — The curvature error. Represents the amount of curvature on the deriv-
ative of the Deconvolved Typecurve. The more linear the derivative curve, the
smaller this value will be.
l Erate — The overall error between the measured and calculated rates.
l ETLS — The total least squares error, which is a summation of all the errors.
Advanced Options
Clicking the Advanced Options button displays the following options:
l Minimization Method — Sets the method to use for fitting (minimization). The
default recommended method is More-Hebdon. Other methods (Line Search and
Double Dog Leg) can be used in this order, if the current method fails to find a
solution.
1. Click the Filter The Data icon on the Control Panel toolbar to open the Filter Well Data dia-
log box (Data Reduction Module). Modify the filtering options as needed to filter the data
accordingly. Any filter settings that are changed are saved and used again in subsequent fil-
ter operations.
Note: Each time you use this filter option, all the flow periods (Select Periods column on the
task bar) are re-selected in order to include them in the filtering process.
2. Select the Filter Data menu option from the popup menu on the Total Test plot. The filter is
applied only to the flow period the mouse hovered over when right-clicking. Note that this fil-
ter option does not remove any points from the dataset, but rather deselects the pressure
points based on the filter, so they are not included when fitting.
Model Manager
Toolbar
The left-side toolbar is located just below the analytical model tab and provides the following options
(described from left-to-right):
l CalcBoth — Calculates a synthetic pressure using measured rate as the input, and
calculates a synthetic rate using measured pressure as an input.
l Derivative Options — Launches a dialog box where you can specify derivative cal-
culation options.
l Copy to / Paste from Clipboard -- "Copy to Clipboard" copies the current values of
model parameters to the clipboard. "Paste from Clipboard" pastes the model para-
meters that have been copied to the clipboard. If model parameters have not been
copied to the clipboard, the Paste from Clipboard option is not available.
Another toolbar is located on the far right. With the first four icons, you can replace one of
the existing plots with a different plot, and with the fifth icon, you can add a new floating view.
When any of these icon are clicked, a menu with available plot types is displayed. There is
also a checkbox for Corrected Pseudo Time.
l The History plot is a Cartesian plot of the measured and synthetic data. It is used to
select the flow and shut-in periods that are displayed on the Specialized and Deriv-
ative plots. Click the bottom section of the plot to select or deselect a region. Multiple
flow and shut-in periods may be selected at the same time.
l The Derivative plot shows the data and the Der (semilog derivative) for all periods
selected on the History plot.
l The Buildup / Falloff plots are the traditional plots used to linearize the data for a
specific flow regime. They only display data from shut-in periods.
This tab has two sections: a Results plot and a Forecast table.
Forecasting
Forecasting in the Advanced Models takes the production history into account and does not com-
mence assuming static reservoir conditions. A forecast may be generated based on your reservoir
model parameters and specified forecasting parameters.
Forecast Parameters
Several additional features exist in the Forecast pane that can be used to customize the forecast.
This pane consists of the following three sections:
1. Forecast Options
2. Forecast Constraints
3. Forecast Results
Note: Each of these sections can be expanded or collapsed by clicking the +/- box.
Forecast Options
The Forecast Options section is where forecast periods are defined according to either durations
(e.g., months), or dates.
l Setting the forecast time method to Duration begins a forecast at a specified start
date. Forecast periods are then defined in the forecast table according to a length of
time (the default is months). For each forecast period, different operating conditions
can be specified.
l Setting the forecast time method to Start Date creates a forecast of a specified length
of time (the default is months). Forecast periods are then defined in the forecast table
according to calendar dates.
If you select Start Date from the Forecast Time Method drop-down menu, enter a value for Duration
(i.e., the total length of the forecast).
The Forecast Flowing Pressure can also be set in the Forecast Options section. If a Gas Condens-
ate well is being analyzed, the Separate Recombined Gas option is available in Forecast Options.
Further details about each forecast period are entered in the forecast options table.
l Time – defines the length of each forecast period and each timestep, and has the fol-
lowing columns:
l Control – defines how the operating conditions change over the forecast period, and
has the following columns:
l Sandface Pressure – If the Control Type is set to Pressure, this column is displayed,
and it is used to set the flowing pressure for the forecast period. If Interpolation is set
to Step, only the initial pressure cell is editable. If Interpolation is set to Ramp, initial
pressure and final pressure are selectable.
l Gas/Oil/Water Rate – If the Control Type is set to Gas, Oil, or Water, this column is
displayed, and it is used to set the rate for the forecast period. If Interpolation is set to
Forecast Constraints
The Forecast Constraints section is where maximum rate conditions and abandonment rate con-
ditions for the forecast can be entered. Available constraints depend on the analysis type, and could
include the following:
l pmin – The minimum allowable sandface flowing pressure during the forecast. In
order to ensure flowing pressure does not go below this value, rates are adjusted.
l (qg)max – Sets a maximum gas rate during the forecast. In order to maintain the max-
imum gas rate constraint, flowing pressure is adjusted.
l (qw)max – Sets a maximum water rate during the forecast. In order to maintain the
maximum water rate constraint, flowing pressure is adjusted.
l (qo)max – Sets a maximum oil rate during the forecast. In order to maintain the max-
imum oil rate constraint, flowing pressure is adjusted.
l (qg)ab – Sets the abandonment gas rate for the forecast. When the abandonment
rate is reached, the forecast ends.
l (qw)ab – Sets the abandonment water rate for the forecast. When the abandonment
rate is reached, the forecast ends.
l (qo)ab – Sets the abandonment oil rate for the forecast. When the abandonment rate
is reached, the forecast ends.
For analytical models, the available constraints include a minimum flowing pressure, as well as max-
imum and abandonment rates corresponding to the fluid being analyzed. For example, a gas ana-
lytical model includes pmin, (qg)max, and (qg)ab. Numerical models include minimum flowing
pressure, as well as maximum and abandonment rates for all fluids in the model.
Forecast Results
The Forecast Results section is where the results of the forecast are summarized. Available results
depend on the analysis type, and could include the following:
Note: For wells with historical data, EUR is defined as the cumulative of the historical data up to
the beginning of the forecast, after which time, the synthetic cumulative calculated from the
model is used until the end of the forecast.
For typecurve and analytical models, EUR and RR of the fluid being analyzed are available. When
analyzing a gas condensate system, EUR and RR of the condensate phase are also available. For
numerical models, EUR and RR of all the fluids in the model are available.
Tables Tab
The Tables tab is a tabular display of the data within each plot in the Plots tab.
Toolbar
The left-side toolbar is located just below the numerical model tab and provides the following options
(described from left-to-right):
l Stop Synthesize — Stops the numerical modeling simulation. This icon is only
enabled during a simulation.
l Copy to / Paste from Clipboard -— "Copy to Clipboard" copies the current values
of model parameters to the clipboard. "Paste from Clipboard" pastes the model para-
meters that have been copied to the clipboard. If model parameters have not been
copied to the clipboard, the Paste from Clipboard option is not available.
Reservoir Tab
The Reservoir tab is where all modeling parameters are entered. The available model parameters
depend on the model type and fluid type, which were selected in the Model Manager when creating
the numerical model. To update the history match and forecast results, update the parameters and
click the Synthesize icon on the toolbar.
l Uniform Grid: Set this option to divide the grid blocks evenly, so that each grid block
is the same size. This results in faster performance, but may not accurately model
transient regions as well.
n #max — Sets the maximum number of divisions to divide the reservoir into.
The fewer the number of divisions, the faster the performance.
l Geometric Grid (Default): Set this option to divide the grid blocks geometrically, so
that there are more grid blocks near the well, and fewer away from the well. This res-
ults in slower performance, but models transient regions more accurately.
n rb — Sets the geometric block ratio, which controls the spacing of the blocks
from the wellbore to the outer portion of the reservoir.
l Fine Wellbore Grid: Turns fine wellbore gridding on / off to further divide the blocks
near the wellbore. This option results in slower performance, but models wellbore
effects such as skin and wellbore storage (afterflow) more accurately.
l Show Visual Output: Enable / disable the gradient results during simulation. If gradi-
ent results are not important, disable this option to increase performance.
Plots Tab
This view is split into four sections, each displaying a schematic, history data, pressure profiles, or
plan views. The item displayed in each section can be changed using the icons on the right-side tool-
bar.
l The schematic is an overhead view of the reservoir with dimensions for size and well
location. Dimensions and well location can be adjusted by clicking-and-dragging over
any of the light blue arrows. Results of pi and OGIP, OOIP, and OWIP are also dis-
played in the upper-right-corner of the view for easy reference.
l The x-y plan views display the pressure or saturation gradients at each time step.
The results can be viewed during, or after the simulation run using the standard video
playback buttons on the toolbar within the view.
This tab has two sections: a Results plot and a Forecast table.
Results Plot
The Results plot shows the production history, pressure data, and the model match, if they exist. In
addition, forecasted data is displayed, which consists of a fluid rate, forecasted pressure, and the
average reservoir pressure. If more than one forecast case is present, the results of each are dis-
played on the Results plot.
Forecasting
Forecasting in the Advanced Models takes the production history into account and does not com-
mence assuming static reservoir conditions. A forecast may be generated based on your reservoir
model parameters and specified forecasting parameters.
Forecast Parameters
Several additional features exist in the Forecast pane that can be used to customize the forecast.
This pane consists of the following three sections:
1. Forecast Options
2. Forecast Constraints
3. Forecast Results
Note: Each of these sections can be expanded or collapsed by clicking the +/- box.
Forecast Options
The Forecast Options section is where forecast periods are defined according to either durations
(e.g., months), or dates.
l Setting the forecast time method to Duration begins a forecast at a specified start
date. Forecast periods are then defined in the forecast table according to a length of
time (the default is months). For each forecast period, different operating conditions
can be specified.
If you select Duration from the Forecast Time Method drop-down menu, the Start Date is auto-
matically populated. (You can change this date later, if needed.)
If you select Start Date from the Forecast Time Method drop-down menu, enter a value for Duration
(i.e., the total length of the forecast).
The Forecast Flowing Pressure can also be set in the Forecast Options section. If a Gas Condens-
ate well is being analyzed, the Separate Recombined Gas option is available in Forecast Options.
Further details about each forecast period are entered in the forecast options table.
l Time – defines the length of each forecast period and each timestep, and has the fol-
lowing columns:
l Control – defines how the operating conditions change over the forecast period, and
has the following columns:
l Gas/Oil/Water Rate – If the Control Type is set to Gas, Oil, or Water, this column is
displayed, and it is used to set the rate for the forecast period. If Interpolation is set to
Step, only the initial rate's cell is editable. If Interpolation is set to Ramp, initial rate
and final rate are editable.
Forecast Constraints
The Forecast Constraints section is where maximum rate conditions and abandonment rate con-
ditions for the forecast can be entered. Available constraints depend on the analysis type, and could
include the following:
l pmin – The minimum allowable sandface flowing pressure during the forecast. In
order to ensure flowing pressure does not go below this value, rates are adjusted.
l (qg)max – Sets a maximum gas rate during the forecast. In order to maintain the max-
imum gas rate constraint, flowing pressure is adjusted.
l (qw)max – Sets a maximum water rate during the forecast. In order to maintain the
maximum water rate constraint, flowing pressure is adjusted.
l (qo)max – Sets a maximum oil rate during the forecast. In order to maintain the max-
imum oil rate constraint, flowing pressure is adjusted.
l (qg)ab – Sets the abandonment gas rate for the forecast. When the abandonment
rate is reached, the forecast ends.
l (qw)ab – Sets the abandonment water rate for the forecast. When the abandonment
rate is reached, the forecast ends.
l (qo)ab – Sets the abandonment oil rate for the forecast. When the abandonment rate
is reached, the forecast ends.
For analytical models, the available constraints include a minimum flowing pressure, as well as max-
imum and abandonment rates corresponding to the fluid being analyzed. For example, a gas ana-
lytical model includes pmin, (qg)max, and (qg)ab. Numerical models include minimum flowing
pressure, as well as maximum and abandonment rates for all fluids in the model.
Note: For wells with historical data, EUR is defined as the cumulative of the historical data up to
the beginning of the forecast, after which time, the synthetic cumulative calculated from the
model is used until the end of the forecast.
For typecurve and analytical models, EUR and RR of the fluid being analyzed are available. When
analyzing a gas condensate system, EUR and RR of the condensate phase are also available. For
numerical models, EUR and RR of all the fluids in the model are available.
Tables Tab
The Tables tab contains a sub-tab for the History Match and another for the Forecast. These sub-
tabs display the results of the simulation in tabular format.
Note: A checkmark indicates the last wizard step that was active when the wizard was closed.
The Wizard is designed to guide you through the software from start to finish. As a result, we recom-
mend that you follow each step in the order presented. To move to the next step in the Wizard, click
the Next> button to continue. Wizard steps may be grayed-out if they don't apply to the current file,
or if information from preceding wizard steps is required.
Click the following links to learn more about specific Wizard steps:
l Or, the pressure adjustment can be calculated based on gauge depth, datum depth,
and fluid gradient.
The pressure adjustment accounts for the hydrostatic head between the gauge depth and datum
depth. As a result, true vertical depth (TVD) should be used for gauge and datum depths when cal-
culating the pressure adjustment. If the gauge depth is below the datum depth, the adjustment will
be negative.
Pressure adjustment = ( Datum Depth (TVD) – Gauge Depth (TVD) ) * Fluid Gradient
Clicking the Show Example button opens the help topic corresponding to that flow regime.
Note: The Analysis Lines option is only enabled in the Wizard menu when the current view is an
analysis or model tab.
If your test type is Minifrac, you will be prompted to select one of the following diagnostic analyses:
In addition to the required inputs, it is often appropriate to modify parameters that have default val-
ues by editing them in the Parameters dialog box.
Saturations
Default values for saturations should be adjusted to reflect your reservoir conditions. These are used
to calculate total compressibility, which is used in all of the analysis line calculations. For the gas
fluid type, gas saturation defaults to 80%, oil saturation defaults to 0%, and water saturation defaults
to 20%. For the oil fluid type, gas saturation defaults to 0%, oil saturation defaults to 80%, and water
saturation defaults to 20%. For the water fluid type, gas saturation defaults to 0%, oil saturation
defaults to 0%, and water saturation defaults to 100%.
Formation Compressibility
The default value for formation compressibility, cf, is calculated using the entered porosity and Hall’s
correlation for consolidated sandstones, and is not suitable for all reservoirs. For example, naturally
fractured, geo-pressured and coalbed methane reservoirs have higher formation compressibility.
Formation compressibility affects the total compressibility, which is used in analysis line calculations.
Wellbore Radius
Wellbore radius (rw) defaults to 0.3 ft (0.0914 m) and should be changed to suit your wellbore dimen-
sions.
By selecting Yes, the Wellbore Calculations dialog box opens. For additional information, see Cal-
culating Sandface Pressures.
l Data from different gauges or datasets contained within a single gauge can by syn-
chronized.
l Data from different gauges (i.e. imported from different files) can be merged .
Any data required for analysis must be transferred to the Production Editor. Click Next to open the
WellTest Wizard: Select data to analyze dialog box.
The effective producing time can be entered manually, when the radio button beside Eff. Prod.
Time is selected. Or, it can be calculated using cumulative production and the rate of the primary
fluid, when the radio button beside Cum. Prod. is selected.
To calculate the effective producing time using this dialog box, select the primary reservoir fluid type
from the Reservoir Fluid Type drop-down list, and enter cumulative production and rate.
After effective producing time has been entered or calculated, it can be applied by clicking Next or
Finish, with the following results:
l The date of this row is back-calculated to honour the effective producing time. The
cumulative time in row 1 is zero, and the cumulative time in row 2 is the effective pro-
ducing time.
l The gas, oil and/or water rates that were entered in the Effective Producing Time dia-
log box will be displayed in row 2. Since WellTest displays end of period rates, pro-
duction at these rates starts at the date specified in row 1, and continues to the date
specified in row 2.
Review the data table to ensure that previous production rates have been applied over the correct
time interval by scaling to display a cumulative time of zero.
The data is segmented into flow and shut-in periods and different filtering parameters may be
applied to each period. It is very important that the flow and shut-in periods have been correctly spe-
cified before filtering. With the checkboxes you can select the filter mode (e.g., logarithmic, arith-
metic, or no filter) for each period. The values in the Logarithmic Points Per Cycle and Arithmetic
Time Step columns are editable. The Current and Remaining Data Points columns keep track of how
many points there originally were, and how many points remain after the filter is applied. The final
row shows the total number of points before and after filtering.
Logarithmic filtering keeps more data at the start of the flow period, when the pressure is changing a
lot with time, and less data at the end, when the change in pressure is more gradual. (This is the
default for all shut-in periods.) Arithmetic filtering takes evenly spaced points (e.g., 1 point per hour)
and is the default for all flowing (or injecting) periods.
The filtering defaults are very conservative and keep way more data than is necessary. For example,
a handful of derivative points could accurately describe the trend within any log cycle, including a
For conventional tests, the data in the History plot is partitioned into flowing (and injecting) and shut-
in periods. Any period can be selected for analysis and the selected period is displayed with a gray
background. To select a period, click the FP button on the History plot’s toolbar.
Then, click the History plot in the region you would like to analyze.
It helps us assess the suitability of the model, characterize the reservoir, and calibrate the model
prior to forecasting.
This dialog box is only used in the Models / Plots tab, and it describes how history matching can be
done.
Select any available model from the list, and click Next to create a new model tab in WellTest.
In many cases, prior production data may be imported and brought into the Production Editor. It is
common for prior production and test production data to come from different ASCII files. In this scen-
ario, you may import both files separately, merge them, and then bring the merged dataset into the
Production Editor. See Loading Data into the Production Editor for additional information.
When prior production data is not available in a format that is easy to import, you can enter it into the
Production Editor manually, or specify an effective producing time.
Using an effective production time will provide reliable analysis results, if production / injection peri-
ods are not interrupted by significant shut-in periods, or major rate fluctuations. Otherwise, the pro-
duction / injection should be entered in sufficient detail to provide a good representation of the rate
history.
Gas gravity impacts the oil correlations and should be entered for oil wells. To enter properties of the
non-primary fluid types, navigate to the Gas Properties, Oil Properties, and Water Properties sub-
tabs within the Properties tab and manually enter the properties (outside of the wizard).
Oil Properties
The oil gravity defaults to 0.88 (30 API), and should be changed to reflect the actual relative density
of your reservoir oil. Saturated or undersaturated conditions must be specified to calculate the appro-
priate fluid properties for your reservoir conditions. The default is undersaturated oil.
Gas Properties
The gas gravity defaults to 0.65, and should be changed to reflect the actual relative density of your
reservoir gas. If there is associated condensate production, the gas equivalent of the condensate
can be estimated to establish the total gas rate and recombined gas gravity.
See Deleting Data in the Production Editor for a procedure on how to remove undesirable data.
Note: For drawdown / buildup tests, fluid types will be grayed-out unless a rate for the fluid type
has been specified in the Production Editor.
Multiphase Options
Numerical models handle multiphase flow more rigorously than analytical models, and do not use
the reservoir fluid type specified in this wizard. The multiphase options in this dialog box, gas con-
densate (two-phase skin), gas condensate (two-phase pseudo-pressure: gas condensate), and oil &
gas (two-phase pseudo-pressure: solution gas drive), are advanced, but approximate methods of
handling multiphase flow analytically, and do not impact numerical models.
Note: Detailed rock and fluid properties are required to use the multiphase fluid types. These prop-
erties may be entered in the Properties / Multiphase Properties tab.
Use the drop-down menus beside each data type to select the datasets that you want to use for your
analysis. The Production Editor table will update according to your selections.
Properly identifying the shut-in point in pressure / rate history data is critical to a proper analysis. If
the shut-in is not placed at the correct time, the analysis plots may not display as expected, and the
data will be more difficult to analyze.
Data is segmented into flow and shut-in periods for filtering, and diagnostic analyses based on the
location of the zero-rates in the production history. Also, the shut-in point is used to identify the final
flowing pressure (pwfo), which is used to calculate total skin (s').
l Can modify imported data. When the test type is Injection / Falloff or Minifrac, impor-
ted rates are converted to negative values.
l Affects default settings. For example, the reservoir pressure defaults to the highest
measured pressure for Drawdown / Buildup tests, and the lowest measured pressure
for Injection / Falloff tests.
l Impacts error checking intended to prevent you from making mistakes. For example,
the cushion pressure (pwfo) must be less than the initial pressure (pi) for Perforation
Inflow tests, and greater than for Perforation Injection tests.
Single-Phase Pseudo-Pres-
Deconvolution Theory
sure
Nomenclature Reporting
In the changing liquid level case, when a producing well is shut-in, the wellbore is not full of liquids
and there exists a gas-liquid interface in the wellbore. As afterflow continues, this liquid level rises in
the case of a buildup test (falls during a fall-off test). The wellbore storage constant in a changing
liquid level case can be expressed as:
(field)
where
Vu = wellbore volumetric capacity, bbl/ft
(metric units)
where
Vu = wellbore volumetric capacity, m3/m
Note: The wellbore volumetric capacity in the above equation is essentially the effective wellbore
area.
If we define:
Where z is the variable used in the least squares minimization process to calculate the derivative
and Pu such that the total least squares error ETLS is minimized.
The error terms used to determine the minimum error are defined as follows:
Where pm is the measured pressure and pc is the pressure calculated from the preceding super-
position equation.
Where qm is the measured rate and qc is calculated by the least squares minimizer. Erate is zero
unless the Adjust Rates option has been turned on.
Where wx is the weighting parameter for each Ex error term. These values can be adjusted to
improve the fit when needed by placing more emphasis on specific error terms.
References
This deconvolution method is based on the work of Thomas Von Schroeter et al. and Michael M.
Levitan et al. for more information on the concepts of deconvolution please refer to the following
papers: SPE 71574, SPE 77688, SPE 84290, SPE 90680.
The following are typical derivative and pressure-time plots with the different time categories
marked:
Wellbore Con-
Early Time Middle Time Transition Late Time
figuration
n Wellbore n Radial Flow n Single No n Pseudo-
Storage Flow Steady
Boundary State
n Linear Flow
Fracture n Linear
Vertical Flow Channel n Steady
Wells Flow State
n Bilinear Flow
Fracture
Flow
n Spherical
Flow
n Wellbore n Horizontal n Linear n Pseudo-
Storage Radial Flow Channel Steady
Flow State
n Vertical Flow
Radial
Horizontal Flow n Steady
Wells State
n Linear Flow
Horizontal
Flow
n Elliptical
Flow
n Wellbore n Early Linear n Pseudo-
Storage Flow Steady
(toward frac- State
n Vertical tures) Flow
Radial
Flow n Early Radial
within the Flow
fractures (around
Multi-frac-
each frac-
tured Hori- n Linear ture prior to
zontal Wells Flow interference
(MFHWs) within the between
fractures fracs)
n Bilinear n Compound
Flow Linear Flow
n Late Radial
Flow
Note: Data affected by wellbore storage contains little or no information about the reservoir.
Wellbore storage is typically controlled by the compressibility of the fluid in the wellbore. For a gas-
filled wellbore, compressibility is high, and wellbore storage effects occur over a longer period of
time. For a liquid-filled wellbore, compressibility is much lower, and wellbore storage effects dis-
sipate more quickly. In some cases, typically in oil wells, both gas and liquid are present within the
wellbore and the liquid level changes after shut-in. In these cases, wellbore storage is also affected
by the changing liquid level, as well as compressibility.
Wellbore storage can be minimized by using a downhole shut-in. This operation reduces the well-
bore volume and consequently the wellbore fills more quickly so you can see reservoir-dominated
flow faster. This kind of operation is typically used to conduct a Closed Chamber Test (CCT) analysis
to reduce the amount of time needed to gather data to see reservoir effects.
This linear flow regime is expected to be the dominant flow regime, as demonstrated from production
analysis results (Nobakht, 2011).
Elliptical Flow
Elliptical Flow occurs when fluid has started to flow from the reservoir at either end of the horizontal
wellbore. It is a transition between linear horizontal flow and horizontal radial flow.
Note: Linear flow within the fractures is not normally observed in data because it is masked by well-
bore storage.
Radial Flow
In this flow regime, flow is in the horizontal radial direction. This type of flow exists in the time period
before the pressure transient has reached the boundaries of the reservoir (infinite-acting time
period).
Radial flow is often observed following other flow regimes. The emergence of radial flow in an infinite
conductivity hydraulically fractured well is shown below. When the radius of investigation (rinv) is
small, linear fracture flow is observed, and as it expands, the flow increasingly becomes radial.
Mathematically, a situation in which a well is next to a sealing fault can be modeled by removing the
fault, and placing an image well with a flow rate equivalent to the producing well as shown in the dia-
gram below.
As shown on the plots below, the pressure response on a semi-log or derivative plot shows a doub-
ling of slope when a single no-flow boundary is present.
Transition Region
The transition region or period occurs between the time the radius of investigation (rinv) reaches the
closest boundary, and the time it reaches the furthest boundary.
Note: This flow regime is not observed in data because it is masked by wellbore storage.
Gas Condensate
If water is being produced during multiphase flow, there are two possible ways to deal with it:
1. If the water flow is caused by "coning" at the well, then that effect is already included as part
of the skin due to damage (sd). There is no way of quantifying the skin effect due to water
production independently.
2. If the water is flowing throughout the formation (and is not just a water cone at the well), then
the skin damage (sd), calculated is indeed reflective of the true damage at the wellbore. The
effects of the water flowing in the reservoir are reflected in the value of kg being a function of
both oil (condensate) and water saturation’s in the bulk formation.
1. If the water flow is caused by "coning" at the well, then that effect is already included as part
of the skin due to damage (sd). There is no way of quantifying skin effect due to water pro-
duction independently.
2. If the water is flowing throughout the formation (and is not just a water cone at the well), then
the skin damage (sd), reflects the true damage at the wellbore. The effects of the water flow-
ing in the reservoir are reflected in the value of ko and kw which are a function of both oil and
water saturation's in the bulk formation.
Analytical and numerical models, with varied MFHW completion details and reservoir properties,
have been used to investigate which of the possible flow regimes can be practically identified in
buildup data (Lougheed, 2013). The following list pertains to both oil and gas reservoirs.
Likelihood of
Flow Regime
Observing
Wellbore storage Likely
Vertical radial flow within fractures Unlikely
Linear flow within fractures Unlikely
Bilinear flow Possible
Early linear flow Likely
Early radial flow (around each fracture, prior to interference
Unlikely
between fractures)
Late (Compound) linear flow Possible
Late radial flow (around MFHW and fracture network) Unlikely
Boundary-dominated flow (PSS flow) Unlikely
As shown below, the buildup derivative shows more character than the drawdown derivative
(Lougheed, 2013). This is caused by superposition. Flow time influences the buildup derivative beha-
vior (Streltsova, 1984).
l The early transition does not occur for high FCD values.
l Larger storage values and shorter flow times prolong this transition region.
l Shorter flow times and larger fracture half lengths make the transition appear flatter.
As learned from well test interpretation of conventional wells, the presence of heterogeneity can
greatly impact the signature of the buildup derivative. Many MFHWs drilled in unconventional reser-
voirs, such as the Eagle Ford, behave as if there is a region around the fractures with an enhanced
permeability. The enhanced permeability region simulates a complex fracture network (branch frac-
tures). The Enhanced Frac Region model (Stalgorova and Mattar (2012)) was used to investigate
buildup trends in the presence of an enhanced permeability region (Lougheed, 2013).
In their study, Dylan Lougheed and Marty Santo found that either an enhanced or a reduced per-
meability around the fractures could result in misleading derivative signatures. A reduction in per-
meability around the fracture is possible if a buildup test is conducted before the well has cleaned
up; because of relative permeability effects in the fluid invaded zone. The following figures are inten-
ded to illustrate examples of buildup derivative signatures that can occur in the presence of het-
erogeneity in MFHWs.
l Gauge malfunction
Note that this differs from the typical definition of the derivative:
By definition, the PPD for any type of reservoir flow regime will always either be a constant or
decreasing. Thus, any non-reservoir or wellbore effects will appear as an increase in the PPD as
shown in the following plot.
The equation for the flow of gas in the reservoir is very similar to that for liquid flow. In well testing,
analytical equations are solved after making certain assumptions. In particular, four assumptions are
very important:
4. Fluid saturations (Sw and Sg) are constant, if Swi is not equal to 0
For liquids, these assumptions are reasonable, since liquid compressibility and viscosity do not vary
significantly with pressure, and the equations can be solved analytically. These analytical solutions
are referred to as the liquid flow solutions, and form the basis of all well test analysis. The result is an
analytical relationship between pressure and time, which, for an infinite-acting reservoir, can be writ-
ten as:
Pressure = Constant * log (time) + ...
For gas, the assumptions listed above are no longer valid, since gas compressibility (cg) and thus
the gas compressibility factor (z) can vary significantly with pressure. Gas viscosity (µg) also varies
with pressure, but not to the same degree. To deal with these changing gas properties, the concept
of pseudo-pressure (ψ) was developed by Al-Hussainy et al. (1966), which is defined as follows:
Note: In some sources, pseudo-pressure is also referred to as the real gas potential (m(p)).
A pressure transient is created when a disturbance such as a change in rate occurs at a well. As time
progresses, the pressure transient advances further and further into the reservoir. This concept is
not theoretically rigorous, but is adequate for practical purposes. Theoretically, when a pressure dis-
turbance is initiated at the well, it has an immediate effect, however minimal, at all points in the reser-
voir. At a certain distance from the well, however, the effect of the disturbance are so small as to be
unmeasurable. The furthest distance at which the effect is detectable is called the radius of invest-
igation, rinv.
The figure below illustrates the basic concept of radius of investigation using a plot of pressure
versus distance into the reservoir.
Note: The constant 948 is one of many that appear in the literature, and in general provides a good
approximation.
Notice that the radius of investigation is only a function of the reservoir properties, and does not
depend on flow rate. Increasing the flow rate results in a greater drawdown (p - pwf), but the radius
of investigation is the same.
The examples below illustrate how radius of investigation can be used. For information on displaying
the radius of investigation in the software, see Displaying Radius of Investigation.
This represents the minimum drainage area before boundary effects are observed. This same applic-
ation can be applied to estimate distances to reservoir heterogeneities.
In this situation, the drainage area of a well is defined by the no- flow boundaries. Using this area,
you could estimate the time it would take to reach the boundaries of the drainage area by rearran-
ging the radius of investigation equation as follows:
1. Tabulate gas viscosity (µg) and compressibility factor (z) at various pressures.
3. Obtain the slope (m) of the semi-log plot of ψ versus log (time function), or from the deriv-
ative or regression analysis.
Buildup:
6. If rock relative permeability data and condensate PVT data are available, krg is calculated as
a function of pressure, as outlined below:
First, kro/ krg is calculated from one of the following equations shown below:
If ρg and ρo are mass densities and Mo and Mg are the molecular weights (g / mol):
If relative permeability data (kro, krg, So) are available, a tabulation of kro/ krg as a function
of oil (condensate) saturation (So) is generated:
Buildup:
This choked zone causes an additional pressure drop that tends to negate some of the benefits of
the fracture. Cinco-Ley and Samaniego (1977 and 1981) have defined this pressure drop as a
choked skin defined as:
References
1. "Effect of Wellbore Storage on the Transient Pressure Behavior of Vertically Fractured
Wells", H. Cinco-Ley and F. Samaniego-V., Paper SPE 6752 presented at 1977 AFTCE,
Denver, CO, October 9 - 12.
l Rate (q)
l Permeability (k)
For gas, the above formula is used combined with pseudo-pressure as shown below:
Note that this equation is a little more complex since it involves converting the pseudo-pressure
term, shown in brackets, to pressure. This can be accomplished by using a pressure versus
pseudo-pressure table.
The skin effect is a dimensionless quantity and is defined as the difference between the actual and
the ideal dimensionless pressure drop in a reservoir, or pressure drop due to skin (∆pskin).
When skin is viewed as an actual pressure drop rather than a dimensionless value, it becomes
easier to determine whether skin is actually a problem, and if steps need to be taken to correct it.
An alternative concept to the skin effect is the effective wellbore radius (reff).
This states that a well with improvement (negative skin) is equivalent to a well with a larger wellbore
radius (rw), while a well with damage (positive skin) has a smaller effective wellbore radius (reff).
For example, in a well that does not penetrate a formation vertical to the bedding plane, the com-
munication (or contact) area with the formation is increased. This reduces the pressure drop
required to obtain a flow rate equal to that of a well that penetrates a formation vertical to the bedding
plane. Therefore, the flow efficiency is improved, which causes a reduction in the apparent or total
skin factor (s'). This reduction in skin factor is referred to as the skin (or pseudo-skin) due to inclin-
ation (sinc). Cinco et al. (1975) have provided the following correlation for calculating sinc:
For a vertical well, θw = 0°, and sinc = 0. Thus, the skin due to inclination is always equal to or less
than 0.
By rearranging the total skin equation, the skin due to damage (sd) can be determined when skin
due to inclination (sinc) and the other skin components are known.
References
"Unsteady-State Pressure Distribution Created by a Directionally Drilled Well", H. Cinco, F.G. Miller
and H.J. Ramey, Jr., JPT (November 1975).
Skin due to partial penetration is always greater than 0, and typically ranges from 0 to 30. Note that
partial penetration greatly magnifies the effect of any skin due to formation damage (sd).
This is more clearly seen in the definition of total skin (s'). By rearranging the total skin (s') equation,
the skin due to damage (sd) can be determined when skin due to partial penetration (spp) and the
other skin components are known.
2. "Numerical Simulations of the Combined Effects of Wellbore Damage and Partial Pen-
etration", R.M. Saidikowski, Paper SPE 8204 presented at 1979 AFTCE, Las Vegas, NV,
September 23 - 26.
By definition, this additional pressure drop or skin is a function of gas flow rate (qg) and the tur-
bulence factor (D) of the system expressed as:
Note that skin due to turbulence is always positive and is one component of the total skin (s'). Thus,
a production test on a stimulated well can still yield a positive total skin (s') value, even if no skin dam-
age (sd) is present, due to the turbulence component (sturb).
By rearranging the total skin (s') equation, the skin due to damage (sd) can be determined when skin
due to turbulence (sturb) and the other skin components are known.
References
"Non-Darcy Flow and Wellbore Storage Effects in Pressure Build-Up and Drawdown of Gas Wells",
H.J. Ramey, Jr., JPT (February 1965) 223 - 233.
By substitution, the effective fracture radius can be defined as an effective wellbore radius as fol-
lows:
As shown, it is always a negative skin since fracturing is a stimulation technique used to reduce total
skin (s') to improve fluid flow.
Even a small skin on the fracture face, sf < 1, can cause a substantial change in the shape of the
derivative curve. Typically, sf ranges from 0 to 1.
References
1. "Effect of Wellbore Storage on the Transient Pressure Behavior of Vertically Fractured
Wells’, H. Cinco-Ley and F. Samaniego-V., Paper SPE 6752 presented at 1977 AFTCE,
Denver, CO, October 9 - 12.
3. "Optimization of the Productivity Index and the Fracture Geometry of a Stimulated Well with
Fracture Face and Choke Skins", D.J. Romero, P.P. Valko and M.J. Economides, SPEPF
(February 2003) 57 - 64.
l Skin due to partial penetration (spp) for a partially penetrated well only
l Skin due to turbulence (sturb) or non-Darcy flow (for gas wells only)
Usually, radial analysis provides the total skin (s') of the system. The value of s' can be positive, neg-
ative, or zero. Sometimes it is important to know what skin components are contributing to the total
skin. To determine this, the relationship between s' and its various contributing components can be
expressed as:
Except for skin due to damage (sd), all other skin components mentioned above are always non-neg-
ative (i.e., are either zero or positive). Note that in the case of a partially-penetrated well, the skin
due to damage (sd) is magnified by a factor of h / hp in the total skin value. Also note that the skin
due to turbulence (sturb) in a gas well is rate sensitive.
Keep in mind that because total skin is a summation of various other skin components, a positive or
negative value may not indicate whether a well should or should not be stimulated. For example, for
an unstimulated well, a positive value of total skin does not necessarily mean that the well is
damaged. It could just be that the sum total of the skin components results in a positive value. Sim-
ilarly, for a stimulated well, a zero or positive value of total skin can result, even if the skin due to dam-
age (sd) is negative. Therefore, the effectiveness of the stimulation is really represented by the skin
due to damage (sd) rather than the total skin.
The skin due to this extra pressure drop of the gas phase is referred to as the two-phase skin. The
method proposed by Raghavan et al. (1995) is used to calculate the s2p value, when pressure,
volume, temperature (PVT) properties and relative-permeability data of the gas condensate system
are available. This method allows you to separate the two-phase skin from the total or effective skin
(s'). The value of s2p can either be zero (when the well pressure is above the dew point pressure,
and no condensate bank exists), or a positive number (when the well pressure is below the dew
point pressure, and a condensate bank exists). If the PVT and relative permeability data are not avail-
able, it is not possible to calculate s2p exclusively. In this case, only the total or effective skin (s')
(which includes s2p) can be calculated. The details of the process of calculating s2p are outlined in
the single phase pseudo-pressure procedure for a gas condensate reservoir.
Typically, two-phase skin (s2p) ranges from 0 to 30. Note that by rearranging the total skin (s') equa-
tion, the skin due to damage (sd) can be determined when two-phase skin (s2p) and the other skin
components are known.
References
"Practical Consideration in the Analysis of Gas Condensate Well Tests", R. Raghavan, W., C. Chu,
and J.R. Jones, Paper SPE 30576 presented at 1995 ATCE, Dallas, TX, October 22 - 25.
Although mathematically rigorous, this variable is totally impractical, as it requires knowledge of the
distribution of saturations in the reservoir, in space and in time. It only applies when a simulation of
the reservoir has been conducted.
A more practical alternative is to define a Sandface Integral (S.I.) which calculates pseudo-pressure
in terms of the saturation at a specified point in the reservoir, usually the sandface, as follows:
This too is impractical to evaluate, but it forms the basis for defining a two-phase pseudo-pressure
that can be evaluated, if appropriate rock and fluid properties are available. Although this definition is
only an approximation to the Reservoir Integral, and is not fully rigorous, it is adequate for analysis
and modeling of solution gas systems. The two-phase pseudo-pressure is defined as:
To calculate ψ2p, the variation of kro and krg with pressure must be known beforehand. In reality,
ψ2p can be defined in terms of kro alone, along with a pressure-saturation relationship. In other
words, the information contained in the krg(p) relationship is already contained in the combination of
kro(p) and pressure-saturation. This is a two-part procedure, the first part is obtaining the relative per-
meability data for the reservoir rock, and the second part is assuming an appropriate model for fluid
flow in the reservoir.
This definition appears to be simpler than the preceding one, which contains both the kro and krg
terms. However, in actual fact, the information contained in this simpler formulation is complete, and
it accounts for, implicitly, both krg and kro, through the pressure-saturation formulation. Because the
relative permeability terms and multiphase effects are built into the definition of two-phase pseudo-
pressure, the analysis results in the determination of the absolute permeability (not the effective per-
meability to oil – k not ko) and the skin due to damage (sd) (independent of gas saturation effects
near the wellbore – gas block).
For analyzing buildup and drawdown tests, the two-phase pseudo-pressure procedure is recom-
mended, provided suitable relative permeability and PVT data are available. The results of the ana-
lysis will include the absolute permeability (k), and the skin due to damage (sd). For drawdown tests,
the producing GOR is used to relate the pressure and saturation terms that occur in the pseudo-pres-
sure function. In the case of buildup tests, the pseudo-pressure function is calculated based on the
producing GOR at the instant of shut-in.
2. From the special core analysis data, tabulate: kro, krg, So, and calculate the relationship
between krg / kro and So.
3. From the laboratory PVT analysis, or from appropriate PVT correlations, tabulate PVT data
for p (psi), Rs (Mscf / stbbl), µo (cp), µg (cp), Bo (Rbbl / stbbl), and Bg (Rbbl / Mscf).
4. For the full range of test pressures, calculate the relationships of shut-in pressure (pws) vs.
Rs, µo, µg, Bo, and Bg.
5. From the producing GOR at the instant of shut-in and PVT data, calculate the relationship
between krg / kro and pws using the following equation:
6. From the relationships of pws vs. krg / kro and krg / kro vs. So obtain pws vs. So.
7. From pws vs. So obtain pws vs. kro using relative permeability data.
8. Calculate two-phase pseudo-pressure m(p) vs. pws using the trapezoidal rule:
9. From the radial semilog plot of shut-in pseudo-pressure (m(p)) vs. shut-in time (t), choose a
straight line which stands for the radial flow, and determine the slope (m) of this line. Cal-
culate absolute permeability (k) and total skin (s') as follows:
Gas Condensate
A more practical alternative is to define a Sandface Integral (S.I.), as shown below, which calculates
a pseudo-pressure in terms of the saturation at a specified point in the reservoir, usually the sand-
face as follows:
To calculate ψ2p, the variation of kro and krg with pressure must be known beforehand. In reality,
ψ2p can be defined in terms of krg alone, along with a pressure-saturation relationship. In other
words, the information contained in the kro(p) relationship is already contained in the combination of
krg(p) and pressure-saturation. This is a two-part procedure, the first part is obtaining the relative per-
meability data for the reservoir rock, and the second part is assuming an appropriate model for fluid
flow in the reservoir.
For a gas condensate system, the above relationship assumes that steady-state flow exists in the
two-phase region. This is a good assumption near the wellbore, where the oil (condensate) sat-
uration is greater than the critical oil saturation and the condensate is mobile. However, it is not valid
in the bulk of the reservoir, where the oil (condensate) saturation is less than the critical oil saturation
and the condensate is not mobile.
Even when the data is available for calculating ψ2p, it is recommended that, for gas condensate sys-
tems, ψ be used as the variable of analysis for calculating s2p and sd, and ψ2p be used to calculate
kg and s'. This recommendation stems from the limitations in the assumptions underlying the ψ2p
definition.
The literature recommends that the two-phase pseudo-pressure, ψ2p, be used for drawdown ana-
lysis, and single phase pseudo-pressure, ψ, be used for buildup analysis. However, drawdown data
is usually very erratic, and in view of the discontinuity that would arise from changing from one
pseudo-pressure to the other, single phase pseudo-pressure, is recommended to model both draw-
down and buildup data. In this case, two-phase skin would be calculated using the single phase
pseudo-pressure procedure. If the two-phase pseudo-pressure procedure is used instead, the skin
calculated represents the total skin (s') at the wellbore, excluding the two-phase skin component
(which is taken into account through the two-phase pseudo-pressure formulation).
If either the rock relative permeability data or the gas condensate PVT data are not available, then
ψ2p cannot be calculated. In such cases, only s2p and sd, can be determined from the pressure
1. First, kro/ krg is calculated from one of the following equations shown below:
If ρg and ρo are mass densities and Mo and Mg are the molecular weights (g / gmol):
If relative permeability data (kro, krg, So) are available, a tabulation of kro/ krg as a function
of oil (condensate) saturation (So) is generated:
6. Drawdown:
7. Buildup:
Gas equations
Oil correlations
Oil equations
Unit conversions
In the modified black oil (MBO) systems, in addition to the above behaviour, oil can condense out of
the gas phase. This process is dominated by pressure and temperature and its functionality is
defined by Vaporized Oil Ratio "Rv" (i.e., Hydrocarbon Liquid content of the gas). Strictly speaking,
Rv is the ratio of barrels of produced stock tank oil per one MMscf of produced separator gas. For
gas condensate and volatile oil system, we need to consider the liquid content in the gas phase.
In P,
1 20.809 -6.7095 0.5136 0 0
psi
2.49E-
2 APId 11.175 -1.2965 0.042311 -0.0005438
06
P Pressure, psi
TR Reservoir temperature, °F
Radial Analysis
Flow Capa-
city
Per-
meability
Pressure
Drop Due
to Skin
Flow Effi-
ciency
Damage
Ratio
Radius of
Invest-
igation
Linear Analysis
Fracture Half-Length
Channel Width
Skin Due to Xf
Vertical Permeability
Horizontal Permeability
Field Units
Radial Analysis
Flow Capa-
city
Total
Skin
Pressure
Drop
due to
Skin
Flow Effi-
ciency
Damage
Ratio
Radius
of Invest-
igation
Linear Analysis
Fracture Half-Length
Channel Width
Skin Due to Xf
Fracture Conductivity
Vertical Permeability
Horizontal Permeability
Where:
Where:
Gas Saturated
Where:
Undersaturated
The oil compressibility (co) used in this equation is obtained from the Vasquez and Beggs cor-
relation.
Reference
"Pressure-Volume-Temperature Correlations for Saudi Crude Oils", M.A. Al-Marhoun, SPE 13718,
1985.
Heavy Oils
A modified Vasquez and Beggs solution gas oil ratio (Rs) correlation is reversed to solve for the
bubble point pressure as follows:
Where:
Extra-Heavy Oils
A modified Standing solution gas oil ratio (Rs) correlation is reversed to solve for the bubble point
pressure as follows:
Heavy Oils
A modified Vasquez and Beggs solution gas oil ratio (Rs) correlation was developed as follows:
Where γg(psp) is calculated the same as above for bubble point pressure.
Gas Saturated
A modified Vasquez and Beggs oil formation volume factor (Bo) correlation was developed as fol-
lows:
Where:
Undersaturated
Oil Compressibility
Gas Saturated
Undersaturated
Modified Vasquez and Beggs oil compressibility (co) correlations were developed as follows:
Heavy Oils
Extra-Heavy Oils
Where γg(psp) is calculated the same as above for bubble point pressure.
Oil Viscosity
Dead Oil
Modified Egbogah-Jack’s oil viscosity (µo) correlations were developed as follows:
Heavy-Oils
Extra-Heavy Oils
Gas Saturated
Modified Kartoatmodjo’s oil viscosity (µo) correlations were developed as follows:
Heavy-Oils
Extra-Heavy Oils
Where:
Undersaturated
Heavy Oils
Extra-Heavy Oils
Reference
"Pressure-Volume-Temperature Correlations for Heavy and Extra Heavy Oils", Giambattista De
Ghetto, Francesco Paone, and Marco Villa, SPE 30316, 1995.
Where:
Where:
Gas Saturated
Where:
Undersaturated
The oil compressibility (co) used in this equation is obtained from the Vasquez and Beggs cor-
relation.
Reference
"Generalized Pressure-Volume-Temperature Correlations", Oistein Glaso, Journal of Petroleum
Technology, 1980.
Gas Saturated
Undersaturated
Oil Density
Note that the density is calculated in metric units (g / cm3).
Undersaturated
Oil Compressibility
Undersaturated
This correlation uses only the oil density (ρob) at the bubble point. Therefore the oil compressibility
(co) is constant for pressures greater than the bubble point.
Gas Saturated
The derivatives dBo/ dRs and dRs / dP were taken from the Vasquez and Beggs correlation.
Oil Viscosity
This correlation calculates the oil viscosity (µo) at any pressure using the corresponding oil density
(ρo).
Where:
Ng and Egbogah
The Ng and Egbogah correlation contains two methods for calculating dead oil viscosity (µod) using
a modified Beggs and Robinson viscosity correlation and a correlation that uses the pour point
temperature. Pour point temperature is the lowest temperature at which the oil is observed to flow
when cooled and examined under conditions prescribed in ASTM D97. The purpose of introducing
the pour point temperature into the correlation is to reflect the chemical composition of crude oil into
the viscosity correlation. To obtain the viscosity for live oils, the dead oil correlations are used with
the Beggs and Robinson viscosity correlation. The data used to derive the correlations was taken
from the Reservoir Fluids Analysis Laboratory of AGAT Engineering Ltd., using a total of 394 oil sys-
tems.
Dead Oils
Note the range for pour point temperature (Tpp) in this equation is from -50 to 15 °C.
Live Oils
Gas-Saturated
Where:
Where:
Reference
"An Improved Temperature-Viscosity Correlation for Crude Oil Systems", J.T.H. Ng and E.O.
Egbogah, Petroleum Society of CIM 83-34-32, 1983.
Where:
Where:
Gas Saturated
Undersaturated
Oil Compressibility
Gas Saturated
The derivatives dBo/ dRs and dRs / dP were taken from the Vasquez and Beggs correlation.
Undersaturated
Note that this equation is valid for an oil compressibility (co) between 2.464x10-5 to 3.507x10-5.
Gas Saturated
Undersaturated
The oil compressibility (co) used in this equation is obtained from the Vasquez and Beggs cor-
relation.
Where:
γo<= γo >
Coef-
30 30
ficient
°API °API
C1 0.036- 0.017-
2 8
C2 1.093- 1.187-
7 0
C3 25.72- 23.93-
40 10
Where:
γo<= γo >
Coef-
30 30
ficient
°API °API
C1 4.677- 4.670-
x10-4 x10-4
C2 1.751- 1.100-
x10-5 x10-5
-
C3 1.377-
1.811-
x10-9
x10-8
Undersaturated
Oil Compressibility
Gas Saturated
The derivatives dBo/ dRs and dRs / dP were taken from the Vasquez and Beggs correlation.
Reference
"Correlations for Fluid Physical Property Prediction", M.E. Vasquez and H.D. Beggs, JPT 968 - 70,
June 1980.
Where:
Where:
The reduced solution gas oil ratio (Rsr) is defined as the solution gas oil ratio (Rs) divided by the solu-
tion gas oil ratio at the bubble point (Rsb) as follows:
Using the above relationship the reduced solution gas oil ratio (Rsr) and the solution gas oil ratio at
the bubble point (Rsb) are used to solve for the actual solution gas oil ratio (Rs) at any pressure
below the bubble point.
Gas Saturated
Where:
After this initial value is known, ρpo is calculated through a 10 step iteration process using the fol-
lowing equations. The values from the ninth and tenth iterations are averaged to yield a final value
for ρpo.
The oil compressibility (co) used in this equation is obtained from the Vasquez and Beggs cor-
relation.
Correlation Limits
Rs Correlation Lim- Pbp Correlation
Variable
its Limits
T 70 to 307 °F 74 to 327 °F
Reference
"Correlation of Black Oil Properties at Pressures Below Bubble Point Pressure – A New Approach",
J. Velarde, T.A. Blasingame and W.D. McCain, Jr., The Petroleum Society 93 - 97, 1997.
Al-Marhoun 1985
0.752 – 1.367 14.3 – 44.6
(Saudi Arabian Oil)
De Ghetto et al
59 – 14.5 –
(Heavy and Extra-Heavy 0.623 – 1.517 6 – 22.3 2.4 – 354.6
177.8 752.2
Oils)
Glaso 0.119 –
0.65 – 1.276 22.3 – 48.1
(North Sea Oil) 106.6
Hanafy et al
0.752 – 1.367 14.3 – 44.6 0.13 – 71
(Egyptian Oil)
Khan et al
(Saudi Arabian Oil)
ρo (g / cm3) ρob (g /
Correlation µod (cp) µos (cp) µob (cp)
cm3)
Al-Marhoun 1985
(Saudi Arabian Oil)
De Ghetto et al 7.7 –
2.1 – 295.9
(Heavy and Extra-Heavy Oils) 1386.9
Glaso
(North Sea Oil)
Hanafy et al 0.648 –
0.428 – 0.939
(Egyptian Oil) 1.071
Khan et al
0.13 – 77.4 0.13 – 17.9
(Saudi Arabian Oil)
Ng and Egbogah
Petrosky and Farshad
(Gulf of Mexico Oil)
Standing
(California Oil)
Vasquez and Beggs
(Generally Applicable)
Velarde et al
(Reduced Variable
Approach)
Radial Analysis
Trans-
missivity
Per-
meability
Pressure
Drop Due
to Skin
Flow Effi-
ciency
Damage
Ratio
Radius of
Invest-
igation
Pro-
ductivity
Index
Linear Analysis
Fracture Half-
length
Channel Width
Skin Due to Xf
Vertical Permeability
Horizontal Permeability
Field Units
Per-
meability
Total Skin
Pressure
Drop Due to
Skin
Flow Effi-
ciency
Damage
Ratio
Radius of
Invest-
igation
Pro-
ductivity
Index
Linear Analysis
Fracture Half-length
Skin Due to Xf
Bilinear Analysis
Fracture Conductivity
Vertical Permeability
Horizontal Permeability
Metric (SI)
T = 15 °C
p = 101.325 kPa
Field
T = 60 °F
p = 14.65 psia
Common Conversions
Metric (SI) To convert from metric to field unit
Field Unit
Unit divide by
103m3 / d MMcf / d 28.17399
mD md 0.9869233
mD m md ft 0.3008142
m ft 0.3048
Pa s cp 1000
K °R 0.5556
Temperature Conversion
°C = (°F - 32) * 5 / 9
°F = °C * 9 / 5 + 32
Minifrac
PAS
Properties
Wellbore
The Alberta Energy Regulator (AER) is a regulatory body of the government. The purpose of the
AER is, in part, to provide appraisal of the reserves and productive capacity of energy resources in
Alberta, prevent waste of those resources, control pollution, ensure safe practices are followed, and
provide useful dissemination of information. The AER requires submission of .PAS files to meet reg-
ulatory requirements for wells.
Note: The AER was previously called the Energy Resources and Conservation Board (ERCB).
Anisotropic permeability is especially important when dealing with horizontal or partially penetrated
wells since flow occurs in both the vertical and horizontal planes.
In general, it is not uncommon to assume that the permeability in the x-direction (kx) is close to that
in the y-direction (ky). However, the permeability in the z-direction (kz) is typically significantly dif-
ferent, usually less, than the horizontal permeabilities (kx and ky).
Thus, vertical radial flow is controlled by ky and kz (where x is the direction of the horizontal well-
bore) and horizontal radial flow is controlled by kx and ky.
l Mead (Simplex)
l Marquardt-Levenberg (Gauss-Jordan)
APE should not be used exclusively since it is possible to find several sets of parameters that yield
an acceptable model match. In other words, the solution can be non-unique. It is always recom-
mended to start with parameters obtained from diagnostic analyses and then fine-tune these para-
meters manually to obtain a close match. APE can then be used for unknown parameters, or
parameters that are not known with confidence. The more parameters that are selected for APE or fit-
ting, the longer the fitting process takes, and the more chance there is of finding non-unique solu-
tions. Also, including more data points in the fit reduces performance.
Marquardt-Levenberg
The Marquardt-Levenberg method is a non-linear regression algorithm used for APE of reservoir
and well parameters when modeling well test data. It is a modified version of LMDIF, a public-
domain non-linear regression algorithm from Argonne National Laboratory. The algorithm requires
partial derivatives of the objective function with respect to each of the parameters (Jacobians). The
derivatives are calculated numerically, using a forward difference approximation. The objective func-
tion is the sum of squares of the difference (residuals) between pressure derivative data and the cor-
responding calculated model data.
The Marquardt-Levenberg routine is generally faster than the Mead (Simplex) method. However, the
requirement of derivative calculations in the routine tends to make it less robust and more com-
putationally intensive when farther away from the solution.
Mead (Simplex)
This a variation of the downhill Simplex method. The Simplex routine is a non-linear regression
algorithm used for APE of reservoir and well parameters when modeling well test data. It requires
only function evaluations of the objective function, and not the derivatives. The objective function is
the sum of squares of the difference (residuals) between observed pressure, or pressure derivative
data and the corresponding calculated model data.
Compared to other non-linear regression methods, this method is not always very efficient because
it can require a large number of function evaluations. This tends to make it extremely slow in some
cases. However, it is straightforward and not encumbered by the requirement of derivatives, and
hence tends to be more robust under any conditions.
1. By measuring the long-term buildup pressure in a bounded reservoir. The buildup pressure
eventually builds up to the average reservoir pressure over a long enough period of time as
shown below.
Note: This time depends on the reservoir size and permeability (k) (i.e., hydraulic diffusivity).
2. Calculating average reservoir pressure from the material balance equation (MBE) – as
described below.
Gas
For gas, the MBE is defined as the relationship between the original gas in place, initial pressure (pi),
cumulative gas production, and the current average reservoir pressure. The basis of the MBE for
gas flow is the volumetric balance of all the fluids at a given time. The following equation proposed
by Ramagost and Farshad (1981) is used to calculate the average reservoir pressure for gas sys-
tems. This equation considers that gas is the only mobile phase in the presence of residual fluid sat-
urations (oil and water) in a compressible formation.
Note that this equation is only valid when the term ce(pi– p) < 1.
For oil:
For water:
where:
This figure shows the wells that are evenly spaced and producing at the same rate, and thus no-flow
boundaries develop (white lines) between the producing wells. In this situation, the drainage area of
each well is defined by the no-flow boundaries around it.
Note: Drainage area can be determined through PSS analysis, radius of investigation, or analysis
of geological data or well spacing.
Quite often, the volume of hydrocarbon stored within the natural fractures is much lower than is
stored in the matrix. In dual porosity systems, the natural fractures have much higher permeability
than the matrix. When the well begins to flow, fluid travels from the high permeability natural frac-
tures to the wellbore, and is rapidly produced. After the natural fractures have been drained, the
large volume of hydrocarbons contained within the bulk of the reservoir (matrix) begins to flow.
These hydrocarbons flow to nearby natural fractures, and virtually all of the fluid is transported to the
wellbore via these fractures.
The signature of dual porosity systems on a semi-log plot is two parallel lines as shown below.
The first semi-log straight line is observed at early time and represents radial flow as the fluid, initially
in the fractures, travels to the wellbore. The second semi-log straight line occurs when the fractures
deliver fluid from the matrix to the wellbore. The transition period between the two semi-log straight
lines occurs when fluid begins to flow from the matrix to the fractures, but has not yet reached a state
of equilibrium.
Naturally fractured reservoirs are characterized by two parameters: the interporosity flow coefficient
(λ) and the storativity ratio (ω).
As shown in the plot below, as the interporosity flow coefficient decreases, the transition between
the two semi-log straight lines is delayed. That is, the larger the fracture permeability is in com-
parison to that of the matrix, the more time the fractures will have to drain before the contribution
from the matrix becomes significant.
Note that there are two types of interporosity flow that are used to model dual porosity systems:
The second parameter, storativity ratio (ω), essentially represents the time separation in log cycles
between the two semi-log straight lines as shown below:
The signature of dual porosity on a derivative plot shows up as two regions of radial flow with the
same conductivity, kh, separated by a transition period. This is often referred to as the dual porosity
dip.
Note that the geometric coefficient (α) accounts for the shape of the matrix blocks. The interporosity
flow coefficient is usually in the range of 10-4 to 10-8 and is always used in conjunction with the stor-
ativity ratio in dual porosity reservoirs.
A dual porosity system gives rise to two parallel straight lines on a semi-log plot. However, often the
first line is obscured by early time effects such as wellbore storage. Moreover, if the interporosity
flow is of the transient mode rather than the pseudo-steady state type, then in the transition between
the two parallel semi-log lines, there exists a semi-log straight line of half of the slope of one of the
parallel lines. In practice, these straight lines are difficult to observe uniquely. When the interporosity
skin is large, then the transient interporosity flow model becomes equivalent to the pseudo-steady
interporosity flow model.
References
"Gas Reservoir Engineering", J. Lee and R. Wattenbarger, Society of Petroleum Engineers Inc.
(1996) Volume 5, 73 - 74, 173 - 181.
"The Behavior of Naturally Fractured Reservoirs", J.E. Warren and P.J. Root, Society of Petroleum
Engineers Inc. (1963) SPEJ 426.
For a model, it is obtained from the continuation of the model line using the specified shut-in time
(tshut-in).
For a well with no depletion, it is equal to initial pressure (pi) as shown in the plot below.
For modeling, it is used as a default anchor point pressure. Thus, it is important to ensure that shut-in
periods are placed correctly to ensure a valid analysis and model.
Note: This is only used in non-conventional analyses and models (i.e., Slug, PITA, CCT).
In the above diagram, if Test 1 is being analyzed, the initial pressure is 2000. If Test 2 is being ana-
lyzed, the initial pressure is 1500, because the reservoir pressure is fully built up before this test, and
the reservoir has forgotten its history prior to Test 2. Thus, it behaves like a new well with an initial
pressure 1500. If Test 3 is being analyzed (which is a combination of Test 2 and some production
prior to Test 2), then the initial pressure is 2000, because it must reflect the reservoir pressure before
the production that is being analyzed during the test.
WellTest files can be imported into Harmony Enterprise, and Harmony Enterprise interchange files
(.hmexp files) can be opened in WellTest. However, scenarios are not imported from .hmexp files.
l Production Editor data — imported to both the Production Editor tab and Data Man-
agement's Prod Data sub-tab.
l Wellbore Parameters — WellTest has one wellbore and uses the Physical Wellbore /
Initial configuration in Harmony Enterprise.
l Analytical and hybrid models — only these models are supported. (CBM and numer-
ical models are not supported.)
l For oil, the liquid gradient is typically between 0.15 - 0.35 psi/ft or 4 - 8 kPa/m.
Introduction
Well testing has been used for decades to determine essential formation properties and to assess
wellbore conditions. There are many different types of tests that can be used to collect this inform-
ation depending on when the test is performed, the well location, the well type, and the formation
type. For the most part, conventional tests (flow / buildup or injection / falloff) have satisfied the major-
ity of our needs. However, under certain conditions, traditional test methods are not feasible for vari-
ous reasons. This is especially true for very low permeability formations that require massive
stimulation to obtain economic production. For these formations, it is extremely important to estab-
lish the formation pressure and permeability prior to the main stimulation. One test that has proved to
be convenient for this purpose is commonly referred to as a “Minifrac” test.
A minifrac test is an injection-falloff diagnostic test performed without proppant before a main frac-
ture stimulation treatment. The intent is to break down the formation to create a short fracture during
the injection period, and then to observe closure of the fracture system during the ensuing falloff
period. Historically, these tests were performed immediately prior to the main fracture treatment to
obtain design parameters (i.e., fracture closure pressure, ISIP (Instantaneous Shut-In Pressure)
gradient, fluid leakoff coefficient, fluid efficiency, formation permeability, and reservoir pressure).
However, since personnel and frac equipment were all waiting on location to perform the main treat-
ment, the falloff period was usually stopped shortly after observing closure, before reliable estimates
of formation pressure and permeability could be obtained. Since these two parameters are critical to
the fracture design and for production / reservoir engineering, it seemed prudent to extend the falloff
period to obtain better estimates, especially since there is little hope of gathering this information
after the main stimulation. Many operators have accomplished this by simply scheduling the minifrac
test well ahead of the main fracture treatment. However, predicting the falloff time required to obtain
meaningful estimates of formation pressure and permeability is difficult, as it depends on having
prior knowledge of the permeability, in addition to knowing the geomechanical properties of the form-
ation. In many cases, the progress of a minifrac test can be assessed with pressure data measured
at the wellhead, eliminating the need for “guessing” when sufficient data has been obtained.
1. Engineers responsible for designing and completing the main hydraulic fracture treatment
prefer to pump the planned fracturing fluid at a high rate / step rate in order to obtain more
representative estimates of ISIP, ISIP gradient, net fracture pressure, fluid efficiency and
fluid loss coefficients. This information is used to help optimize pad volume, select the best
fluid-loss additives for the main treatment, and design the pumping schedule that would cre-
ate the optimum fracture from a productivity point of view (SPE 140136). Step rate tests can
give additional information about matrix leakoff (at low rates) and frac extension (after break-
down at higher rates) (SPE 78173). The higher rates give a bigger minifrac representation
and distinct leakoff characteristics. Unfortunately, a bigger minifrac means it will take longer
to close. In general, these tests are designed to see closure. In tight / shale formations, such
as Eagle Ford and Bakkan, we are not likely to collect enough after-closure data to get good
estimates of permeability, or initial pressure based on this type of test. In unusual cir-
cumstances where we can monitor falloff data for a few weeks, we can see reservoir- dom-
inated flow (linear / radial flow) for these tests.
Theoretically, the key parameters can be obtained from either method, but in practice the test object-
ives must be weighed, and the test designed to meet those objectives. Alternatively, two tests can
be performed back-to-back. First, a test can be performed with the injected volume minimized to
obtain virgin-rock breakdown pressure, permeability, initial pressure, and initial estimates of leakoff
characteristics. Then another test can be performed with the fluid and rate more representative of
the main treatment.
l ISIP gradient
l Fluid efficiency
The convergence of the flow lines near the wellbore cause an additional pressure drop near the well-
bore, an effect similar to that caused by wellbore damage. In some cases, when the penetrated
region is very small compared to the formation thickness or net pay (h), spherical flow may be
observed. The effects of partial penetration are accounted for by treating it as a skin effect called
skin due to partial penetration (spp). This skin is always positive and typically varies from 0 to 30. It is
a function of the height of the perforated interval (hp), the distance from the top of the zone to the top
of the perforations (htop), and the horizontal to vertical permeability ratio (kh/kv).
The perforated interval has a maximum value dependent on the net pay (h) and the net wellbore
inclination (σ).
When a formation is only partially penetrated, the location of the perforated interval has an effect on
the skin due to partial penetration (spp). Thus the diagrams below show the same net pay (h) and
For example, in a well with partial penetration the fluid has to travel vertically because the whole of
the net pay (h) is not open to the wellbore as shown below.
A large horizontal-to-vertical permeability ratio implies a relatively low vertical permeability, which
creates a larger pressure drop near the wellbore due to the vertical component of flow. Thus, this
increase in pressure drop near the wellbore is represented as an increase in the skin due to partial
penetration (spp).
TRG, GRD and PRD PAS files can be read and generated by IHS applications to meet the AER’s
requirements for initial pressure, deliverability, and production. IHS applications do not read or gen-
erate fluid analysis PAS files (GAN, OAN, or WAN) For more information about these PAS files, try
contacting your fluid sampling provider (e.g., Core Labs, AGAT Laboratories, etc.).
If you do not deal with wells in Alberta, Canada, you can hide the PAS tab.
Important Notes
There are exceptions to the requirements described below (e.g., non-control wells in a Development
Entity, or wells producing under Self-declared (SD) Commingling). Please refer to the AER website
re: Directive 40 and 65 to learn more about these exceptions.
Exemptions on Pressure Survey requirements can be requested. These have been granted for low
permeability, heavy oil, small reserve, declining productivity pools, etc. Please refer to the AER web-
site re: Directive 40, section 4.6.
l gas wells: within the first three months of production (one well per pool per section)
l oil wells: prior to any sales production (one well per pool per quarter section)
l Analyze a static gradient and generate and submit a GRD.PAS file using IHS
ValiData.
To meet the initial pressure requirement with a static gradient, one of the following criteria
must be met:
l the pressure must not build more than 2 kPa/hr over a period greater than 2 hours
while the gauges are on bottom.
l OR, the gauge was on bottom (with well shut-in) for a minimum of 14 days.
For more information, see Meeting the AER Initial Pressure Requirements.
l oil wells: not required, however, if an IPR is available, the test data and interpretation
should be submitted to the AER.
For more information, see Meeting the AER Initial Deliverability Requirements.
l oil wells: not required, however, AER rules state that if you have cleanup or flow data
in PRD.PAS format it shall be submitted.
l If the well is producing inline, the TRG-PAS file generated to meet the initial deliv-
erability requirements will also meet the initial production requirement.
In a petroleum reservoir, the rock is usually not fully saturated with a single phase fluid. Generally,
saturations in the reservoir rock consist of different amounts of gas, oil, and water. These saturations
change the effective permeability of the rock.
Permeability can be measured in a laboratory from core analysis. While this is sometimes done with
the core in its native or restored state, the more common method is to clean and dry the core, and
measure its absolute permeability (usually to air or nitrogen, but the same value would be obtained if
water were used instead). Permeability can also be determined by Pressure Transient Analysis
(PTA). However, it must be remembered that the permeability determined from analysis is the in-situ
effective permeability of the primary reservoir fluid type, and not the absolute permeability.
The in-situ effective permeability is usually significantly less than the absolute core-derived per-
meability by a factor ranging from 2 to 200, depending on the reservoir. In addition, when a per-
meability is determined from a pressure transient test, it reflects the average permeability of the
reservoir within the radius of investigation of the test (often several hundred feet). This is in contrast
to a core measurement that represents only a few inches of the reservoir.
In PTA, effective permeability (in-situ) can be determined either by semi-log (radial) analysis, or mod-
eling (matching) pressure data. The results of these two techniques should be consistent; any incon-
sistencies should be accounted for by reviewing and modifying the information and data provided,
and/or modifying the analysis/model accordingly.
When gas, oil, and water are being produced during a test, the effective permeability for each phase
can be estimated by assuming that only that fluid phase was flowing.
Note: In most of the pressure transient equations, the permeability term occurs as a mobility (k / µ)
or a transmissivity (kh /µ) term.
For example, if a well flows at a 1,000 stbbls/d with a flowing (sandface) pressure of 1500 psi and at
an average reservoir pressure of 2,000 psi, then the productivity index is:
PI = 1000 / (2000 - 1500) = 2 stbbls/d / psi
After the PI is known, the equation can be rearranged to determine the deliverability rate as follows:
This equation is only valid for a well in an undersaturated reservoir. For wells in saturated reservoirs,
or for gas wells, the relationship is not as straightforward, and the simple relationship described
above does not apply. In these situations, either an IPR (for oil) or AOF (for gas) analysis should be
performed.
The recombination calculation takes the volume of condensate, vaporizes it, and adds it the gas
volume to obtain a single phase rich gas as it exists in the reservoir. A recombined gas gravity is also
calculated for use in the analytical models.
Given the specific gravity of measured separator gas, the CGR, and the specific gravity of the con-
densate, the recombined gas gravity can be calculated using the following equation:
where:
The recombined gas rate factor (RGRF) is defined as recombined gas rate divided by the measured
gas rate:
In a single-phase system such as a dry gas or an under-saturated oil reservoir, the effective per-
meability of flow of the mobile fluid through the reservoir varies little during production because the
fluid saturations do not change. However, when more than one phase is mobile, the effective per-
meability to each mobile phase changes as the saturations of the fluids change in the reservoir.
In a two-phase system, the fluids might consist of oil and water or oil and gas. In a three-phase sys-
tem, all three fluid phases occur. Each fluid, as it flows through the porous media, interferes with the
fluids because capillary forces exist that reduce the flow rate of each individual phase in a non-linear
fashion. Consequently, the sum of the relative permeability of each phase is always less than one.
Wettability
Wettability in a reservoir is a measurement of the ability of a fluid to coat the rock surface. Wettability
and heterogeneity have a significant impact on the shape of the relative permeability curves. The
wetting fluid relative permeability curve is concave upwards whereas the non-wetting fluid has an “s”
shape. In the case where there is no interfacial tension between the fluid phases, the relative per-
meability curves simplify to straight lines between the endpoints.
l Swc is the connate or irreducible water saturation. This is the water saturation
below which water is not mobile because of capillary forces. The relative per-
meability of water at water saturations below Swc is zero.
l Sorw is the residual oil saturation or critical oil saturation. This is the oil saturation
below which the oil is immobile, that is, its relative permeability is zero.
l Sgc is the critical gas saturation. This is the minimum saturation for gas to
become mobile.
l Sorg is the residual oil saturation to gas. This is the immobile oil when gas is the
displacing fluid.
l krgc is the relative permeability of the gas at the residual oil saturation.
Corey
This model assumes the wetting and non-wetting phase-relative permeabilities to be independent of
the saturations of the other phases and requires only a single suite of gas/oil-relative permeability
data.
Honarpour
Developed using data from oil and gas fields in the continental US, Alaska, Canada, Libya, Iran,
Argentina and the United Arab Republic.
Sandstone
Water Wet
Intermediate Wet
Limestone
Water Wet
Intermediate Wet
Any Wettability
Three-Phase Correlations
Three-phase relative permeability can be generated from the two-phase relative permeability curves
of the oil-water system and the relative permeability curves of the gas-oil system. The two-phase
curves represent the end curves when either the gas saturation or water saturation equals zero.
Stone I is widely used in the industry as the benchmark for oil simulation. It is a better predictor than
Stone 2 in low oil saturation regions, is more appropriate for water-wet systems, and is not suited for
intermediate wet systems.
Stone II
Stone's Model II is a modified version of Stone I. It is a better predictor than Stone 1 in high-oil sat-
uration regions. It is more appropriate for water-wet systems and is not suited for intermediate wet
systems.
Time to Stabilization is defined as the time it takes for the radius of investigation to reach all of the
reservoir boundaries.
References
"Theory and Practice of the Testing of Gas Wells, Third Edition", Energy Resources Conservation
Board (ERCB), pp. 3-34 to 3-36, 1975.
The forecast gives the analyst an understanding of the reservoir and the economics of the well. For
gas or oil reservoirs, it can help illustrate and quantify the value of adding suction, improving skin, or
changing the design of the gathering system. Note that for oil reservoirs, the forecast assumes con-
stant compressibility, constant fluid saturations, and single-phase flow; therefore, it is only valid
above the bubble point pressure, and long-term forecasts may not be reliable.
The forecast below is a plot of rate and cumulative production versus time for a specified back pres-
sure. The calculations are done on a daily basis for the first month, and then on a monthly basis up to
the specified forecast time, and consist of sequential material balance and deliverability calculations.
Multiple forecast cases (up to three) can be performed as shown in the plot below.
A contract rate can also be specified to determine how long that rate can be maintained at the spe-
cified back-pressure as shown in the plot below.
l Start forecast at beginning of test — starts the forecast at the initial average
reservoir pressure calculated from the production data and model.
l Start forecast at end of test — starts the forecast at the final average reservoir
pressure calculated from the production data and model.
l Start forecast at t = — starts the forecast at the average reservoir pressure cal-
culated at the specified time in months from the beginning of the test using the
production data and model.
l qc — specifies the contract rate to use in the forecast. This rate is maintained as
long as possible at the specified flowing pressure. Rate decline occurs when this
rate can no longer be supported.
l Time — specifies the total time or duration of the forecast. The default is 12
months.
l Calculate — click this button to generate the forecast based on the current fore-
cast parameters.
l Wellbore — specifies that the flowing pressures used are at wellhead con-
ditions. When using this option, the wellbore parameters also need to be spe-
cified to convert the wellhead flowing pressure to sandface conditions. Note that
this option requires more computational time as an iterative process is required
to convert the wellhead pressure to sandface. This is due to the fact that the con-
version is dependent on the rate being forecast, and the rate being forecast is
dependent on the sandface pressure. Thus, an iterative process is required to
solve this scenario.
l Sandface — specifies that the flowing pressures used are at sandface con-
ditions.
l Edit Wellbore Parameters — click this button to open the Wellbore Properties
dialog box and modify the current wellbore configuration. This option is only avail-
able when the Wellhead option (see above) is active.
l pwf1 — specifies the case 1 flowing pressure to run the forecast at. This value
must be less than the reservoir pressure (pR).
l s1 — specifies the case 1 skin to run the forecast at. The default value is the skin
specified in the model.
l pwf2 — specifies the case 2 flowing pressure to run the forecast at. This value
must be less than the reservoir pressure (pR).
l s2 — specifies the case 2 skin to run the forecast at. The default value is the skin
specified in the model.
l pwf3 — specifies the case 3 flowing pressure to run the forecast at. This value
must be less than the reservoir pressure (pR).
l s3 — specifies the case 3 skin to run the forecast at. The default value is the skin
specified in the model.
l t2 — specifies the time in months to make either a flowing pressure (add com-
pression), or skin (stimulate well) change. This can be applied to all three cases.
l t3 — specifies the time in months to make either a flowing pressure (add com-
pression), or skin (stimulate well) change. This can be applied to all three cases.
For a multi-rate test (AOF LIT analysis), s' and D have different effects and can be determined inde-
pendently.
Note: The turbulence factor cannot be a negative value. If it is, it reflects errors in the data or ana-
lysis and it should be set to zero.
Distance to Boundary
The distance to boundary is the distance between the well and the boundary. The boundary may be
a no-flow or constant-pressure boundary, and may be circular or linear as shown below.
If a well is located near a single boundary, the reservoir is considered to be infinite acting in all the
other directions.
No-Flow Boundary
A no-flow boundary is a boundary that does not allow flow through it. This kind of boundary usually
occurs in reservoirs with sealing faults, or is created between producing wells that are equally
spaced and producing at the same rate. Pseudo-steady state flow signifies that all no-flow bound-
aries have been reached.
Observation Well
In modeling, pressure can also be calculated at a different location or observation point away from
the active well. This point is known as an observation well. The location of the observation well is
measured with respect to the active well location (e.g. Xw and Yw) and has the coordinates ∆Xo, ∆Yo,
Zo for rectangular reservoirs or ro for cylindrical reservoirs. This coordinate system for rectangular
reservoirs is shown in the diagram below.
Well Location
The location of the active (producing or injecting) well for rectangular- shaped reservoirs is specified
using an x, y, z coordinate system, as outlined below. The distance is represented from the origin in
the lower-left hand corner (Xw, Yw and Zw).
Note: The Zw distance is only applicable for a horizontal well and represents the distance or height
from the origin (base of the formation) to the center of the horizontal leg of the wellbore. For
cylindrical-shaped reservoirs, the well is always assumed to be in the center.
n Formation properties
n Fracture properties
n Gas properties
n Oil properties
n Reservoir temperature
n Water properties
It is of the same order of magnitude as the oil compressibility (co) or the water compressibility (cw),
approximately 10-6 1 / psi.
Thus, in the absence of gas saturation (Sg), the formation compressibility is relatively significant and
must not be ignored.
However, when there is gas present in the pores, the gas compressibility (cg) is one or two orders of
magnitude higher, and the formation compressibility is relatively small and may be ignored.
The formation compressibility may be measured in the laboratory, but it is usually derived from cor-
relations. It is obvious that the nature of the formation, its degree of consolidation, its stress field, etc.
all affect the formation compressibility, yet the correlations that exist do not take all of these factors
into account. They are very simplified and relate the formation compressibility to the total porosity
(φ). The correlation used in this software is applicable to sandstones and carbonates. It was derived,
from laboratory measurements, by Hall and published in the Transactions, AIME (1953) 198, p. 309.
If the gas reservoir is part of a gas cap in a saturated oil reservoir, the space may contain gas, oil,
and water. Thus the gas saturation is rarely 100% but varies from 30% to 90%. Gas saturation is
used directly in the calculation of the reserves and also in the calculation of total compressibility
(ct). Note that the gas saturation has a significant effect on the total compressibility (ct) since the gas
compressibility is much higher than the other fluid compressibilities (oil (co) and water (cw)).
From this definition, the total compressibility takes into account the compressibility and saturation of
gas (Sg), oil (So), and water (Sw) as well as the (rock) formation compressibility (cf).
Note that the gas saturation (Sg) has a significant effect on the total compressibility (ct) since the gas
compressibility (cg) is much higher than the other fluid compressibilities (oil (co) and water (cw)).
Thus, if there is any gas present in the reservoir, the total compressibility is dominated by the gas
compressibility (cg) component.
The total compressibility has only a small effect during pressure transient flow analysis (it occurs
within the log term in the skin calculation equation for radial flow), but it has a significant effect during
pseudo-steady state (PSS) flow as it directly affects the total energy content of the reservoir.
In pressure transient analysis, the total compressibility is assumed constant, except for the case of
gas pseudo-time in which the total compressibility is allowed to vary with pressure. Note that when
forecasting production rates in a gas reservoir, if the forecast duration is much longer than the time
to pseudo-steady state, significant errors in the late-time forecast may result (the rate will be under-
predicted) when assuming a constant total compressibility.
When using this parameter for pressure transient analysis, it is assumed that the pore space is filled
with a single-phase fluid (either gas, oil, or water), such that the total porosity equals the porosity of
that specific fluid. In reality, this is rarely the case, and the total porosity is usually filled with a mix-
ture of gas, oil, and water.
Total porosity is used to calculate formation compressibility (cf) and thus affects the total com-
pressibility (ct). The total porosity has only a small effect during pressure transient flow analysis (it
occurs within the log term in the skin calculation equation for radial flow), but it has a significant
effect during PSS flow as it directly affects the reserves contained in the reservoir.
Water saturation may range from 10% to 50% for an oil or gas reservoir, and it will be 100% for an
aquifer. The water saturation obviously affects the oil and gas reserves, and is also used in the cal-
culation of total compressibility (ct).
A large value of fracture flow capacity (>10,000 md ft) represents an infinite conductivity fracture,
and yields a linear fracture flow response on the derivative. A small value of the fracture flow capa-
city (<10,000 md ft) represents a finite conductivity fracture and may yield a bilinear fracture flow
response on the derivative. When the value of fracture flow capacity is divided by the product of form-
ation permeability (k) and fracture half-length (Xf), the result is known as the dimensionless fracture
conductivity (FCD) defined as:
Note: This dimensionless form is a more common measure of fracture conductivity that is found in
the literature.
The fracture half-length depends on the size of the fracture treatment and varies from a few feet to a
few hundred feet. In pressure transient analysis, it is estimated from the linear fracture flow analysis.
Based on the gas properties, Condensate Gas Ratio (CGR), Condensate Gravity (gcnd), Separator
Temperature (Tsep), and Separator Pressure (psep), a Recombined Gas Rate Factor (RGRF) and
Recombined Gas Gravity (gr) can be calculated. The measured gas rates are multiplied by the
RGRF to get the total sandface rate. The gas at sandface is richer than the gas at surface, so the
higher Recombined Gas Gravity is used for analysis. In WellTest, you can either use a constant
RGRF calculated in Gas Properties, or calculate a different RGRF for each rate based on the current
CGR on the Production Editor tab.
Note: The condensate is separated from the gas in the separator and measured in the stock tank.
However, the fluid in the stock tank has been flashed and will no longer contain components
< C5.
The gas compressibility is a very strong function of pressure, and increases as the pressure
decreases. Mathematically, it can be expressed as:
Where p is the specified pressure and z is the gas compressibility factor (z) at that pressure. Thus,
the magnitude of gas compressibility (cg) is of the order of 1 / p. Moreover, when the pressure draw-
down at the wellbore is large, the difference in compressibility values at the initial pressure (pi) and at
Note: It is a very strong function of pressure, and a weak function of temperature and gas com-
position.
The gas gravity affects the calculations of gas viscosity (µg), compressibility (cg), compressibility
factor (z), and solution gas oil ratio (Rs). The effect of gas gravity on these properties is not sig-
nificant, so it is not detrimental if the gas gravity is not precisely known, using the default or an aver-
age value is sufficient in most cases.
For sour gases, this correlation is preferred to the Lee, Gonzalez, and Eakin formulation which does
not account for H2S, CO2, and N2. Gas viscosity is used in numerous equations, most notably in the
definitions of pseudo-pressure (ψ) and pseudo-time (ta). Typically, gas viscosity is in the range of
0.015 to 0.03 cp or 15 to 30 micro-Pa.s.
In the analysis of gas wells, the reservoir temperature is used directly in the pressure transient ana-
lysis calculations, as well as calculating all the gas fluid properties such as gas formation volume
factor (Bg), gas viscosity (µg), and gas compressibility (cg). For gas wells, the reservoir temperature
has only a weak effect on the pressure transient analysis calculations.
In the analysis of oil wells, the reservoir temperature does not enter into the pressure transient ana-
lysis equations, directly, but rather It is used in calculating the oil fluid properties such as oil form-
ation volume factor (Bo), oil viscosity (µo), oil compressibility (co), and solution gas-oil ratio (Rs). It
also enters into the calculation of the total reservoir fluid rate (qtBt) through the gas formation
volume factor (Bg), and in the calculation of total compressibility (ct) through the calculation of gas
compressibility (cg).
The value of water compressibility can be obtained from laboratory Pressure, Volume, Temperature
(PVT) measurements, or determined from correlations. The magnitude is approximately between
1.0 x 10-6 and 9.0 x 10-6. It is a weak function of pressure, temperature, and salinity.
Water formation volume factor can be measured in the laboratory, or determined from correlations.
Under most conditions it has a value of approximately 1.0. It is a very weak function of pressure, tem-
perature, and salinity.
Water specific gravity has no effect on calculated properties such as water compressibility (cw), form-
ation volume factor (Bw), and viscosity (µw). It is used, however, in the wellbore pressure drop cal-
culations when converting pressures from wellhead to sandface.
n Number of Segments
l Simplicity
l Cost effectiveness – not having to pull pumps and rods in order to land sub-
surface recorders
l Less disturbance to the formation – not having to kill the well, or damage the form-
ation with kill fluid
Because of these significant advantages, AWS measurements are widely used throughout the oil
industry as a means of obtaining subsurface pressures.
Quality
The best way to obtain a subsurface pressure is to measure it directly at the sandface. Subsurface
pressure gauges usually yield acceptable results, except in the case of a total failure of the gauge.
By contrast, you can use AWS as an indirect way of determining the pressure at the subsurface, and
the quality of the results can be quite variable. In many cases, the results are quite acceptable
(within the objectives of the test). However, in many other cases, the results from AWS are grossly in
error for a number of reasons.
In other words, unless special precautions are taken, it is very easy to obtain the wrong answer with
AWS.
Fekete’s experience and expertise in AWS surveys is summarized in this abstract. It provides an out-
line of procedures to assist you in your AWS sandface pressure surveys and calculations, and con-
siders both practical and theoretical requirements for a successful survey.
One of the best ways of achieving this goal is to have a survey coordinator who is responsible for
planning, scheduling, and coordinating personnel and equipment needed to conduct the field work
and analysis of acquired data. The survey coordinator must have a thorough understanding of the
field procedures involved with conducting an AWS survey, as well as test interpretation skills to
In wells where wax problems occur, a hot oil (or treated water) program should be performed long
enough in advance to allow the well to return to normal pumping conditions prior to the survey. Also
at this time, the wellhead should be inspected for any leaks, and casing valves serviced.
One day prior to conducting a foam depression test, flow rate measurements should be conducted
and compared to historical data. Experience suggests that two measurements, of 12 hours each, be
conducted to verify measurements. Reid measuring devices should be recently calibrated to ensure
accuracy.
The most recent and complete tubing tally and wellbore configuration should be obtained from the
well files. Remember, the collars of the tubing are what we are counting. The degree of error here is
reflected in the accuracy of the sandface pressure, even when employing the Acoustic Velocity
method, which must be calibrated using the count of the collars.
In order to determine proper parameters for the analysis of the data, the most recent Pressure,
Volume, Temperature (PVT), petrophysics, and geology data should be acquired for the pool.
Sampling and analysis of the oil, water, and gas should be conducted by qualified personnel only.
l Where possible, sampling should be conducted from a separator, noting the tem-
perature and pressure at which the sampling was conducted.
l Otherwise, stock tank oil samples, casing gas samples, and representative water
samples should be obtained and analyzed.
l All inconsistencies should be resolved prior to the pressure data being analyzed.
Required personnel should be notified well in advance of the survey. If a field survey is being con-
ducted, ensure enough personnel are available to conduct the survey to its completion. Don’t forget
support staff (i.e., field operators, chart readers, etc.). Remember, pressure surveys may require 24-
hour, 7 days a week monitoring to return survey wells to production as soon as possible.
Field Supervision
Supervision of the data gathering in the field is extremely important. As the saying goes "garbage in -
garbage out". Incorrect field-data-gathering techniques result in inaccuracies of the final pressure
calculations.
Wellbore Dynamics
To better understand and quantify the results of an AWS survey, experience has shown that well-
bore dynamics must be considered during a pressure buildup monitored by AWS.
Where
pcs is the shut-in casing pressure measured by either a deadweight, or a digital pressure
transducer. A normal pressure gauge is not recommended due to normal "wear and tear" of
field operations.
Note: Prior to commencing and throughout the survey, ensure that there are no leaks at the well-
head by inspecting it visually and listening for small leaks.
∆pGAS is the pressure exerted by the gas column in the annulus. This is easily calculated
using the Cullender and Smith method with very acceptable accuracy.
∆pLIQUID is the pressure exerted by the liquid column in the annulus. This is where the dif-
ficulty comes and consideration of wellbore dynamics is essential.
l Dual Channel – These machines produce a strip chart indicating tubing collars
on the bottom and liquid level on the top. This provides easy verification of the
liquid level at a glance.
l Single Channel – The single channel fluid level machine has the capability of
providing tubing collar, or fluid level strip charts individually. In order to compare
both charts, they must be placed side-by-side and verified.
l Acoustic Velocity – This method of determining fluid levels is based on the velo-
city of sound in the gas, and the time it takes for the sound wave to be reflected
off the gas-liquid interface. Acoustic velocities are sensitive to gas composition,
temperature and pressure, and thus, may vary during the test.
All of these methods may be affected by the presence of a "foam" column, as the acoustic reflection
indicates the top of the "foam" column rather than the top of the "liquid" column.
l The lower the liquid level is, the smaller is the error in the calculated pwf.
However, care should be taken not to depress the fluid lower than one joint
above the pump, so as not to have gas bypass occur in the pump.
A foam depression test can vary in time to complete; therefore a review of past tests assists you in
planning your AWS survey. The foam depression test is usually performed immediately preceding
the buildup test.
l PVT. Correlations
l Static Gradients
Note: In many properly supervised AWS tests, the margin of error calculated is within very accept-
able limits.
Conclusions
In order to have a successful survey, you need:
l Proper design
l Knowledgeable supervision
A large value of dimensionless wellbore storage is associated with a large wellbore volume (Vw). In
practice, a range of values from 500 to 10,000 has been observed. Occasionally, a value much
greater than 10,000 is observed when matching actual test data against a type curve. In such cases,
it is considered that some part of the formation (high permeability or fracture) is in communication
with the wellbore, and is acting as an extension to it.
In a horizontal or slant well, the dimensionless wellbore storage constant (CD) is sometimes
expressed in terms of the effective wellbore half-length (Le/ 2). In such a case, the value for wellbore
radius (rw) in the above equation is replaced by the effective wellbore half-length, which makes CD a
very small number.
Wellbore storage is normally assumed to be constant during a test, and in practice, this assumption
is often reasonable. However, there are numerous situations where wellbore storage is not
constant. This changing wellbore storage may be caused by a changing wellbore fluid com-
pressibility, by phase redistribution, or by a change in the type of storage from a changing liquid level
to a liquid filled wellbore. Changing wellbore storage is accounted for using a modified form of the
dimensionless wellbore storage definition, which includes additional dimensionless parameters,
apparent dimensionless storage (CaD), and a storage pressure parameter (CpD).
The wellbore storage constant and dimensionless wellbore storage constant are calculated from the
early time data, by performing an afterflow or wellbore storage analysis.
Note: A value of 100 segments is recommended and handles most situations adequately.
When the wellbore fluid is highly compressible (e.g., gas), the wellbore can store more fluid with
increasing pressure. Subsequently, wellbore storage is higher, and the effects are seen for a longer
period of time.
Theoretical development of the wellbore storage equations assumes a single-phase fluid in the well-
bore with a constant compressibility. The value of the compressibility used corresponds to that of the
primary, or dominant fluid type in the wellbore. In reality, this is not the case, since the wellbore often
contains a mixture of oil, gas, and water. Thus, the fluid compressibility is not constant, but is a func-
tion of pressure. As a consequence, any value used for wellbore fluid compressibility is only an estim-
ate.
This radius is used to determine the sandface area (2rwh), which represents the area through which
all the produced reservoir fluids must flow. The sandface area has a significant effect on the flowing
wellbore pressure, and thus, the wellbore radius is found in many of the flow equations used in Pres-
sure Transient Analysis (PTA).
In real life, the area of contact between the wellbore and the formation is rarely cylindrical as seen
above in Figure II. It depends on the perforations (density, phasing, effectiveness, etc.), and is also
affected by the type of perforating gun, casing, cement, etc. Thus, a true wellbore radius does not
exist (except for open hole completions), and the wellbore radius used in the PTA equations is an
approximation at best. A reasonable value to use can be the drill bit radius, or the outside diameter
of the casing. The default value used in the software is 0.3 ft (0.091 m).
For non-conventional analysis and modeling (Slug / PITA / CCT), wellbore volume is used in cal-
culating the influx rate, or the afterflow rate into the wellbore at the sandface.
l Area of the annulus when the tubing is closed off, either by a valve or hydrostatic
head in the tubing
l Area of the tubing and the annulus when the liquid level is rising in both
The value of the volumetric capacity typically varies from approximately 0.02 to 1 ft2 (0.002 to 0.1
m2), and is determined by the tubular configuration of the wellbore. It is used in non-conventional
analyses (e.g., Slug, PITA, CCT) when there is a changing liquid level in the wellbore.
Vertical Model
l Phase redistribution
l Change in the type of storage from a changing liquid level to a liquid filled well-
bore
The phenomenon of phase redistribution occurs in a well that is shut in at the surface with gas and
liquid flowing simultaneously into the tubing. In these situations, the gravity effects cause the liquid
to fall and the gas to rise to the surface. Because of the very low compressibility of liquid and no
extra room for gas to expand in a closed chamber, the redistribution of phases causes a net increase
in the wellbore pressure. When the phenomenon is present in a buildup test, the extra pressure
surge in the wellbore is relieved through the formation. Eventually, equilibrium will be attained
between the wellbore pressure and the formation pressure adjacent to the wellbore. However, at
early times, the wellbore pressure may exceed the formation pressure causing an anomalous hump
in the buildup pressure, which cannot be analyzed in a conventional way with only the dimensionless
wellbore storage constant (CD). In order to deal with phase redistribution, two models have been pro-
posed by Fair (1981) and Hegeman et al. (1993), which introduces two additional dimensionless well-
bore constants, apparent storage (CaD) and the pressure parameter (CpD).
Fair considered the exponential expression for the dimensionless anomalous pressure (ppD) rise
as:
Later, Hegeman et al. showed that the negative CpD values in the Fair model can be used for
buildup data that has an anomalous pressure decrease. Therefore, for these wells, they argued that
using an error function to model the anomalous pressure may allow for better modeling of field data
with increasing or decreasing storage. Thus, Hegeman et al. proposed that:
The following plots illustrate the effects of the three dimensionless storage parameters on the dimen-
sionless type curves.
"The Primary Pressure Derivative (PPD) – A New Diagnostic Tool in Well Test Interpretation", L. Mat-
tar and K. Zaoral, JCPT (April 1992) Volume 31, No. 4, 63 - 70.
"Pressure Buildup Analysis with Wellbore Phase Redistribution", W.B. Fair Jr., SPEJ (April 1981)
259 - 270.
"Well-Test Analysis with Changing Wellbore Storage", P.S. Hegeman, D.L. Hallford, and J.A.
Joseph, SPEFE (September 1993) 201 - 207.
One way to account for a variable permeability over time is to modify the definition of pseudo-pres-
sure and pseudo-time.
where
The modified PSS equation, which incorporates a variable permeability with pressure, is as follows:
1. Depletion Part:
(1)
2. Inflow Part:
(2)
The equation, 1. Depletion Part, is a material-balance equation, and thus, has no dependence on
the mobility terms k and m. The mobility terms only enter into the inflow, 2nd part of the PSS equa-
tion. When the PSS equation is coupled, however, the k/m term must be included in pseudo-time.
It helps us assess the suitability of the model, characterize the reservoir, and calibrate the model
prior to forecasting.
Analytical Models
Analytical models have relatively short calculation times, so you can history match with these mod-
els quickly. More than that, you can even match the model automatically using APE.
The main disadvantage for these models is that they can only simulate single-phase flow – either
liquid or gas, but not gas and liquid flowing together. Additionally, analytical models do not fully
account for changing fluid properties with pressure. In case of liquid flow, fluid properties are
assumed to be constant, and in case of gas flow, the change of fluid properties with pressure is
accounted for by using pseudo-pressure and pseudo-time.
Unfortunately, pseudo-time is not an exact transformation. As a result, changing gas properties are
only accounted for to a certain extent. Therefore, in cases where pressure varies significantly across
the reservoir (e.g., pays with low permeability), long-term forecasts for gas analytical models may be
inaccurate. There is no known way to fully account for the change in fluid properties when using ana-
lytical models.
Numerical Models
With numerical models, you can model multiphase flow, and account for changing properties of each
phase, and the interaction between phases.
The main disadvantage for these models is their long computation time. With the numerical mod-
el,the reservoir is divided into a number of cells, and then the flow is modeled simultaneously in all
cells. To be able to assume constant pressure and constant fluid properties within each cell, the size
of each cell needs to be relatively small. Therefore, the number of cells required for an accurate solu-
tion becomes large, which results in a long computation time. Note that decreasing the number of
cells to reduce computation time may significantly affect calculation results.
Performing a manual history match for numerical models is time consuming. In addition, an auto-
matic history match for numerical models is complicated, and thus not widely available in com-
mercial software.
Equation 1
where:
φ = porosity
ρ = density
Mass flow from cell i to the adjacent cell j during timestep n+1 can be calculated as:
Equation 2
In equation 2:
μ = viscosity
Mass produced at the well during timestep n+1 (if well penetrates cell i) can be cal-
culated as:
Equation 3
WI = wellbore index (as per Peaceman 1978 or Babu and Oden, 1989)
Equation 4
Using equations 1, 2, and 3, and bringing all the terms to the left, we can rewrite equation 4
as:
Equation 5
Assuming that the pressure distribution has been calculated for timesteps 1, 2, ... n, to calculate pres-
sure distribution at timestep n+1, we formulate material balance (equation 5) for each cell, and solve
all these equations simultaneously. Unknowns of the system are for each cell i; therefore, the
number of equations is the same as the number of unknowns, and the system can be solved.
Equation 6
where:
If the modeled fluid is liquid, its properties do not change much with pressure; therefore, pressure in
any cell can be used.
In the case of gas flow, properties significantly vary with pressure; therefore, selecting the correct
way of estimating properties becomes important.
In classical numerical simulation (including numerical models in Harmony), properties are estimated
at the pressure of the upstream cell (i.e., the cell with the higher pressure). For such a simulation to
be accurate, you should ensure that properties in the adjacent cells are close. This could be
achieved by having small grid cells. Unfortunately, having smaller grid cells results in longer com-
putation times.
With the hybrid model, we use the pseudo-pressure formulation to estimate fluid properties.
Note: The definition of hybrid model pseudo-pressure is very similar to a traditional pseudo-pres-
Equation 8
Therefore, the hybrid model pseudo-pressure and traditional pseudo-pressure are different by a con-
stant factor.
Equation 9
Equation 10
Integrating both parts from 0 to L with respect to x, and using the left-side of the equation is a con-
stant with respect to x, we get the following equation:
Equation 11
Note: We used the definition given in equation 7 for the last transformation.
To summarize, when fluid properties are changing with pressure, the mass flow rate through the
cross-section can be calculated as:
Equation 12
Therefore, equation 2 can be modified to account for changing properties. The mass flow from cell i
to cell j during timestep n+1 is calculated as:
Equation 13
Equation 3 can be modified in a similar fashion. Mass produced at the well during timestep n+1 is cal-
culated as:
As a result, the material balance equation for each cell (equation 5) is modified to:
Assuming that the pressure distribution has been calculated for timesteps 1, 2, ... n, to calculate pres-
sure distribution at timestep n+1, we formulate a modified material balance (equation 14) for each
cell, and solve all these equations simultaneously. Unknowns of the system are for each cell i;
therefore, the number of equations is the same as the number of unknowns, and the system can be
solved.
Note: Well constraints are treated in a similar way to general numerical model formulations.
1. Smaller number of grid cells: In classical numerical simulation, grid cells have to be small
enough to consider constant fluid properties within each cell. Pseudo-pressure formulation
accounts for a variation of fluid properties between the centers of two adjacent cells; there-
fore, it is possible to have larger grid cells. By having a smaller number of grid cells, cal-
culation speed increases.
2. Faster solution for a non-linear system: While performing numerical modeling, equation
5 is solved at each timestep. This system of equations is non-linear; therefore, it is solved
iteratively, using the Newton-Raphson method. This method involves calculating derivatives
of each matrix element with respect to each unknown. Calculating these derivatives is faster
for pseudo-pressure formulation (equation 14), because we have to deal with one fluid prop-
erty function ( ) as opposed to a combination of three functions ( )
for traditional formulation. More than that, due to the integral nature of , its derivative
is calculated easily.
3. Symmetry: When it comes to solving systems of equations, there are faster algorithms avail-
able for the case when the matrix of the system is symmetrical. Therefore, another advant-
age of using the formulation given in equation 14 over equation 5, is that equation 14 forms
a symmetrical matrix, while estimating properties at the upstream cell for equation 5 results
in an asymmetrical matrix.
The underlying assumption of the analytical models for production data analysis is single phase flow
in the reservoir. In order to accommodate multiple flowing phases, the model must be able to handle
changing fluid saturations and relative permeabilities. Since these phenomena are highly non-linear,
analytical solutions are very difficult to obtain and use. Thus, numerical models are generally used to
provide solutions for the multiphase flow problem. (The numerical engine used in the software is
based on a general purpose black-oil simulator.) Numerical models can be created with less sim-
plifying assumptions for reservoir properties than analytical models. The reservoir heterogeneity,
mass transfer between phases, and the flow mechanisms can be incorporated rigorously.
Numerical models solve the nonlinear partial-differential equations (PDEs) describing fluid flow
through porous media with numerical methods. Numerical methods are the process of discretizing
the PDEs into algebraic equations, and solving those algebraic equations to obtain the
solutions. These solutions that represent the reservoir behaviour are the values of pressure and
phase saturation at discrete points in the reservoir and at discrete times.
The advantages of the numerical method approach are that the reservoir heterogeneity, mass trans-
fer between phases, and forces / mechanisms responsible for flow can be adequately taken into con-
sideration. For instance, multiphase flow, capillary and gravity forces, spatial variations of rock
properties, fluid properties, and relative permeability characteristics can be represented accurately
in a numerical model. In general, analytical methods provide exact solutions to simplified problems,
while numerical methods yield approximate solutions to the exact problems. One consequence of
this is that the level of detail and time required to define a numerical model is more than its equi-
valent analytical model.
Gas-condensate and volatile oil systems contain gas that may have non-negligible amounts of vapor-
ized liquid hydrocarbons, and this may have a significant impact on fluid behaviour.
Common black oil numerical formulations do not consider changes in the liquid hydrocarbon content
of the gas phase. In the modified black oil model, two new properties are added to account for the
liquid contained in gas:
To explain these new properties, consider the figure below, which illustrates the distribution between
the phases at given reservoir conditions. Excluding water, in both black oil and modified black oil
modeling, it is assumed that there are two components (separator gas (G) and stock tank oil (N)),
and two phases (gas phase (g) and oil phase (o)); and each of these components can exist in either
phase. The amount of the produced gas at the separator is the summation of gas components that
come from the gas phase (Gg), and the gas component that comes from the oil phase (dissolved
gas, Go). Likewise, the stock tank oil (N) comes from both the oil phase (No) and the gas phase
(vaporized oil, Ng):
Ng is neglected in the conventional black oil numerical formulation, but is honoured in the modified
black oil numerical formulation.
Dry gas formation volume factor (Bgd) is another distinct modified black oil property. The wet gas
formation volume factor (Bg) is defined as the ratio of the gas phase volume at a given pressure and
temperature (Vg) to the equivalent volume of the whole gas phase at standard conditions:
The dry gas formation volume factor (Bgd) is the ratio between the volume of the gas phase at given
conditions (Vg) to the volume of its gas component (Gg) at standard conditions:
In conventional black oil applications, Ng is zero, which causes Bgdand Bgto become the same.
However, for modified black oil applications, we need to differentiate between Bgand Bgd. Although
where:
The Gas Condensate and Volatile Oil Numerical models use modified black oil properties; therefore
the condensate drop-out is modeled correctly.
Gas Condensate Numerical models use the Gas Rate and Condensate Rate columns from the Pro-
duction Editor tab as an input for production data. Condensate properties are set in the Properties
tab’s Oil section. The Gas Condensate Numerical model reports dry gas rates and condensate
rates.
Volatile Oil Numerical models use the Gas Rate and Oil Rate columns from the Production Editor tab
as an input for production data. The Volatile Oil Numerical model reports dry gas rates and oil rates.
Before a Gas Condensate Numerical model or Volatile Oil Numerical model starts a simulation,
WellTest verifies that the properties entered into the simulators are physically meaningful. The
model performs multiple property consistency checks, and if these checks do not meet the estab-
lished criteria, an error or warning message is displayed. The checking criteria are as follows:
co > 0
cg > 0
ρo ≥ ρ g
1 / Rv ≥ Rs
Bgd / Rv ≥ Bo
Bo / Rs ≥ Bgd
μo ≥ μg
At p = pbp or pdew:
Simulation can stop: Although the numerical engine is robust, there are several cases that cause
the execution of the model to fail. (This is true of any gridded simulator.) This usually happens when
the model shortens its time steps too much, such as when there are rapid operational changes in
rate or pressure. These issues can usually be overcome by smoothing the data, and eliminating dra-
matic shifts in rate and/or pressure.
Numerical Errors: There are several types of errors associated with the numerical solution. The
first type is called truncation error, which is caused by a truncated Taylor series expansion replacing
the spatial derivative and time derivative. The order of truncation error is proportional to ∆x (grid size)
There is another kind of approximation that can result in one more type of error that is caused by the
well model incorporated in the numerical model. In order to obtain the accurate wellbore pressure or
flow rate, very fine grids around the wellbore are required. This is especially critical for compressible
fluids and low permeability reservoirs. The best gridding method around the well is cylindrical grids,
but they are only applicable for single-well modeling.
Gridding: Fine gridding is required to obtain the pressure or flow rate at the wellbore. This is espe-
cially true for compressible fluids and low permeability reservoirs. Although cylindrical grids provide
the best solution around the wellbore, they can't be used for multi-well modeling. Cartesian gridding
is the most common grid type used in the industry. The Peaceman Wellbore Index model for
Cartesian grids is used by most general numerical simulators, but can overestimate the wellbore
pressure, if rate is the constraint, and overestimate rate, if pressure is the constraint. It also produces
artificial wellbore storage (afterflow) effects during the early transient period. Our numerical models
use the Peaceman Wellbore Index model for vertical wells, and the Babu and Odeh Wellbore Index
model for horizontal wells.
Shut-In Time
Horner Time
Superposition Time
Equivalent Time
Root Time
Complex time functions like superposition linear equivalent time and many others are simply derived
from the above equations. Subscripts not shown above, but used for linear and bilinear functions are
L and B, respectively.
The constant rate solution for bilinear fracture flow data is:
For gas:
Based on these definitions, bilinear fracture flow data appears as a straight line on a plot of pressure
versus quad root time (p vs. ∆t0.25). For a complex sequence of flow rates, the superposition bilin-
ear time function is used instead.
In the conventional definition of pseudo-time, the com-pressibility and viscosity terms are evaluated
at average reservoir pressure conditions. Clearly, the average res-ervoir pressure is only a function
of original-gas-in-place (OGIP) and cu-mulative production. During the transient flow period, before
any boundary effects are observed, the flow behaviour of two different-sized reservoirs should be
similar, and independent of the OGIP – they are both infinite-acting reservoirs. However using the
conventional definition of pseudo-time, the result would be that the producing rates would be dif-
ferent because their average reservoir pressures are different due to the different OGIPs.
If a well is producing under boundary-dominated conditions, the average reservoir pressure is a very
reasonable datum at which to establish fluid properties such as cg. However, if the well production is
still in transient flow and no reservoir boundaries have been observed, the average reservoir pres-
sure based on total reservoir volume is not an appropriate datum to use. Consequently, pseudo-time
can cause anoma-lous model responses under certain conditions.
Anderson and Mattar (2007) proposed that the average reservoir pressure used in the pseudo-time
calculation during the transient flow period should be calculated based on the gas-in-place of the
investigated volume at that time (as determined by the radius of investigation). This way, during tran-
sient flow, pseudo-time is independent of OGIP. As soon as the reservoir enters boundary-dom-
inated flow, the conventional definition of pseudo-time is automatically resumed, because at that
time the region of investigation is the whole reservoir.
This volume is adjusted for the effect of the reservoir boundaries, and the coalescence of regions of
influence caused by interference between wells.
For analysis, it should only be used for constant rate drawdown tests, or for buildup tests which have
been preceded by a long duration flow (greater than time to stabilization).
Note that this can also be defined as the flow time (tc) divided by the Horner time as follows:
This definition is derived for radial flow and thus works very well for analyzing wellbore storage and
radial flow regimes, but breaks down in heterogeneous reservoirs, or when boundaries come into
play. When analyzing bilinear or linear flow regimes, corresponding definitions of bilinear or linear
equivalent time are used. Equivalent time has little effect on early time data, but compresses late
time shut-in data significantly. The maximum value of equivalent time is flow time (tc) as shut-in time
(∆t) approaches infinity. Thus, no matter how long the well is shut-in, equivalent time cannot be
greater than tc. Thus, the behaviour of a buildup depends heavily on the flow time.
The preceding definition of equivalent time assumes a constant rate before shut-in. When the rate
has been variable, a complex superposition equivalent time function should be used.
Effective producing time is only valid if the final flow rate before shut-in is constant. If this is not the
case, superposition time, which is more rigorous but also much more computationally intensive,
must be used.
Effective producing time has been called a number of different names (some of them misnomers) in
literature. Some common aliases for effective producing time include:
l Flow Time
Superposition in time of the radial flow equation for the case of a buildup test is as follows:
The time function in this equation is called Horner time defined as:
By plotting pws versus Horner time on a semi-log (radial) plot, radial flow buildup data appears as a
straight line and can thus be analyzed to determine permeability (k) and apparent or total skin
(s'). Note that based on the above definition, Horner time decreases as shut-in time (delta time)
increases, and as shut-in time approaches infinity, Horner time approaches the limit of one. Due to
this inverse relationship with shut-in time, the time axis on the semi-log (radial) plot is plotted in
reverse so that increasing shut-in time is still represented from left-to-right.
Since Horner time is based on the radial flow equation, it should only be used for analyzing radial
flow. For linear or bilinear flow, linear or bilinear time functions should be used instead for analysis.
Horner time is valid only when the reservoir is infinite-acting and the rate prior to shut-in was con-
stant. When the rate has not been constant, superposition radial time should be used instead.
The constant rate solution for linear fracture flow data is:
For gas:
The constant rate solution for linear channel flow data is:
For gas:
Based on these definitions, linear flow data is displayed as a straight line on a plot of pressure
versus square root time (p vs. ∆t0.5). For a complex sequence of flow rates, the superposition linear
time function is used instead.
Note: Both time functions are bounded by 1.0 (at early shut-in time) signifying the time of closure,
and 0.0 (at late time) representing an infinite shut-in time.
It has a minimum value of 0.0 (at late time) representing an infinite shut-in time. P* is the extra-
polation of the radial line to FR = 0, which corresponds to an infinite shut-in time.
In many cases, injection rates on the order of 1 to 2 bpm (1440 – 2880 bbl/d) have been observed to
quickly achieve formation breakdown. With average injection times of 3 – 5 minutes, fracture closure
was observed within 1-to-24 hours of falloff. In some cases, radial flow developed quickly after clos-
ure, while in others, bilinear and linear flow was observed. When radial flow developed soon after
closure, a falloff period on the order of 1 – 3 days was sufficient to obtain reliable estimates of reser-
voir pressure and formation permeability.
When bilinear or linear flow is observed after the main closure event, the development of radial flow
is delayed, and may not be observed in a reasonable time-frame. Although this introduces greater
error in estimating the reservoir pressure and formation permeability, extending the falloff period
helps reduce the error. For these situations, the required falloff duration is case sensitive, and
depends on the importance of the information. Even if radial flow is not achieved, “upper limits” of
reservoir pressure and formation permeability can usually be determined, which are very useful.
An example of a minifrac test analyzed using WellTest is shown in the figures below. The test was
conducted on a vertical well at a formation depth of 10,000 ft. Pressures were monitored at the well-
head, and converted to sandface values for analysis. The total test duration was about 24 hours. Fig-
ure 1 shows the pressure profile during the 18-minute injection period performed at 1 bpm (1440
bbl/d). The sudden drop in pressure shortly after injection commenced indicates that formation
breakdown occurred very quickly.
Figure 1
Figure 2
Figure 3
Figure 4 represents the falloff data plotted with the radial time function for the Soliman / Craig solu-
tion.
Figure 4
Figure 5
Note: It has a minimum value of 0.0 (at late time) representing an infinite shut-in time.
Note: It has a minimum value of 0.0 (at late time) representing an infinite shut-in time.
Note: It has a minimum value of 0.0 (at late time) representing an infinite shut-in time. P* is the
extrapolation of the radial line to FR = 0, which corresponds to an infinite shut-in time.
For liquids, these assumptions are reasonable, since liquid compressibility and viscosity do not vary
significantly with pressure, and the equations can be solved analytically. These analytical solutions
are referred to as the liquid flow solutions, and form the basis of all well test analysis. The result is an
analytical relationship between pressure and time, which, for an infinite-acting reservoir, can be writ-
ten as
Pressure = Constant * log (time) + ...
For gas, most of the assumptions listed above are no longer valid. Gas compressibility (cg) varies
significantly with pressure. Gas viscosity (µg) and Gas Compressibility Factor also vary with pres-
sure but not to the same degree. Pseudo-pressure and pseudo-time (ta) are used to deal with these
changing properties and linearize the flow equations for gas. With the introduction of pseudo-pres-
sure and pseudo-time, the gas flow equation can be written in a manner similar to the liquid
equation. Therefore, the liquid flow solution can be used for gas well test analysis and forecasting
provided pressure is replaced by pseudo-pressure, and time is replaced by pseudo-time as follows:
Pseudo-Pressure = Constant * log (Pseudo-time) + …
l Depleting reservoirs
The first plot shows that type curves can be inaccurate if pseudo-time is not used. The second plot
shows that the buildup pressure does not reach up to the average reservoir pressure, when pseudo-
time is not used. Thus, it is imperative that pseudo-pressure and pseudo-time be used for accurate
analysis, modeling, and forecasting of gas production.
Development of Pseudo-Time
It should be noted that the concept of pseudo-time is not amenable to a completely rigorous solution,
as is the case for pseudo-pressure, because the gas properties change with pressure, not time. The
form of the diffusivity equation dictates which properties are grouped with pressure and which are
grouped with time.
In the 90s, when the gas flow equations were being used for analyzing or forecasting data affected
by reservoir depletion, it was realized that the Agarwal definition of pseudo-time for buildups, was
inappropriate for boundary-dominated flow (depleting systems). Moreover, the Agarwal pseudo-time
definition did not solve the problem, because it was using a simplified version of the total system
compressibility (ct).
Blasingame et al. introduced a new definition of pseudo-time to account for the depletion
effects. Instead of defining the pseudo-time transformation in terms of wellbore conditions like Agar-
wal did, they defined it in terms of the average reservoir pressure. This pseudo-time is best
described as material balance pseudo-time, which is appropriate for boundary-dominated flow.
Fekete's Rahman, Mattar, and Zaoral showed that a rigorous definition of compressibility, which
uses material balance principles, is required because fluid saturations can change over time.
Fekete's Anderson and Mattar found that, in a reservoir with significant transient flow, it was more
appropriate to define pseudo-time in terms of the average pressure within the region of investigation,
rather than the average reservoir pressure.
Two different pseudo-time formulations are needed: one for buildup and a different one for draw-
down. The equations look similar, but they are radically different from each other.
Definition of Pseudo-Time
Pseudo-time (ta) is defined as:
Note that the units for this definition of pseudo-time are hr/cp psi-1. In order to convert this unit to hr,
the integral is multiplied by µgi cti to obtain normalized pseudo-time as follows:
Buildup Pseudo-Time
Pseudo-time for use in buildup analysis, is defined in terms of pressure at the wellbore:
If the wellbore has significant storage, a difference in the early time behaviour will be noted, since
pressure and thus compressibility and viscosity are changing most near the wellbore. Late time
behavior is not typically affected by buildup pseudo-time since the pressure, viscosity, and com-
pressibility become constant at the wellbore and in the reservoir.
Drawdown Pseudo-Time
Pseudo-time for use in drawdown analysis, is defined in terms of average reservoir pressure:
Note that the gas in place (GIP) must be known in order to perform the material balance calculations
necessary to determine the average reservoir pressure.
TM
The analytical models in the Advanced Models tab in WellTest have an option to use Corrected
pseudo-time. This option uses the definition of drawdown pseudo-time that is based on the average
pressure within the region of investigation, rather than the average reservoir pressure. This is the
default setting, which we recommend you keep as is.
Agarwal suggested that the total system compressibility (ct) be defined as:
Where the compressibilities and saturations of residual fluids (oil (So) and water (Sw)) are con-
sidered constant. As shown by Rahman, Mattar, and Zaoral, the above definition of total system
compressibility (ct) does not use the material balance equation principles in most scenarios of deple-
tion. Therefore, a rigorous definition of total system compressibility (ct) for the purpose of computing
drawdown pseudo-time has been proposed as:
By applying drawdown pseudo-time, the shape of the derivative curve in late time changes depend-
ing on the size of the reservoir. As the reservoir depletes the average reservoir pressure changes,
which in turn affects the compressibility and viscosity in the late time region. Conversely, as the GIP
is increased, the effect of pseudo-time becomes less noticeable (approaches real time), since aver-
age reservoir pressure approaches initial pressure and compressibility and viscosity become con-
stant throughout the late time region.
Geomechanical effects can be modeled using pressure-dependent permeability and formation com-
pressibility. Further adjustment to the definition of pseudo-pressure and pseudo-time are used to
model these effects.
References
1. "A New Method for Computing Pseudo-Time from the Material Balance Equation for Real
Gas Flow", N. M. Anisur Rahman, Louis Mattar, and Karel Zaoral, Paper CIPC 2004 - 182
presented at the CIPC, Calgary, June 8 - 10, 2004.
2. "Theory and Practice of the Testing of Gas Wells, Third Edition", Energy Resources Con-
servation Board (ERCB), pp. 2 -16 to 2 - 30, 1975.
3. "Real Gas Pseudo-Time - A New Function for Pressure Buildup Analysis of MHF Gas
Wells", R.G. Agarwal, Paper SPE 8279 presented at AFTCE, Las Vegas, NV, September 23
- 26, 1979.
4. "The Flow of Real Gases through Porous Media", R. Al - Hussainy, H.J. Ramey, Jr., and P.B.
Crawford, JPT, pp. 625 - 636, May 1966.
5. "An Improved Pseudo-Time for Gas Reservoirs with Significant Transient Flow", D. Ander-
son, L. Mattar, JCPT 07-07-05, July 2007.
The Minifrac Models section contains a brief discussion about the placement of the start of the injec-
tion line. The Start of Injection line is usually placed at breakdown; however, some people prefer to
place it after the wellbore fill-up. The placement of the Start of Injection line impacts the time func-
tions for Pre-Closure and the Nolte After-Closure Analyses. In addition, the injected volume is cal-
culated as the cumulative injection from the Start of Injection line to shut-in. The injected volume is
used in the Nolte and Soliman / Craig After-Closure Analyses and the Soliman / Craig models.
Superposition time is also required when multi-rate analyses are used to approximate situations in
which the rate is slowly varying.
Note: The superposition time functions defined here and used in the software are different from
those published in the literature. Superposition radial time is defined as the anti-log of that
published in literature. Using this definition, data plotted on a semi-log plot of p vs. super-
position radial time can be compared to a semi-log plot of p vs. Horner time.
The following equations are the generalized forms of superposition time for each flow regime. They
can handle any number of step changes in rate.
Radial Time
Linear Time
Bilinear
Time
Spherical
Time
Pseudo-
Steady
State
For a variable rate drawdown test, the following values are provided:
t pwf q
(hrs) (psia) (stbbl / d)
0 5000 0
0.20 4890 80
2. For the second rate, a well is flowing at –50 stbbl / d from t = 0.1 hrs to the time point t = 0.2
hrs.
4. To get the superposition radial time, we add up all of the above times and take the anti-log of
the result:
1. For the first rate, a well is flowing at 200 stbbl / d from the start of the test to the time point t =
0.2 hrs.
2. For the second rate, a well is flowing at –50 stbbl / d from t = 0.1 hrs to the time point t = 0.2
hours.
4. To get the superposition linear time at t = 0.2 hrs, we add up all of the above times and take
the square of the result:
Bgi Gas formation volume factor at initial pressure (Rbbl / scf or Rbbl / MMscf)
Boi Oil formation volume factor at initial pressure (pi) (Rbbl / stbbl)
Bob Oil formation volume factor at the bubble point (Rbbl / stbbl)
D Turbulence factor
Der Derivative
DR Damage ratio
FCD
Dimensionless fault conductivity
(fault)
FE Flow efficiency
k Permeability (md)
m Semi-log slope
Inverse slope of the stabilized deliverability line (varies between 1.0 for completely lam-
n
inar flow and 0.5 for fully turbulent flow)
Nw Number of segments
p Pressure (psia)
pD Dimensionless pressure
pr Reduced pressure
rwc Corrected wellbore radius (ft) for htop > 0 (rwc = rw for htop = 0)
Rsb Solution gas oil ratio at the bubble point (scf / stbbl or MMscf /stbbl)
s Skin
t Time (hr)
ta Pseudo-time (hr)
tD Dimensionless time
Vsur-
Fluid volume at surface conditions (scf or stbbls)
face
ws Depth of fluid-loss damage into the formation and normal to the fracture face (ft)
zm Distance between the top of the productive zone and the middle of the open interval (ft)
γo Oil gravity
θr Relative temperature
Angle of inclination of the well, measured in degrees from a normal to the formation bed-
θw
ding plane (angle of slant)(degrees)
3 3
ρ Liquid density (lbm / ft or kg / m )
Pseudo-liquid density at reservoir pressure and standard temperature (60 °F) (lbm /
ρbs
ft3 or kg / m3)
3 3
ρg Gas density (lbm / ft or kg / m )
ρo Undersaturated oil density (lbm / ft3 or kg / m3)
ω Storativity ratio
2. Lougheed, D., Santo, M. and Ewens S. “Is That Radial Flow? What Can Be Learned From
Buildup Analysis of Multiply-Fractured Horizontal Wells in Unconventional Reservoirs” SPE
Paper 164525 prepared for presentation at the Unconventional Resources Conference-
USA, Woodlands, Texas, USA. April 2013.
3. Nobakht, M., Clarkson C., Kaviani, D. "New Type Curves for Analyzing Horizontal Well With
Multiple Fractures in Shale Gas Reservoirs" SPE 149397, Canadian Unconventional
Resources Conference, November 2011.
Major References
1. "Pressure Transient Testing", John Lee, John B. Rollins and John P. Spivey, SPE Textbook
Series, Vol. 9, Henry L. Doherty Memorial Fund of AIME, Society of Petroleum Engineers,
Richardson, Texas, 2003.
2. "Well Test Analysis: The Use of Advanced Interpretation Models", Handbook of Petroleum
Exploration and Production, Vol.3, Dominique Bourdet, Elsevier, Amsterdam, 2002.
3. "Fundamental and Applied Pressure Analysis", T.S. Daltaban and C.G. Wall, Reservoir
Modeling Division, Imperial College Press, London, 1998.
4. "Gas Reservoir Engineering", John Lee and Robert A. Wattenberger, SPE Textbook Series,
Vol. 5, Henry L. Doherty Memorial Fund of AIME, Society of Petroleum Engineers, Richard-
son, Texas, 1996.
5. "Well Test Analysis", Rajagopal Raghavan, Prentice Hall Petroleum Engineering Series,
New Jersey, 1993.
7. Pressure Transient Analysis, J.F. Stanislav and C.S. Kabir, Prentice-Hall, Inc., 1990.
10. Gas Well Testing - Theory and Practice, E.R.C.B. Fourth Edition, 1979 – Guide G-3.
11. Advances in Well Test Analysis, R.C. Earlougher, Jr., SPE Monograph, Vol. V, 1977.
12. Pressure Buildup and Flow Tests in Wells, C.S. Matthews and D.G. Russell, SPE Mono-
graph, Vol. I, 1967.
Minifrac References
1. “Holistic Fracture Diagnostics”, R. D. Barree, SPE, and V. L. Barree, Barree & Associates,
and Craig, SPE, Halliburton, Paper SPE 107877, Presented at the Rocky Mountain Oil &
Gas Technology Symposium held in Denver, Colorado, USA, 16-18 April 2007.
3. “After-Closure Analysis of Fracture Calibration Tests”, Nolte, K.G., Maniere, J.L., and
Owens, K.A., Paper SPE 38676, Presented at the 1997 SPE Annual Technical Conference
and Exhibition, San Antonio, Texas, 5-8 October, 1997.
4. "Background for After-Closure Analysis of Fracture Calibration Tests", Nolte, K. G., Paper
SPE 39407, Unsolicited companion paper to SPE 38676 Jul, 1997.
7. "New Method for Determination of Formation Permeability, Reservoir Pressure, and Frac-
ture Properties from a Minifrac Test", Soliman, M.Y., Craig D., Barko, K., Rahim Z., Ansah
J., and Adams D., Paper ARMA/USRMS 05-658, 2005.
8. "Design and Appraisal of Hydraulic Fractures", Jones, J.R. and Britt, L.K., SPE book, 2009.
2. Gringarten, A. C. and Ramey, H. J. Jr. “The Use of Source and Green’s Functions in Solving
Unsteady-Flow Problems in Reservoirs.” SPE Paper 3818. SPE Journal 13(5): pp. 285-296,
October 1973.
4. Thompson, L.G., Manrique, J.L., U. and Jelmert, T.A. “Efficient Algorithms for Computing
the Bounded Reservoir Horizontal Well Pressure Response.” SPE Paper 21827 prepared
for Presentation at the Low Permeability Reservoirs Symposium, Denver, Colorado. April
1991.
Gas Properties
"Gas Well Testing - Theory and Practice Appendix A", ERCB GUIDE G-3, Fourth Edition, 1979.
Oil Properties
1. "Estimating the Viscosity of Crude Oil Systems", Beggs, H.D. and Robinson, J.R., JPT
Forum, September 1975, pp. 1140 - 1141.
2. "Correlations for Fluid Physical Property Prediction", Vazquez, M.E. and Beggs, H.D.,
Journal of Petroleum Technology, June 1980, pp. 968 - 970.
Ng and Egbogah
"An Improved Temperature-Viscosity Correlation for Crude Oil Systems", J.T.H. Ng and E.O.
Egbogah, Petroleum Society of CIM 83-34-32, 1983.
Water Properties
1. "Applied Petroleum Reservoir Engineering", Craft, B.C. and M.F. Hawkins, Prentice-Hall,
1959, p. 131.
2. "A Correlation for Water Compressibility", Meehan, D.N., Petroleum Engineer, November
1980, pp. 125 - 126.
4. "Correlations for Physical Properties of Petroleum Reservoir Brines", Numbere, D., W.E.
Brigham and M.B. Standing, Petroleum Research Institute, Stanford University, November
1977, pp. 8 - 16.
Rock Properties
"Applied Petroleum Reservoir Engineering", Craft, B.C. and M.F. Hawkins, Prentice-Hall, 1959, p.
132.
Specific References
AOF
1. "Theory and Practice of the Testing of Gas Wells", L. Mattar, G. Brar, and M. Mumby, Energy
Resources Conservation Board (1978) Third Edition, Chapter 3.
2. "The Primary Pressure Derivative (PPD) – A New Diagnostic Tool in Well Test Inter-
pretation", L. Mattar and K. Zaoral, JCPT (April 1992) Volume 31, No. 4, 63 - 70.
3. "Pressure Buildup Analysis with Wellbore Phase Redistribution Effects", Fair, W.B., Jr, Soc.
Pet. Eng. J. (April 1981) pp. 259 - 270.
4. "Well-Test Analysis with Changing Wellbore Storage", P.S. Hegeman, D.L. Hallford, and
J.A. Joseph, SPEFE (September 1993) 201 - 207.
3. "Optimization of the Productivity Index and the Fracture Geometry of a Stimulated Well with
Fracture Face and Choke Skins", D.J. Romero, P.P. Valko and M.J. Economides, SPEPF
(February 2003) 57 - 64.
Composite Model
1. "An Analytical Model for Composite Reservoirs Produced at Either Constant Bottomhole
Pressure or Constant Rate", Olarewaju, J.S. and Lee, W.J., SPE 16763 (1987).
Deconvolution
Thomas Von Schroeter et al. and Michael M. Levitan et al.: SPE 71574, SPE 77688, SPE 84290,
SPE 90680.
Derivative Analysis
1. "Derivative Analysis Without Type Curves", Louis Mattar, JCPT Special Edition (1999) vol.
38.
2. "Use of Pressure Derivative in Well-Test Interpretation", D. Bourdet, J.A. Ayoub, and Y.M.
Pirard, SPE Formation Evaluation, (June 1989).
Dual Porosity
1. "Gas Reservoir Engineering", J. Lee and R. Wattenbarger, Society of Petroleum Engineers
Inc. (1996) Volume 5, 73 - 74, 173 - 181.
2. "The Behavior of Naturally Fractured Reservoirs", J.E. Warren and P.J. Root, Society of Pet-
roleum Engineers Inc. (1963) SPEJ 426.
2. "Pressure Analysis for Fractured Wells", A.C. Gringarten, H.J. Ramey, Jr., and R.
Raghavan, Paper SPE 4051 presented at the 1972 AFM, San Antonio, TX, October 8 - 11.
Horizontal Model
1. "Efficient Algorithms for Computing the Bounded Reservoir Horizontal Well Pressure
Responses", L.G. Thompson, J.L. Manrique and T.A. Jelmert, Paper SPE 21827 presented
at 1991 Rocky Mountain Regional Meeting and Low-Permeability Reservoirs Symposium,
Denver, CO, April 15 - 17.
IPR
"Inflow Performance Relationship for Solution-Gas Drive Wells", Vogel, J. V., JPT, Jan. 1968.
2. "The Application of the Laplace Transformation to Flow Problems in Reservoirs", A.F. van
Everdingen and W. Hurst, Trans. AIME (1949) vol. 186, 305 - 324.
3. "A Study of the Behaviour of Bounded Reservoir Composed of Stratified Layers", H.C.
Lefkovits, P. Hazebroek, E.E. Allen and C.S. Matthews, SPEJ (March 1961) 43 - 58.
4. "Transient Pressure Behaviour of Commingled Reservoirs", F.J. Kuchuk and D.J. Wilkinson,
SPEFE (1991).
2. "A Study of the Behaviour of Bounded Reservoir Composed of Stratified Layers", H.C.
Lefkovits, P. Hazebroek, E.E. Allen and C.S. Matthews, SPEJ (March 1961) 43 - 58.
3. "Wells Producing Commingled Zones with Unequal Initial Pressures and Reservoir Prop-
erties", L. Larsen, Paper SPE 10325 for Presentation in AFTCE, San Antonio, TX, October 5
- 7, 1981.
4. "An Efficient Algorithm for Computation of Well Responses in Commingled Reservoirs", J.B.
Spath, E. Ozkan and R. Raghavan, SPEFE (June 1994) 115 - 121.
2. "The Use of Source and Green’s Functions in Solving Unsteady-Flow Problems in Reser-
voirs", A.C. Gringarten and H.J. Ramey, Jr., SPEJ (October 1973) 285 - 296.
3. "Transient Pressure Behaviour of Commingled Reservoirs", F.J. Kuchuk and D.J. Wilkinson,
SPEFE (1991).
PITA Analysis
1. "Use of PITA for Estimating Key Reservoir Parameters", N.M. Anisur Rahman, Mehran
Pooladi-Darvish, Martin S. Santo and Louis Mattar, Paper CIPC 2006 - 172 presented at
2006 Canadian International Petroleum Conference, Calgary, Alberta, Canada, June 13 -
15.
3. "Perforation Inflow Test Analysis (PITA)", N.M. Anisur Rahman, Mehran Pooladi-Darvish
and Louis Mattar, Paper CIPC 2005 - 031 presented at 2005 Canadian International Pet-
roleum Conference, Calgary, Alberta, Canada, June 7 - 9.
PITA Model
1. "Use of PITA for Estimating Key Reservoir Parameters", N. M. Anisur Rahman, Mehran
Pooladi-Darvish, Martin S. Santo and Louis Mattar, Paper CIPC 2006 - 172, presented at
7th Canadian International Petroleum Conference, Calgary, AB, June 13 - 15, 2006.
2. "Development of Equations and Procedure for Perforation Inflow Test Analysis (PITA)", N.
M. Anisur Rahman, Mehran Pooladi-Darvish and Louis Mattar, Paper SPE 95510, presen-
ted at 80th Annual Technical Conference and Exhibition of the SPE, Dallas, TX, October 9 -
12, 2005.
Pseudo-Pressure
1. "Theory and Practice of the Testing of Gas Wells, Third Edition", Energy Resources Con-
servation Board (ERCB), pp. 2 -16 to 2 - 30, 1975.
2. "The Flow of Real Gases through Porous Media", R. Al - Hussainy, H.J. Ramey, Jr., and P.B.
Crawford, JPT, pp. 625 - 636, May 1966.
Pseudo-Time
1. "A New Method for Computing Pseudo-Time from the Material Balance Equation for Real
Gas Flow", N. M. Anisur Rahman, Louis Mattar, and Karel Zaoral, Paper CIPC 2004 - 182
presented at the CIPC, Calgary, June 8 - 10, 2004.
2. "Theory and Practice of the Testing of Gas Wells, Third Edition", Energy Resources Con-
servation Board (ERCB), pp. 2 -16 to 2 - 30, 1975.
4. "The Flow of Real Gases through Porous Media", R. Al - Hussainy, H.J. Ramey, Jr., and P.B.
Crawford, JPT, pp. 625 - 636, May 1966.
2. "Numerical Simulations of the Combined Effects of Wellbore Damage and Partial Pen-
etration", R.M. Saidikowski, Paper SPE 8204 presented at 1979 AFTCE, Las Vegas, NV,
September 23 - 26.
3. "Optimization of the Productivity Index and the Fracture Geometry of a Stimulated Well with
Fracture Face and Choke Skins", D.J. Romero, P.P. Valko and M.J. Economides, SPEPF
(February 2003) 57 - 64.
Slug Model
1. "Analysis of Slug Test or DST Flow Period Data", H.J. Ramey, Jr., R.G. Agarwal and I.
Martin, JCPT (July - September 1975) 37 - 47.
2. "Annulus Unloading Rates as Influenced by Wellbore Storage and Skin Effect", H.J. Ramey,
Jr. and R.G. Agarwal, SPEJ (October 1972) 453 - 462.
3. "Analysis of Slug and Drillstem Tests", A.M.M. Peres, Ph.D. Thesis, The University of Tulsa
(1989).
Two-Phase Skin
"Practical Consideration in the Analysis of Gas Condensate Well Tests", R. Raghavan, W., C. Chu,
and J.R. Jones, Paper SPE 30576 presented at 1995 ATCE, Dallas, TX, October 22 - 25.
Vertical Model
1. "The Application of the Laplace Transformation to Flow Problems in Reservoirs" A.F. van
Everdingen and W. Hurst, Trans. AIME (1949) vol. 186, 305 - 324.
2. "The Use of Source and Green’s Functions in Solving Unsteady-Flow Problems in Reser-
voirs", A.C. Gringarten and H.J. Ramey, Jr., SPEJ (October 1973) 285 - 296.
2. "Computing Pressure Distribution in Wedges," C.-C. Chen and R. Raghavan, SPEJ (March
1997) 24 - 32.
Support for our products is available Monday to Friday, 8:00 a.m. – 4:30 p.m. mountain time. Limited
support is available on North American holidays. We will respond to all inquiries within 24 hours, or
by the next business day.
For WellTest-specific information, see WellTest Software Training. For general information on all of
our offerings, see Global Upstream Energy Education & Training.
Resources
We have technical videos, published papers, and unconventional energy maps for your use. To view
these resources, see our Engineering Technical Resources.
Customized training
In addition to our regularly scheduled courses, we offer private training for our software to groups.
Allow us to customize your training by using your data to create workflows to suit your specific
assets.