SPE 75223 Using Statistics To Evaluate A History Match
SPE 75223 Using Statistics To Evaluate A History Match
SPE 75223 Using Statistics To Evaluate A History Match
David Uldrich, SPE, BP Kuwait ; Saad Matar, Kuwait Oil Company; Hugh Miller, ChevronTexaco Kuwait
Copyright 2002, Society of Petroleum Engineers Inc.
This paper was prepared for presentation at the SPE/DOE Improved Oil Recovery Symposium
held in Tulsa, Oklahoma, 1317 April 2002.
This paper was selected for presentation by an SPE Program Committee following review of
information contained in an abstract submitted by the author(s). Contents of the paper, as
presented, have not been reviewed by the Society of Petroleum Engineers and are subject to
correction by the author(s). The material, as presented, does not necessarily reflect any
position of the Society of Petroleum Engineers, its officers, or members. Papers presented at
SPE meetings are subject to publication review by Editorial Committees of the Society of
Petroleum Engineers. Electronic reproduction, distribution, or storage of any part of this paper
for commercial purposes without the written consent of the Society of Petroleum Engineers is
prohibited. Permission to reproduce in print is restricted to an abstract of not more than 300
words; illustrations may not be copied. The abstract must contain conspicuous
acknowledgment of where and by whom the paper was presented. Write Librarian, SPE, P.O.
Box 833836, Richardson, TX 75083-3836, U.S.A., fax 01-972-952-9435.
Abstract
Three methods are proposed for quickly evaluating the history
match of a numerical simulation to actual reservoir
performance. All of the methods rely on computing a set of
deviation values, each of which is defined to be a calculated
simulator result minus the corresponding surveillance
measurement value.
For any particular type of surveillance data, such as rates,
watercuts, or gas-oil ratios, the deviation values can be
grouped by well, by area, or combining all measurements in
the database. The first two proposed methods rely on simple
graphical presentations of each group of deviation values to
show how well the simulation results match the surveillance
data. Plotting together the results from more than one
simulation run allows a quick comparison of the match for
each run, which is useful during the history match process.
The third method converts each deviation value to a
quantity called Match Factor, which is a relative measure of
the confidence that the simulator actually reproduced the
particular reservoir performance at the time the surveillance
measurement was made. Weighted-average Match Factors can
reveal the degree of match by well, by area, and by data type.
These techniques are especially valuable when matching
reservoirs with a large volume of surveillance data. They can
help focus the history matching process by identifying areas
less well matched. They can identify when the history
matching process is not significantly improving the match and
can stop.
Introduction
One of the more challenging aspects of history matching large
numerical simulators is assessing how well the simulator
results match observed field behavior. Unfortunately, there are
few objective measures of the degree of match readily
available.
The traditional approach is to plot the observed data values
versus time, along with the corresponding simulator output,
and visually assess how well the simulator reproduced the
measured values1. The quality of this evaluation can vary,
depending on the experience and judgment of the simulation
engineer. It can also be very time consuming, especially for
fields with many wells and years of surveillance data.
This paper describes a two-stage approach, developed by
BP Kuwait for Kuwait Oil Company, and validated on large
models of giant reservoirs in Kuwait. In the first stage,
observed surveillance data values are directly compared to the
corresponding predicted values extracted from simulator
output. Each pair of observed and predicted values defines a
deviation value. Groups of deviation values are presented in
two graphical ways, showing the overall degree of match for
the type of data in that group for that run. Plotting together
results from two or more simulation runs can be used to
quickly compare the matches for the runs.
In the second stage of the analysis, each deviation is
converted to a Match Factor value, which represents the
confidence level that the simulated result matches the actual
field behavior represented by the observed surveillance value.
This conversion takes into account the inherent uncertainty of
the field measurement technique, and the limits of the
simulator calculation. Plotting on a map the Match Factors
averaged by well quickly shows where a model is better
matched compared with other areas. This helps guide where
changes in the reservoir description should be made in
subsequent history match runs.
Averaging the Match Factors for all values in one data
type quickly provides a measure of the overall degree of
match for that data type for that simulation run. Comparing
this average to that from other runs, shows whether the history
match is improving or deteriorating. If the average Match
Factor value has stabilized, the useful end of the history match
process may have been reached (unless a significant alteration
in the reservoir description is made).
SPE 75223
SPE 75223
selecting the case with the median value and the overall curve
shape closest to zero deviation.
As the examples of Cases A, B and C show, by plotting
results from two simulation runs in a single Deviation
Distribution Plot (such as successive cases in the history
match process, or different sensitivity cases after history
matching), it quickly becomes clear whether one case is better
matched than the other. If there is no clearly superior case
during history matching, it may be an indication that the
process has reached its useful end. Finding more than one
reservoir description that matches past performance equally
well (such as when running sensitivity cases after history
matching), allows the simulation engineer to assess the impact
of reservoir description uncertainty on future predictions.
Future predictions of rates and recoveries from simulators
equally matched to past performance should be equally
probable, which is useful when defining the uncertainty range
for those future predictions.
History Match Evaluation Using The Deviation Band Plot
Finding positive or negative bias in the match is important in
comparing two history match cases and deciding what changes
in the reservoir description to make next in the history match.
But evaluating the total amount of deviation, regardless of
whether it is positive or negative, can be useful as well. This is
equivalent to repeating the method used to build the Deviation
Distribution Plot, this time using the absolute value of the
deviations instead of the actual values.
A convenient way to visualize this concept is to picture
each of the observation values plotted in Figure 1 with a
vertical bar defining the same deviation range above and
below each observed surveillance value. At each step in the
deviation band analysis, the length of the bars defines a
deviation band above and below the observed data, in the
way illustrated by Figure 4. For any particular deviation band
width (length of the bars), a certain number of the predicted
values from the simulation will fall within the deviation band.
If the width of the deviation band is zero, only those predicted
values which perfectly match the observed values will fall
within the deviation band. As the width of the deviation band
increases, a larger and larger percentage of the predicted
values fall within the band. When the band width equals the
maximum of the absolute values of all of the deviations, the
band will enclose 100% of all of the predicted values.
A Deviation Band Plot presents the relationship between
the width of the deviation band and the cumulative percentage
of the predicted values which fall within the band. Case A in
Figure 5 shows the same static BHP deviation data from
Figure 3, replotted in the Deviation Band Plot format.
Although usually described in a different way from that
presented here, the concept behind the deviation band has
often been used by simulation engineers to specify a
quantifiable requirement for declaring a simulator sufficiently
history matched. The requirement is usually expressed as a
target: to match a certain percentage of observations within a
certain range (positive or negative). In this example, 82% of
the predicted values are within 100 psi of the observed
SPE 75223
SPE 75223
SPE 75223
SPE 75223
rate more than half of the time (Match Factor less than 50).
Plotting these average match factors on a bubble map will
show whether the problem is well-specific (problem
wells widely scattered amongst other wells), or regional
(problem wells grouped into one or more regions). Finding a
regional trend would suggest there may be a problem with the
reservoir description, causing insufficient reservoir pressure
around the producers in some areas. In this way, Match Factor
analysis of the well control rate performance can provide
valuable data to assess a history match and guide reservoir
description modifications.
Example 5: Match Factor Analysis Of The Overall Match
(Comparing New To Previous Full-Field Simulation)
Because Match Factors are a relative measure of the degree of
history match, they can be compared with Match Factor values
from other data types directly. A bar chart is a convenient way
to summarize the history match evaluation, by presenting the
average Match Factors for each of the data types in the
surveillance dataset. When comparing the match for two
simulator runs, bars representing the average Match Factors
for each data type are plotted side by side.
For the large Kuwait oil field being evaluated here, the
surveillance dataset (for both the new and previous
simulations) included Total Liquid Production Rate (TLPR),
watercut (WC), gas-oil ratio (GOR) static bottomhole pressure
(SBHP), and RFT pressure (RFT). For each data type, Match
Factor curves were developed, using the concepts described
above. Since both simulator versions used essentially the
same surveillance dataset, the same Match Factor curves
were used to evaluate both. Weighting factors were not used in
this example.
The resulting averaged Match Factors for the previous and
new simulator runs (after each were history matched) are
plotted in Figure 13. The TLPR match is very good, as would
be expected for the data type used to control well rates in the
simulation. The watercut is also well matched in both cases,
reflecting the time and effort devoted to it during both history
matching efforts. The GOR match is slightly improved in the
new simulation, a result of improvements in the description of
fluid properties.
The SBHP Match Factor value is somewhat lower in the
new simulator than in the previous version, although
considerable effort was spent on matching that data type
(because of the importance of reproducing the proper level of
aquifer support on predicting future recovery). Although this
drop in average Match Factor is not a large enough drop to
invalidate the new simulation, it indicates the need to review
the methods used to build the reservoir description before this
simulator is updated again.
Although the average RFT Match Factor is very consistent
for both versions of the simulation, it has the lowest value of
the five data types. Part of the reason for this is that the RFT
pressure measurement is one of the most accurate and least
ambiguous surveillance measurements made. The depth is
defined by wireline techniques to within a few feet. The
reservoir pressure is typically measured by a very accurate
SPE 75223
SPE 75223
Conclusions
In this paper, we have proposed using three simple statistical
methods for evaluating the degree to which a simulation is
history matched to past performance, as defined by the
surveillance data collected in the field. The goal is to quickly
extract simulator results and generate graphical displays which
can be used to 1) guide subsequent history matching iterations,
2) determine when further history matching no longer adds
value, and 3) compare how closely two or more simulations
match the surveillance data.
All three methods rely on extracting specific results from a
simulation run and calculating the deviation of that value from
the corresponding surveillance observation value. The
Deviation Distribution Plot shows whether the simulator is
over or under-predicting surveillance data of a particular type,
and the general degree of that match. The Deviation Band Plot
also demonstrates the degree of match and provides a
quantitative measurement of the match: the apparent standard
deviation (ASD). In Match Factor analysis, each deviation is
converted to a relative measure of the confidence that the
simulator actually reproduced the particular reservoir
performance at the time the surveillance measurement was
made. Arithmetically averaging Match Factors by well
allows bubble maps to be created that quickly show which
areas are better matched. Bar charts of Match Factors
averaged by data type summarize the entire history match
evaluation for one or more simulation runs. Weighting factors
are used to emphasize the impact of key Match Factor values
on the average.
Because of the volume of data involved in history
matching simulators with large surveillance databases, a
simulator-specific computer program is required in order to
extract the required data from the simulators output files.
However, once the extraction routine is constructed, the
analyses proposed here can be carried out in a matter of
minutes following the completion of a simulation run. That
makes it well suited to the iterative nature of history matching.
In addition, the ability to evaluate the match well-by-well,
focuses the process on those areas which need it the most.
Acknowledgements
The authors wish to thank the management of Kuwait Oil
Company for their support during the development and
implementation of the techniques described in this paper.
Specifically, our thanks go to Mr. Ali Al-Shammari, the
General Superintendent of the West Kuwait group in the
Fields Development Department of the Exploration And
Development Division in the Kuwait Oil Company (KOC).
Also, our thanks go to Mr. Jaber Al-Hunaif and Mr. Fahad AlMehdadi, the Superintendents responsible for the original
development work and implementation.
In addition, the authors also wish to thank the management
of BP Kuwait and Chevron Kuwait for their support. Finally,
from among the numerous colleagues in KOC, BP Kuwait and
Chevron Kuwait who contributed their input and suggestions,
we wish to draw special attention to the support provided by
Dr. Robert Merrill and Alan Clark, both KOC secondees from
References
1.
2.
3.
4.
5.
10
SPE 75223
SPE 75223
11
12
SPE 75223
Observed
Values
Predicted
Values
Deviations
Time
25
20
15
10
5
0
-300
-200
-100
100
200
300
Cumulative Frequency
(%)
SP3E 7522
100
90
Case C
80
70
60
50
40
30
20
10
0
-800 -600 -400 -200
13
Case B
Case A
Observed
Values
Predicted
Values
Deviation
Band Bars
TIME
Figure 4: Example History Match Plot Showing Deviation Band
Cumulative Frequency
(%)
14
100
90
80
70
60
50
40
30
20
10
0
SPE 75223
Case A
Case D
0
100
200
300
400
500
600
700
Match Factor
100
90
80
Maximum
Value
70
60
50
40
30
20
10
0
-0.3
Zero
Confidence
Limits
-0.2
-0.1
0.0
0.1
0.2
0.3
15
Observed
Values
SP3E 7522
Predicted
Values
Match
Factor
Bars
TIME
Cumulative Frequency
(%)
100
90
80
70
60
50
40
30
20
10
0
-600 -400 -200
Well 01
Well 02
Well 03
Well 04
Well 05
Well 06
Well 07
Well 08
Well 09
Well 10
Well 11
Well 12
Cumulative Frequency
(%)
16
100
90
80
70
60
50
40
30
20
10
0
SPE 75223
40% Of Data At
Zero Deviation
0.0
0.2
0.4
0.6
0.8
1.0
Well Location
Problem Well
MF 95-100
Well
Matched
Area
MF 90-95
MF 80-90
MF 60-80
MF 30-60
Poorly
Matched
Area
MF 0-30
Figure 10: Example Bubble Map Showing Average Watercut Match Factors By Well
SP3E 7522
17
Match Factor
100
90
80
70
60
50
40
30
20
10
0
-4
-3
-2
-1
Figure 11: Match Factor Curve For Total Liquid Production Rate Example
100
90
80
70
60
50
40
30
20
10
0
75% Of Wells
At Maximum
Match Factor
5% Of Wells
At Very Low
Match Factor
10
20
30
40
50
60
70
80
90 100
18
SPE 75223
98 99
100
83 85
75
80
81
66
60
53
52 51
40
Previous
New
20
0
TLPR
WC
GOR
SBHP
RFT
Data Type
Figure 13: Average Match Factor Summary Bar Chart
Completely
Interval
Completely
Unswept:
Unswept:
Swmax
Sorw
TDT
SwSw
Swi
Match Factor
Match
Factor
Model Layer
100
Swi
Block Sw
Sorw
Swmax
SP3E 7522
Completely
Interval
Completely
Swept: Swept:
Swmax
Sorw
TDT
SwSw
Swi
Model Layer
Match Factor
Match
Factor
100
Swi
Block Sw
Sorw
Swmax
Contact
In Layer:
Water
Contact
Within Interval :
Swmax
Sorw
TDT
SwSw
Average Sw
From Log
Swi
New OWC
Match Factor
Match
Factor
Model Layer
100
Swi
Block Sw
Sorw
Swmax
19
20
SPE 75223
TransitionWater
Complex
Zone:Saturation Distribution:
Swmax
Sorw
TDT
SwSw
Average Sw
From Log
Swi
Match Factor
Match
Factor
Model Layer
100
Transition
Zone
New OWC
Swi
Block Sw
Sorw
Swmax
T1
T2
T3
T4
T5
T6
5
4
1
PreDrill
Period
6
Shut-In
Period
Producing Time 1
Producing Time 2
Figure B 1: Weighting Factor Example #1 For Data Types Defined Only When Wells Active
Time
SP3E 7522
T1
T2
T3
T4
T5
2
1
T6
4
PreDrill
Period
21
5
Shut-In
Period
Time
Figure B 2: Weighting Factor Example #2 For Data Types Defined Throughout History Match Period