RC Manual V6
RC Manual V6
A M ULTI -D ISCIPLINARY
T EAM A PPROACH
PetroSkills, LLC.
2930 South Yale Avenue
Tulsa, Oklahoma 74114-6252
United States
1.918.828.2500--Office
1.800.821.5933--USA
1.918.828.2580--Fax
www.petroskills.com
The information contained herein and/or these workshop/seminar proceedings (WORK) was prepared by
or contributed to by various parties in support of professional continuing education.
For purposes of this Disclaimer, “Company Group” is defined as PetroSkills, LLC.; OGCI Training, Inc.;
John M. Campbell and Company; its and their parent, subsidiaries and affiliated companies; and, its and
their co-lessees, partners, joint ventures, co-owners, shareholders, agents, officers, directors, employees,
representatives, instructors, and contractors.
Company Group takes no position as to whether any method, apparatus or product mentioned herein is
or will be covered by a patent or other intellectual property. Furthermore, the information contained
herein does not grant the right, by implication or otherwise, to manufacture, sell, offer for sale or use any
method, apparatus or product covered by a patent or other intellectual property right; nor does it insure
anyone against liability for infringement of same.
Company Group does not guarantee results. All interpretations using the WORK, and all
recommendations based upon such interpretations, are opinion based on inferences from measurements
and empirical relationships, and on assumptions, which inferences and assumptions are not infallible, and
with respect to which competent specialists may differ. In addition, such interpretations, recommendations
and descriptions may involve the opinion and judgment of the USER. USER has full responsibility for all
interpretations, recommendations and descriptions utilizing the WORK. Company Group cannot and does
not warrant the accuracy, correctness or completeness of any interpretation, recommendation or
description. Under no circumstances should any interpretation, recommendation or description be relied
upon as the basis for any drilling, completion, well treatment, production or other financial decision, or any
procedure involving any risk to the safety of any drilling venture, drilling rig or its crew or any other
individual. USER has full responsibility for all such decisions concerning other procedures relating to the
drilling or production operations. Except as expressly otherwise stated herein, USER agrees that
COMPANY GROUP SHALL HAVE NO LIABILITY TO USER OR TO ANY THIRD PARTY FOR ANY
ORDINARY, SPECIAL, OR CONSEQUENTIAL DAMAGES OR LOSSES WHICH MIGHT ARISE
DIRECTLY OR INDIRECTLY BY REASON OF USER’S USE OF WORK. USER shall protect, indemnify,
hold harmless and defend Company Group of and from any loss, cost, damage, or expense, including
attorneys’ fees, arising from any claim asserted against Company Group that is in any way associated
with the matters set forth in this Disclaimer.
No part of the text or original figures of the WORK may be reproduced in any form without the express
written consent of PetroSkills, LLC. Use of the WORK as a reference or manual for adult training
programs is specifically reserved for PetroSkills, LLC. All rights to the WORK, including translation rights,
are reserved.
PROFESSIONAL EXPERIENCE
Teaching Experience
Instructor/coach at Shell’s introductory course for graduates teaching Geology,
Geophysics, Petrophysics, Reservoir Engineering and economics (Rijswijk, 2006-
present day)
Mathematics teacher at Stanislas college (Delft, 1999-2005)
Physics teacher at Technical High School (Amsterdam, 1967)
Shell Experience
Liaison Officer between Shell in the Netherlands and Shell in Nigeria (1996-
1999)
Chief petrophysicist Shell Group (The Hague, 1994-1996)
Procurement and Logistics Manager (Nigeria, 1993-1994)
Chief petrophysicist (Nigeria, 1989-1993)
Manager petrophysical research and development (Rijswijk, 1987-1989)
Section head rock mechanics (Rijswijk, 1983-1987)
Section head reservoir engineering (Brunei, 1981-1983)
Reservoir Engineer (Den Haag, 1979-1981)
Research petrophysicist (Rijswijk, 1978-1979)
Operational petrophysicist (Aberdeen, 1976-1978)
Research petrophysicist (Rijswijk, 1974-1976)
EDUCATION
Gymnasium-β (1965)
Physics (University of Amsterdam, 1969)
PhD in Physics (University of Amsterdam, 1972)
C H A P T E R 2 D A T A F O R C H A R A C T E R I Z I N G F I E L D S ................................................. 2-1
OBJECTIVES ............................................................................................. 2-1
DEVELOP DETAILED CHARACTERIZATION PLAN ............................................ 2-2
PROBABILITY DISTRIBUTIONS ..................................................................... 2-3
MULTIPLE REALIZATIONS OF PERMEABILITY ............................................... 2-10
DEPTH OF INVESTIGATION AND RESOLUTION OF LOGGING TOOLS ................ 2-16
SCALES OF MEASUREMENT ...................................................................... 2-17
QUALITY/ERROR MINIMIZATION ................................................................ 2-23
DATA MANAGEMENT ................................................................................ 2-24
EXTRAPOLATION TO OTHER AREAS........................................................... 2-29
HOW CAN WE INTEGRATE DATA ................................................................. 2-31
DATA - SUMMARY .................................................................................... 2-31
SUMMARY OF CHAPTER 2 ......................................................................... 2-32
C H A P T E R 4 G E O S T A T I S T I C A L M O D EL I N G .............................................................. 4-1
INTEGRATED MODELING ............................................................................. 4-2
MODELING CONCEPTS ............................................................................... 4-3
GEOSTATISTICAL APPLICATIONS ................................................................. 4-4
OBSERVATIONS ABOUT RESERVOIRS .......................................................... 4-4
GEOSTATISTICAL OVERVIEW ...................................................................... 4-5
METHODS TO ASSIGN VALUES .................................................................... 4-6
UNIVARIATE ANALYSIS ............................................................................... 4-8
BIVARIATE ANALYSIS ................................................................................. 4-8
GEOSTATISTICAL METHODS ..................................................................... 4-10
GEOSTATISTICAL TOOLS .......................................................................... 4-11
MAPS AND MAP MAKING (CHAMBERS ET AL) .............................................. 4-13
BUILDING A VARIOGRAM........................................................................... 4-13
BOOLEAN OR OBJECT BASED SIMULATION (SRIVASTAVA 1994) .................. 4-29
SIMULATED ANEALING (SRIVASTAVA 1994) ............................................... 4-36
GEOSTATISTICS CASE STUDIES ................................................................ 4-40
THE TRUTH ............................................................................................. 4-44
SUMMARY OF CHAPTER 4 ......................................................................... 4-49
GLOSSARY OF VARIOGRAM TERMS, KUPFERSBERG AND DEUTSCH, 1999 .... 4-50
C H A P T E R 5 S E Q U E N C E S T R A T I G R A P H Y ................................................................. 5-1
SCENARIO TREES ...................................................................................... 5-2
RESERVOIR CHARACTERISATION WORKFLOWS ............................................ 5-3
STATIC RESERVOIR MODELS ...................................................................... 5-4
FACIES ..................................................................................................... 5-5
STRATIGRAPHY ......................................................................................... 5-5
BASIC IDEAS & DEFINITIONS ....................................................................... 5-6
POWER OF SEQUENCE STRATIGRAPHY ....................................................... 5-8
SEQUENCE STRATIGRAPHIC CYCLES ......................................................... 5-11
BIOSTRATIGRAPHY .................................................................................. 5-21
LITHO- AND CHRONOSTRATIGRAPHY ......................................................... 5-22
FLOW CHART FOR MODELING ................................................................... 5-26
SUMMARY OF CHAPTER .......................................................................... 5-26
C H A P T E R 6 L I T H O L O G I C M O D E L S ........................................................................... 6-1
MONTE CARLO SIMULATION ........................................................................ 6-3
TRAPS ...................................................................................................... 6-4
C H A P T E R 10 E C O N O M I C S A N D R I S K I N G ................................................................ 10-1
RISKING / ECONOMICS CHALLENGE ........................................................... 10-1
OBJECTIVES FOR ECONOMICS .................................................................. 10-2
RISKED BUSINESS VENTURES ................................................................... 10-8
NET PRESENT VALUE (NPV) ................................................................... 10-10
CASH SURPLUS AND ACCOUNTING PROFIT ............................................... 10-11
VALUING RESERVOIR STUDIES / DATA (VOI) ............................................. 10-12
RC EXAMPLE ........................................................................................ 10-14
SUMMARY OF CHAPTER 10 ..................................................................... 10-15
C H A P T E R 11 O R G A N I Z A T I O N A L S T R U C T U R E ........................................................ 11-1
TEAM STYLES ......................................................................................... 11-1
ORGANIZATIONAL MANAGEMENT .............................................................. 11-3
WHITNEY CANYON FIELD TEAM - AMOCO - HISTORY ................................. 11-11
TEAM LEADERSHIP SKILLS ...................................................................... 11-12
TEAM BUILDING ELEMENTS .................................................................... 11-18
TEAM COMMUNICATION.......................................................................... 11-20
SUMMARY OF CHAPTER 11 ..................................................................... 11-24
C H A P T E R 12 P R O J E C T M A N A G E M E N T .................................................................... 12-1
PROJECT PLANNING AND IMPLEMENTATION ............................................... 12-3
SETTING RESERVOIR CHARACTERIZATION TEAM GOALS ........................... 12-12
SUMMARY OF CHAPTER 12 ..................................................................... 12-15
C H A P T E R 13 P R O B L E M S , P I T F A L L S , A N D T H E F U T U R E O F R E S E R V O I R ........ 13-1
C H A R A C T E R I Z AT I O N ............................................................................ 13-1
PROBLEMS, PITFALLS .............................................................................. 13-1
THEMES FOR SUCCESSFUL RESERVOIR CHARACTERISATION ....................... 13-3
THE FUTURE OF RESERVOIR CHARACTERIZATION ...................................... 13-3
C H A P T E R 14 R E F E R E N C E S ....................................................................................... 14-1
C H A P T E R 15 E X E R C I S E M A T E R I A L ....................................................... 15-1
G A S P O O L E X E R C I S E - P A R T 1 ........................................... 15-1
WELL CORRELATION ................................................................................ 15-3
VARIOGRAM EXERCISE KINGFISH FIELD ..................................................... 15-5
LITHO- VERSUS CHRONO-STRATIGRAPHIC CORRELATION ............................ 15-8
G A S P O O L E X E R C I S E - P A R T 1 ........................................................... 15-9
PRESSURE VERSUS DEPTH FOR THE GRASSHOPPER FIELD ........................ 15-11
SEISMIC INTERPRETATION ...................................................................... 15-17
MATERIAL BALANCE EXERCISE ................................................................ 15-18
WHAT IS THE RIGHT PERMEABILITY? ........................................................ 15-20
G U I D E L I N E S F O R D EL I V E R I N G A N E F F E C T I V E P R E S E N T A T I O N ...... 15-24
It must be economic
All activities - to produce the evidence - in line with the corporate vision
The models developed during the course will be based on the application of state-of-
the-art technical applications within the framework of a multi-disciplinary team approach.
The integration of data required for modern reservoir characterization encourages
changes in the organizational structure of companies. Modern multi-disciplinary team
approaches relieve many of the conflicts that result from traditional reporting
relationships.
There are several key elements in successfully characterizing reservoirs. These can be
categorized into related areas of learning; the course will treat these categories in depth
in the following sections; case examples will be included in the discussions to clarify the
concepts and to demonstrate practical methods of solving these problems.
Sources
Scale of the data/extrapolation to other areas
Cross-disciplinary applications/integration
Quality/error minimization and measurement
Data management/planning
Understanding what is represented by data
Geostatistical Applications
Reservoir models
Sequence stratigraphic
Lithologic
Geophysical
Reservoir engineering
Economics and risking
Volumetrics
Probability of success
Financial returns of project
Organizational structure
Team styles
Organizational management
Team communications
Assessment and evaluation
Problems, pitfalls and the future of reservoir characterization
3. Data gathering to describe the reservoir normally begins with seismic and continues
with well information and then actual production results.
4. An initial model of the reservoir may be built and refined on a continual basis as
further information is obtained so as to understand the movement of fluids and to
develop a plan of production with the ultimate goal of maximizing the hydrocarbon
recovery while achieving a positive cash flow.
RC PRINCIPLES
It is business
MODELING OBJECTIVES
Assessment of reserves
Production of the field development program (FDP)
Identification of additional exploitation opportunities
Particularly, the third objective can have a significant impact on the economics of
projects. If an expensive pipeline is to be laid, the economics will improve if additional
reservoirs can be found along the track of the pipeline that may feed into it. If satellite
fields can be found in the vicinity of an offshore platform, the lifetime of the platform may
be extended improving the economics significantly.
The hydrocarbon life cycle shows that RC plays a role from delineation op to tertiary
production. Throughout, creativity in RC is required. In this creativity, divergent and
convergent thinking is important. For any project, different stages are identified. The
early stage is a phase with divergent thinking, during which all possible options for the
project are being considered. At the end of this stage a feasibility report is presented.
In the second and third stages, the thinking is convergent in order to select and define
the project. The second stage produces the FDP (field development plan) and the final
stage a project execution plan including a specific analysis of required investments.
After the third stage, management can take the decision to go ahead with the project
and should make the necessary funds available.
Depending on the nature of the study, typical elements that may be included in the
study report, are:
Uncertainty analysis
Volumetrics
Flow modeling
Clearly, a seismic project will have different elements from thos of a water injection
project. However, the three staged approach with divergent and convergent thinking is
applicable to all projects.
and when a field or reservoir has been on production for some time. Different study
methods are applied in these phases, namely volumetric methods and methods based
on production performance. The range of uncertainties will reduce over time, but even
at the time of decommissioning uncertainties remain. Yes, the actual recovery may be
well known, but STOIIP or GIIP may still be uncertain. With time, when more data and
knowledge become available, the relative risk will reduce.
Each project is special depending on the nature of the project and the techniques and
technologies that may be applied to provide the evidence for the business decisions.
The slide on the next page shows a range of techniques and technologies such as
using the MDT (modular dynamic tool), PVT (pressure, volume and temperature
analysis), SCAL (special core analysis) and core-calibrated FAL (formation attribute
log). When applying these tools, the purpose is to reduce the uncertainties and increase
the reserves.
Uncertainties cannot be reduced to zero, but are reduced from an initial level to an
acceptable level. Acceptable means that there are contingency plans available that can
deal with the uncertainties. This implies that any decisions being taken, are taken under
risk.
In this context, the concept of reservoir management must be highlighted. The purpose
of reservoir management is to manage the reservoir in line with corporate objectives.
There may be different corporate strategies such as optimizing production or optimizing
ultimate recovery. These are often conflicting strategies and cannot always be pursued
at the same time.
The purpose of any hydrocarbon project is to realize an ultimate cash surplus upon
completion:
Ultimate
Cash Surplus
Time
Payout
Maximum Time
Exposure
First
Oil
During the life cycle of a project, different stages are followed, namely:
Exploration
Delineation
Development
Maturity
In each of these phases, there are opportunities to optimize the project, which may lead
to acceleration of the project. In the following pages these acceleration opportunities are
being discussed.
STAGE 1: EXPLORATION
The first stage in the life of a reservoir may span a period of two years. At this stage, the
reservoir and its contents are just beginning to be defined and understood. Early
interpretation of geologic and seismic data is necessary to confirm the shape and extent
of the reservoir and its contents. Such interpretations and analyses lead to reduction of
uncertainty and risk in the early stage of field delineation and development. More
specifically, they lead to optimal placement of future wells.
Acceleration Options
Drill fewer wells
– Understand reservoir better
Better earth model
3D data
Early, conceptual simulation
– Horizontal/multi-lateral wells
STAGE 2: DELINEATION.
The delineation stage of a reservoir can range from two to five years. Initial wells are
completed in the producing zone with the aim of delineating the size and extent of the
reservoir. Direct and indirect information gathered from these wells and from seismic
surveys significantly improves our understanding of the reservoir structure, productivity
and reserves. During this stage, field development plans are created and a variety of
wells-vertical, horizontal and multilaterals-may be designed and drilled to build up
production and deplete the reservoir efficiently.
Acceleration Options
Earlier production
– Contemporaneous drilling/
production
STAGE 3: DEVELOPMENT
Reservoir development can take place over a period of five to fifteen years depending
on the size of the reservoir. This is the phase of maximum activity in the field, which
further improves reservoir understanding.
Acceleration Options
Build alternative plans (workovers/
recompletions) at time of initial
completion
Empower team to make operational
business choices at fairly high levels
Alliances with service companies to
update technical level and their local
knowledge
Involve field and plant personnel in
planning/expenditures
Visit adjacent facilities for best practices
STAGE 4: MATURITY
During the mature stage, which ranges from 15 years to abandonment, the reservoir is
well understood. Oil depletion by production wells causes the gas cap and aquifer to
expand and make the oil column thinner.
Wells are stimulated by fracturing, acidizing or other means to improve productivity and
prolong their economic life. Water or gas injection and artificial lift continue. At the end
of this stage, the wells are plugged and abandoned as their production falls below
economic levels.
----------------------------------------------------------------------------------------------------------------
No Match Results
Negative
Collect New YES Additional
Data Information Information Acceptable
Available? Results
NO
Proceed to
Identify New Not Next RM
Cost-Effective
Data to Cost-Effective Plan Step
to Collect
Improve Models to Collect
Revise RM
Plan Step
Team Review
Petrophysical data
Normalize/calibrate/edit against cores
Identify cutoffs for mapping
Core data
Fully describe/interpret cores
Correlate with logs
Calibrate with test/production data
Seismic data
Develop structural interpretation
Seismic stratigraphic correlation
Establish correlations with reservoir parameters
Define facies
Well History
Completed interval, recompletions
Stimulation history
Production data/injection
Fluid data
Geologic structure
Stratigraphic model
Modeling techniques - geostatistics, attribute assessment and modeling
Fluid distribution - vertical and horizontal
Production history
Facilities
RC team members
2 3D surveys
1987 - 12,000’ cable, 1990 - 19,000’ cable
INTEGRATION
To achieve the best results in RC, integration is of vital importance. Usually, RC teams
are set up with participation of the traditional subsurface disciplines. However,
depending on the scope of a RC project it is also important to have less conventional
disciplines on the team, such as legal, environmental and other.
Disciplines Data
Geoscience Geologic
Reservoir engineering Geophysical
Production engineering Petrophysical
Drilling Core
Negotiating/legal Pressure
Financial Production
Fiscalists Well test
Field personnel Financial
Facilities engineering Simulation
Management
Technology Tools
PSDM Geostatistics
Horizontal wells Seismic
Multi-lateral wells Logs
Production system simulators Reservoir simulators
Tomography Work stations
Etc. Analogues
SUMMARY OF CHAPTER 1
The purpose of Reservoir Characterisation (RC) is to provide evidence for business
decisions. Therefore, critical elements are being studied following a six staged RC
process. During the third step of this process, agreement is reached between the RC
team and their sponsors regarding the terms of reference for the RC project.
In the technical arena the static model needs to be studied in terms of its architecture
and internal structure as well as the dynamic model producing production forecasts and
estimates of reserves. The static and dynamic models are strongly interlinked
representing in fact one single model. Such a model is only complete when in addition
all the uncertainties have been specified.
In RC different stages are recognized, namely an early stage in which the feasibility of
an opportunity is assessed. This stage is characterized by divergent thinking. Stages
two and three are characterized by convergent thinking, when the opportunity is
specified in detail and the way forward is defined including an investment plan.
During projects, there may be opportunities to accelerate them. Value drivers such as
‘advance first production’ or ‘increase plateau rate’ will improve the economics of
projects.
For optimum success, integration of data, people, technology and tools is of vital
importance. With respect to integration of people, disciplines - such as for instance
legal, HSE and public affairs - different from the usual subsurface need to be included
when necessary.
The first step in data analysis involves creating an inventory of the data available and
the quality of those data. This preliminary review is focused upon a scoping of the data
- not interpretation. Key questions are whether the data are viable and adequate for the
purpose of the study.
OBJECTIVES
DATA CONSIDERATIONS
Sources
Scale of the data
Cross-disciplinary applications/integration
Quality/error minimization
Data management
Do data represent larger area?
Extrapolation to other areas
PROBABILITY DISTRIBUTIONS
An Important Concept - Even though we describe potential values with a range, the
truth is that a single value represents the prospect or play value under consideration.
For example, if we predict that a prospect has most likely gross reservoir thickness of
300 ft (90 m), a minimum gross reservoir thickness of 50 ft (15 m) and a maximum
gross reservoir thickness of 550 ft (165 m), we are describing our assessment of the
geologic conditions based on our best knowledge. The actual gross reservoir thickness
is a single value that – hopefully - falls within the predicted range - ideally near the most
likely expected value.
Exceedance/Cumulative*
Normal (gaussian or bell-shaped)
Lognormal
Histogram
Equal
Rectangular
Triangular*
Log-triangular*
Symmetrical triangular distributions are safest distributions to use, because the mean
value or "weight" of the variable in the final assessment is closer to the Most Likely
value (the assessor's Best Interpretation) than it is if the triangle is skewed. The mean
values of both normal and log triangles can be calculated as shown below.
Symmetrical distributions.
CUMULATIVE
In cumulative curve displays, small probabilities values indicate small numbers
(example P10 is a thin reservoir).
20
0
9 18 27 36 45 54 63 72 81 90
Cumulative curve.
EXCEEDANCE
In exceedance curve displays, small probabilities values indicate large numbers
(example P10 is a thick reservoir). Companies use both displays and care needs to be
taken to understand which is displayed. Usually the context will indicate the distribution.
Exceedance curves will be used as the default display of population in this course.
Some cumulative presentations from other authors may be used.
120
Economic Threshold
- 40'
100
Exceedance
80 distribution
A chart that shows the
60 number or proportion (or
percentage) of values
greater than or equal to
40 a given amount.
20
0
9 18 27 36 45 54 63 72 81 90
Exceedance curve.
NORMAL DISTRIBUTION
Describes many natural phenomena (IQ's, people's heights, the inflation rate, or
errors in measurements)
Continuous probability distribution
Parameters are:
− Mean
− Standard deviation
Some value is the most likely (the mean of the distribution)
The unknown variable could as likely be above or below the mean (symmetrical
about the mean)
The unknown variable is more likely to be close to the mean than far away
− Approximately 68% are within 1 standard deviation of the mean
Some equations:
68% of
Data
95% of
Data
99.7%of
Data
–4 –3 –2 –1 0 1 2 3 4
Normal distributions.
LOGNORMAL DISTRIBUTION
Widely used in situations where values are positively skewed (where most of the
values occur near the minimum value)
− Financial analysis for security valuation
− Real estate for property valuation
− Distribution of reserves in a play
Continuous probability distribution
Financial analysts have observed that the stock prices are usually positively
skewed.
− Stock prices exhibit this trend because the stock price cannot fall below the
lower limit of zero but may increase to any price without limit.
The parameters for the lognormal distribution
− Mean
− Standard deviation
Three conditions underlying a lognormal distribution are:
1. The unknown variable can increase without bound, but is confined to a finite
value at the lower limit.
2. The unknown variable exhibits a positively skewed distribution.
3. The natural logarithm of the unknown variable will yield a normal curve.
HISTOGRAMS
Histograms demonstrate the results of many trial of a given set of geologic conditions.
The resulting plots show the frequency of answers within given ranges. This is a
desirable approach to test the models; unfortunately, it is not often possible to generate
many solutions. What do we do in this case?
A frequent answer is to define the potential outcomes of the models by estimating
minimum and maximum values and a central value. There is a great deal of
controversy around how to do this. The most important judgment is to apply sound
geoscience reasoning to the outcome and to assure that whatever technique is used,
that the results are realistic. Some of these approaches are described below and
demonstrated here.
RECTANGULAR DISTRIBUTION
All values between the minimum and maximum are equally likely to occur.
Continuous probability distribution.
The parameters for the uniform distribution are minimum and maximum.
Three conditions:
1. The minimum value is fixed.
2. The maximum value is fixed.
3. All values between the minimum and maximum are equally likely to occur.
TRIANGULAR DISTRIBUTION
The normal triangular distribution, as well as other distributions, can be transformed into
a cumulative distribution. The triangle is ideal for assessment input because it
emphasizes a Most Likely (Best Interpretation) value, the Minimum and Maximum
limits are fixed (there are no surprises), and its mathematics are relatively simple. The
third leg of any symmetrical triangle can be calculated if the other two are known.
Shows number of successes when you know the minimum, maximum, and most
likely values.
Continuous probability distribution.
The parameters for the triangular distribution are minimum, maximum, and
likeliest.
− For example, you could describe the number of cars sold per week when past
sales show the minimum, maximum, and most likely number of cars sold.
Three conditions:
1. The minimum number is fixed.
2. The maximum number is fixed.
3. The most likely number falls between the minimum and maximum values,
forming a triangular shaped distribution, which shows that values near the
minimum and maximum are less likely to occur than those near the most
likely value.
NORMAL TRIANGLE
(e.g., 2 - 4 - 6)
MAXIMUM = 2 ML - MIN = 2 x 4 - 2 = 6
Normal triangle.
LOG TRIANGLE
(e.g., 2 - 4 - 8)
Log triangle.
LOG-TRIANGULAR DISTRIBUTION
The log triangular distribution is also useful, especially for describing those variables
that seem to be lognormally distributed in nature. The normal and log triangles can
closely approximate the bell-shaped normal and lognormal distributions. These figures
also show the significant difference between normal and log triangular distributions that
have the same Minimum and Maximum endpoints.
Figure 2-1. An asymmetrical data set. The three measures Figure 2-2. A symmetrical data set. The three measure of
of central tendency are different. central tendency are identical.
WELL THREE
MEAN 10.32
0.16
COEF. OF VAR 0.39
MAXIMUM 24.12
MEDIAN 9.08
FREQUENCY
0.12
MINIMUM 4.72
0.08
0.06
0.00
0 5.0 10.0 15.0 20.0
WELL THREE
DATA LIMITATIONS
Impact of each parameter can’t be predicted
(e.g. the impact of the permeability distribution on the production forecast)
Usually large inter-well regions without any data apart from seismic data
Planning for acquisition is often ignored to the detriment of the overall plan. A multi-
disciplinary approach to planning for data acquisition and usage will maximize the
benefit of the data over time and across disciplines and will minimize expense
associated with re-acquiring data at a later time of acquiring unnecessary data.
Scale of investigation
(Sech et al, 2009)
1. For the list of terms below, indicate the scale of investigation (in feet or meters), the
purpose of the tool or technique and whether the data type is measured or
interpreted.
4. You may also wish to consider the different scales of investigation laterally and
vertically.
Core
Outcrop
Seismic - 2D
DST
Well test
Production
Seismic methods trade coverage of the reservoir for resolution because it is impractical
today to achieve both high resolution and high coverage. Crosswell methods fill a
resolution “gap” between sonic log measurements and vertical seismic profiles.
In a tool either resolution or depth of investigation has to sacrificed to gain the other.
– Small sensor means better resolution and poor depth of investigation and vice
versa.
– Large sensor spacing means poorer resolution but good depth of investigation.
The physics of measurement also dictate the resolution and depth of investigation
for some tools (specially radioactivity based tools).
SCALES OF MEASUREMENT
For the following plot, show the scale of processes to be investigated in a reservoir
characterization through the use of a bar scanning the appropriate ranges. Include in
your plots the following processes:
Diagenesis
Bedding units
Facies
Fracturing
Flow units
Oil/gas
Production
SOURCE OF DATA
Most professionals are familiar with the type of data required for reservoir
characterization. There are many sources for this data, but at times, the required
information is not available from measured sources. If it is necessary to estimate a
value or range of values, the professional is well advised to review the model
(geological, geophysical, engineering, etc. - each discussed in a later section). With the
model in mind, better interpretations can be made in selecting the values. There are a
great number of published data sources that may also be used. Satter and Thakur
(1994) provided a table that categorizes data (see below):
Production and Oil, Water and Gas Production Production and Production and
Injection Rates and Cumulative Injection Reservoir
Productions, Gas and Water Engineers
Injection Rates and Cumulative
Injections and Injection and
Production profiles
Data Considerations - 2
Distinguish between measured and interpreted data
Understand importance of distribution of range of values for reservoir parameters
Evaluate the error range in each data set
Evaluate the uncertainty in each data element
Communicate with data utility and limitations clearly indicated
SCALE OF DATA
Professionals from across the disciplines frequently use common terms to describe the
data sets that are generated within their discipline. Unfortunately, there is often a
communication problem in that the professionals are not describing the same data. For
example, geologists and geophysicists may each describe a formation boundary. The
geologist may well be thinking about a depositionally-defined unit with a scale on the
order of feet while the geophysicist is thinking of a unit defined by reflection
characteristics with a scale of tens of feet. Worthington demonstrates the scales of
investigation and scales of heterogeneity as shown in the two illustrations below:
SCALES OF MEASUREMENT
GAMMA RAY/OUTCROP/SEISMIC
SCALES OF HETEROGENEITY
Data Considerations -3
Field scale variations controlled by viscous forces - capillary forces negligent
Small scale - viscous and capillary forces interact with sedimentary structures
Viscous forces drive oil and along permeable laminae
Capillary forces force the wetting phase from hi to low K zones
Orientation of laminae very important
An additional aspect of dealing with scales is that of upscaling. This critical element in
reservoir characterization refers to the process of taking data at one scale and making it
applicable at another. This subject is treated in the section on geostatistics, but needs to
QUALITY/ERROR MINIMIZATION
To be most useful, the data used for reservoir characterization should fully represent
that for which it was measured. Without an understanding of the restrictions on the data
due to the many factors that influence it, that data set is subject to considerable risk -
particularly when those data are to be combined and massaged through several steps
in the reservoir characterization process. Field data is subject to many errors
(sampling, systematic, processing, random, interference etc.). Error and quality can be
understood better when questioned and subjected to testing and validation. Cross-
validation against other disciplines and (where possible) is a useful technique to validate
the data. Sprunt offers the following observations:
Inaccuracy can be removed by calibration. If data are accurate, but imprecise the
solution is to increase the number of data.
Care should be taken to distinguish between measured data and interpreted data.
Measured data are those that are directly observable and require no interpretations to
derive (ex. depth of a core sample); interpreted data require some sets of assumptions
before deriving the value (ex.. SW saturation based on log interpretation). Those data
that are interpreted are subject to potential error if the assumptions behind their
derivation are flawed or mis-applied.
A final important point is the manner in which values are chosen to represent reservoir
parameters. Typically, a single value is chosen to represent the entire reservoir or a
significant portion of it. The manner in which values are chosen is critical to the
reservoir characterization process. This topic will be addressed later in the section on
lithology, but is worth emphasizing here. Data forms provided as a part of the course
will assist in making these choices.
DATA MANAGEMENT
Access to data and availability of data sets across the disciplines are critical
characteristics of successful reservoir characterization projects. Part of the problem lies
in the data itself, as addressed above; in addition the framework on which
interpretations and manipulation of data are conducted have an impact. The
geosciences have typically worked on a UNIX platform and the engineers on a PC
platform. Bringing these two, platforms together or networking them--does not eliminate
gaps. To manage petroleum reservoirs more profitably, more efficient methods of
integration between PC and UNIX applications will become essential. Valusek offers
this observation:
Databases to integrate:
Base map information
Geological, geophysical
Petrophysical
Engineering, production
Within this hierarchical framework there are many tasks, sub-tasks, and iterations.
Fundamentally however, the process involves nothing more than a series of information
flows (i.e., the results of one task are used in the next). This process is highly non-linear
and iterative, and information frequently flows in loops, often more than once. The
discussion that follows is meant to be a general guide. The specific steps that one
should follow are highly dependent on the nature of the study.
LOADING PROBLEMS
Well Name
Surface Name
Logs Readings
Log Orders
Log Correction
Deviation Surveys
DATA TYPES
Core
Geophysical
Logs - all types (digital and paper)
Fluid
Perforations
Completions and workovers
Production and well test
Production
Petrophysical data
Normalize/calibrate/edit against cores
Identify cutoffs for mapping
Core data
Fully describe/interpret cores
Correlate with logs
Calibrate with test/production data
Seismic data
Develop structural interpretation
Seismic stratigraphic correlation
Establish correlations with reservoir parameters
LOG DATA
1. Must have good correlations
2. Edit data
Depth edit
Baseline shifts
3. Normalize data
Identify by type, contractor, vintage
Compare to core data
Calibrate to common response in same lithology
Normalize to field average and standard deviation
CORE DATA
1. Fully describe/interpret cores
Depth calibrate
Note location of special analyses
2. Correlate with logs
3. Develop understanding of ranges in parameter values
4. Calibrate with test/production data
Figure 2-3.
Permeability profiles on logarithmic (a) and linear (b) scales from plug and probe data
over the interval of DST-3 (2874.5–2882.5 m) in well 30/6-19. Note that the patterns of
permeability seen in the plug and probe data are comparable. Running average (13-
point) probe data are also compared with the core plug data (c) and show the probe
data to be slightly higher.
Permeability profiles on logarithmic (a) and linear (b) scales from plug and probe data
over the interval of DST-3 (2800.0–2809.0 m) in well 30/6-23. Note that the patterns of
permeability seen in the plug and probe data are comparable. Running average (13-
point) probe data are also compared with the core plug data (c) and show the probe and
plug data are similar.
Many reservoir characterization factors are common among several different fields; if
this were not the case, we would have to develop a unique investigation for every
property without the benefit of analogues. The question then follows: What is and is not
transferable from one field setting to another? The transferable items are the principles
and resulting processes that control each of the factors. Specific products (for example
- absolute permeability of an eolian sand) of the processes may or may not transfer, but
should be examined with great care and not applied directly without analysis.
Pressure Temp
Volumetrics
Fluid Types
Saturations
Continuity
Abs Perm
Rel Perm
Porosity
DEPOSITIONAL CHARACTERISTICS
External Geometry O O
Internal Geometry O O O O O
Pattern of Variation
O O O O O
(Quantitative Res. Properties)
Absolute Values
X X X X X X X X
(Quantitative Res. Properties)
DIAGENETIC CHARACTERISTICS
Compaction (including Geopressure) O O O
Cementation O O O O O
Dissolution O O O O O
Recrystallization O O O O O
TECTONIC CHARACTERISTICS
Folding O O O
Faulting O O
Fracturing O O O
FLUID CHARACTERISTICS
Hydrocarbons O
Formation Water O O
Injected Fluids O O
DATA - SUMMARY
Upscaling
SUMMARY OF CHAPTER 2
Data form the basis for the evidence that is produced by the RC process. The evidence
supports the decisions to be taken. In order to take the best decisions the critical data
elements should be acquired and in addition, the uncertainties in such data must be
specified.
Data must be fit for purpose and the costs of acquisition may not be greater than the
value of that data. Hence, certain data may not be acquired resulting in greater
uncertainties. The uncertainties in the data lead to the risk that is present when taking
decisions.
IN RC, the main question to be asked is what is the evidence that can be produced
based on the available data. Often there is no time or money to acquire additional data
to reduce uncertainties.
Some observations:
Solution:
Notes:
Size of solution space such that the least likely models on the low side (10%) and
the least likely models on the high side (10%) are excluded.
Reduce the size of the solution space (if necessary) by appraisal (data
acquisition)
The key challenges for organizations that want to characterize reservoirs accurately
involve several aspects:
Implied in the discussion above is not only a method to maximize the technical
assessment of a reservoir, but a requirement to communicate across disciplines to
create multi-disciplinary approaches that solve complex problems. Rather than work
sequentially within a specific discipline and then pass the final interpretation to the next
discipline, this new approach succeeds when an asset management team works
through the problem jointly. For example, geologists will talk to reservoir engineers in
order to develop a better understanding of reservoir dynamics from the existing test and
production data and then use that data to improve the depositional model; the
geophysicist may work with the petrophysicist to derive reservoir parameter indicators
from the geophysical data. To summarize, the work of developing models is integrated
and supported by feedback rather than sequentially passing the interpretations to the
next discipline.
Exploration concept
Research to confirm (geologic concepts, geophysical data gathering)
Drill test well(s)
Petrophysical analysis
Reservoir engineering analysis
Create development plan
Drill development wells
Monitor production
Each of the disciplines is involved throughout the project. For example, in planning the
3D seismic program the reservoir engineer may be consulted about the potential for
imaging reservoir compartments, the petrophysicist may offer suggestions that will
enhance attribute acquisition to be converted to porosity and permeability
measurements and the geologist could offer suggestions for design that would
maximize resolution of individual genetic units. This multidisciplinary approach to
solving the entire problem yields better answers that are more specifically focused on
the issues at hand. Strong potential exists for significant cost savings by focusing on
the specific needs of the other disciplines rather than gathering all of the data possible.
The following sections will discuss data interpretation techniques that have been found
effective in oil and gas development around the world. Even though each discipline is
discussed separately, the inference still exists that multidisciplinary approaches are
used throughout where applicable. This does not mean that all decisions or
interpretations need to be made in a team or group setting. The test for whether to use
a synergistic approach is simple - Does the involvement of other disciplines give a
better (cheaper or more comprehensive) answer?
MODEL TYPES
Static
Architecture of the field
Giga, Mega, Macro and Micro scale
Dynamic
All fluid flow within the reservoir
Production, Pressures, Saturations
Fluid contacts, Pore system model
Define facies
Establish depositional model
Define reservoir architecture
Well History
Production data/injection
Fluid data
One of the PetroSkills Alliance partners with the assistance of PetroSkills consulting
services has mapped the requisite skills for Reservoir characterization. These skills are
applied to each RC project and are a good checklist for workflows in other projects. The
skill groups and specific skills are shown below:
Reservoir Modeling
– Geocellular Model Development
– Data Analysis
– Geocellular Model Framework
– Facies Modeling
– Attribute Modeling
– Fracture Modeling
– Data Integration Methods
– Upscaling the Geocellular Model
– Model Validation
– Reservoir Simulation Support
– Integrate Dynamic Data
General Characteristics:
Lithology Porosity amt. Porosity Permeability
type
Interp. Method:
Net/Gross Net/Gross Water sol.
Interp. Method:
Texture Fabric
Interp. Method:
Reservoir Dimensions:
Modern Min ML Max Ancient Min ML Max
Length Length
Width Width
Thickness Thickness
Heterogeneities:
Stratigraphic Modern Min ML Max Ancient Min ML Max
Length Length
Width Width
Thickness Thickness
Control
Production Profile:
Analogue Expected
Critical Elements:
Analogue Expected
RESERVOIR PROPERTIES
(example only - repeat for each factor modeled)
Lowest value
ML Second mode
Highest value
Lowest value
ML Second mode
Highest value
SUMMARY OF CHAPTER 3
Each reservoir is unique, but analogues (other reservoirs, outcrops) are used to
understand a reservoir under consideration. The truth about a reservoir may never be
known. Therefore, it is necessary to model multiple realisations of the reservoir to obtain
a complete understanding of the possibilities that may be present. The realisations are
within the solution space, which covers 80% of all the models possible. The solution
space is defined such that reality is also encompassed by it. Development plans for the
modeled realisations will then be meaningful for development of the reservoir under
consideration.
In his excellent primer on “Geostatistics and Reservoir Geology”, Journel offers the
following description of geostatistics:
Modeling allows the geoscientist to predict reservoir occurrence in much more powerful
ways than previously possible. The capability to offer realizations of the reservoir with
many options built into a three dimensional presentation, bring capabilities that did not
exist previously. Geostatistics lends itself to these multiple realizations allowing effective
economic analysis. Geostatistics also serves as an excellent mechanism for the
integration of geological and geophysical data to generate a three dimensional model of
a reservoir that incorporates all available data; using the modeling process brings the
disciplines together to consider the possibilities that exist within the data sets and ways
to integrate the divergent data into a cohesive model.
The following figure shows one slice of a geostatistical model of a large Middle East
field. A huge amount of computation is needed to build such models. Efforts are now
underway to speed up the geostatistical modeling process to enable numerous model
scenarios to be tested.
INTEGRATED MODELING
MODELING CONCEPTS
MODELING STEPS
Definition of Layering Scheme
GEOSTATISTICAL APPLICATIONS
Data distributions
– Representing uncertainties
Geometry of parameters
Testing hypotheses
GEOSTATISTICAL OVERVIEW
The traditional approach to geologic modeling was define the reservoir through the use
of horizon maps, facies distribution maps and top and bottom structural maps to define
reservoir thickness and geometry; sand/shale ratios, porosity, permeability and other
properties are typically interpolated from well to well and displayed in probability
distributions (in the more sophisticated approaches). This approach works reasonably
well with fairly isotropic reservoir, but is lacking in sufficient detail and variability for
reservoirs that are more complex. Geostatistics provides a way for reservoir scientists to
more closely approximate the heterogeneous nature of the earth. It will model both
structural information and reservoir property data into a three-dimensional model which
may he used as a basis for a reservoir simulation study.
Object based
Spatial distribution of “geo-objects”
Conditioned to well observations and sequence stratigraphy
Sequential simulation (useful for selecting values for individual parameters such
as porosity or reservoir thickness)
Simulation in which the reservoir model is constructed by iterative trial and error
Probability field simulation that uses original well data to construct controlled
distributions of data
As indicated in the Bashore et al. study (see sequence stratigraphy section), selection
of the model can have significant implications on the modeling approach that is used.
When using a geostatistical approach to developing a reservoir model, geoscientists
distinguish between hard and soft data. Hard data are those data that are reproduced
exactly by all methods, while soft data are reproduced with some degree of tolerance
(Journel). Examples of hard data are well data including thickness of a reservoir,
porosity, and measured permeability; soft data are those that may replicate hard data
such as a seismic attribute that has a relationship to porosity - not a direct reading from
a specific location. Hard data are very desirable, but are much less abundant than soft
data. It is important, then, to establish relationships between data sets so that the more
pervasive soft data can be used to fill in for the less frequent - but more reliable - hard
data. Therefore, selection of the most appropriate model is extremely important to the
success of reservoir characterization.
Srivastava (1994) offers the following warning about selection of a geologic model for
the reservoir:
Several advantages to using detailed geological modeling are mentioned in Tyler et al.
Among them are the following:
Several realizations are available for studies of future predictions, contrary to the
traditional approach of a single simulation with the one probable model for
sensitivity studies.
Approach
Organize, present, and summarize a dataset of a single variable
Purpose
Quality control
Better understanding of the data
DESCRIPTIONS
Graphical
Histograms
Cumulative Histograms
Statistical
Measures of location
Measures of spread
Measures of shape
BIVARIATE ANALYSIS
Approach
Comparison of distributions for two different variables
Purpose
Investigate the existence and degree of statistical relationship between the
two variables
DESCRIPTIONS
Graphical
Crossplots
q-q Plots
Statistical
Correlation coefficient
Rank correlation coefficient
Linear regression
Conditional expectation
Graphical Analysis
Crossplots - Qualitative feel to how two variables are related; draws attention to
extreme (possibly aberrant) data pairs.
100
Y Data Value
Y Data Ranks
0 100
-100
-200
0
0 10 20 0 100 200
X Data Value X Data Ranks
q-q Plot
Straighter the line, the more similar the distributions; tails are of particular
interest
q-q Plot
200
100
Y Data Quantiles
-100
-200
0 1000 2000
X Data Quantiles
GEOSTATISTICAL METHODS
Sequential simulation
Variograms
Kriging
Cokriging
KED
GEOSTATISTICAL TOOLS
Uncertainty
VARIOGRAMS
Variograms - also called semi-variograms - are useful tools to assign values from points
of control to points between wells. The method can be used for many attributes and is
often used to develop permeability distribution throughout the field area. Stalkup offers
an excellent explanation of how to calculate a Variogram (emphasis is added):
Variogram Principle
As the distance between two points decreases, the correlation between data at those
points becomes larger.
Variogram Definition
The variogram is a function that measures one half the mean square error between data
at two sample locations a distance h apart.
At any distance larger than the correlation length, the co-variance becomes equal
to the overall variance of the data
Variogram
BUILDING A VARIOGRAM
Consider 4 different wells with coordinates (X,Y) and with a value of a particular quantity
(Z). Z may be net footage in the well, the porosity of a zone or any other quantity
relevant in the well. The distance between two wells can easily be calculated by
sqrt((X2-X1)2 + (Y2-Y1)2) and the semi-variance between two wells by ½ (Z2-Z1)2. In a
variogram, the semi-variances are plotted on the vertical axis against the distances on
the horizontal axis:
To remove the scatter (graph on the left), the data are binned (graph in the middle).
Then a variogram model (graph on the right) can be drawn.
The variogram model should extrapolate back to the origin. Sometimes, already over
short distances the variogram jumps to a finite value near the origin, which is called the
nugget effect. In addition, the distance from the origin over which the variogram reaches
its maximum is the range of influence or correlation length. For distances beyond this
correlation length, correlation between the wells is not meaningful anymore.
Furthermore, the variance in data may differ in different horizontal directions, which is
often the case for permeability (permeability anisotropy):
In this example, it is seen that the correlation length in the N-S direction is somewhat
greater than in the E-W direction.
KRIGING
Kriging is a mapping technique that estimates values for a parameter and is consistent
with the Variogram model. According to Wolf et al., kriging answers the question: “What
value should the grid point be assigned that best fits the Variogram model?” Kriging is
developed by mathematically calculating grid values as a weighted average of the
surrounding control points. Hence, the Kriging estimate is a linear combination of the
values at the surrounding control points. The weights become less when the distance
between the point to be Kriged and a control point is larger.
KRIGING PROCEDURE
CROSS-VARIOGRAMS
Cross-variograms are used to integrate “soft” data sets with “hard” data sets. For
example the data projections into the unknown areas can be improved over a simple
Variogram approach by using a seismic data set that correlates with the well properties.
Care must be exercised that the data sets have a correlative relationship. Cross-
variograms quantify the spatial relationships between the two data sets.
DEFINITION OF SEMI-CO-VARIANCE
COKRIGING
An improvement on the kriging process involves the integration of a second, denser
(“soft”) data set with the data to be gridded. Typically, this is a seismic data set that
correlates with the original data set. The second data set (the guide data) offers
weighting for the calculated relationships and honors the variograms, and cross-
variograms. Wolf et al. present the excellent graphical depiction of this process shown
below:
In this approach, a relationship is sought between the guide data (the seismic data) and
the well data. For instance, if the well data are net footages (isochores), seismic
amplitude at the well locations may be used to derive a relationship between peak
amplitude and sand thickness. The model used in the KED consists then of the
variogram of the well data and the relationship between amplitude and footage.
Practical applications have shown little difference between cokriging and KED. This
process is illustrated by Wolf et al. in the diagram on the next page:
CONDITIONAL SIMULATION
Rather than developing a single "best fit" answer as in the kriging or cokriging process,
conditional simulation offers a number of plausible solutions - each of which fits the
conditions described; each of these equi-probable solutions is a potential answer to the
conditions. Typically, a great number of solutions are derived, and statistical analyses
are performed. From these analyses, the probability of any occurrence can be
determined. For example it might be determined that 32% of the solutions to the sand
thickness at a specific location might be greater than 40 feet.
Wolfe et al. demonstrates a geostatistical approach to assessment. In their method, Realizations are generated for
thickness of sand at a prospect location, and, then ranked from low to high.
MAPPING TECHNIQUES
The seismic horizon equivalent to the sandstone – note amplitude changes along
the interval.
TWO KILOMETERS
80 FEET
40 FEET
-
0 FEET
The thicknesses were arranged in ascending order and plotted on a graph (see
following illustration). This is a compilation for sand thickness estimates at the location.
In the prospect risk assessment, there were many trials run to predict the thickness of
reservoir at the prospect drill site. These predictions are shown here in a histogram.
Thirty-five percent of the outcomes were below the economic minimum of 40 feet of
sandstone (this economic minimum was derived separately in a parallel evaluation
process). The risk for reservoir thickness in this prospect is therefore 0.35; this means
that if you drill a prospect like this, there is a 35% probability of failure (less than 40 feet
of sand). The adequacy for reservoir at this prospect location is 0.65 (1 – risk or 1.0 -
0.35). Keep in mind that when this prospect is drilled, only one thickness of sandstone
will be found!
S a n d D is trib u tio n
30
Economic
25 Threshold 40'
20
15
10
35% 65%
5
0
8
0
9
The following map simply displays the data generated in the last few slides at all points
on the map and computes the possibility or the risk of having less than 40 feet of sand.
One could also display the probability of having more than the economic minimum (a
more positive view). This probability is called the adequacy.
Random numbers are used to generate the bodies in the reservoir model. The shape
and the size of the bodies are determined by the rules, whereas the selected bodies
must be such that the well data (i.e. the control data) are being obeyed.
Application of this technique is demonstrated in the following pages for the Statford field
in the North Sea. In this field five different ‘valey fills’ can be distinguished, which were
modeled using geostatistical conditional simulation:
464
SOME EXAMPLES
N. Sea
one well with water contact, 3D survey
Stochastic depth maps predicted net/gross, SW’s, oil volumes
Oklahoma Morrow
Risk maps to predict presence of economic sand thickness
N. Cowden Ranch
7, 14, 55 well cases
Each accurately anticipated porosity trends
9.1
9.9
7.8 9.4
7.5
7.1
8.3
THE TRUTH
VARIOGRAM – 17 WELLS
400
Binned data points
300
Variance
200
Variogram model
100
0
0 3 9 12 18 27 36 45
30
15
Elevation
-15
-30
0.0 Distance 100.0
KRIGING SOLUTION
“The Truth”
30
Kriged Solution
15
Elevation
Oil/Water
0
Contact
-15
-30
0.0 Distance 100.0
Kriged Solution
30 Average of 100 conditional simulations
15
Elevation
-15
-30
0.0 Distance 100.0
SIX SIMULATIONS
ANDRE JOURNEL
“Different equiprobable numeric reservoir models can be built all honoring the
same hard and soft data”
Areas to improve
- model size,
- how to select simulated models from a set of equiprobable models for flow
simulation and development planning.
SUMMARY OF CHAPTER 4
Geostatistics is a technology which produces results that are in line with all the data
used in the analysis (all well data are obeyed, seismic data may be used to guide the
evaluation). As such results have become independent of the interpreter.
Geostatistics starts with analyzing the data in order to define variograms. Variograms
show the correlation length. Correlation beyond this correlation length is not meaningful.
Linear estimation techniques (such as Kriging, Co-kriging and KED) produce one single
solution, whereas simulation techniques (such as conditional simulation, Boolean
simulation and simulated annealing) produce many equiprobable solutions.
The application of geostatistics has grown over recent years and is nowadays widely
used.
Anisotropy Ratio: The ratio of the range in the direction of maximum continuity and in
the direction of minimum continuity.
Correlogram: A chart of correlation coefficients vs. lag separation distance (see “lag”
below). In general, the correlation between petrophysical properties separated by
some distance decreases as the separation distance increases. Given some
assumptions, the correlogram is 1.0 minus the standardized or relative variogram
(see below).
Correlation Scale: The distance scale at which the correlogram correlation is zero.
From a geostatistical modeling perspective, well data have no influence on modeling
beyond the correlation scale. The correlation scale is commonly called the range of
correlation.
Covariance Function: The covariance is a measure of correlation. Recall that the
correlation coefficient is the standardized coviariance, that is, the covariance divided
by the square root of the product of the variance of each variable. So the covariance
function is like the correlogram (see “correlogram”), except that it has not been
standardized or normalized to a correlation coefficient between -1 and 1.
Distance Classes: The variogram is calculated as the average squared difference in
the petrophysical property for a given distance. In practice, well data are not on a
regular grid and the calculation of an experimental variogram must consider some
tolerance (say, 100 m ±20 in). Such a range of distance is called a distance class.
Experimental Variogram: The ultimate goal of a variogram study is to arrive at a licit
three-dimensional (3-D) variogram model that must obey specific mathematical
properties to be useful in further studies. Experimental variogram points for specific
distance classes (see “distance classes”) are calculated from the available well data
before fitting such 3-D variogram models. They are called experimental in the sense
that they precede the determination of an ultimate variogram model.
Facies: In this paper, facies simply may be considered as different rock types or
groupings of data that share certain properties. For example, limestone and dolomite
would be considered two different facies. A variogram study would have to be
conducted within each significant facies.
γ (h) or gamma (h): See the text for a formal definition of the variogram. We should
note that γ (h) is technically the semi-variogram, whereas 2 γ (h) is the variogram.
Gaussian Variogram: Depending on the context there can be two meanings: (1) a
continuous petrophysical property such as porosity or permeability commonly is
transformed to a Gaussian or normal distribution prior to variogram calculation; the
experimental variogram in this case is sometimes referred to as a Gaussian
variogram or (2) more commonly, there is an analytical function that is used to fit
experimental variograms known as a Gaussian variogram model. This analytical
function is called Gaussian because it has an ‘exponential-squared” term, which is
the same as the Gaussian probability distribution.
Lag Increments: When the experimental variogram is calculated, the data pairs are
grouped into classes of similar distances. Any of these classes is also called a lag,
and the difference between the average distance of the data pairs in one class and
the average distance of the data pairs in the subsequent class is referred to as the
lag increment.
Lag Separation Vector: The vector that describes the change of direction and average
distance between one lag and the subsequent lag is called a lag separation vector.
Lag Tolerance: To rigorously define a distance class, the two parameters lag increment
and lag tolerance are used. If the distance of a data pair falls into the distance range
described by the lag increment plus or minus the lag tolerance, the data pair is
sorted into this distance class. The most common choice for the lag tolerance is one-
half the lag increment.
Nested Structure: The experimental variogram that is defined only at discrete
distances is modeled by an analytical function called the theoretical variogram. The
theoretical variogram can be composed of a sequence of individual analytical
functions each of which describes only the dissimilarity between data pairs within a
particular distance interval. In this case, the theoretical variogram is said to have a
nested structure. The term “structure” refers to each analytical subcomponents.
Normal Score Transform: The transformation of the given data so that their cumulative
distribution function corresponds to a standard Gaussian cumulative distribution
function with zero mean and unit variance.
Positive Definiteness: Each analytical function that can be used to model the
experimental variogram has to be positive definite. Positive definiteness is a
mathematical condition that ensures (1) existence of the solution of the kriging
matrix, (2) uniqueness of the solution, and (3) that the variance of any linear
combination of the data values will be positive.
Relative Variogram: The relative variogram measure is computed as one-half of the
squared difference between two data standardized by the squared mean of the data
used for the lag. This standardization distinguishes the relative from the traditional
variogram (see equation 1 in the text).
Sensitivity Study: If a model consists of several parameters that need to be
determined, it is useful to figure out those parameters that influence the model result
the most. Typically, the parameter values are varied, and the model response is
monitored. The model is most sensitive to parameters where a little change causes
a significant different model result. This procedure is called a sensitivity study.
Stratification: In this paper, stratification refers to the sequence of different fades in the
vertical direction. Inter-facies variation, in this context, is considered as the variation
of facies properties from one facies to the next fades.
Vertical Trend: A variable is said to have a vertical trend if the variation of data within a
neighborhood (vertical direction) can be described as a smoothly varying function of
the coordinates.
The recent geologic literature is full of examples of sequence stratigraphic models that
incorporate geologic models with the high quality of recent seismic data. Even though
emphasis has focused on sequence stratigraphic approaches, the basics were
developed in the last century through the work of Wagner and many others that
followed. It is not necessary to develop new and burdensome terminology to explain
the facies distribution in a sedimentary sequence.
SCENARIO TREES
Data
Fully
Interpretation
Correlatable Meandering Distributary
Channels Channels
Realization (facies)
Realization (K, )
(Facies related distributions)
Flow Simulator
2. Define facies
FACIES
STRATIGRAPHY
OBJECTIVES OF STRATIGRAPHY
Define distribution of reservoir and non-reservoir facies
Establish correlations
We can recognize sequences based on Outcrops, Seismic, Cores and Well data
EMPHASIS
Facies Analysis - emphasizes internal character of strata.
Unconformity
Correlative Conformity
Unconformity
Correlative Conformity
WALTHER’S LAW
HISTORY
Subsidence/Uplift
Accommodation
Sediment Supply
Base Level
Transit Cycles
Volumetric Partitioning
Facies Differentiation
Our Tools
1. Seismic interpretation of stratigraphy and depositional systems – 2d and 3d
2. Core and well-log interpretation of stratigraphy and environments – 1d stacking
patterns, 2d-3d correlation and mapping
3. Modern analogs – depositional systems, sediment patterns, bathymetry, physical
processes, object models, 1d-3d
4. Ancient analogs (outcrops) - stacking patterns, stratal patterns, object modeling,
1d-3d
What More Do We Need?
1. Better understanding of depositional systems (processes and products) –
modern/ancient
2. Improved understanding of geobody shapes, dimensions, composition, etc. –
catalog
3. Improved understanding of stacking patterns
4. Improved recognition of lithofacies from indirect detection – logs (image), seismic
5. Improved confidence in correlation strategies (refine our understanding of
sequence stratigraphy)
6. Higher frequency seismic data to improve resolution
E+T=S+W
Eustasy minus Tectonic Subsidence equals Sediments plus Water depth.
DEFINITIONS
Hiatus:
Fall in relative sea level or base level or
Change in base level rise or
Climate or tectonic processes in source area influencing erosion and clastic
influx
Parasequence: Basic building block of sequences (Van Wagoner et al., 1988)
Characterized by aggradation, progradation, and retrogradation
Affected by sediment supply and minor base level fluctuations
Accommodation space:
Space available for potential sediment accumulation (Jervey, 1988); space
available below base level (Shanley and McCabe, 1994)
Provided by tectonic subsidence and eustasy
PARASEQUENCE-STACKING PATTERNS
Summary diagram showing for each third-order sequence the characteristics of the
sedimentary system, the sequence stratigraphic model, and the typical reservoir, source
rock, and seal facies. The differences between these sequences are largely due to the
dominant controlling factor: eustatism for sequence I, environmental change through
clay influx for sequence II, and tectonism for sequence III. HST = high-stand systems
tract; TST = transgressive systems tract; LST = low-stand systems tract.
1. Construct a tentative correlation for this sequence that is prograding to the right.
Hint: Draw time lines that might result from the depositional processes prior to
interpreting lithologic distribution.
Typically, support for the interpretation of the sequence stratigraphic model comes from
the geological depositional models, study of analog environments (including outcrops
and modern environments), seismic studies and reservoir engineering data (principally
well test data). The distribution of the reservoir and non-reservoir facies provide the
basis for many of the interpretations, but it is essential to include an interpretation of the
time distribution of the facies in order to develop a complete reservoir characterization
model. Numerous interpretations of the physical data are possible, but an integrated
approach including data from other disciplines and a chronostratigraphic approach will
give a truer representation of the reservoir.
Fluid flow interpretations results indicate that the geologic model had little affect
on the breakthrough time or oil recovery, but fluid displacement patterns may
be substantially different. Correlation strategies significantly affect how the
seismic (soft) data are applied within the model and therefore have significant
affects on developing continuity relationships when interpreting porosity and
permeability distributions. They indicate that an integrated model may prove
crucial for solving swept and unwept areas, developing infill drilling programs,
conformance control and designing tertiary recovery processes.
BIOSTRATIGRAPHY
GR
GR GR
GR
Flooding shale - datum
GR GR
Flooding shale - datum
The following figure is a schematic longitudinal section along axis of the Gironde incised
valley illustrating facies and stratigraphic patterns of the transgressive fill.
SUMMARY OF CHAPTER 5
Sequence stratigraphy provides a time framework for correlating and mapping
sedimentary strata. In this work, correlation lengths needs to be considered. Since in
any stratigraphic evaluation, there are uncertainties, it is essential to maintain multiple
working hypotheses and build multiple different models.
An understanding of the difference between time units and rock units assists in the
definition of flow units.
Lithology
Net/gross ratio
Porosity
Hydrocarbon saturation
Hydrocarbon types
Diagenetic alteration
Facies distribution
When coupled with the analyses of the facies derived from sequence stratigraphic
studies, analyses of the presence of barriers to flow, interpretation of the physical
characteristics of the reservoir facies and consideration of structural implications, flow
units can be defined. These flow units are the basic building blocks of the model which
will predict reservoir behavior. Cross describes the process well:
Cross in the same paper discusses the usage of analog data in the application of facies
models form other areas. He makes several significant points:
Like all processes in reservoir characterization, this process shows a feed-back loop. A
reservoir characterization is never finished. It is always important to test concepts and
modify them where necessary, e.g. when new data become available.
The steps in the process are vital steps required to define reserves. An elegant way to
arrive at reserves resulting from static evaluations is applying Monte Carlo simulation.
With Monte Carlo simulation volumes in place and reserves can be calculated in a
probabilistic way taking into account the uncertainties in all contributing variables
captured by probability density functions.
The Monte Carlo method is highlighted in the figure on the next page:
Monte Carlo simulation (probabilistic modeling) showing inputs with their uncertainties, the various stages
in the 3D modeling process and the corresponding outputs. Final result is the exceedance probability
curve of reserves.
TRAPS
The following figure illustrates several types of traps that might be encountered.
LITHOLOGIC COMPONENTS
Lithology
Net/gross ratio
Porosity
Hydrocarbon saturation
Hydrocarbon types
Diagenetic alteration
Facies distribution
DEPOSITIONAL ENVIRONMENTS
OUTCROPS
Outcrops are essential for obtaining statistical data on sand bodies. In the following
picture the channel sand body is clearly visible:
Reynolds (AAPG) produced an excellent paper providing many statistical data for sand
bodies.
A log-log plot of width vs. thickness for all available data with sandstone
bodies differentiated by sand-body type. Sand-body types show clear
clustering and only limited overlap of dimensions. Two lines recording
thickness to width ratios of 1:100 and 1:1000 are shown for reference.
(l) Valleys are much wider than channels. Crude averages show that channels
(distributary channels, crevasse channels, etc.) are narrow, less than 1 km in width, with
aspect ratios of 1:100. Whereas valleys average around In km in width and have aspect
ratios on the order of 1:1000.
(2) Shoreline-shelf sands are huge sheets tens to hundreds of kilometers in length, with
mean widths that range from 7 to 25 km and vary according to systems tract. Shoreline-
shelf sands deposited in highstand systems tracts are, on average, twice as wide as
those deposited in transgressive systems tracts. Thicknesses of shore-line-shelf sands
relate strongly to their position within a sequence set.
(3) Valleys and shoreline-shelf sands have areal extents comparable to giant oil fields.
Flood tidal deltas and mouth bars are comparable in area to small fields.
(4) Distributary channels crevasse channels, tidal flats, and crevasse splays are small
areally.
(5) Some systematic trends occur (e.g., in flood tidal deltas), allowing the prediction of
width and length from thickness. Other data sets show a high degree of scatter, but
maximum, mean, and minimum values can be determined for width, length, and
thickness. There is a clear partitioning of certain sand-body types into certain systems
tracts.
Typical shapes of sand bodies are shown in the next two diagrams:
ALLUVIAL FANS
The relationship between particle size and velocity of the water can be used to
understand the deposition of particles in the various environments.
An alluvial fan deposit occurs on the edge between the mountains and the continental
plain. At the transition, a sudden change in velocity happens:
Using the relationship, a range of particles that can not be carried any more by the
lower velocity water will be deposited. Hence, the sorting in this deposit is relatively poor
implying a reduced permeability.
RIVERS
Two types of rivers occur, namely braided rivers and meandering rivers. The braided
river is usually somewhat higher in the mountains, whereas the meandering river
meanders across the continental plain.
The river banks of the braided river are relatively strong and well defined causing the
river to flow along a relatively well-defined path. The meandering river can meander
over an area of several kilometers wide. Consequently, the deposition in the braided
river is localized whereas deposition in the meandering river can be found spread-out
over the pain as visualized in the following diagram:
The velocity of the water in the braided river is fairly constant, whereas for the
meandering river the velocity in the inner bends of the river where the deposition of
material takes place, is reducing with time. As a result, the sorting in a braided river
deposit is very good, whereas in the meandering river deposit a fining-upward sequence
is found.
For a braided river, the distance of transport of the particles has been relatively short so
that the size of the particles is still relatively large. This leads to high permeabilities,
which may as high as 2 D. In the meandering river deposit the quality of the deposits is
reducing when going up implying deteriorating permeability when moving up (fining-
upward sequence).
– River-dominated
– Tide-dominated
– Wave-dominated
When a river flows into the sea or ocean, the velocity of the water is gradually reduced.
As a result, the finest particles are transported further away from the river mouth so that
larger particles are stacked on top of smaller particles. In this way, a coarsening-upward
sequence is deposited. Also when particles are picked up by the waves and transported
e.g. along shore (see Galveston picture), first the finest particles are picker up and
transported over the largest distance leading again to coarsening-upward sequences.
Consequently in coastal deposits, the best developed reservoirs are found near the top
of deposits.
In summary:
Eolian deposits are found in deserts. The material is picked up from beaches by the
wind and transported over sometimes very large distances to the place of deposition.
Cap rocks of Eolian deposits are often transgressive shale in combination with salt
deposits.
TURBIDITES
Some material is carried across the continental shelf and dropped into deep water
forming turbidites.
The material has been transported over a long distance. Hence, turbidites consist of fine
often shaly material and have a thinly-laminated character. Nevertheless, very good
production can be obtained from turbidites. Bonga offshore Nigeria, the first turbidite
development offshore Africa, produces 240,000 bbls per day.
Turbidites are often found in deep water and therefore, often developed using FPSOs.
An excellent paper was presented by Weber and van Geuns (JPT), who proposed a
categorization of reservoirs in three different types so as to distinguish their complexity.
These types are:
Labyrinth reservoirs
LAYER
–CAKE RESERVOIR TYPE
In summary:
Jigsaw K-contrast 2 to 4
Labyrinth Inter- 10 to 30
connectivity
Tools that can be used for definition of the different types are:
Seismic
- sequence stratigraphy
- 3D seismic, VSP
- well-to-well seismic (possibly)
Wells
- logs and cores
- borehole imaging tools
- RFT / MDT
Production
- production tests
- well interference tests, pulse tests
- tracer tests
Outcrops
- shale length
- san body characteristics and size
- depositional models
Recovery from clastic reservoirs may vary between 5 and 80% depending on many
variables. In the next diagram, some statistics are shown depending on depositional
environment and drive mechanism.
The following figure shows three channel belts present in part of the Colville River flood
plain, Alaska. Note the juxtaposition of the narrow channel belt of the Kogosukruk River
(left) formed by a single, sinuous channel and the multiple, sinuous-to-braided channels
in the much wider Colville River channel belt (center). Active and abandoned channels
and bars are easily discernible. To the right is another abandoned channel belt. Each of
these channel belts could be identified in the subsurface using core and wireline-logs,
and accurate predictions of their dimensions could be made. Photograph from July 1979
in the National Petroleum Reserve Alaska, approximately 40 km (25 mi) northeast of
Umiat.
And the next figure shows simplified plan forms, cross profiles in various positions and
orientations, and large-scale stratal geometries of channel belts for braided and
unbraided sinuous rivers. Braid bars and point bars migrated downstream and
expanded laterally. Large-scale inclined strata are shown schematically on cross
profiles, and upper bar and lower bar deposits are distinguished. The total thickness of
bar deposits varies laterally and is generally less than the maximum bankfull depth.
(A) Two-dimensional stratigraphic simulation of fluvial channel belts (stippled) in the Travis
Peak Formation, North Appleby field, East Texas basin (zone I of Tye, 1991! Davies et al).
(B) Wireline logs and channel-belt sandstone bodies of zone I as originally correlated by Tye
(1991). (C) Revised well-to-well correlations using same data as in (B), but constrained by
recalculated channel-belt widths. Note that limiting channel-belt widths decreased sandstone
body continuity in two dimensions.
Note that width/thickness ratios of sand bodies of some of the meandering systems can
in places be relatively small, because of the thick packages of cohesive sediments in
which these systems develop. Apparently lateral accretion in meandering systems can
be restricted under such conditions, just as in anastomosing systems. Note also that
multistory sand bodies have a significantly different w/t trend than single-story sand
bodies.
After Allen
SEDIMENTARY
FACIES
ENVIRONMENT
A. Fine-grained sandstones Fluvial channels
A2: Course - V. coarse sandstones Braided Stream
B. Rooted siltstones Paleosol
C. Very course - to fine - grained sandstones with clay drapes Upper Estuarine channels
D. Parallel laminated black shales with climbing ripples Central Estuarine bay
E. Flaser and wavy bedded sandstones and siltstones Restricted tidal flag
F. Inclined deformed sandstones and siltstones Tidal channel
G. Laminated calcareous mudstones Lower estuary
H. Poorly to moderately fossiliferous planar cross-bedded Estuarine mouths
sandstones and pebble conglomerates
I. Highly fossiliferous planar cross-bedded very coarse- to medium- Upper shoreface
grained sandstones and pebble conglomerates
J. Rarely to moderately burrowed planar cross-bedded medium- to Proximal middle shoreface
fine- grained sandstones
K. Moderately to thoroughly burrowed rippled fine-grained Distal middle shoreface
sandstones
L. Thoroughly burrowed fine- to very fine-grained silty sandstones Lower shoreface
and fading ripples
M. Thoroughly burrowed very fine-grained silty sandstones and Offshore transition
sandstones with interbedded normally graded sandstones
N. Thoroughly burrowed siltstones Offshore
O. Parallel laminated black shale Shelf
LAGOONAL LITHOFACIES
Typical vertical sequences of lithofacies from different parts of sandy channel bars and channel fills.
Idealized gamma-ray logs also given.
kH,max
Measurements of paralic sandstone bodies. (A) Sandstone bodies were measured at the parasequence scale.
An individual parasequence may yield a number of measurements in this example three channel
measurements (1 to 3) and a shoreline-shelf sandstone body measurement (4). (B) In plan view in each case,
width is the minimum sand-body dimension and Length is the maximum sand-body dimension. (C) Sand-body
thickness data have been collected by measuring maximum thicknesses. A randomly placed well penetration is
likely to encounter a reduced thickness, closer to the mean thickness. By contract, in sheet sands a randomly
placed well is more likely to penetrate a thickness close to the maximum. (D) Multistory and multilateral channel
sandstone bodies, respectively, record numerous stacked and adjacent channel positions. At a parasequence
scale, multistory and multilateral channels appear to be rare in paralic successions, and as a result, it was
possible to collect data from individual channel sandstones.
Models of major depositional environments. The curve on the left shows the SP or gamma ray response and
the curve on the right shows the relative grain size profile. The size of the dots next to the vertical profile
indicates the relative magnitude of permeability expected in such a sequence. (Parts c and d are from Walker,
1984, and parts f, h, and I are from Galloway and Hobday.)
The following case study is a study using the Paria river outcrop. The Paria River is a
contributory to the Colorado River in the US. The study simulated a water flood under
various permeability distributions.
The study objective was to understand the effect of permeability heterogeneity on the flood
front.
(a) The exposure of the Paria River a sandstone incised valley fill used in the reservoir simulations. (b) A
lithological interpretation is made initially, before this interpretation is corrected for any perspective distortion
developed when acquiring the outcrop photomosaic. (c) The corrected panel is then gridded into 661 × 81
rectangular simulation cells of 0.74 × 0.61 cm. Corrected interpretation is used as input for Eclipse reservoir
simulations. The Paria River photomosaic is a section through the A sandstone incised valley fill unit,
perpendicular to paleoflow, that is, paleocurrent was into the picture.
Water saturation plots at 0.25 pore volumes injected for four realizations of the full model with an injector-
producer pair injecting to obtain an interstitial velocity of 12 cm/day where sand is (a) homogeneous or has
a stochastic distribution with a volumetric average permeability of (b) 238 md, (c) 280 md, or (d) 182 md.
Arrows indicate the location of wells and direction of flow within.
Water saturation plots for the same permeability distribution plotted at 0.25 pore volumes injected where
(a) heterolithics have been replaced by shale and where sand flow properties have been assigned to (b)
heterolithics, (c) shales, or (d) both heterolithics and shales. Arrows indicate the location of wells and
direction of flow within.
The results of this study lead to a very important conclusion that needs to be realized
when upscaling models for reservoir simulation purposes
A second case study involves data from the Lower Jurassic Kayenta outcrop. In this
study, the objective was to evaluate the contribution of very low permeability (~ 1 mD)
conglomerates on the recovery efficiency.
Ephemeral-fluvial sandstones form the main reservoirs in many oil and gas fields.
Production histories and development difficulties suggest these deposits are more
complex than is commonly supposed in the petroleum industry. There is inadequate
published information to account for this observation. Relevant sedimentological
literature on such deposits is qualitative and unsuitable for detailed reservoir
characterization studies. Therefore, integrated sedimentological, petrophysical and fluid-
flow simulation studies have been conducted on outcrop analogs in the southwestern
United States. Our key objectives were to define flow units and to find areas of
sensitivity inherent both in the sediments and in the analytical and statistical procedures
used to describe and model such systems.
Model 1
Each cell assigned mean permeability value
Model 2 - Assess influence of mudclasts on flow
All conglomerates impermeable (Cgu, Cgc)
Model 3A - assess permeability range in facies
Stochastic permeability model - resembled #1
Model 3C - A.A., but set conglomerates to 1 md
Model 4 - Utmost realism
Added spatial variability to sandstones
Note under-running in basal sequence
DIAGENESIS
It is easy to see that diagenesis in carbonates can take place in the depositional
environment; it is not necessary to bury it or expose it to start this process.
Type/Description
I. Rocks with pores capable of producing gas without natural or artificial fracturing.
II. Rocks with pores capable of producing gas with natural or artificial fracturing
and/or interbedded with type I rocks.
III. Rocks too tight to produce at commercial rates even with natural or artificial
fracturing.
Step/Action
1. Estimate grain size and sorting using standard size-sorting comparators, thin
section and SEM photomicrographs, and rock photographs.
Lateral Heterogeneity
Low Moderate High
Wave dominated Delta-front Meander belt,
delta, Barrier mouth bar, Fluvially
core, Barrier Proximal delta dominated
Low
shore face. front, Tidal delta, Back
Sand-rich strand deposits, mud- barrier
plain rich strand
Eolian, Wave- Shelf bars. Braided
modified delta Alluvial fan, Fan stream, Tide –
Vertical
Moderate (distal) delta. dominated
Heterogeneity
Lacustrine delta
delta, Distal
Basin-flooring Coarse-grained Back barrier,
turbidites meander belt, Fluvially-
Braided delta dominated
High
delta, Fine-
grained
meander belt
After Galloway et al.
CARBONATE RESERVOIRS
CARBONATES SUMMARY
Clastics Carbonates
Deposition after transport Deposition in situ
Climate no constraint Sediments occur in tropical areas
Worldwide at all depths Shallow water
Both marine and terrestrial Marine
Grain size related to distance of transport Grain size reflect size of skeletons
DEPOSITIONAL MODEL
Under low-stand conditions and in a humid environment, a rimed platform may undergo
karstification due to widespread subaerial exposure. This may lead to subterranean
drainage and the development of caves near the water table. At the shelf margin, low-
stand reefs and/or grain shoals may form. The shallow-marine carbonate factory is
greatly reduced in size. Although slope failure and the generation of sediment gravity
flows can occur during any sea-level position, some triggering mechanisms may favor
low-stand conditions. Debris-flow deposits may form lobes or aprons at the toe-of-
slope, whereas turbidites may extend further out onto the basin floor as apron or fan
lobes. Handford and Loucks.
Generalized depositional models for Clear Fork and lowermost Glorieta formations, in North Robertson unit,
Gaines County, Texas.
Fossils
Pellets
Sand-size, spherical or ellipsoidal particles
Formed by deposit-feeding organisms that eat mud.
These organisms digest organic matter and excrete the non-digested mud as
pellets.
Ooids
Composed of Aragonite
Nuclei in the center
Concentric or radial
Well Sorted
Intraclasts
Large rounded grains
Angular sometimes
Internal structure
Reworked carbonates
Limestone
Dolomite
LIMESTONE
While most carbonate grains are easy to recognize in hand specimen when seen in thin
section they can be difficult to identify. The problem is, as Majewski (1969) remarked,
that a "variety of shapes are produced by random cuts through a single geometric
pattern, and shapes can be duplicated in such cuts by different designs; also,
characteristic features may be obscured and others may become apparent" when the
grains are exposed in a specific plane. Never the less the carbonate grains of this
collection are separated from one another on the basis of their shape, size and internal
structure.
Carbonate grains can be separated from one another on the basis of their shape, size
and internal structure. Because the grains commonly collect near their site of origin,
they can be used, in conjunction with other rock characteristics including vertical and
lateral facies relationship and sedimentary structures, to determine the depositional of
the rocks they occur in. Information about grain types and the manner in which they
occur in rocks can be communicated by means of limestone classifications.
Carbonate sedimentary particles may be subdivided into micrite (lime mud) and sand-
sized grains. The grains can be separated on the basis of their shape and internal
structure and can be subdivided into two major groups: skeletal and non-skeletal.
It is still used to describe the fabrics in core today. Note some clear limitations in
reservoir description for petrophysical responses; limestone versus dolomitized fabrics;
and challenges in characterisation of interpartical porosity present in packstone (Lucia
1999). In this context, the figure indicates where REEF fits in addition to more familiar
Clastic split of SAND and MUD - grain size separator for carbonate sands and muds -
so muds BOTH very shallow in lagoon and very deep in the basin/slope
The figure below shows the interaction of geological processes - deposition and
diagenesis - in determining the microgeometry of the present-day samples studied for
the Khuff and L-III samples. The column on the left shows the initial sediments, their
fabrics are differentiated by the four small cartoons: ranging from grainstone (top)
through packstone, wackestone and mudstone (bottom).
After deposition the sediments are altered by the physical and chemical processes of
diagenesis. For these cases diagenesis proceeded along one of four paths. The fabric
which results (take a close look at the cartoons) may be very different from the fabric at
the time of deposition and, for a given rock, depends on the depositional class and on
the diagenetic path.
There are many more paths than this. There are even subpaths and variations within
the paths within the case studies, but the simplification represented here is sufficient to
make the point that diagenesis plays a major role in determining microgeometry and
fabric.
In path 2 the dolomitization process simply replaced the original material retaining the
depositional fabric, including pore space both within and between grains. Compaction
and cementation significantly reduced the porosity.
In path 3 the primary porosity was reduced by cementation. There was a later
dissolution phase, which preferentially dissolved some of the original grains and
resulted in the formation of moldic porosity. Compaction and cementation significantly
reduced the porosity.
In path 4 the limestone was leached by a process that did not follow the depositional
texture. The leached channels which result cross grains and matrix.
Geological Processes.
The microgeometry classes which we have recognized in these case studies are shown
here in the context of their depositional and diagenetic history.
Most geologic schemes for carbonate classification emphasize the solid phase of the
rock – types of fossils and mineralogy – seeking to predict the reservoir scale
distribution of rock textures through an understanding of the depositional and diagenetic
processes. In many cases depositional class impacts diagenesis, but this is by no
means consistent.
Unfortunately existing geologic schemes are not optimal for the purposes of
petrophysics. Knowing only the depositional fabric does not specify the present day
microgeometry. Knowing in addition the diagenesis still does not (in general) specify the
microgeometry well enough.
The following figure shows cartoons of idealized end member fabrics recognized in case
study samples. Each of these fabrics is associated with a particular microgeometry
class in our case studies. However there may be other fabrics, which we did not
encounter, corresponding to the same microgeometry classes in other formations.
Petrographically we found these fabrics to be dominant, however other fabrics do occur
on the thin section scale. The cartoons are drawn to help visualize the NMR and flow
properties of the samples associated with the microgeometry classes in our case
studies and to help geologists relate to microgeometry classes.
DOLOMITE
Dolomite = Dolostone
Magnesium replace Calcite
Dolomite Texture
The distribution of the size of crystals may be:
Recrystallization in Carbonates
Anhedral
Crystal is not bounded by crystal faces
Fabric is named Xenotopic
Subhedral
Crystal is partly bounded by crystal faces
Fabric is named Hipidiotopic
Euhedral
Crystal is wholly bounded by crystal faces
Fabric is named Idiotopic
Idiotopic dolostone is one of the most important carbonates reservoirs for oil and gas. It
is named as Sucrosic Dolostone.
From: Dehghani, et
al, 1999, AAPG
A comparison between historical field production data (hist.) and the results from a scaled-up permeability model
not considering vuggy zones (Model). Water-cut (bbl of water per day/total of oil and water) and oil rate (stock tank
bbl/day) are plotted for 60 yr of production history. The model results do not match the water-cut and oil rate
histories.
A comparison between historical field production data (hist.) and the results from a scaled-up permeability model
with superimposed high-permeability vuggy zones; in this case 4 d of permeability was added to the secondary
porosity cube using a longer correlation length of 400 ft (122 m) (4SP2). Water cut (bbl of water per day/total of oil
and water) and oil rate (stock tank bbl/day) are plotted for 60 yr of production history. The model closely matches
the oil rate, as well as the water cut.
SUMMARY OF CHAPTER 6
Since data are often sparse, modeling of reservoirs is significantly improved if
knowledge about the environment of deposition is known. It is important in this modeling
to test concepts and to modify models where necessary. Particularly, when new data
becomes available re-modeling becomes necessary.
Reservoir complexity can be described using concepts such as layer cake, jigsaw
puzzle and labyrinth. In this order, reservoir development requires more wells and a
denser well spacing.
After deposition, diagenesis starts to play a role causing either porosity increases or
porosity decreases.
Aims of the petrophysical analysis are to develop models that can be used to project
data across the study area. The data that are available for critical analysis is
concentrated in a very small area when considering the breadth of the reservoir. It is
incumbent, then, to be able to project data from the well bores to the large areas
between the wells.
WHAT IS PETROPHYSICS
Petrophysics is the science of measuring rock properties and the study of their
relationships to each other.
Petrophysical data
Normalize/calibrate/edit against cores
Identify cutoffs for mapping
Core data
Seismic data
Develop structural interpretation
Seismic stratigraphic correlation
Establish correlations with reservoir parameters
RC PETROPHYSICAL PROGRAM
Open Hole Log Evaluation
NMR Logging
Reservoir Surveillance
Define facies
Do I have a reservoir?
Do I have hydrocarbons?
What is the type of hydrocarbons?
Can we produce it (permeability)?
How much is there (porosity and saturation)?
3D reservoir description
Calibrate with production tests, pressure transient tests, decline curves
Product: Input into simulation model
It is clear in this approach that integration is vital in petrophysics. ‘Define the rock types’
must be done in close cooperation with geologists and sedimentologists, ‘Integrate with
formation evaluation’ is an area where petrophysicists, geologist and reservoir
engineers meet, and for “Calibrate with seismic and/or geostatistics’ petrophysicist,
seismologist, seismic intepreters and geologists need to come together.
Petrophysics is involved in data acquisition. A major part of this data acquisition is done
in borehole by logging (while drilling or after drilling). The issue is that measurements
are done in an environment (the borehole) that is different from the virgin formation, for
which data are required. In addition, borehole fluid invasion complicates the situation:
The properties are measured in the borehole, whereas the conditions in the uninvaded
zone need to be known. The challenge is:
Steps:
Wireline Logs
PETROPHYSICAL PERSPECTIVES
Borehole Seismic
The scale of petrophysical data varies from the micrometer scale up to as much as the
100 of metres scale:
At all stages the measurements assist to define the static and dynamic models Whereas
often the static and dynamic models are referred to as separate models, it must be
emphasized that they form one single reservoir model. Changes in the static model
impact on the dynamic model and vice versa. The following diagram provides a
summary:
The combination of density and neutron logs can be used to distinguish between the
major facies. In addition, the combination assists in distinguishing between oil and gas
zones. When the neutron and density logs are overlain while the scale on the neutron
log is reversed, a large separation is shown between the logs in a gas zone, whereas
this separation is small in a oil zone in fact being marginally larger than in a water zone
(see the diagram on the next page).
It should be noted that in gas condensate reservoirs, the density of the fluid in situ is
such that hardly any separation is visible. There is a risk that in such a reservoir the
hydrocarbon is interpreted as oil. A production test will of course resolve this issue.
However, if under circumstances no production test can be done (say the rainy season
started early), there is a log evaluation technique which will distinguish between the
different options of oil or gas. This technique will be demonstrated below.
The density and neutron logs are a combination of two logs. More logs can be used in
cluster analysis, which is in fact a ‘crossplot in more than two dimensions’.
Fluid densities can also be evaluated using pressure versus depth plots.
CLUSTER ANALYSIS
In the following case study (Tucker et al, AAPG) six different log responses are used to
distinguish between different lithologies, called clusters. The logs used are:
Gamma Ray
Density
Neutron
Sonic
Resistivity
Photo-electric effect
As six logs are being used also six cluster can be distinguished. The cluster need to be
defined using core data. The geologist is needed to describe the core and the
petrophysicist is needed to interpret the logs and to add the fluid data, again a very
integrated job.
Once the clusters have been defined and the model is calibrated using the core data,
the model can be used to simulate a core in other wells in a reservoir. In the
development phase when many wells (several hundreds of wells is easily possible) this
technology provides a means to evaluate lithologies and facies in detail.
Note that permeability can be related to grain size distributions and sorting.
Sneider’s procedure for this purpose is given below.
Qualitative:
Is there a mud cake? Mud cake is being built up against permeable zones.
From micro-resistivity logs (they measure porosity in the invaded zone, and
hence permeability; low resistivity implies high permeability).
Remember:
We will look at SEM (Scanning Electron Microscopy) pictures of rock samples that have
approximately the same porosity, but significantly different permeabilities. The pictures
are followed by Sneider’s pore classification system and a cross-plot of the K - φ data.
Sneider’s procedure to estimate permeability from grain size distributions and sorting
Sneider and King developed a simple method of classifying pore types from cuttings.
The classification of clastic rock pore types from cuttings is made by comparing pore
types with production tests and log analysis. The pore types are as follows:
Type/Description
I. Rocks with pores capable of producing gas without natural or artificial fracturing.
II. Rocks with pores capable of producing gas with natural or artificial fracturing
and/or interbedded with type I rocks.
III. Rocks too tight to produce at commercial rates even with natural or artificial
fracturing.
Step/Action
1. Estimate grain size and sorting using standard size-sorting comparators, thin
section and SEM photomicrographs, and rock photographs.
Ghawar
Discovered in 1949, Oil column 1300 ft, Np = 80 billion bbls (2011), Potential = 5-6
million bbls/d, pressure maintenance by water injection, simulated with a 1.032 billion
grid block black oil simulator.
In the picture, the colors vary from red (oil) to blue (water)
ACTIVATION LOGGING
The generated gamma-rays, which have an energy typical for the elements they
interacted with, are counted by the tool at least those that reach the detector. Tools can
be designed to measure specific gamma-ray energies or a complete gamma-ray energy
spectrum. The investigation depth is a couple of inches. The tool can be run both in
open and in cased hole.
The responsibility for the data acquisition using this technology lies in Petrophysics, but
the technology has applications in many areas of the business. Hence, close
cooperation between disciplines is again required to make optimum use of this
technology.
Production logging
Activation of silicon and aluminum. The detector is placed below the neutron source and
the log is run slowly so that oxygen and other activated elements have decayed before
the detector passes the activated interval. A high count rate indicates a high quantity of
silicon in a sand pack or aluminum in a bauxite pack.
A somewhat fancier application is oxygen activation to detect and quantify the flow of
water in and around the borehole. Oxygen can be activated to become a nitrogen
isotope, which upon decay (half-life only 7.1 sec) produce gamma-rays of 6.13 MeV.
There are detectors above and below the source to monitor flow upwards and
downwards.
formation. The combination of induced and natural spectroscopy logs is often referred to
as the geochemical log used in formation evaluation.
Gamma-ray spectroscopy is applied in gas- and oil-shales to determine the TOC (total
organic content) of the formation.
Evaluation method
• Method 1
- total porosity from density log
- mineralogy from spectroscopy log
- carbon content from mineralogy
- excess carbon content = carbon content minus carbon associated with minerals
- TOC is proportional to excess carbon content
• Method 2
- effective porosity from NMR log
- total grain density (ρgr) from density logh and NMR porosity
- matrix density (ρm) from spectroscopy log
- TOC volume fraction of matrix = (ρm-ρgr)/(ρm-ρTOC)
- Toc weight fraction = (ρTOC/ρgr) x TOC volume fration
Barnett Shale
Well 1 Well 2
FORMATION MICRO-SCANNING
On the images sinusoidal features are recognized, which highlight planes that intersect
the borehole. The higher the amplitude of sinuses the greater the relative dip between
the plane and the borehole.
Applications are:
In the following two displays FMS images are compared with layering seen in outcrops
showing that indeed the FMS provides similar information as what can be seen in
outcrops.
FRACTURES
One of the most difficult aspects of characterizing reservoirs is how to deal with fractures.
Some reservoirs are totally controlled by fractures, others are unaffected and most are
influenced to some degree. Recognition and prediction of these fractures is difficult at best –
although modern imaging logging techniques are major improvements. Antonellini and
Molemna have published an important paper that analyzes fracture systems in carbonates
at each level of occurrence – regional, fault-based, local fractures, brecciation, micro-
fractures. RC teams might find a similar approach useful in their reservoirs.
FMI image, 3.3 m (10 ft) in length, shows fractures and bedding planes. Image log with interpretations is shown as
folded image shown as borehole cylinder in 3-D space with fracture planes (a) and unfolded image with sinusoid
fracture curves (b). Basic interpretations include fracture dip, dip azimuth, and depth.
The above example is characterized by the different fracture spacings in the different
limestone layers (which one would expect for different layer thicknesses). This fracture
spacing is a key parameter to model along with fracture orientation.
Borehole image logs provide important sub-seismic images. Current data analysis
methods, however, do not fully use the information contained in these images. Instead
of using borehole images to extract a three-dimensional (3-D) visualization of fracture
networks, most interpretations are limited to rose diagrams, stereonets, and histograms,
thus ignoring the rich 3-D information contained in the image logs. The ignored
information can be used for representation of
There are other borehole image logs in addition to the FMS. They have been widely
used for detecting fractures and have become the most important tool for imaging sub-
seismic features. These logs include electric, acoustic, radar, and ultrasonic methods.
Oriented cores also provide firsthand, reliable fracture data, and efforts have been
made to integrate data from core and different imaging tools to validate image
interpretation. Core, however, is commonly of limited extent compared to borehole
images, and is not oriented in most cases. Therefore, fracture orientation is not always
measurable in cores. A new logging while drilling (LWD) density log ROSI provides a
complete 360 degrees image, but its resolution (16 pixels/circumference and 2-10
pixels/ft) is still not yet adequate for most fractures.
In the wellbore
Production logging
Corrosion inspection
In the formation
Reservoir monitoring
CORES
So far, we saw how cores were used to obtain mini-permeameter logs, which provide a
continuous profile of permeability versus depth. Such profiles need to be calibrated with
routine permeability measurements on one inch core samples.
In addition, cores were used to calibrate the model in cluster analysis. After such
calibration the cluster model can be applied in other wells to simulate cores in uncored
wells.
In general, core data are used to calibrate logs. Porosities from density logs are
calibrated using core porosities. Resistivities can only be translated into hydrocarbon
saturations after the saturation exponent (n) and lithological exponent (m) have been
measured on core samples. Such measurements represent in fact the calibration.
Obviously, for good calibration, cores need to be correlated with logs. This is done by
correlating gamma ray logs run during the logging operations and those run when the
cores were taken.
Systematic observations from cores starts with completing a list such as:
Cores are also used (in combination with logs run in boreholes) in studies to determine
depositional environments. In the geology session, we saw already how the gamma ray
log is used to identify fining-upward and coarsening-upward sequences representative
for meandering fluvial deposits and coastal deposits. In braided river environments we
were expecting fairly constant gamma ray logs showing little radio-activity.
In this chart supplementary or special core analysis is mentioned. Here, three of them
are discussed in more detail, namely relative permeability, wettability and capillary
pressures.
Realtive permeability
Wettability
Cores brought to the surface are made dirty by the drilling fluid. In the laboratory, core
samples cut from cores are cleaned by cooking them in a solvent, e.g. toluene. After
this treatment, the core samples are fully water wet. Then, they are saturated with
formation brine, and subsequently they are aged by flooding the hydrocarbon from the
reservoir through the cores until residual oil saturation is reached.
The question is whether by this treatment the in situ wettability will be restored? And in
addition, is it possible that in situ the rock was oil wet and that the original wettability can
be restored.
In the literature there is an ongoing debate. Some authors claim that they encountered
reservoirs that are oil wet, whereas others believe that all anomalous reservoir
behaviour can be explained without having to introduce oil wetness.
Oil wet reservoirs can be recognised on logs if high resistivities leading to oil saturations
of over 98% are observed, or during the production phase, if unexpected early
breakthrough of water and anomalously high watercuts are observed.
Capillary pressure
Capillary pressure data are used to determine hydrocarbon saturations as a function of
height above the free water level (FWL). This application is particularly useful in thinly
layered reservoir for which the height of the thin layers is less than the resolution of the
logging tools. Clearly, core is required for this application.
SUMMARY OF CHAPTER 7
Petrophysics plays a central role in data acquisition across the entire hydrocarbon life
cycle. As such, petrophysics is a key player in integration.
When data are acquired, the logging environment must be taken into account.
Corrections have to be made or the impact of the logging environment on the data
needs to be specified.
Fluid determination can be done using logs and pressure versus depth measurements.
Formation micro-scanning may replace cores after the FMS images have been
calibrated with core. Application are in thinly bedded formation and in fracture detection.
Petrophysics also plays an important role in both routine and special core analysis.
Special core analysis include saturation en lithological exponents required for resistivity
interpretation, relative permeability and wettability, and capillary pressure curves.
1. THE PROBLEM
– Declining liquid rates
– High operating costs
– Increasing injection gas costs
3. Assignments
1. Designing an approach to solve problem
2. Developing common data base with consistent measurements
3. Adequately representing inter-well heterogeneities
Suggestions:
Consider the support material
Clarify the problem to use as a basis for planning the study
Choose the style for the team
Discovery
The Anschutz Ranch East field was discovered in 1979 by Amoco Production Co. with
the completion of the No. 1 Bountiful Livestock well. The well spudded 29 March 1979
on a farmout lease from BWAB, Inc. and Tom Brown, Inc. and was completed 31
December 1979. The discovery well penetrated 940 ft (286.5 m) of gross pay in the
Jurassic Nugget Sandstone in what is now the West Lobe of the field. The East Lobe
was discovered by The Anschutz Corporation, Anschutz Ranch East No. 12-26W (26-
13~121W, Uinta County, Wyoming, completed 23 December 1980. The two lobes
produce gas and gas condensate from the Nugget Sandstone and have separate gas-
water contacts. The West Lobe discovery well is perforated from 12,816 ft (3909 m) to
13,515 ft (4122 m) and had an initial flowing rate of 1054 BCPD, 4053 MCFD, and 28
BWPD. The East Lobe discovery well is perforated from 14,750 ft (4499 m) to 14,910 ft
(4548 m) and flowed 141 BCPD, 750 MCFD, and 96 BWPD (Lelek, 1983).
The No. 1 Bountiful Livestock well was drilled to test for the possibility of an anticlinal
fold in the Nugget Sandstone as interpreted from seismic data and inferred from
shallower folding in the Champlin Petroleum Co. No. 404 "A" No. 1(35-13N-121W, Uinta
County, Wyoming).
Discovery Method
Approximately 25 mi. (40 km) of conventional 2D seismic and shallow correlations to a
nearby Cretaceous Bear River Sandstone test (Champlin
Petroleum Co. No. 404, Amoco "A" No. 1) were interpreted to define the Anschutz
Ranch East prospect and locate the discovery well. The Champlin No. 404 well was
used to tie shallow seismic reflections, and it revealed an overturned fold in the Jurassic
Preuss Formation. The presence of folding in the Preuss, stratigraphically above the
Nugget-Twin Creek interval, and the seismic interpretation suggested the possibility of
an anticlinal fold at the deeper Nugget level (Lelek, 1983). A 3D seismic survey
including 40 mi. (102 km) of 12-fold data was acquired in 1981 to help detail the
structure for further field development.
Post-Discovery
Shortly after the field's discovery it was determined that the West Lobe of Anschutz
Ranch East field contained a rich retrograde gas-condensate with a dew point only 230
Psia (1586 kPa) below the original reservoir pressure. Pressure depletion below the
dew point would have resulted in a condensate liquid dropout of up to 40% of the
hydrocarbon pore volume (HPV). Such a dropout could have led to a 75% loss of the
reserves in place. Because of this, a specialized plan-of-depletion was designed to
optimize production while maintaining the reservoir pressure. This plan called for
nitrogen generation and injection to replace the hydrocarbons produced from the
reservoir.
The original injection plan in the West Lobe utilized an inverted nine-spot pattern and
proposed a 10% HPV "slug" of injected buffer mix containing lean residual field gas and
no greater than 35% nitrogen. Subsequent pressure maintenance has resulted in the
total injection of approximately 50% of the hydrocarbon pore volume with similar
nitrogen and field gas mixtures. In 1986, the injector layout was modified to an
"imperfect line-drive pattern" to take advantage of permeability trends inherent in the
eolian reservoir sandstones. The continuing challenge is to maintain a constant sweep
of the reservoir to minimize hydrocarbon loss and to maintain field-wide pressure above
the dew point. Extensive coring (approximately 12,000 ft [3658 m. of core), core
analyses, and reservoir modeling programs have been undertaken to better predict the
field's response to various depletion schemes.
Currently, the West Lobe of Anschutz Ranch East field contains 30 producing wells and
18 nitrogen-gas injection wells. Production from the West Lobe is expected to continue
into the early part of the twenty-first century.
In the East Lobe, sufficient differences exist between the original and dew point
pressures to tolerate partial pressure maintenance (cycling) and pressure depletion
alternatives. Furthermore, the amount of reserves contained in the East Lobe cannot
economically justify full pressure maintenance. The East Lobe contains seven
producing wells and is entering into the pressure depletion phase of production.
Significant drilling problems were encountered during the development of the field. Salt
in the Jurassic Preuss Formation often caused hole problems and casing failures. Early
in the field's development, a method was devised to drill and simultaneously under-ream
the salt to minimize salt encroachment around the drill pipe (Holt and Johnson, 1985).
More recently, oil-based drilling fluids have proved effective for controlling the thicker
sections of salt. As a result of these improvements and of lower day-rate rig and
subcontractor costs, drilling expenses are now half of those in the early1980s.
Rock Properties
The Nugget reservoir is anisotropic and inhomogeneous with respect to fluid flow. Many
of the reservoir properties are tied to the original depositional architecture and only
slightly modified by subsurface physical and chemical diagenesis (Lindquist, 1983).
The Nugget reservoir in the West Lobe can be effectively divided into three
Petrophysical zones. The upper two-thirds of the formation contains large-scale eolian
dune sandstones and is the best reservoir rock. The lower third of the interval is
characterized by smaller-scale eolian sandstones, some with water-influenced
depositional textures. This lower zone comprises the poorest Nugget reservoir rock in
the field.
The average matrix porosity of the Nugget reservoir rock is 12% with the highest
observed porosities approaching 22%. The average reservoir matrix permeability is
10.9 MD with some permeability’s ranging as high as 400 md. Correlation quality of the
Kmax plot (reservoir porosity plotted against maximum permeability parallel to the
bedding plane) and the Kh90 plot (reservoir porosity vs. permeability perpendicular to the
bedding plane) for this data set are 31% and 51%, respectively. This high variability
makes reservoir modeling predictions difficult.
Lindquist (1983) gives a detailed description of the Nugget eolian sandstones and their
relationship to reservoir petrophysics. Zone 1 is characterized by laminated sandstones
that are interpreted as stacked climbing wind-ripples and grainflow deposits in relatively
small cross-bed sets. Zone 2 is the best reservoir interval with generally higher
porosities, thicker cross-bed sets, and more grainflow deposits. Zone 3 is the poorest
Fractures probably influence the directional permeability within the field. Core
descriptions commonly identify gouge-filled fractures in Zone 1, only minimal gouge-
filled fractures in Zone 2, and numerous open fractures in Zone 3. The distribution of
fracture morphology is similar to that discussed by Nelson (1985), in which higher
porosity sandstones commonly shear and create a low-permeability gouge. Fracturing
in the lower-porosity sandstones more often creates open fractures that can enhance
reservoir permeability. The fracture sets in the Nugget Sandstone at Anschutz Ranch
East field are interpreted to be of dominantly tectonic, fold-related origin with generally
predictable senses of motion or offset, as set forth by Nelson (1985). There is also a set
of open fractures oriented at low angles to the structural axis that are thought to be
related to the period of post-thrust extension.
Partners
The field is owned by six companies with differing interests throughout the field.
Regulation
The field straddles the Utah-Wyoming line. Each has jurisdiction for a portion of the
field. The field is located in a fairly sensitive environmental area.
Oil and gas companies are using an increased number of three-dimensional (3D) seismic
surveys in reservoir management because the detail needed to develop reservoirs is
increasing greatly. The amount of data requisite to find an accumulation is insufficient to
develop that same accumulation. Robertson (1991) describes the reservoir
characterization process and the application of 3D seismic to the resolution:
The first thing to note is that this definition is not geophysical or geological;
in fact, it's not even an engineering definition. Reservoir management really
is the economic process of raising the worth of a property to the highest
possible level….
He continues:
A 3D seismic survey is one of the tools in the evaluation tool kit. An initial
interpretation of the survey impacts the original development plan. As
subsequent events like the drilling development wells occur, the added
information is used to revise and refine the original interpretation. Often, as
time passes and the database builds, elements of the 3D data that were
initially ambiguous begin to make sense, and the interpretation becomes
more detailed and sophisticated.
WAVE PRINCIPLES
For the lithological rock column, an acoustic impedance (ρ x v) trace is made. At layer
boundaries, where ρ x v changes, acoustic waves are reflected and transmitted depending
on the contrast in ρ x v. Reflection coefficients are spikes at the layer boundaries. This
reflection trace is convolved with the seismic wavelet yielding the seismic trace. Because the
wavelet is not a spike but rather an acoustic disturbance of finite duration, averaging takes
place, which becomes visible in the seismic trace.
ACQUISITION
In the old days, seismic lines were acquired onshore using dynamite as the energy
source. Today, most onshore seismic data is acquired using vibroseis, where the energy
source is a giant vibrator truck in which a pad in the middle of the truck rests directly on
the ground to raise the truck so that the pad supports the full weight of the vehicle. The
truck then vibrates to put energy into the ground.
This technology was developed to acquire seismic e.g. in housing areas, where dynamite
was of course prohibited. However, because the energy transmission into the subsurface
is very controlled with vibroseis, the technique is nowadays also being used in remote
areas such as deserts.
Vibroseis truck
The main advantage of 3D surveys is that “out-of-plane” reflections can be recorded and
treated properly. In 2D surveys it is assumed that all the sound from the subsurface
comes from a vertical plane underneath the line of receivers. This assumption is OK
when dips are relatively small. In faulted reservoirs when reflections against the faults are
also recorded, the assumption breaks down completely. 3D seismic solves this problem.
In the next two pictures a 2D and 3D seismic survey are being compared.
Modern seismic acquisition technology includes three component seismic including shear
waves, wide azimuth seismic allowing for instance viewing underneath salt packages,
OBC (Ocean bottom cables) avoiding the negative impact of the sea water layer,
wireless acquisition offshore using ROVs (Remotely operated vehicles) placing acoustic
receivers on the sea floor, and wireless acquisition on land e.g. in rain forests minimizing
the operational footprint in pristine environments.
PROCESSING
In the two examples above, it is clearly seen that the Pre-stack Depth Migration or
Imaging is by far to be preferred.
In PSTM the migration is done in time, the domain in which seismic is recorded.
Basically, this means that a constant velocity is being assumed for the entire subsurface,
which is not the case. In PoSDM, the migration is done after stacking. Indeed stacking
removes noise, but also part of the signal. In PSDM, every trace recorded is considered
separately, which removes the negative effects of the stacking. The term ‘imaging’ is
used when a forward modeling approach is applied while iterating until the constructed
seismic image resembles the recorded image sufficiently.
GEOPHYSICAL APPLICATIONS
INTERPRETATION
An excellent study was carried out by Gartner and Schlager (AAPG), who constructed
seismic sections for a synthetic example using different frequencies:
Nowadays seismic interpretation is done on workstations that allow full integration with
the business of other disciplines. The example shown is a display from Landmark’s
‘Decision Space 3D’.
This integration on workstations allows inclusion of all data in a consistent and uniform
way including comparison of seismic with geological models, assessment of drained
areas, well planning, etc.
STRUCTURAL INTERPRETATION
SEISMIC CUBE
After processing of 3D seismic, the ‘seismic cube’ is obtained. The seismic cube is an
areal array of seismic traces with time in the third dimension. These traces may be
compared with synthetic seismograms obtained from well data, if at a particular areal
location a well was drilled.
The cube can be sliced in any vertical dimension to create vertical sections or sliced in
horizontal dimensions to create time slices, which represent constant time. In addition,
sections through the seismic cube may be created following a particular loop in the
seismic, providing a horizon slice possibly representative for say the top of a reservoir.
Robertson offers the following warning: “minor character changes in 3D seismic data
don’t need to correlate with real geologic changes”. It is, therefore, critical that a multi-
disciplinary approach to reservoir characterization includes interpretation of the seismic
data. The greater the understanding of the reservoir model characteristics such as the
expected Lithology, stratigraphic sequences, continuity and structure the better the
resulting interpretation. This requires detailed communication across the disciplines.
The next example compare a 2D and a 3D interpretation. The 2D is smooth, whereas the 3D
shows a lot more detail. This is not necessarily noise, but a result of the much denser
acquisition of the 3D survey. In the 2D interpretation much more averaging is required.
The next example shows the interpretation of a reef. At various times, time slices are cut
from the 3D cube displaying the variation of the reef with time, or after a time-to-depth
conversion with depth.
The last example shows a subsea structure map of the top of a reservoir sandstone. The
map was created by multiplying the time structure map, as interpreted from the 3D
seismic data set, by the velocity gradient map:
The color transition from white to green represents the pool’s original gas/water contact.
The size of the grid is 400 x 400 ft.
SEISMIC COHERENCE
For each point in the 3D seismic volume, the wafeform of adjacent traces (e.g. red trace
compared to the blue traces) over a short vertical time window:
SEISMIC ATTRIBUTES
In addition to seismic coherence, there are many other seismic attributes, such as
amplitudes, phase changes, frequency variations, dip and azimuth data, interval travel
times between events, velocity ratios (vp/vs). Typical information to be derived from these
attributes include the presence of gas (bright spots) from amplitudes, lateral continuity
and bedding configuration from phase information, bed thicknesses from frequency data
using spectral decomposition, and as we discussed, faults and fractures and other
stratigraphic information from coherence.
Integration of the seismically derived parameters with other data is often a difficult
procedure - especially if the desire is to be able to incorporate hard data rather than soft
data. The geostatistical section above offers some suggestions on methods to integrate
data sets. Chambers et al. provide the following observations:
The integration of secondary data sets, such as acoustic impedance data from the
inversion of 3D seismic data, can significantly reduce interwell estimation uncertainty.
The soft data also add small-scale variability not sampled from well data alone. Extra
care is required in the data processing and inversions steps if seismic data are to be
quantitatively integrated into a reservoir description, especially if the final product is a
fluid flow model. To realize the full potential of the secondary data set, data types must
be properly correlated and scaled.
Loren et al. offer methodology that is applicable to areas of minimal log control. The lack
of control may be due to a lack of enough wells drilled in the area or due to old logs with
inadequate information. The method they suggest is:
Araktingi et al. offer a comprehensive review of methods of integrating seismic and well
logs. This presentation offers a primer for those who want to apply seismic data sets in
the reservoir characterization process. Included in their paper is a compilation of
limitations, problems and pitfalls. Significant points include:
As a result of these limitations, extreme care and in-depth understanding of the geologic
models are very important to deriving appropriate data from the seismic data set.
processing.”
In the following a number of examples are presented:
As a result, when entering a gas zone there is a strong reflection yielding high amplitudes
compared to a much smaller reflection down structure in the water zone.
In the seismic section, this effect may be observed. Such a reduction in amplitude when
moving off-structure is often accompanied bij a flat reflector indicative for a gas-water
contact. Schematically, this is shown in the next picture:
The next areal map is a field case from the Imo River M1.0 reservoir in Nigeria (courtesy
Shell). The red color represents the high amplitudes. Noticeable in this map is that the
amplitude anomalies conform very well with the fault pattern in the reservoir.
The next example compares a synthetic seismogram derived from well data with the
seismic section. Using the well data a stratigraphic interpretation of the section can be
made:
The next picture shows the Spectral Decomposition of a 50 ms window centered on the
"D" sand of an example reservoir. The 30 Hz slice appears to be the best at imaging the
"D" sand.
To obtain the best results from attribute analysis, a combination of attributes should be
used. The following field case from the C5 Project (Crave Norte area, Colombia,
September 2011) is an excellent example using the horizon slice itself (Picture A),
seismic coherence (Picture B) and spectral decomposition (Picture C).
The next example compares a stratigraphic section with a geologic impedance model. In
the impedance picture some of the internal structure of the reef van be observed.
4D SEISMIC MONITORING
Time-lapse reservoir monitoring
Fluids alter impedance in reservoir
Change in fluid content changes impedance
Differences may result from fluid changes
Technical challenges:
Matching datasets
Compatible processing
Positioning
Quantification
Good reservoir model
Good analysis of the Value of Information (VOI)
Staging 3D seismic surveys over time to monitor a producing reservoir can determine the
position of fluids being moved by production and can predict future fluid movement more
accurately. Monitoring the reservoir in this way helps:
Baseline survey across a field prior to production is the first step in the 4D process.
Operators in mature areas may use seismic reservoir monitoring to prevent premature
abandonment. The technical concept behind 4D is to relate changes in seismic
attributes with time to changes in the reservoir due to production. 4D will not work for all
fields: it depends on the rocks, the pore fluids, reservoir thickness, structure and depth,
and data quality and consistency. 4D seismic is often combined with other forms of time-
dependent information. Other 4D datasets can include pressure histories, temperature
monitoring from producing wells, cased-hole worker logs such (pulsed neutron) and
reservoir simulations of the expected vs. actual behavior of the field. Companies now
expect to remotely sense the drainage pattern in subsurface reservoirs. Time-lapse
images of producing reservoirs are acquired through the repeated acquisition of 4D
seismic surveys to track hydrocarbon movement with time. The promise of identifying
bypassed and poorly drained reserves are very encouraging. 4D seismic truly
represents an integrated approach that brings together disciplines that were separated
from each other almost entirely in approaches that are more traditional.
After the 3D seismic survey of 1992 was acquired, there were only three
wells active in the fault block: new wells B-5ST, B-6ST in Block 330 and
the A-12 in Block 338. A new 3D seismic survey was acquired in 1994,
and it clearly shows drainage, or dimmed amplitudes, caused by
production in the intervening years. About 2.4 MMBOE were produced
from these wells, and about 600 acres were dimmed. Assuming a 100 foot
pay sand with 30% porosity and a drop in Sw, of 30%, that would require
that each acre-foot produce about 400 bbl, quite reasonable for clean
Pleistocene sands in the GOM.
The bypassed oil remaining after 1994 will require the placement of yet
another well into the fault block, and therein lies the power of 4D seismic
monitoring. Each new well recovers an additional percentage of oil-in-
place, increasing the ultimate recovery efficiency. In addition, the 4D
seismic changes can be modeled within a reservoir simulator to risk the
new well. Amplitudes drained by the new well are predicted to account for
an additional 2 MMbbl for a well placed horizontally along the 330/338
property line to the four corners. In another two years, that well will water
out, and additional amplitudes are predicted to remain after that. Then a
well projected more north-south across the property line will be required to
recover an additional 2 MMbbl, and so on into the future.”
4D SUCCESS FACTORS
Best applied in soft unconsolidated reservoirs
Works best when gas phase is appearing or disappearing
Oil / water transitions more difficult to image
This study was carried out by S. Njerve and I. Mæland of Statoil, Norway. The study was
presented at the 74th EAGE Conference in Copenhagen, Denmark, 4th to 7th June 2012.
SUMMARY
The decision of sub-sea redevelopment of the Visund North field that was closed in since
2006, was made. GRM was to be applied to update existing models. By thorough cross-
disciplinary use of production data and geological and geophysical interpretation of 4D
seismic, the Visund Nord dynamic reservoir model was updated. Based on the new
model, two new oil producers east in the area and an additional exploration pilot were
planned. Drilling started in Q3 2012. It is likely that another two wells will be drilled from
the sub-sea template in the future.
- Closed-in in 2006.
As the reservoir simulation led to conservative recommendations (new wells too close to
existing wells), the project remained marginal and it was decided to do a full GRM study
in order to determine whether integration with the seismic data would make the project
feasible.
The seismic database consisted of four different surveys in 1990/1991, 2004, 2007 and
2010. An example of the quality of the seismic can be seen from the following dip line
(blue line in map view):
In the difference picture, the anomaly caused by the water flooding is clearly observed.
Since the oil production started in 1999, the flooding indicated by the seismic has
developed from north to south. Between 1991 and 2004, the strongest flooding has been
mapped in the north, while between 2004 and 2007, the anomaly has moved further
south. This flodding pattern was confirmed by production data from the two previous
wells, of which the most northern one observed high watercut quite early.
Since the height of the flooded part of the oil rim is far below the seismic resolution, 1D-
modeling of different flooded cases was used to understand the tuned 4D signal and to
quantify the movement of the OWC. The principle of this 4D tuning works as follows:
The two stratigraphic columns are convolved with a wavelet to produce two synthetic
seismograms, of which the difference is taken (note the polarity difference in the two
synthetics).
With this principle synthetics for different levels of the OWC were constructed:
Using these synthetics, relationships were derived between the OWC movement (m) and
time difference (ms), and secondly, between the OWC movement and the amplitude ratio
(amplitude at time t divided by amplitude of base) difference:
Detailed analysis of the various time-lapse seismic surveys led to maps showing the
flooding development as derived from the seismic. These maps were the main input from
the geophysicist to the reservoir engineer:
The objective of the history matching was to match the produced volumes, the observed
reservoir pressures and the flooding maps obtained from seismic. The match was
achieved by varying the size of the gas cap and of the aquifer, by adjusting the
transmissibility across the faults in the model and by varying the vertical communication.
The results of this reservoir simulation are shown in the next picture:
The base case for redevelopment of Visund Nord is to drill new oil producers in the north
and the south away from the identified 4D seismic anomalies in the area. Although the
4D affected areas are not considered to be completely flooded, it is expected that
production east of the previous wells will drain the remaining reserves without too high a
watercut.
In addition, an exploration well will be drilled on the eastern side, where hardly any 4D
anomalies have been observed.
There are two more slots on the sub-sea template, which may be used in future for two
additional development wells.
In this process, iteration through the various steps may be required, which makes a good
integration between geophysics, geology and reservoir engineering more than vital.
PetroSkills is most grateful that Statoil made the above study available for the benefit of a
broader audience. The main authors, S. Njerve and I. Mæland, are commended for an
excellent piece of work. Their integrated approach to 4D reservoir modeling
demonstrates best practice, which PetroSkills wishes to recommend to all in the
integrated reservoir characterization arena.
New acquisition techniques allowed denser lateral sampling and better vertical
resolution
– Better resolution from the newer data images layers 40 feet instead of 80 feet
– Lateral spacing of 50 feet instead of 100 feet
Structural or stratigraphic details that explain an amplitude remaining isolated from
producing wells
In the picture on the next page, the original oil water contact (OOWC) is shown. The
brighter colors represent areas with hydrocarbons. It is seen that the flooding of the
reservoir can be well observed. In the lower left corner of the picture, remaining
hydrocarbons are observed:
WELL SEISMIC
CROSS-WELL SEISMIC
Crosswell data are collected by placing a seismic source in one well and a receiver string
in a nearby well. Energy which propagates directly between wells without being scattered
serves as the basis for constructing velocity images (tomograms). Energy which is
reflected is used to construct reflection images.
Cross-well seismology has the capability to define reservoirs with much more detail.
Paulsson demonstrates this capability using computer tomography modeling. The
process can add detail not available previously and has implications in identifying non-
continuous reservoirs, location of pitchouts and details of the interwell structure. The
following illustration is a model of such an application:
The cross-well reflection and velocity images (center) fill the resolution gap between
modern surface seismic data (left), and the sonic log data (right center) and core
measurements (far right). These data are from a West Texas carbonate reservoir.
A crosswell acquisition system developed by Texaco EPTD using an air gun source
produces greater low-frequency output than the GRI/TomoSeis RCPTM piezoelectric
source used in high-speed, commercial operations. The air gun and RCP (resonant
cavity piezoceramic) sources are similar in many ways: both are fluid coupled and can be
operated "on-the-fly," firing while the source is moving. As a result, GRI, Texaco EPTD
and TomoSeis worked together to identify a gas reservoir application in which both
sources could be operated at well spacings greater than 40 acres
Example - McElroy Field, W. Texas –Tucker, et al, 1998 AAPG. S-Wave images between
wells 1068 and 1202.
This example is from the Clark Wilson sandstone in the Mist field. This reservoir was
selected for a gas storage project.
The following figure is a vertical seismic section from 3-D survey parallel to 2-D seismic
section of figure 8-44. The top of the Clark and Wilson Sandstone (C&W) has a strong
trough (red) amplitude anomaly. There is a “relict” amplitude anomaly near the downdip
edge of the original gas-water contact (GW). There is no amplitude anomaly associated
with the “new” gas-water contact (NG/W).
Cross section through a Mist gas pool after primary production. Note the position of the
“new” gas-water contact. Discovery well on left; injection/withdrawal well on right. Vertical
exaggeration 1.7-to-1.
The next figure is a vertical seismic section parallel to the path of a "horizontal” well.
Looking at many sections like this, and at time slices and flattened time slices, is
invaluable when planning a well path. The amplitude anomaly continuity gives a "feel” for
a variability of internal reservoir stratigraphy or structure.
PROJECT ECONOMICS
The NPV is US$ 10.8 million. OK, a small project but profitable.
The following is a flowmap of modeling procedures showing the feedback loops with 3-D
seismic processing (via 3-D velocity model) and 3-D seismic interpretation (volume
analysis, dine-depth conversion). As a final step, 3-D synthetic seismic can be used for
quality control of the reservoir model This process also allows a quantitative assessment
of the likely probabilities to be made for the different modeling scenarios.
Input Modelling
1-D Well Data
The above figure is a schematic drawing showing a composite section that illustrates the
tectonic and sedimentary history of the Malampaya/Camagopla~formations. The
Oligocene part of the Nido Limestone overgrows a preexisting rug geomorphology for
nearby top-ended clastics deposits prior to the Nido Limestone. Note that the Eocene nor
mal faults along the western margin of the Malampaya~ buildup were reactivated during
the middle Miocene transgressional movements, causing tectonic inversion Some
reverse fault movements also are observed along the eastern margin of the platform. The
figure is not to scale; distances between the wells are given for reference.
DIAGENESIS
This figure shows the subsidence and relative sea level history of the
Malampaya/Camago platform as derived from Sr isotope stratigraphy (absolute ages)
and carbon isotope measurements (suggesting repeated exposure). High resolution of Sr
isotope dating in the Oligocene—early Miocene allows recognition of time gaps of 0.5—
2.0 m.y. during the growth history of the platform. Constrained data points are indicated
as circles on the curve. Note the change in tectonic subsidence during the late Rupehan
with rates changing from about 10 m/m.y. to 2O m/m.y. (not compensated for absolute
eustatic sea level change).
The above figure shows thin-section photomicrographs illustrating the reservoir rock types
(RRTs) for the Nido Limestone in the Malampaya/Camago area. Five RRTs have been
defined based on thin section analysis, prosily, permeability, and saturation-height
functions. Note that RRTs do not necessarily follow the classical texture-based carbonate
classification schemes, but rather combine the effects of depositional texture and
diagenetic overprinting into classes of common petrophysical properties. This approach
was necessary because the diagenetic overprinting is not controlled by fabric or lithofacies.
Average properties for porosity and permeability and main pore types are listed for each
RRT Well name and sample depth are given below each photomicrograph.
TIME-DEPTH CONVERSION
The following figure shows time-to-depth conversion in carbonates is predominantly a
function of the 3-D porosity distribution within the carbonate platform. Reservoir modeling
was first performed in the time domain (A). This approach allowed construction of a 3-D
Nido limestone velocity model that honors the actual 3-D reservoir architecture. Based on
this velocity model, the porosity dependent 3-D depth conversion of horizons was
performed, matching in all scenarios the well data (B). For this approach, porosity-
velocity functions for gas- and oil-bearing zones and the gas-oil contact (GOC) were
provided to the model (steps 1 and 2 in the box). Using this information, depth (z) was
calculated on a cell-by-cell basis (step 3 in the box).
SUMMARY OF CHAPTER 8
The seismic sequence consists of acquisition, processing and interpretation. The three
major processing techniques are deconvolution to improve seismic resolution, stacking to
remove noise and migration to improve the seismic image.
PSDM (pre-stack depth migration) processing enhances the seismic record significantly.
The latest development in seismic processing is Full Waveform Inversion (FWI), which
has become possible as a result of ever growing computer power. Clearly, with FWI,
more data are considered and a better image of the subsurface may be expected.
Seismic visualization has made major steps forward thanks to the rapid improvement in
seismic work stations. Workstations allow fully integrated reservoir characterization,
because all data is being dealt with in a consistent way.
4D seismic has become a powerful technology to monitor fluid flow in reservoirs. The
Geophysical Reservoir Modeling (GRM) process consists of three steps, namely 1)
interpretation of seismic differences in terms of flooding patterns; 2) history matching of
the flooding patterns derived from seismic together with history matching of produced
volumes, pressures, fault sealing patterns, etc; and 3) construction of synthetic seismic
using the history matched reservoir model and comparing this synthetic seismic with the
recorded seismic.
Well-to-well seismic is powerful in cases where well distances are relatively small.
Seismic while drilling may be a means to view ahead of the bit, e.g. for overpressure
prediction.
Assembly of the pieces of the reservoir characterization puzzle often occurs within the
reservoir engineering process. The traditional approach was to assemble the pieces
from the other disciplines and put together a model that explains as much as possible.
Typically, little feedback was given to the other disciplines for some very good reasons:
The software used to process the models was based on different platforms (pc
environment for engineers, Unix for geoscientists)
Prediction of the future performance under the prevailing reservoir conditions (i.e.
estimation of the reserves and recovery factors)
Estimation of the recoverable reserves and recovery factors under various other
producing methods that are presently known or projected.
Periodic updating of the model as the quality and quantity of the data improve
during the life of the reservoir
However, such responsibilities can only be carried out successfully in close cooperation
with other disciplines. Ultimately, there is one single reservoir model with its
uncertainties (leading to multiple realisations) for which input of all disciplines is
required. Every time that new data is obtained, the reservoir model needs to be updated
using the input from the different disciplines.
Weber showed that input from all disciplines is required to enable integrated reservoir
studies as displayed in the next diagram:
Gravity (buoyancy)
- segregates fluids
Capillary forces
- Increase wetting phase saturation
(effect more in tighter rocks)
Viscous forces
- impedance to flow
DRIVE MECHANISMS
Critical data appropriate for these responsibilities include the analysis of the drive
mechanism active in the reservoir. Variations in reservoir pressure, production rates,
gas-oil ratios, water influx and gas cap expansion are determinants in describing the
drive mechanism present. Important contributions in analyzing the drive mechanism are
enhanced by integrating lithologic, geophysical and sequence stratigraphic data. Drive
mechanisms recognized include:
Water drive
Gas cap
Dissolved gas
Compaction drive
Gravity drainage
Gas depletion
Combinations
The capacity to produce oil from a reservoir is dependent upon the reservoir fluid
pressure available and the use of this pressure in moving fluids through the reservoir
pore system into the well bore. Production of reservoir fluids is dependent upon a
pressure draw down; therefore, the pressure drop that is produced between the static
reservoir pressure and the well bore flowing pressure. The relation between reservoir
draw-down pressure, flow rate, and permeability will be discussed in the next chapter
when individual oil well performance is considered. The use of pressure to move fluids
into the well bore, however, does not explain reservoir performance, and, hence, the
natural energy system used to produce oil during the primary production period.
The natural energy or reservoir drive that is used during primary production operations
can be visualized by considering that each unit volume of oil produced must be replaced
by something in the reservoir since a vacuum cannot exist. In describing the various
drive mechanisms that do exist, therefore, each mechanism will be discussed
concerning the material that is replacing the oil being produced. As a consequence,
three major primary drive mechanisms will be discussed, which are: (1) the solution
gas-drive reservoir; (2) the gas-cap drive reservoir; and (3) the water-drive reservoir.
Additionally, a combination-drive reservoir will be discussed that is, in effect, a
combination of all three primary drive mechanisms. Finally, an additional energy source
called gravity segregation (potential energy) will be discussed although it is not
considered a primary drive mechanism, but a source of additive energy that can be
available under certain circumstances. Each of these drive mechanisms will be
discussed individually, and it should be remembered that for all practical purposes, the
reservoir is a constant temperature tank.
WATER DRIVE
A water-drive type reservoir is one in which the oil column is associated with a very large
aquifer. The oil column can be either an undersaturated oil or a saturated oil having a gas
cap. For the system to be specifically a water-drive reservoir, however, the gas cap does
not play a part in the energy drive mechanism. For a saturated oil reservoir to be
performing as a true water drive and not a combination drive, therefore, the pressure
would have to be maintained such that the gas cap does not expand. For this to be the
case, it would be necessary for an equal volume of water to encroach and replace the
voidage created by oil production. As can be observed, a constant pressure head exists
in the outcropping aquifer. With maximum access area available for water encroachment
into the oil column (resulting from the bottom-water drive) as oil is produced, water could
readily replace the voidage created. This would result in an optimal reservoir performance
during primary recovery operations and should lead to good primary recovery depending
upon the degree of reservoir homogeneity and oil mobility. It is easy, of course, to
observe the energy in an artesian-type water-drive reservoir since a pressure head exists.
However, for the non-artesian type of associated aquifer the question arises as to why
water would influx into the oil column as oil is produced.
The principal reasons for water encroachment during oil production would be due to: (1)
the cumulative result of aquifer water expansion resulting from the initial pressure draw
down and (2) aquifer pore volume compression in response to any pressure reduction.
Of course, both water expansion and pore volume compression in response to each psi
of pressure draw down would be small unless magnified many times. The volumetric
change in water and the volumetric change for rocks would require an associated
aquifer that is very large, hence, humongous in relation to the oil reservoir itself. A third
factor that could play a part in forcing water to encroach into the oil system would be the
evolution of gas from the aquifer water once the pressure in the aquifer is lowered
below the saturation pressure of the water. This condition may not occur, however,
unless the other two mechanisms are not sufficient to maintain pressure in the system.
Ultimate recovery by this type of primary production drive can most commonly be
expected to range between 40 and 60%, although higher recoveries have been
observed. A good rule of thumb for primary recovery efficiency by a water-drive type
reservoir would be 50%. The many factors that could influence the recovery efficiency
will be discussed in later chapters. Since recovery is usually quite high as a result of this
drive mechanism it would be a poor risk when considering enhanced recovery
operations. Since the oil has already been displaced by water, a more exotic and
expensive operation, as compared to waterflooding will have to be initiated.
Pressure production performance of some water-drive pools. (After Elliott, courtesy of AIME.)
Characteristics Trend
1. Reservoir pressure Falls slowly and continuously
2. Surface gas-oil ratio Rises continuously in up-structure wells
3. Water production Absent or negligible
4. Well behavior Long flowing life depending upon size of gas cap
5. Expected oil recovery 20 to 40 percent
this free gas volume that exists initially in the reservoir at initial reservoir pressure and
temperature conditions that substantially alters the performance behavior of this system
during the primary producing life of this reservoir. Wells are drilled and completed in the
oil column, and, as oil is produced and conceptual voidage created, the contained
reservoir fluids must expand to fill that voidage. This will be, as one might expect,
associated with a reservoir pressure decline. The dominant fluid expansion, however, is
associated with the free gas present in the gas cap. In fact, since the reservoir oil is a
saturated oil, the remaining oil volume will shrink with production, and gas will evolve
within the reservoir from this remaining oil volume. The voidage conceptually created by
production is, therefore, replaced by gas-cap expansion and solution-gas evolution
within the oil column. If the gas cap is of any consequential size, the dominant influence
in voidage replacement would be due to the gas-cap expansion in effect overshadowing
the influence of solution gas in the replacement of the voidage created by oil production.
It should also be pointed out that if any associated water exists, it is a very limited
amount and not a factor in this producing mechanism nor would be the influence of pore
volume compression.
Typical gas-cap-drive performance, Goldsmith San Andres Dolomite Pool in West Texas. (After Craze, courtesy of
AIME.)
Characteristics Trend
1. Reservoir pressure Declines rapidly and continuously
2. Surface gas-oil ratio First low, then rises to maximum and then drops
3. Water production None
4. Well behavior Requires pumping at early stage
5. Expected oil recovery 5 to 30 percent of original oil in place
Once the bubble-point pressure is reached in the reservoir, however, the performance
characteristics change drastically. The voidage created by production must still be filled
by expansion of the remaining reservoir fluid, but now the associated pressure drop is
below the bubble point and the oil will shrink. Of course, associated with this shrinkage
will be evolution of solution gas within the reservoir, which readily expands with
pressure reduction and readily fills the voidage created by both production and
shrinkage of the oil. The evolution of solution gas in the reservoir strongly influences the
performance characteristics. It can be noted that rapid initial decline in pressure to the
bubble point is somewhat arrested due to the expansion capacity of the gas, which
tends to support the reservoir pressure and cushion its decline. One would also expect
that the producing gas-oil ratio would reflect this gas evolution in the reservoir with a
reflected decline in the producing gas-oil ratio since not all the initial gas in solution is
retained in the oil entering the well bore. This lower gas-oil ratio will exist as long as the
gas being evolved in the reservoir is below the critical saturation that is the saturation of
gas required for gas to flow. Once this reservoir gas saturation builds up to the point
where gas starts to flow as a free gas phase into the well bore, a rapid rise in the
producing gas-oil ratio will be observed. This rising gas-oil ratio will continue over time
as the gas saturation in the reservoir continues to increase. It will not continue to rise
throughout the primary production period, however, as a function of three combined
influences. These influences are: (1) the reservoir pressure is declining, and, hence, the
capacity of the gas to expand will decline; (2) the oil production rate will become more
stabilized later in the primary producing life of the reservoir; and (3) the oil that is being
produced later in time contains less gas in solution. This gas-oil ratio behavior becomes
a fingerprint for a solution-gas drive reservoir. Finally, the oil production behavior is also
quite sizably altered once the bubble- point pressure in the reservoir is reached. The oil
production rate will undergo a rapid decline caused by the loss of mobility resulting from
the combined loss of effective oil permeability and the increase in oil viscosity. This
rapid loss of the capacity to produce oil should be expected as a result of the loss in
mobility and pressure decline. In time, the rate of loss will diminish as the production
history goes through a transitional period and ultimately reaches a somewhat stabilized
but low oil producing rate. If this stabilized oil production rate is still economic, it may be
projected to continue for many years.
This drive mechanism, however, is not efficient with expected primary performance to
result in recoveries ranging from as low as 5% to 20%. A good rule of thumb to apply as
a generalized recovery factor would be about 15% of the initial oil in place recovered
during primary production operations. Higher recovery efficiencies might be
encountered if certain conditions exist, as will be discussed in later chapters. One mode
of improved primary recovery performance would be that of pressure maintenance,
hence the injection of water and/or gas into the reservoir to maintain pressure. Pressure
maintenance would keep gas in solution, thereby optimizing oil mobility as well as
maintaining pressure and ultimately improving primary recovery. The concept of
pressure maintenance will be elaborated upon in a later chapter. It should also be
obvious that with such low expected primary recovery efficiency that this type of
reservoir is a prime candidate for post-primary injection operations and in particular, a
waterflood, since it is the most economic of the post-primary injection operations. These
operations will be discussed in a later chapter.
Performance of the Slaughter San Andres Limestone Pool, West Texas, under predominantly solution gas drive
(After Sessions, courtesy of AIMT).
COMBINATION DRIVE
COMPACTION DRIVE
When the reservoir pressure is declining during production, the effective stress (being
the difference between the overburden stress and the pore pressure) is increasing. This
effective stress is acting on the rock frame as such decreasing the pore space. This
phenomenon is called compaction. Effectively, the hydrocarbons are squeezed out of
the pore space, which is called compaction drive.
GRAVITY DRAINAGE
GAS RESERVOIRS
Gas reservoir fluid systems were described in Chapter 4 and consist of the gas
reservoir (whether defined as a wet or a dry-gas system) and the retrograde gas
condensate reservoir. These reservoirs are the non-associated type, do not exist in
equilibrium with an oil column and, therefore, are not gas-cap systems. The
non-associated gas reservoirs, however, may coexist with an aquifer that can provide
some degree of water drive.
When producing a gas reservoir, the reservoir engineer's goal is to obtain maximum
recovery just as it is when producing an oil reservoir. This optimal recovery includes
maximum recovery of both gas and condensable liquid hydrocarbons. For the reservoir
engineer, this presents several problems: (1) the loss of hydrocarbon liquid in the
reservoir when producing a retrograde gas condensate reservoir by pressure depletion,
and (2) low recovery efficiencies (particularly sweep and conformance efficiencies)
when producing, in particular, a retrograde gas condensate reservoir by water drive.
There are several methods available for obviating this liquid loss in a
pressure-depleting, volumetric-type reservoir, which are: (1) gas cycling, and (2) water
injection pressure maintenance. The most desirable of these would be gas cycling since
it will usually result in higher recovery efficiency when compared to water injection
pressure maintenance. This process shows how the retrograde gas condensate fluid is
produced from the reservoir and replaced with a dry gas. The original reservoir
retrograde phase envelope is one in which, when pressure depleting Reservoir A at
constant reservoir temperature, the dew point would be crossed and liquid condensate
can be expected to occur within the reservoir. By producing this system while it is at
high reservoir pressure conditions, stripping the condensable liquids in the surface
processing system, and then re-injecting the dry gas the problem would be eliminated.
The re-injected dry gas miscibly displaces the retrograde gas. This results in a reservoir
gas system that can be pressure depleted, with no liquid condensation occurring in the
reservoir, once cycling is completed. This is because the phase envelope for the dry
gas system has a cricondentherm less than the reservoir temperature.
The injection of water to maintain reservoir pressure above the dew point pressure for
this volumetric system having no water drive would lead to lower recoveries due to the
inevitable bypassing and entrapment of pockets of gas as the system is depleted. Initial
Some equations:
PV nzRT
z factor introduced to compensate for deviations from ideal gas law. z is a function of
pressure and temperature.
G p GIIP GIP
gas produced equals gas initially in place minus gas remaining in place.
P Pi Gp
1
z zi GIIP
FIELD EXAMPLE
If the primary reservoir energy is insufficient, then energy may be added to the reservoir
by injection of water or gas into the reservoir.
When injecting gas usually hydrocarbon gas is being utilised. However, particularly in
retrograde gas condensate reservoirs, also nitrogen injection is useful, When injecting
dry gas or nitrogen into the reservoir, the composition of the hydrocarbon is shifting to
lighter components. As a result, when at the end of the producing life the reservoir is
blown down, liquid saturation in the reservoir will be negligible or at least limited to a
minimum.
WATERFLOOD PATTERNS
PROBLEM ADDRESSED
Prior to the onset of a waterflood, there must be an accurate evaluation and prediction of reservoir
condition. This includes an analysis of the architecture, structure, presence and orientation of
fractures, as well as reservoir limits, rock, and fluid property and pressure distributions.
TECHNOLOGY OVERVIEW
Since waterfloods are typically implemented in mature developed properties, it might be assumed
that the reservoir is well defined, with few potential surprises. Experience indicates otherwise.
However, operators who are aware of common pitfalls can avoid mistakes that will reduce recovery
and profitability when implementing waterfloods.
The growth of a waterflood prospect, from conception to abandonment, is divided into two stages.
The first stage, which is the topic of this workshop, deals with the initial analysis up to the
development of the reservoir model. It is especially critical that engineers and geologists work
together during this initial conceptual stage. The second stage, which is more engineering- oriented,
deals with development of the reservoir model, implementation, and operation until project
abandonment.
Typically, an operator’s first step in evaluating a waterflood candidate should be to assemble and
evaluate preliminary reservoir data. Information about the gas:oil ratio (GOR) and water production
provide key information about whether a project is an attractive candidate. Operators should examine
GOR data, both on the basis of the well’s initial potential and cumulative production basis. Oil
reservoirs exhibiting cumulative GORs that are less than 3000/1 are generally attractive candidates.
Wells with a GOR less than 3000/1 are generally attractive candidates. Gas reservoirs do not make
favorable waterflood candidates because of the lack of sufficient oil to bank. The operator should,
however, be aware that a well’s high GOR can be caused by comingling of fluids from different
stratigraphic layers. Detailed stratigraphic correlation should be able to distinguish the types of fluids
coming from various zones.
Operators should take special care looking for indications of a gas-oil contact or gas cap. If present, a
gas cap will distort GOR data and its presence must be considered when designing a waterflood. It
also is critical to note water-oil contacts, calculate fluid saturations, and identify sources of water
production. In analyzing gas and water data, operators should be aware that they are not always
accurately reported, so considerable judgment may be required.
A two-step analysis of reservoir data is used to evaluate a waterflood candidate. First, the operator
reviews initial production potential, gas-to-oil ratios (GORs), amount and type of water production,
pressure levels, and core information. If this information meets the requirements, a second, more
thorough review of data of the reservoir’s fluids pressure, volume, and temperature (PVT) properties
should be initiated in cooperation with the company’s engineering department. Ii also is important to
understand the source of formation water within the reservoir and calculate its saturation. Formation
water production can be caused by many processes, such as water saturation above irreducible water
saturation, water entering the wellbore accidentally from another reservoir (e.g. hole in casing), water
introduced mechanically (e.g. an injection well), or natural encroachment due to migration of the oil-
water contact.
Production history should be analyzed, too. An analysis of primary recovery will determine not only
if sufficient oil saturation remains for successful water-flooding but also if there are potential
problems.
Determining net pay for a waterflood is important because the calculations affect original -oil-in-
place (OOIP) values, reserve estimates, economics, and performance characteristics. Properties that
can have an impact on net pay include porosity-permeability correlations, average porosity and
thickness, impermeable parts of the reservoir, and fluid contacts. Mapping net pay may be the
primary contribution of geologists to evaluating waterflood candidates.
Although it’s often the least-evaluated aspect of a waterflood project, the reservoir’s structure can
have major influences, such as those caused by faulting (and fracture); degree of dip; and
determination of overall homoclinal, anticlinal, or synclinal nature of the reservoir.
The proper techniques for core analysis begin even before drilling and continue through laboratory
examination. Cores can be used to significantly increase data available about lithology, depositional
environment, diagenesis, and rock type. Core samples are required to define pore geometry changes
with depth, areal distribution, permeability, and porosity.
But fractures can be a positive benefit for wells under primary production because of increased
permeability accessibility. This holds true whether the fractures are
natural or induced, and whether they are produced from stress caused by local, regional or basinal
factors.
Fractures also can be a detriment to hydrocarbon production because they may connect to downdip
water. Because fracture permeability is often several orders of magnitude higher than the matrix
permeability, injected water will generally travel through fractures. This can be disastrous for a
waterflood if water is injected into a fracture that is in communication with a producing well.
Economic Implications:
Economic analysis, with all its incorporated data, predicts the success of the waterflood. As a result,
neither the geologist nor the engineer alone is responsible for the prediction. A realistic price
structure for the economic evaluation is often the deciding factor for economic success of
waterfloods.
Long-range price trends are important, since the future value of a waterflood project is based on an
extrapolated price of oil. Because the future price of oil is unknown, project economic analysis
should incorporate a conservative price structure. Further, prices must be adjusted, taking inflation
into account, as part of the economic analysis.
Based on a workshop sponsored by P7TC’s South Midcontinent Region and the Oklahoma Geological Survey on
June 10-11, 1998, in Oklahoma City, OK
reservoir. Such Enhanced Oil Recovery schemes are aiming at increasing the mobility,
i.e. the ratio of permeability over viscosity, of the oil in the reservoir.
Steam injection
By injecting steam the temperature in the reservoir is being increased.
That means that the viscosity of the oil is being decreased and
consequently that the oil mobility is being increased.
Polymer flooding
It often is the case in a water drive that the water mobility is greater than
the oil mobility (unfavourable mobility ratio). In this situation water will
finger through the oil and eventually bypass the oil leaving remaining oil
behind in the reservoir.
In this situation will the water mobilty clearly be greater than the polymer
mobility. Hence, the water will finger through the polymer and eventually
bypass the polymer. The polymer slug has, however, been designed such
that breakthrough of the water through the polymer will only occur when
the head of the polymer slug has arrived at the production well.
Surfactant flooding
Introduction of surfactants in the reservoir will reduce the interfacial
tension between oil and water. In this way the relative permeability of the
oil with respect to water will increase resulting in a higher oil mobility.
CO2 flooding
CO2 is injected into a reservoir at residual oil saturation. Two effects occur,
namely the saturation of the oil containing the CO2 increases above
residual saturation and becomes mobile and secondly, the viscosity of the
mixture is reduced. Both effect contribute to improved oil recovery.
In-situ combustion
When oxygen is being injected into the reservoir, combustion of the oil
around the injector well may occur. Such in-situ combustion increases the
temperature (hence reducing oil viscosity) and the pressure (hence adding
Microbial EOR
Microbes are injected in the formation together with nutrients such as
sugars, phosphates or nitrates. The microbes ferment the hydrocarbons
and also produce surfactants and CO2 which lead to processes as
described above.
UNCONVENTIONAL RESOURCES
Gas
- shale gas
- tight sandstones
- coal bed methane
- coal gassification
- deep HP gas
- methane hydrates
Oil
- tar sands
- oil shales
- bio-fuels
- GtL, Gas to Liquids
Volumes of hydrocarbons in place commonly are derived through techniques that are
based on:
Volumetric calculations
Decline curves
Material balance
Mathematical simulation
Each of these is well known, and details of their implementation are outside the scope
of the course. It is worth pointing out, however, that the geostatistical techniques
described in a preceding section can be used to provide multiple realizations
(simulations) based on the expected ranges of components; these realizations can
provide the basis for a Monte Carlo type estimation of reservoir volumes (see section on
economics and risking).
One of the most important functions of the reservoir engineer is to determine the reservoir
oil and gas in place and predict the recovery anticipated under the prevailing reservoir
mechanism. This determination of the initial hydrocarbon-in-place (oil and/or gas) and the
projected recoverable reserves are essential for economic planning. These evaluations, of
course, are a continuing process throughout the life of the reservoir. Hydrocarbon volumes
in place and reserve estimates are made initially to optimally develop the reservoir system.
Continued monitoring of initial hydrocarbons-in-place and reserve predictions are necessary
to confirm earlier estimates and for long-term planning pertaining to maximizing oil and/or
gas recoveries.
To accomplish these reservoir evaluations, the engineer requires an array of data falling in
four broad categories. These categories include: (1) reservoir rock properties, (2) reservoir
fluid properties, (3) a reservoir pressure history, and (4) the reservoir production history.
Previous chapters have discussed the acquisition and to some degree the application of this
data. This chapter will further emphasize the need for accurate and representative data as
required in reservoir evaluation.
The material balance is a tool the reservoir engineer can use to predict initial hydrocarbons-
in-place using production history or reserves using initial hydrocarbon-in-place. It can be
expressed as a balance between the voidages created by production being equal to the
change in volume of all the reservoir fluids, plus the volume of all the extraneous fluids
which enter the reservoir. As such, it can be visualized as the reservoir engineer's
accounting system.
The material balance is a tool the reservoir engineer has used to simulate reservoir
performance. The problem with the material balance approach is, however, two fold. First
allowances cannot be made for the variation of rock and fluid properties with location in the
reservoir and secondly, the dynamics of fluid movement within the system are overlooked.9
Material balance applications are, of course, still used since they are technically correct and
economically justifiable.
In order to overcome these limitations analog models were used. This approach, however,
also has many drawbacks some of which were that they were expensive, difficult to make
and operate and lacked adaptability.
The recent application of 4-D seismic for reservoir monitoring has greatly enhanced
reservoir management and will enable even greater recovery to be achieved. 4-D
seismic monitoring is the use of repeated 3-D seismic surveys obtained over time.
Differences between 3-D seismic surveys can be evaluated to reflect changes in fluid
distribution over time. To do this competently, however, all factors affecting the seismic
data must be considered. When monitoring fluid movement in a reservoir, therefore, it is
necessary to integrate the engineering process and geologic model with the seismic
survey. 4-D seismic provides an insight into the interwell area that cannot be obtained
from any other source. This insight enables the engineer to "see" the reservoir as it
responds to production operations. Reservoir drainage patterns can be observed,
bypassed oil and gas defined, fluid influx more optimally controlled, and additional wells
(or laterals) drilled to optimize recovery.
RESERVOIR MODELING
Additionally, characterization of the reservoir as individual flow unit or - most likely - a
series of interrelated flow units (zones within the reservoir with similar fluid flow) is a
key step in understanding the reservoir and developing simulations. Flow units are best
described with a thorough understanding of the geologic model developed through
collaboration by the reservoir characterization team.
RESERVES
Working separately, the Society of Petroleum Engineers (SPE) and the World
Petroleum Congresses (WPC) produced definitions for known accumulations –
adopted in 1996
CONCEPTS
“Proved reserves should be based on current economic conditions, including all
factors affecting the viability of the projects …. the term is general and not
restricted to costs and price only”
ESTIMATION OF RESERVES
Method of estimation is called:
Deterministic if a single best estimate of reserves is made based on known
geological, engineering, and economic data
Probabilistic when the known geological, engineering, and economic data are
used to generate a range of estimates and their associated probabilities
VOLUMETRIC DETERMINATION
NR = GBV N/GF(1-Sw)
Where:
NR = hydrocarbons in place at reservoir conditions
GBV = Gross Bulk Volume of reservoir
N/G = Net to Gross ratio
F = Porosity, fraction
Sw = Water saturation, fraction
MORE ON RESERVES
Estimates will generally be revised as additional data becomes available or
economic conditions change
PROVED RESERVES
“Those quantities of petroleum which, by analysis of geological and engineering
data, can be estimated with reasonable certainty to be commercially recoverable,
from a given date forward, from known reservoirs and under current economic
conditions, operating methods, and government regulations”
In certain cases, proved reserves may be assigned on the basis of well logs
and/or core analysis analogous to other proved reservoirs
UNPROVED RESERVES
Similar to proved reserves, but
– Technical, contractual, economic, or regulatory uncertainties preclude such
reserves being classified as proved
PROBABLE RESERVES
“Those unproved reserves which analysis of geological and engineering data suggests
are more likely than not to be recoverable”
Types
Formations appear to be productive based on well log characteristics - lack
definitive data/tests
Area of the formation separated from the proved area and structurally higher than
the proved area
POSSIBLE RESERVES
“Those unproved reserves which analysis of geological and engineering data suggests
are less likely to be recoverable than probable reserves”
When probabilistic methods are used, there should be at least 10% probability
that the quantities actually recovered will equal or exceed the sum of estimated
proved plus probable plus possible reserves
RESERVES SUMMARY
Proven reserves, 0 to P90
(developed or undeveloped
FLOW UNITS
Major types of geological and petrophysical data applied to flow unit zonation of a well.
As seen in this Figure, the flow units are delineated on the basis of permeability
contrasts due to lithofacies changes and to the presence of a laterally continuous barrier
to vertical permeability (flow unit 4). According to the definition of flow unit, it is
permissible to define a flow unit that exhibits only weak flow or no flow through it. This
property of flow units makes it possible to use a single numbering system for identifying
both obvious flow units and probable permeability barriers that can be mapped at the
same scale as reservoir quality flow units.
Some examples of lithofacies and flow unit subdivisions of clastic and carbonate reservoirs, (a) Lithofacies and (b)
flow unit subdivision of the Shannon Sandstone body in the Hartzog Draw field, Powder River basin, Wyoming.
Modified from Hearn et al., 1984.) (c) Lithofacies and (d) reservoir facies (flow unit) subdivision of the Rainbow Lake
reef reservoir (“A” Pool), Alberta, Canada. (Modified from Langston and Chin.)
HETEROGENEITY
Baffles to flow
Permeability heterogeneity
Contrasts in vertical and lateral permeability
IMPEDIMENTS TO FLOW
Barrier
Baffle
Hi K contrast
Low K contrast
HETEROGENEITY CAUSES
Causes:
– Depositional and stratigraphic controls
– Diagenesis
– Pore system connectivity – diagenetic or depositional
– Intersecting structural and stratigraphic patterns
– Hydrodynamics and overpressure
Detailed study of faults and clay smearing in Denver Basin (K) – Shale Gouge
Ratio (SGR) used to calculate connectivity
– Supporting data – BHP bombs, pressure transient tests, GOR, production
history
Diagenetic and oil saturation differences in Red River in Williston basin identified
two trends
– (NE-SW pinchout of porosity and NW plunging nose)
– best wells at intersection of trends
Heterogeneities can occur at all scales within the reservoir. To appropriately simulate
the reservoir, an adequate understanding of these heterogeneities and their affects are
needed. Lassater et al. discusses reservoir heterogeneity in detail; in their paper, they
classify the heterogeneities into small (centimeters in scale), medium scale
(depositional units) and large scale (many depositional units). They summarize their
investigation as follows:
At the medium scale, gravitational forces are more important and capillary
forces less important. Relative to the small scale, fewer pore volumes
sweep a medium-scale region. Heterogeneities tend therefore to have a
greater effect. In comparing the channel sand and bar sand examples, we
showed that while the channel appears to have a less favorable saturation
profile, behind the front there is probably little difference in sweep
efficiency. Additional cases need to be examined in order to reach a more
general conclusion.
CATEGORIES OF COMPARTMENTS
By this point in the course, the student should understand several methods to identify
each of these categories. The producing and depleted reservoirs seem obvious; the
question here remains the size and distribution of the reservoir (or flow unit).
Stratigraphic modeling techniques allow a prediction of the reservoir dimension. These
dimensions can be verified with reservoir simulation techniques. This is an important
practice to develop confidence in the predictive capability of the reservoir
characterization techniques being used.
SCALES OF HETEROGENEITY
FLUVIAL FACIES
DELTA-PLAIN FACIES
The next figure shows a reservoir model based on the proximal delta-front fades
association. The association is characterized by alternation between laterally extensive
delta-front siltstones, and mouth-bar and channel sandstones. The result is a stacked
succession of reservoirs vertically sealed by the delta-front siltstones. Internally, the
ands display low heterogeneity. Coarsening-upward and fining-upward profiles occur
and are likely to be reflected in systematic permeability trends. In addition, individual
reservoirs are likely to be compartmentalized by sealing faults.
The previous figure shows a reservoir model based on the distal delta-front facies
association. The association is characterized by tubular sandstone beds interpreted to
be the product of hyperpycnal underflows separated by fine-grained siltstones and
mudstones. This type of reservoir Is likely to have zero effective vertical permeability
and to have horizontal permeability strongly controlled by vertical faults.
RESERVOIR MANAGEMENT
Maximize recovery
Treatment (fracturing, acidizing = inflow improvement)
Secondary recovery techniques, rates, patterns
EOR
Tools
Reservoir Characterisation (static model, flow units)
Simulation
Reservoir monitoring
When optimizing the business, there is a choice of three options as shown in the
following diagram:
VIR = Net present value of the cash flow divided by the net present value of the investments
FIELD SURVEILLANCE
Pressure
Tools: MDT, permanent downhole pressure gauges
Rates
Tool: PLT
O/W contact
Tools: c/o logging, RST
HORIZONTAL WELLS
PI hor Lhor kv
PIF
PI vert Lvert kh
where
PI productivity index
Lhor net length of the reservoir drained in horizontal well
Lvert net perforated interval in the vertical well
kv vertical permeability
kh horizontal permeability
After: Hurst and Bong Poh Yuk, SPE IPTC 15199 (Nov 2011)
UPSCALING
Reservoir simulation models consist of the assembly of individual blocks or cells that
represent discrete portions of the reservoir. Reservoir properties and fluid movement
dynamics then are calculated for each block over time. These complex calculations are
time consuming and computer intensive, therefore, in the interest of minimizing time and
cost, the total number of blocks often is reduced. Upscaling is the term that applies to
averaging the detail from a data source - usually for reasons of efficiency. The need to
diminish the model size is balanced against the need to adequately represent the detail
available and that required to accurately characterize the reservoir. Achieving balance
requires tough decisions and significant effort within the reservoir characterization team.
Ideally, all of the detail available to the reservoir characterization should be used in an
actual simulation model, but modern reservoir simulation software running on computer
workstations has a practical limit of about 100,000 blocks (models may consist of a
million or more blocks.). This limit is set by current software products and the enormous
memory requirements needed by most simulation software. Therefore, without extreme
care being exercised, much of the detail generated by the reservoir characterization
team could be eliminated from the simulation model. The conclusions to be drawn are:
Planning prior to data gathering - with the needs of the simulator in mind - is time
well spent
This upscaling process has the most potential for creating inaccuracies in large
heterogeneous reservoirs. The challenge is to develop a balance between providing
details required to depict the reservoir versus minimizing computational requirements. It
also is important to be able to develop multiple realizations of the simulation in order to
look at a wide range of possibilities and to reach a good history match. Lassater et al.
demonstrate upscaling graphically:
UPSCALING – CELLS
THE STRUCTURES
reflect the physical processes of sediment transport and deposition and include planar
lamination, cross lamination and cross bedding ( A), which can be described as 2D or
3D, depending upon the shape of the original sedimentary bedform ( B).
The transport mechanisms responsible for the different structures can lead to variation
in sediment texture both within and between laminae, and between adjacent structures,
with porosity and permeability controlled by grain size, packing and sorting. Permeability
can vary by an order of magnitude or more between adjacent foreset laminae. In recent
years, the impact of such small-scale permeability variation on hydrocarbon recovery
has been demonstrated, and upscaling procedures such as the Geopseudo method
have been developed to take advantage of the hierarchy of geological length scales.
The Geopseudo Atlas library of lamina, lamina set and bed upscaling models is
intended to be a comprehensive representation of the types of sedimentary structures
likely to be encountered in reservoir rocks. Dependent upon the 3D geometry of the
sedimentary structures, 3D ( C) and 2D ( D) models of sedimentary structures can be
constructed, by use of geologically relevant input.
UPSCALING CHOICES
RESERVOIR SIMULATIONS
One of the goals of reservoir characterization is to develop data needed for input into a
computer-based simulation or characterization of the reservoir. Reservoir simulators
are generally classified as black oil, compositional, thermal and chemical, depending
upon fluid flow, mass and heat transport behavior as discussed below (from Satter and
Thakur1994):
Thermal simulators account for both fluid now and heat transport and chemical
reactions. They are used for simulating steam-flood and in-situ combustion
processes.
Chemical simulators account for fluid flow and mass transport due to
dispersion, absorption, partitioning, and complex phase behavior. They are used
for surfactant, polymer, and alkaline flooding.
HISTORY MATCHING
Starting point: The static and the dynamic model are one integrated model.
Changes in the simulation model must be checked against what is possible in the
petrophysical geological model.
It is still an integrated effort
SIMULATORS - OBSERVATIONS
Running a simulator requires special expertise.
Simulators require significant computing power.
The up-scaling may have simplified the model.
SIMULATORS - METHODS
Differential equations (hydraulic diffusivity equation) relate pressure, volume and
time at every location in the reservoir. The hydraulic diffusivity equation conbines
three basic principles, namely the conservation of mass, the equation of state
and Darcy’s law.
No-flow boundaries may be recognized that subdivide the reservoir into sectors.
MECHANICS OF SIMULATION
Reservoir divided into grid cells using a X, Y, & Z coordinate system
Computations are carried out for all phases in each cell at discrete time steps
SIMULATORS – INPUT
History matching
Geologic modeling
Porosity, thickness, saturation, permeability
Upscaling
Technology integration
4 D seismic
Logs, tests, performance data
Reservoir data
Parallel simulation use new, parallel architecture computers that break the simulation
model into smaller pieces and distribute them among multiple processors. In this way,
each processor effectively solves a much smaller model. As a result, the operator can:
Parallel simulation allows the operator to perform full physics simulation on a geologic
scale and receive the results in a reasonable timeframe. An initial model for the geologic
scale (200,000 grid blocks), broken into 16 segments and simulated with a multinode
IBM 590 processor is performed in less than 12 hours.
CONDITIONING DATA
Fit data to model
Integrate all data sets
Fit to local conditions
Normalize throughout field
Carefully correlate data sets
Determine sensitivities to known variables
Petrophysics
Permeability*
Porosity*
Connate water saturation*
Relative perm.s*
Capillary pressure curves
Variability in porosity and perm.
Reports, documents
Fluid properties
PVT analysis
Black oil table*
Fluid composition*
Equation of state
Reports, documents
STREAMLINE SIMULATION
Streamline simulation is most effective when there are clear velocity direction, i.e. in
convective processes. If the velocity field is not so clear, as in diffusive processes such
as initial depletion generating a secondary gas cap or primary gas cap expansion.
The last two graphs in the picture concentrate on the situation between the injector and
well 4. The streamlines are coloured blue signifying the shortest Time-of-Flight between
the injector and the producer and pink the longest Time-of-Flight. The watercut graph
is also coloured with blue showing the early watercut and pink the late watercut. With
these two graphs the watercut performance in years is related to a particular set of
streamlines. Let say that the watercut performance in the yellow period does not fit the
history, then changes in those parts of the reservoir model need to be made that are
penetrated by the yellow streamlines.
Disadvantages are:
The assumption that streamlines are fully independent is not fully correct.
It can not easily capture physics that is transverse to the flow direction (e.g.
gravity, capillarity, diffusion, compressibility, transverse thermal effect), although
this disadvantage has to some extent be solved by ‘operator splitting’.
Thiele, M.R., Batycky, R.P. and Fenwick, D.H., Streamline Simulation for modern
reservoir engineering workflows, JPT, January 2010
RESERVOIR MANAGEMENT
The amount of oil and gas likely to be recovered from a new field depends partly on
where the wells are drilled, the type of wells used (vertical or horizontal), the type of
hydrocarbons being produced, and whether gas or water is injected to maintain
pressure. ExxonMobil's reservoir engineers are known for their ability to get the most
value from oil and gas resources. That is an important advantage in the high-cost
business of deepwater development.
Ubit Field, Nigeria - increased the amount of recoverable oil by 500 million
barrels through the application of 3-D seismic, advanced reservoir modeling
techniques, horizontal wells, and by upgrading surface processing equipment
Cold Lake, Canada - increased the recovery of heavy (tar-like) oil from 10
percent to more than 25 percent by injecting steam into the formation to heat the
oil and make it flow toward the producing wells.
Arun Field, Indonesia - recovering about 94 percent of the gas and 67 percent
of the condensate from a mature, lean gas condensate reservoir using
compositional reservoir modeling to balance gas injection and maximize
condensate recovery.
Brent Field, North Sea - increased the amount of recoverable oil and gas by
400 million oil-equivalent barrels through the early reduction of reservoir pressure
to liberate gas from residual and bypassed oil. The use of advanced reservoir
modeling, reworking existing wells, drilling new wells and upgrading surface
equipment have contributed to the increase in reserves.
To meet the world's growing demand, energy companies must be able to find oil and
gas from increasingly difficult locations and bring them to market at a competitive price.
Thanks to 3-D seismic imaging and advanced reservoir modeling, we're doing that now
and will continue to do so well into the future.
Complex interplay of facies in Aramco Shuaiba (Cretaceous) - Rudist barrier (orange), lagoonal (blue), and basinal
(green).
The depositional model is a prograding carbonate ramp with local patch reef
accumulations in the shallower waters. The reefs have lower permeability than the
grainstones surrounding them and are "baffles" to vertical communication within the
reservoir. The unique placement of these reefs caused the lagged oil to form where it
did. The simulation model results, which are shown in the diagram as changes in oil
saturation in 5-year increments, reproduced the lagged oil correctly. This was
considered a positive feedback on the geology of the model, which continues to be
predictive many years later.
HANIFA RESERVOIR
The Hanifa reservoir in Berri field has been on production under peripheral water flood
since the mid-1970's. During that time, a significant quantity of the reserves has lagged
the water encroachment. In order to evaluate alternative strategies to recover this
lagged oil, using reservoir simulation modeling, a detailed geological study and
geocellular modeling effort was launched in the late 1980's. Previous simulation efforts
had been unable to model the lagged oil in the field. This time, a sequence stratigraphic
framework was developed which changed the layering and facies distribution within the
geological model. This was the first sequence stratigraphic framework model
constructed in Saudi Aramco and the subsequent successful reservoir simulation
history-match and predictions have confirmed its power as a predictive geological tool.
Structural growth has raised more of the low quality facies into the oil column, leaving
more of the high permeability grainstones in the water let. Facies were distributed in the
model using a stochastic algorithm constrained by well control and parasequence
isopachs.
Porosity was distributed in the model using a stochastic algorithm constrained by facies.
As a result, the porosity model distribution by facies matches that of the well data.
Laboratory experiments
RESERVOIR PROPERTIES
Porosity – 16% 3800 acres
Bubblepoint – 625#
WATERFLOOD STUDY
Feasibility for waterflood at Horse Creek
Not a good candidate due to poor water/oil relative permeability
Rejected as possibility
ARC TESTS
Assessed oil oxidation parameters under reservoir conditions
Qualitative tests indicating oil reactions with air
Indicated oil would react at low temperature conditions
EXERCISE – CONTRIBUTIONS TO RC
1.
Establish Stratigraphic Framework
2.
Establish Facies, Flow Units and
Geologic Prototype Areas in Cored Wells
3.
Develop Digital Log-Based Statistical Models for
Facies, Flow Units and Permeability in Cored Wells
4.
Test Statistical Models Against Core Data
5.
Map Facies, Flow Units and Geologic
Areas Field-Wide with Log Data
SUMMARY OF CHAPTER 9
Reservoir studies require a great deal of integration between disciplines. However,
there is also a lot of work to be done that is typically carried out by disciplines.
Therefore, study teams have to decouple, do their specific discipline-type of work, and
then recouple to share information.
Understanding of the drive mechanisms that play a role in the reservoir is necessary to
enable accurate estimates of ultimate recoveries and reserves.
Unconventional resources both in oil and in gas have become a major part of
hydrocarbon development.
In (finite difference) reservoir simulation, it is not always necessary to model the full
field. No-flow boundaries may subdivide the reservoir into sectors that can be modeled
separately. In sector models concepts (such as size of drainage areas per well, water
and gas coning, mobility ratio effects) can be studied providing guidelines for optimum
reservoir development.
Horizontal wells have become bread and butter in hydrocarbon resource development.
The objective of a reservoir study is to build a reservoir model including all geology
(geomodel), all petrophysics (petromodel) and all dynamic properties (dynamodel) in an
integrated way while the model is only complete when all the uncertainties are also
specified. In addition to all technical work, the economics (econmodel) of the reservoir
development need to be considered. The close integration of these four models into the
reservoir model, will lead to the best understanding of the reservoir and consequently,
to the best development of the reservoir.
RISKING/ECONOMICS CHALLENGE
A = Pool
0.30
0.80 0.24 A Pool, B Dry
B = Dry
0.20
A = Dry
0.14 A Dry, B Pool
B = Pool
0.70
Objectives for economic and risking portions of the reservoir characterization process
are to:
CHANCE OF ADEQUACY
(1- RISK)
Used for individual factors or combined for entire accumulation
Adequacy of prospect or play equals product of individual control factors:
Source, Reservoir, Trap, Seal, Recoverability
Adequacy minimum economic accumulation not presence or absence
RISKING GUIDELINES
1.0 Certain
0.7 Likely
0.5 Equally likely/unlikely
0.3 Unlikely
0 Impossible
The first and most crucial step in assessment is to have a firm understanding of the
geologic factors controlling oil and gas occurrence. These control factors become the
geologic risks. List below all the factors that you can think of that are essential to the
geologic creation of an oil or gas field having significant potentially recoverable
reserves. You may list more than 10 if you want to.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
RC Hint:
Always understand the economic impact
of your work or recommendation
– Act as a business person not just a technical expert
Economics and risking are critical to every oil and gas project. Each oil and gas
company has competing demands for its limited funds, and those that are most
successful are the ones that evaluate their business decision most accurately over time.
Underestimating the value - or overestimating the risk - of a project means that the
project may be funded at the expense of a more suitable investment. The converse is
true as well. An ideal assessment scheme is one in which the sum of the risked means
of all project assessments equals the sum of the reserves found.
The group decided that risk could best be defined as "the difference between
the value of an acceptable outcome and the value of the possible range of
outcomes at any given time. The definition itself gives clues as to the
elements involved in estimating risk: there must be an estimate of the
acceptable outcome followed by an estimate of the range of possible
outcomes. Unfortunately, neither task is especially easy. An estimate of the
acceptable outcome entails a reasonable understanding of the economic
climate under which the project is being undertaken. This factor is highly
process specific inasmuch as the requirements of acceptable rate of return
(short-term benefit) must be balanced against the needs of reserve additions
(long-term benefit). Unfortunately, the acceptable outcome estimate is itself
uncertain owing to variances in economics, politics and costs.
Estimating the range of uncertainty involves using the model developed above
to translate the effects of individual uncertainties into the cumulative
uncertainty of the model output. How to do this efficiently is the subject of
active research, as evidenced by the number of presentations dealing with
mathematical modeling and stochastic assignments at this conference. The
cost of many of these techniques is large; however, this cost, like those
associated with fathering and analyzing data, must be factored into the notion
of a acceptable outcome as discussed above.
As a final point, the group noted that the idea of risk manifests itself in numerous ways in
current practice. The focus of this discussion has been on economic risk, but there is now
a risk, incompletely manifest in economics, that is associated with environmental factors
and a risk associated with spending too much money on data gathering and analyses.
EVALUATION PROBLEMS
Critical to the evaluation process is understanding the models: modeling the right
concepts with the appropriate methodology. The problems come when:
trap volume
net/gross ratios
porosity
permeability
HC saturation
oil and gas fractions
HC fill
oil shrinkage factors
free gas volume factors
recovery factors
Care should be taken when assigning values to each of the assessment factors.
Assistance from other disciplines can provide better answers to such factors as whether
the factors are dependent or independent or what an informed estimate of a value might
be. Even with the great advances in geostatistical approaches and other new methods it
is important to remember the old computer mantra “garbage in, garbage out”.
Cumulative probability distribution of sand thickness at the proposed location. Each asterisk represents sand
thickness for one map. Vertical line intersects the distribution at the mean (45 ft) which equals the KED solution.
In this example, the median sand thickness expected is 465 feet and there is and 80%
probability of having 32 feet or more of sand at the modeled location. This type of
analysis can be conducted for any location on the map or for any number of
parameters. Another possibility with the data set is to plot the economic threshold for a
specific location and then directly read the probability of exceeding that value. The next
illustration demonstrates that approach:
This presentation indicates that there is a 35% probability that the sand thickness at the
proposed location is less than the economic threshold (40 feet). A risk map could
demonstrate the probability of reaching the economic threshold for sand at all of the
locations on a map. This map could be the contoured probability of finding less that the
threshold. In the case Wolfe et al. presented there appeared to be a lot of error in sand
thickness map and offered suggestions of means to increase confidence in the
interpretation:
To quote them:
The narrower the probability density function, the less error or uncertainty.
Primary factors that control the narrowness of the probability density
function or geostatistical estimate are:
Correlation coefficient between the seismic and well data, where the
higher the correlation coefficient, the narrower the probability density
function (PDF)
Variogram model range of the well data, where the longer the range, the
narrower the PDF
Distance of the proposed location to the nearest well control, where the
lesser the distance, the narrower the PDF (if it is at an existing well, the
PDF is a spike)
.021 21
.014 14
.007 7
.000 0
29.33 38.75 48.16 57.58 67.00
BCF
The net present value is the summation of the present values (income and expenses)
over all years of the project. Note that always a reference date needs to be specified.
In this analysis, concepts such as money of the day (simply the money one can produce
from his/her purse), real term money (money that maintains its spending power and
therefore requires a reference date) and inflation (devaluation of money over time) are
relevant.
25% Hi Case 80
EMV 1 Dev’mnt 50% Mid Case 50
No info (mid
case) 25% Low -20
Case
EMV 1 = 0.25 x 80 + 0.5 x 50 + 0.25 x -20= 40
Concepts:
Acquiring knowledge reduces uncertainty (requires investment/delay)
Increases value
IMPACT OF APPRAISAL
RC EXAMPLE
FIELD ASSUMPTIONS
Discount factor 12
Recoverable reserves (MBO) 1000
Appraisal wells (cost $10M/well) 4
Reservoir studies (2 year expenditure - $M) 11.5
3D seismic (2 year expenditure - $M ) 26
Development wells (cost $8M/well) (cost $10M/well) 16
Expected well productivity (BOPD) 20,000
Development drilling program – years 4
Facilities (2 year expenditure - $M ) 250
Oil price margin (after tax $) 3
SCENARIO RESULTS:
1. Appraisal only $1057M
2. Appraisal + reservoir studies $1097M
3. 3D + appraisal + reservoir studies $1235M
4. Perfect knowledge $1605M
SUMMARY OF CHAPTER 10
A Reservoir Characterisation study is only finished when the economics and risks of the
recommendations have been fully assessed. In other words, economic analysis is an
integral part of the Reservoir Characterisation study.
Cash surplus and accounting profit are different quantities. In cash surplus CAPEX
appears, whereas in accounting profit depreciation plays a role.
The uncertainties in the evidence on which basis the decisions are taken, is usually
significant giving rise to risks in the decisions. Data acquisition or appraisal may reduce
the uncertainties and hence the risks. However, data acquisition requires expenditure
and causes delay of the development, which must be justified by evaluating the Value of
Information (VOI).
TEAM STYLES
TEAM BUILDING:
MANAGING QUALITY
Traditional approach
Quality through inspection
Team-based approach
Quality through adherence to high standards
Lower overhead
High profitability
OBSERVATIONS
Performance challenges create teams
Permanent teams: continue for long periods and take on a continuing series of
assignments.
Most reservoir characterization teams fit into one of two relationship categories:
ORGANIZATIONAL MANAGEMENT
The organization may take a number of approaches to forming teams. The methods
will vary with the company's culture, the setting and the nature of the field to be
described. Some of the approaches require reducing the controls that the company has
over its own destiny through the inclusion in a partnership or an alliance. The benefits
of any approach should outweigh the costs of doing the work in the new way (the costs
include both the financial and intangible costs associated with the change process).
Approaches that companies have taken include:
The bottom line is forcing many oil companies to redefine their core activities. Oil and
gas companies are reviewing costly research and development programs and retaining
only proprietary technology that gives them a competitive advantage. The rest is being
outsourced to service companies or university consortia. Many companies see less
advantage in owning technology, more in being able to access it quickly and to use it to
improve efficiency.
Outsourcing has the added advantage of eliminating duplication and raising the
industry's overall efficiency. The trend is toward closer relationships with the service
industry. This provides dramatic opportunities for improved communication. The oil
company can guide and follow service industry developments from beginning to end
and, conversely. the service industry gets the chance to understand a client need from
the inception of each project. Partnerships greatly facilitate success in the new world of
cooperation.
A natural outgrowth of this trend that has been growing for a long time is the formation of
formal agreements to jointly conduct reservoir characterization projects through a strategic
alliance. These alliances typically consist of the participating partners, service companies,
suppliers and possibly regulatory agencies. For these approaches to be successful, a
number of elements need to exist. A strong alliance or partnership must have:
It is apparent from this list of criteria that companies cannot simply form a strategic
alliance on an ad hoc basis; the formation of successful alliances requires a significant
amount of planning and commitment and development of skills that have not been
required for development projects in the past.
On their way to becoming mature and fully functioning, reservoir characterization teams
go through a predictable five-stage transition process that often lasts from three to five
years.
Stage 1: KICKOFF
Stage 2: STORMING
Stage 3: NORMING
Stage 4: FORMING
Stage 5: INTEGRATED
SUPER TEAM
Management Shares team Overwhelmed by Pulls away from Conflict with team, Coordinates,
model confusion and needs, daily op.s, offers focus on org. vision facilitates cross team
threatened guidelines comm.
Support Skeptical Feel left out of process, Relat. improve as High level of coop. w/ Active in offering
Groups about teams confused auth. is defined team integrated sol.s
Organization Traditional Ident. systems that Modifies: perf. Monitors team Expands team
approaches need to be changed eval., auth., prof. activity and adjusts approach throughout
growth, rewards, as needed organization
training, discipline
systems
CHALLENGES OF STAGE 1
Overall focus - Sharing info about:
model
expectations
identifying issues
CHALLENGES OF STAGE 2
Overall focus - identifying confusion and building plans for:
managers letting go
team empowerment
identifying systems to modify
CHALLENGES OF STAGE 3
Overall focus - becoming a unit
leadership exertion
defining authorities with management, support, other teams
trust building
CHALLENGES OF STAGE 4
Overall focus - effective team, fully integrated, rotational task leadership
working within the organization
sharing results with other teams
quality management
CHALLENGES OF STAGE 5
Overall focus - developing new organizational norms
sharing ideas up, across and down
participation in strategic initiatives
TRADITIONAL ORGANISATION
MATRIX ORGANISATION
MATRIX TEAMS
Work is conducted through team
DEDICATED TEAMS
Reporting is solely through team structure
Gas plant
operating at 1/2 capacity
considering merger with Chevron facility
Successful project
Production increases, faster turnarounds, better drilling/completions,
employee buy-in, faster hookups
Kick-off
Understanding organization and team model, high trust level
Confusion
Conflict resolution, commitment to team and project
Dependency
Multi-disc. Skills, planning, decision-making
Team loyalty
Sharing of leadership, org. Integration, team learning
Integrated
“Big picture” orientation
Girgis et al. presented a matrix that defines roles for an AMT (asset management team)
Leader:
THE LEADER
One member of the team may be formally designated as leader either by a
management or by the team itself. Unlike a manager, a leader may have no power to
discipline, punish or reward; the leader is usually empowered to do only three things:
preside over meetings, see to it that the team completes its assignment on time and
report the team's conclusions to higher management. Nevertheless, even though the
leader's role is procedural, not hierarchical, the leader can gain the strong respect of the
team and emerge as a genuine inspirational leader.
The process: A team without a designated leader or an informal leader can still
have leadership. The leadership is provided on a situational basis by the use of
process skills, goal setting and problem skills by all members.
Leaders' responsibilities
Envision opportunity - identify and anticipate change
Portray the vision and mission both internally and externally and articulate new
strategic directions
Inspire total belief and commitment to the new directions to make the team
culture responsive
Strive to develop and lead the team
Achieve measurable results
Communicate with the rest of the organization
The role of the leadership and the management of the organization is to:
MANAGEMENT'S JOB
Prepare the environment, discipline to maintain direction
Select and empower talented people
Create an inventory of challenging problems
Create the vision for the future
CANDID FEEDBACK
Clean feedback:
- Specific
- Focus on behaviour
- Descriptive
- Designed to help the receiver of the feedback
Dirty feedback
- General
- Guessing motives of receiver
- Prescriptive
- Designed to help the giver
RECEIVING FEEDBACK
Clean:
- Listen
- Be silent
- Ask questions for clarification
Dirty
- Argue
- Justify
- Attack / defend
- Rationalise (“rational lies”)
- Withdraw
- Seek revenge
SERVING LEADERSHIP
Vision owned by the leader
Cannot be disputed!
Bottom LINE
Embracing computer technologies can increase efficiencies and provide a positive economic
impact through integrated reservoir characterization studies.
PROBLEM ADDRESSED
Many oil and gas fields in mature basins have not benefited from interdisciplinary char-
acterization studies. These reservoirs provide an opportunity for innovative independent
operators to maximize the economic potential through integrated teams using low-cost PC-based
software packages.
TECHNOLOGY OVERVIEW
Reservoir heterogeneity, often due to compartmentalization, has a profound impact on oil and
gas recovery. Heterogeneities occur at the wellbore, interwell, and on field-wide scales. Types of
compartrnentalization can include: stratigraphic sedimentologic, structural, petrophysical/fluid,
diagenetic, pressure, and combinations. The use of integrated, multidisciplinary teams has
proven to be an effective way to increase recovery and profits. Experience confirms that putting
the right people with the right skills (and tools) in the right projects—with firm management
support—is the key to success.
using scanners; vector data is input using a digitizing tablet or stylus. The average cost of
digitizing well data in-house is $42 to $57 per well.
FIELD RESULTS
Field studies have demonstrated several ways to use this digitized data. Examples of integrated
reservoir characterization are described below for two fields:
Sorento Field in Colorado produces primarily from flu-vial sandstones deposited in incised
valley fills cut into marine mudstones of the Morrow Formation (Pennsylvanian) on the eastern
flank of the Denver Basin. To understand compartrnentalization of the field, a reservoir
characterization study was conducted integrating geology, geophysics, and petroleum
engineering information.
Sequence stratigraphy was used to define the development of reservoir architecture. The fluvial
architecture is highly heterogeneous and compounded by several minor fluctuations of sea level
that stacked and partially eroded individual depositional stages. High-frequency 3-D
compressional seismic data were used to distinguish the incised valley outline and lateral barriers
to flow. Low-permeability carbonate-cemented zones and floodplain deposits are associated with
stacked channel sequences within the valley. These deposits created baffles and barriers to fluid
flow that were mapped using high-frequency 3-D compressional data.
When the distribution of diagenetic features was combined with sandstone geometry, the
reservoir compartments became apparent. Production data were used to corroborate geological
and geophysical interpretations. It was found that the reservoir is comprised of four flow units
that display separate fluid contacts but maintain a baffled pressure relationship. This implies that,
in developing similar fields, drilling a location structurally lower than a known oil-water contact
may be justified. A large 3-D seismic survey may pay for itself in a complex reservoir setting if
it can decrease the number of dry holes. Based on a dry hole cost of $115,000 per well and a
seismic cost of $50,000 per square mile, a 3-D survey of five square miles that prevents two dry
holes would cost the same as drilling them. If more than two dry holes were prevented, which is
likely in a large field, the 3-D survey is even more economic. Also, the operator gains
understanding that should result in more efficient development and management of the field.
Slattery Field in ~4~oming contains approximately 65 wells and produces primarily from
aeolian sandstones of the Minnelusa Formation (Permian) on the eastern flank of the Powder
River Basin. The study of South Slattery sought to: maximize oil recovery, demonstrate an inte-
grated approach to reservoir characterization, and create a “cookbook” that independent
operators could use to develop their own methodologies for reservoir characterization.
Geological work for this study was done using a desktop computer loaded with GeoGraphix GES
(well data base manager, mapping and cross section software) and GeoQuest QLA2 and RIS
LESA (log analysis software programs). Digital log data was used, even though it tends to be
expensive to acquire.
The geologist and engineer were physically located in facilities 1,000 miles apart and relied upon
e-mail for communication and file transfer. This arrangement illustrates that operators can use
remote expertise without having to hire additional staff to complete reservoir characterization
work.
Most of the time was spent mapping, correlating, and working on cross-sections. While the
geologist had a working knowledge of the software at the beginning of the project, nearly 25% of
the time was spent learning even more about the software. If the same software were used in
future projects, the learning curve should decrease.
If digitized log data had been purchased, about 20% of the geologist’s time could have been
conserved. The amount of time the geologist spent preparing for simulation (grid preparation)
was small, suggesting that simulation would not add a significant time commitment.
Based on this study, a realistic estimate for the time a geologist needs to characterize a reservoir
is approximately one well per day. Using the mean salary for a petroleum geologist ($300 per
day according to the American Geological Institute), the cost of a geologist’s time for reservoir
characterization is $300 per well. The cost of software (amortized over the first year. with one
well per work day and 250 work days/year) is about $82 per well. The cost of data works out to
be $27 per well. Thus, the total cost of developing the geological model for reservoir
characterization is approximately $400 per well.
The highest cost consideration most likely is the geologist’s time, which is three times that of the
data and software combined. Efforts to save money should, therefore, be directed at equipment,
software, and other methods to make the time spent by the geologist more productive. Simple
suggestions for increased geological efficiency include purchasing data and using geotechnical
assistants.
Based on a workshop sponsored by PTTC’s Rocky Mountain Region and the American
Association of Petroleum Geologists, in Denver, CO.
Performance appraisals*
Reward systems*
Simplifying systems*
Effective communications
Clear roles /responsibilities / assignments
Fostering innovation
Asking for appropriate help
Need-based training
TRAINING CATEGORIES
Interpersonal skills training
Without the conventional channels of authority for giving and receiving direction, work
team members need to talk with, explain to, agree with, disagree with, decide, listen to
and convince more people more frequently than ever before. In other words, they need
to be skilled communicators - both one-on-one and in groups. Interpersonal skills
training starts with the basics, like listening, giving feedback and problem solving. As
the teams evolve, training moves on to more advanced skills, such as influencing
others, handling diversity, resolving conflicts and peer reviews.
TEAM COMMUNICATION
Define the problem (by identifying the symptoms and the probable cause); list
the criteria of a good solution
Generate trial solution
Evaluate each trial solution in the light of pertinent data
Choose the best of the trial solution in the light of pertinent data
Develop an implementation plan and arrange follow-up
Brainstorming
A frequently used term in analyzing issues is "brainstorming". Frequently the term is used
informally to indicate a general discussion of an issue. A more rigorous application of the
technique is shown below. This very useful technique can yield a great many useful
solutions. It is recommended that teams use the more formalized approach until they
become more skilled at brainstorming. Frequently individuals are assigned additional
responsibilities (recorder, facilitator, and spokesman). The steps are:
Have each team member individually write down as many ideas as possible that
might solve the problem or answer the question.
Record and post all ideas without discussing or evaluating any them.
Encourage people to build on other ideas and continue recording and posting all
ideas.
Clarify the ideas. Give everyone a chance to look at the ideas and ask questions
about the meaning of any idea. Do not evaluate any of ideas at this time.
Roles of Facilitator
1. Help keep the group on task
2. Watch the time to assure that the group completes its task
3. Assure that everyone is able to participate - no one dominating, no one
excluded
4. Remind people to listen as others are talking
5. Encourage people to respect and use their different perspectives and views
Roles of Recorder
1. Listen for key words; do not edit - use exact words
2. Capture the basic ideas
3. Write rapidly, legibly (1 to 1-1/2 inches high)
4. Number each sheet; reference topic, group name, etc.
Selection grid
A method for selecting one option out of several possibilities. It involves deciding what
criteria are important and using them as a basis for reaching a acceptable decision.
Impact analysis
A procedure to discover what impact a situation has on the organization. It involves
getting specific impacts to the organization and other available information about the
decision.
Letting team members share their viewpoints and come to agreement about what
exactly is happening
Suggesting how things should be
An illusion of invulnerability
A tendency to discount warnings
Unquestioning belief in the team
Stereotyped views of other groups
Pressure on any member who contradicts the team consensus
Self-censorship by individuals
An illusion of unanimity
SUMMARY OF CHAPTER 11
The dedicated team structure is very suitable for the Reservoir Characterisation
process, because it facilitates the necessary inter-disciplinary sharing that is so
important for optimum RC. Both team members and management need to have the
appropriate skills to make dedicated teams successful. Suitable training is required for
this purpose. Three training categories need to be addressed, namely interpersonal,
administrative and technical training.
Project management skills are much needed in the reservoir characterization process.
Much of the work previously done by management personnel will be passed along to
the reservoir characterization team. The skills include traditional project management
skills plus the ability to understand the whole of the project within the context of the
company’s overall strategies and tactics. The process consists of a number of elements
including:
There is computer software available to facilitate these activities, but the task confronting
a reservoir characterization team charged with a large project can be quite daunting. Bob
Sneider demonstrates the activities that a reservoir characterization team needs to
undertake to accomplish the reservoir characterization goals:
– How do the goals, roles, procedures and relationships affect people's work?
With most reservoir characterization projects, one of the most difficult activities will be to
bring together a team that consists of representatives from several partner companies,
service companies, suppliers and the government representatives (Refer to the section
above that outlines the methods for establishing an effective strategic alliance.)
Appropriate milestone along the way to completion of the project and rewards
appropriate to the contributions expected are significant success factors. Paramount
among success factors are the focus on the team product and the adherence to the
shared vision for success of the project.
Tyraskis, et al. have charted the RC process for Saudi Aramco. Their conclusions
provide insights to the project planning process:
Improving reservoir characterization has a direct positive impact on the company’s core
mission — the more efficient production of hydrocarbons. The perceived benefits to the
Saudi Aramco E&P organization include the following:
Better planning model for the acquisition and allocation of key resources,
such as computer hardware/software and professional expertise.
A detailed knowledge of the time, physical resources, and people necessary for these
projects also provides a better platform from which to make decisions concerning
outsourcing of related work.
The existing model was constructed in the most general way possible to cover the
widest range of contingencies. The methodology employed in defining the IRCPM is
easily adaptable to other organizations and processes. Consequently, it does not apply
in detail to any current project, although it does come closer to some than others. The
next step in the implementation of a more rigorous planning methodology is to apply the
existing framework to a specific project and to predict project milestones, resources and
costs (Tsingas et al.).
The chart shown below illustrates an outline for the processes of a typical reservoir
characterization project. Larger and more complex projects or those with multiple
partners will be more complex. The basic outline consists of definition of the problem
usually by a management group and more senior discipline leaders followed by
development of a preliminary plan. Resources and constraints help to define the plan
that is the next step. An important and frequently forgotten step is to identify all of the
stakeholders in the reservoir characterization process; stakeholders are any parties
that are affected by the project. Typically include are internal stakeholders
(management, professionals, support staff), partners, service companies, suppliers,
environmental groups and the government agencies. Successful plans take into
account the needs of the company plus the needs of these stakeholder groups.
The plan for the reservoir characterization project should always be considered a
working document subject to change based on new data or conditions. It is very
important for the team members to have a significant hand in developing the plan which
should include not only activities, but should also define working relationships on the
reservoir characterization team and its relationship to management. The most
successful reservoir characterization teams are those with a clearly-defined and shared
vision for the project; the team members should be empowered to make routine daily
decisions and more substantial decisions in their own discipline - in other words to act
as the experts that they are. The shared vision for the project is the driving force behind
the work that is done; it is the unifying factor that all decisions should revolve around.
After development of the plan, the more technical activity of reservoir characterization
begins. Each discipline member should be charged with developing the best model
tested by the most recent techniques within the constraints of the project (time,
resources, data, and financial constraints). As has been stated many time before in this
manual, cross discipline coordination will enhance the final product immensely and
will lead to significant efficiencies. Each team member will need to develop and
exercise new ways of doing business.
Finally, the team will be given the opportunity to present its findings and
recommendation; these findings should come as little surprise since perceptive
managers will already have sampled the findings on an ongoing basis throughout the
project. The successful completion of the reservoir characterization project is only the
first step in a long multi-disciplinary approach to producing the most reserves at the
least cost.
VALUE CREATION
The planning process described above seems long and cumbersome, but can save a
great deal of time later in the project. Even though there are a great many steps to the
planning process, these are accomplished rather quickly –m particularly if experienced
professionals are accessed in a collaborative manner. The steps that follow are the
heart of the RC process and comprise the bulk of the time in characterizing reservoirs.
The process follows from development of a static model to development of a dynamic
model to test the model. Many of the steps in the process are contemporaneous and
not sequential. Note that there are many iterative loops built into the RC process so
that interpretations can be enhanced based upon other data sets or subsequent
interpretations.
Petrophysical data
Normalize/calibrate/edit against cores
Identify cutoffs for mapping
Core data
Fully describe/interpret cores
Correlate with logs
Calibrate with test/production data
Seismic data
Develop structural interpretation
Seismic stratigraphic correlation
Establish correlations with reservoir parameters
Define facies
General Characteristics:
Lithology Porosity amt. Porosity type Permeability
Interp. Method:
Net/Gross Water sol.
Interp. Method:
Texture Fabric
Interp. Method:
Analogues:
Modern
Stratigraphy Observed (lithology, source, seal, reservoir?):
Overlying zone:
Undelying zone:
Ancient
Stratigraphy Observed (lithology, source, seal, reservoir?):
Overlying zone:
Undelying zone:
Reservoir Dimensions
Modern Min ML Max Ancient Min ML Max
Length Length
Width Width
Thickness Thickness
Heterogeneities
Stratigraphic Modern Min ML Max Ancient Min ML Max
Length Length
Width Width
Thickness Thickness
Production Profile
Analogue Expected
Critical Elements
Analogue Expected
Production data/injection
Fluid data
Develop flow unit correlations
Who - which people, resources - are required in order to accomplish those tasks?
(roles)
How do the goals, roles, procedures and relationships affect people's work?
(interpersonal)
GOAL SETTING
Only after the team has considered its vision and mission statement, identified its
needs, identified available and required resources, should work begin on goal setting.
Any move to set goals without the appropriate background is premature.
CHARACTERISTICS OF GOALS
Specific
Measurable
Achievable
Realistic
Timely
Realistic
Requires optimal effort
Specific
Comprehensive
Understandable
Management sets overall goals for the organization compatible with its shared
vision and available assets (technical, financial, people)
The reservoir characterization team writes goals for the project based on the
management guidelines
The management may outline goals independently of the reservoir
characterization team
During a working session, the reservoir characterization team presents a plan
that is accepted or modified based on information from the management about
priorities, budgets, resources, etc.
The reservoir characterization team and management agree to a performance
contract with specific action plans to achieve the goals and periodic review
sessions (including specific dates and times during which the goals will be
modified as needed or reviewed against performance milestones)
At the end of the period the reservoir characterization team evaluates itself (the
management will also conduct an evaluation, but will not share that evaluation
except where disagreements occur)
During a forward-looking meeting, the management and reservoir
characterization team will review the project, analyze opportunities and threats,
strengths and weaknesses - this is a good time to reaffirm the shared vision and
to discuss professional growth opportunities for individuals in a constructive,
collaborative manner
This process works equally well with individuals as well as with teams; it focuses the
accountability to perform directly on the individual (or team) and builds an ownership in
performance and output driven work.
TYPES OF GOALS
There are many types of goals that can be used with teams and individuals:
Continuing challenge goals - (example - bring field recovery rate to 45%, etc.)
Moving target goals - something to aim for, but is probably not achieved
ACTION PLANS
In order to attain goals, specific action plans are needed. Action plans should include
timetables, personnel responsible for each task, tasks to be accomplished and the order
in which the tasks are to be completed. These plans will also provide a method for
evaluating progress toward the goals. Project planning software can assist the process
significantly.
SUMMARY OF CHAPTER 12
Successful projects create value for the organization. Most value is created in the early
stages when projects are envisaged and designed.
Projects are carried out for the benefit of the organization. Therefore, clarity about the
projects carried out must be visible in the organization.
PROBLEMS, PITFALLS
Some of the major pitfalls and problems that may be encountered are discussed below.
Dependence on technology
Good geologic and engineering reasoning are more important than the dependence
upon the latest technology without a thorough understanding of the technology -
there are no “black box” solutions.
Understanding scales
As indicated above, the reservoir characterization team must be using common
terms with common definitions. Care in understanding the scale represented by a
data set will increase understanding and avoid misunderstandings.
Fractures
The influence of fractures in reservoirs is just beginning to be understood, but the
knowledge is not widespread or well understood.
Diagenesis
The influence of diagenesis is profound is many reservoirs, but many professionals
simply map primary depositional environments and then apply interpretations based
on those primary environments.
Dependency of risk
Not understanding the significance of dependency and independence in evaluating
risks can yield significant errors.
Flow physics
Another area that is yet not adequately understood. Significant potential exists to
improve interpretations with concentration on identifying and understanding flow units.
Neglecting priorities
The tendency to solve every problem to the ultimate degree is wasteful of time and
effort. Some activities require a 95% solution others need only an 80% effort;
identifying priorities is a critical project management skill.
Reservoir characterization will continue to bean important tool for developing oil and gas
accumulations. The future will probably contain a great number of technical
advancements - particularly in the manipulation, integration and visualization of data -
as well as organizational efficiencies. Advances in our understanding of:
Organizational advances have bees underway for many years; continuation of the trend
to recognize the value, contribution and impact of the working professionals will
enhance the reservoir characterization effort. Outsourcing and the formation strategic
alliances will expand the cross-disciplinary communication demands on working
professionals. A new set of skills will have to be developed at all levels in the
organization. In summary, a combination of procedural, organizational and technical
advances will drive the oil and gas companies of the next century.
Future growth was summarized in the training manual Oil and Gas Teams: How to
make them work (Treesh):
the efficiency of the individuals who are performing the critical tasks.
Traditionally there has been a great deal of attention paid to the upgrading
of technical applications in the search for, development of and production
of oil and gas reserves. In addition, in recent years there has been a
dedicated effort in most companies to assure that the most "efficient"
processes are applied to each phase of exploration/production. The area
that has received little attention is that of the application of human
resources; team building through its associated human efficiency
elements is an important method through which to significantly improve
the human element. …concurrent improvement in each of these area can
yield significant magnitudes of improvement in the productivity and
profitability of oil industry companies.
- The triangle signifies that in the business decisions have to be taken under risk
based on evidence produced by activities in the Reservoir Charaterisation (RC)
process.
- The P-model says that in the business a profit needs to be made, that the business
must have a purpose, that in the business three aspects need to be addressed,
namely the technical plan, the people issues and the organization with processes,
and finally that the business has to be sustainable.
- The steps represent the six steps in the RC process. Note that in step three the
agreement of the RC plan is concluded.
- In the RC process, critical elements are studied that support the decisions in an
optimum way as shown by the exclamation mark.
- The solution space is represented by the circle, in which multiple realisations (M1
etc.) are shown. If reality (R) is also within the solution space, then the RC modeling
will be meaningful.
- The inter-connected circles manifest that an integrated effort is required with
scientists and engineers from different disciplines as required for the particular RC
study. Note that not only the traditional subsurface disciplines should be present,
but also people from disciplines like legal, environmental etc. as required by the
study.
- The RC team should work as a super-team as shown by the graph with the different
phases in team formation.
- Last but not least, everything in RC must be SMART: specific, measurable,
achievable and agreed, realistic and time-bound.
SUMMARY OF CHAPTER 13
Reservoir Characterisation is continuously developing in order to meet the challenges of
the modern time and the future. As resources become more difficult to develop and
more marginal, continuous improvement in technology will be required.
GENERAL REFERENCES
REFERENCES CITED
Bashore, W.M., Araktingi, U.G., Levy, M., Schweller, W.J., 1994, Importance of a
Geological Framework and Seismic Data Integration for Reservoir Modeling and
Subsequent Fluid-Flow Predictions pp. 159-176, in Stochastic Models And
Geostatistics: Principles, Methods, And Case Studies, edited by Yarus, J.M. and
Chambers, R.L, edited by Yarus, J.M. and Chambers, R.L.
Begg, S.H., A. Kay, E.R. Gustason and P.F. Angert, 1996, Characterization of a
Complex
Fluvial-Deltaic Reservoir for Simulation, SPE Formation Evaluation, Sep. 1996, p. 147-153.
Bracco, Guido L Gartner and Wolfgang Schlager, 1999, Discriminating Between Onlap
and Lithologic Interfingering in Seismic Models of Outcrops, AAPG Bulletin, V 83.
No 6 (June 1999), P. 952-971.
Chambers, R.L., Zinger, M.A., and Kelly, M.C., 1994, Constraining Geostatistical
Reservoir Descriptions with 3-D Seismic Data to Reduce Uncertainty, pp.143-158,
in Stochastic Models And Geostatistics: Principles, Methods, And Case Studies,
edited by Yarus, J.M. and Chambers, R.L.
Cross, T.A., Field-Scale Reservoir Characterization. 1986, Pp. 493-496, in Reservoir
Characterization II, edited by Lake, L.W., Carroll, Jr., H.B., and Wesson, T.C.,
Academic Press.
Dehghani, Kaveh, Paul M. Harris, Kelly Edwards, and William T. Dees, 1999, Modeling
a Vuggy Carbonate Reservoir, McElroy Field, West Texas, AAPG Bulletin, V 83.
No 1 (January 1999). P 19-42
Deutsch, C.V., and Journel, A.G., 1994, Integrating Well Test-Derived Effective
Absolute Permeabilities in Geostatistical Reservoir Modeling, pp. 131-142, in
Stochastic Models And Geostatistics: Principles, Methods, And Case Studies,
edited by Yarus, J.M. and Chambers, R.L.
Dubrule, Olivier, 1998,Geostatistics in Petroleum Geology; AAPG Continuing Education
Course Note Series #38Eschard, R, P. Lemouzy, C. Bacchiana, G. Desaubliaux, J.
Parpant and B. Smart, Combining Sequence Stratigraphy, Geostatistical
Simulations, and Production Data for Modeling a Fluvial Reservoir in the Chaunoy
Field (Triassic, France), AAPG Bulletin. V 82, No 4 (April 1998).
Girgis, John, Sneider, R.M., and Thomas, B., 1995, Multi-disciplinary Teams - What is
the “Right” Structure? Based on Ten Years of MTS’s in Indonesia, pp. 545-564,
Indonesian Petroleum Association, Proceedings.
Grotsch Jurgen and Christophe Mercadier, 1999, Integrated 3-D Reservoir Modeling
Based on 3-D Seismic: The Tertiary Malampaya and Camago Buildups, Offshore
Palawan, Philippines; AAPG Bulletin, V.83, No. 11 (November 1999), P. 1703—
1728.
Haldorsen, H. H. and. Damsleth, E, 1990, Stochastic Modeling, pp. 404—112, Journal
of Petroleum Technology, (April 1990).
Haldorsen, H.H., and Chang, D.M., 1986, Notes on Stochastic Shales; From Outcrop to
Simulation Model, pp. 445-486, in Reservoir Characterization edited by Lake, L.W.
and Carroll, Jr., H.B., Academic Press.
Hardage, Bob A,. Carr, David L Robert, Finley, J., Lancaster, David E., Elphick, Robert
and Ballard, James R, 1995, Secondary Natural Gas Recovery: Targeted
Applications For Infield Reserve Growth In Midcontinent Reservoirs, Boonsville
Field, Fort Worth Basin, Texas; GRI-95/0454.2, VOLUME II: APPENDIX
Hardage, B.A., and Fertl, W.H., 1991, Synergism in Formation Evaluation Enhances
Detailed Reservoir Description, pp. 123-129, in The Integration of Geology,
Geophysics, Petrophysics and Petroleum Engineering in Reservoir Description and
Management, Proceedings of the First Archie Conference, Oct 22-25, 1990,
published by AAPG.
Hart, Bruce S, 1999, Geology plays key role in seismic attribute studies, Oil and Gas
Journal, July 12, 1999, P 76-80.
Hart, Bruce S., David M. Sibley, and Peter B. Flemings, 1997, Seismic Stratigraphy,
Facies Architecture, and Reservoir Character of a Pleistocene Shelf-Margin Delta
Complex, Eugene Island Block 330 Field, Offshore Louisiana, AAPG Bulletin, V
81.No 3 (March 1997). P 380-397
Hohn, Michael Ed., Ronald R. McDowell, David L. Matchen, and Ana G. Vargo, 1997,
Heterogeneity of Fluvial-Deltaic Reservoirs in the Appalachian Basin: A Case
Study from a Lower Mississippian Oil Field in Central West Virginia, AAPG Bulletin.
V 81. No. 6 (June 1997). P. 918-936
Holland, D.S., Leedy, J.B., and Lammlein, D.R., 1990, Eugene Island Block 330 Field -
USA, Offshore Louisiana, pp. 103-144, in Structural Traps III, Treatise of
Petroleum Geology, Atlas of Oil and Gas Fields, published by AAPG.
Hoover, Andrew R., Burkhart, Tucker, and Flemings, Peter B.; 1999, Reservoir and
Production Analysis of the K40 Sand, South Timbalier 295, Offshore Louisiana,
with Comparison to Time-Lapse (4-D) Seismic Results; AAPG Bulletin, V. 83, No.
10 (October 1999), P. 1824—1641.
Jones, Alistair, McGill, Colin, and Williams, John, 1993, Renormalization: a Multilevel
Methodology for Upscaling, pp. 823-828, in Reservoir Characterization III, edited
by Linville, Bill, Pennwell Books.
Journel, A.G.,1994, Geostatistics and Reservoir Geology, pp. 19-20, in Stochastic
Models And Geostatistics: Principles, Methods, And Case Studies, edited by
Yarus, J.M. and Chambers, R.L.
Kuhn, Tibor, and Kolosi, Zsolt, 1993, How to Estimate the Reliability of Reserves? A
Quick and Practical Solution, pp. 769-780, in Reservoir Characterization III, edited
by Linville, Bill, Pennwell Books.
Kupfersbergerz, H. and C. V. Deutsch, 1999, Methodology for Integrating Analog
Geologic Data in 3-D Variogram Modeling, AAPG Bulletin. V. 83. No 8 (August
1999) P 1262- 1278
Lasseter, T. J., Waggoner, J. R. and Lake, L. W., 1986, Reservoir Heterogeneities and
Their Influence on Ultimate Recovery, in Reservoir Characterization edited by
Lake, L.W. and Carroll, Jr., H.B., Academic Press.
Lia, Oddvar, Henning Omre, Hakon Tjelmeland, Lars Holden, and There Egeland, 1997,
Uncertainties in Reservoir Production Forecasts, AAPG Bulletin. V 81.No 5 (May
1997). P 775-802
Loren J.D, Kuhla, J.T., Renbarger, K.S., and Sneider, R.M., MacDonald, A. C. and
Aasen, T.O., 1994, A Prototype procedure for Stochastic Modeling of Facies Tract
Distribution in Shoreface Reservoirs, pp.91-108, in Stochastic Models And
Geostatistics: Principles, Methods, And Case Studies, edited by Yarus, J.M. and
Chambers, R.L.
Sailer, Arthur F, Dickson, J.A.D. and Matsuda, Fumiaki, 1999, Evolution and Distribution
of Porosity Associated with Subaerial Exposure in Upper Paleozoic Platform
Limestones, West Texas; AAPG Bulletin, V.83, No. 11 (November 1999), P.
1835—1854
Satter, Abdus, and Thakur, Ganesh, 1994, Integrated Petroleum Reservoir
Management: A Team Approach, Pennwell Books, 335p.
Sneider, R. M., 1996, Multidisciplinary Teams: How and Why They Make Your Money,
AAPG Continuing Education Course Notes #37, 140 p.
Sneider, R. M., and Johnson, H.D.,1996, Development Geology, . Training manual
published by OGCI.
Sprunt, E.S., 1993, Quality Control of Special Core Analysis Measurements, pp. 254-
269, in Reservoir Characterization III, edited by Linville, Bill, Pennwell Books.
Srivastava, R. Mohan, 1994, An Overview of Stochastic Methods for Reservoir
Characterization, pp. 3-16, in Stochastic Models And Geostatistics: Principles,
Methods, And Case Studies, edited by Yarus, J.M. and Chambers, R.L.
Stalkup, F.I., 1986, Permeability Variations Observed at the Faces of Crossbedded
Sandstone Outcrops, pp. 141-180. in Reservoir Characterization edited by Lake,
L.W. and Carroll, Jr., H.B., Academic Press.
Tillman, R.W. and Pitman, E.D., 1991, Reservir Heterogeneity In Valley-Fill Candstone
Reservoirs, Southwest Stockholm Field, Kansas, pp. 51-105, in Reservoir
Characterization III, edited by Linville, Bill, Pennwell Books.
Tinker, C.N., Chambers, L.D., Ritch, H.J., McRae, C.D., and Keen, M.A., 1991, The
Benefits of a Synergistic Approach to Reservoir Characterization and Proration,
Rose City Prairie du Chien Gas Field, Ogemaw County, Michigan, pp.- 360-372, in
The Integration of Geology, Geophysics, Petrophysics and Petroleum Engineering
in Reservoir Description and Management, Proceedings of the First Archie
Conference, Oct 22-25, 1990, published by AAPG.
Tobin, R.C., 1997, Porosity prediction in frontier basins: a systematic approach to
estimating sub- surface reservoir quality from outcrop samples, in JA. Kupecz, J.
Gluyas, and S. Bloch, Reservoir quality prediction in sandstones and carbonates:
AAPG Memoir 69, p. 1-18
Tucker, Karla E., Paul M. Harris, and Richard C. Nolen-Hoeksemn, 1998, Geologic
Investigation of Cross-Web Seismic Response in a Carbonate Reservoir, McElroy
Field, West Texas, AAPG Bulletin, V 82, No 8 (August 1998), P. 1463-1503
Treesh, M. I., 1995 revised 1999, Oil and Gas Teams: How to Make Them Work.
Training manual published by OGCI.
Tyler, K., Henriquez, A., and Svanes, T., 1994, Modeling Heterogeneities in Fluvial
Domains: A Review of the Influence on Production Profiles, pp. 77-90, in
Stochastic Models And Geostatistics: Principles, Methods, And Case Studies,
edited by Yarus, J.M. and Chambers, R.L.
Tyler, K., Svanes, T, and Omdal, S., 1993, Faster History Matching And Uncertainty In
Predicted Production Profiles With Stochastic Modeling, SPE Paper 26420, SPE
Annual Technical Conference and Exhibition, Proceedings, p. 3143.
Tyraskis, P.A., C. Tsingas, B. B. Toelle, T. H. Keith, B.L. Nebrija and G. C. Robinson,
1996, Integrated Reservoir Characterization Process Model: a planning and