Opus Terra
Optimization & Uncertainty Solutions
Terra 3E SAS
Dominique Gurillot
Jrmie Bruyelle
Outline
Opus Terra toolbox
Example of Petrel* workflows
History Matching
Optimization
Uncertainties
PUNQ-S3
Presentation
Geological description
Dynamic data
Geological modeling
History Matching
Prediction
Conclusion
* Mark of Schlumberger
Opus Terra : Optimization & Uncertainty Solutions
Opus Terra Toolbox
Toolbox contains plug-ins for Petrel*
Glhis : Global History Matching (CMA-ES)
CMA-ES has been recognized as one of the most powerful continuous
optimization algorithms on benchmark problems (Hansen et al., 2010)
and real-world problems
Sirenn : Simulator Reservoir Neural Network
Neural networks have been developed to reproduce complex physical
phenomena
Neural networks are very well adapted to represent nonlinear
phenomena
These tools are fully integrated in Petrel*
* Mark of Schlumberger
Opus Terra : Optimization & Uncertainty Solutions
Example of Petrel Workflows History Matching
Opus Terra : Optimization & Uncertainty Solutions
Example of Petrel Workflows Optimization
Opus Terra : Optimization & Uncertainty Solutions
Example of Petrel Workflows Uncertainties
Opus Terra : Optimization & Uncertainty Solutions
PUNQ-S3
The PUNQ-S3 case has been
taken from a reservoir
engineering study on a real
field performed Elf
Exploration Production.
It was qualified as a smallsize industrial reservoir
engineering model.
The geological description is
based on knowledge of the
regional geology.
Opus Terra : Optimization & Uncertainty Solutions
PUNQ-S3 Geological Description
Layers 1, 3, and 5 have linear streaks of high-porous sands (phi > 20 %),
with an azimuth somewhere between 110 and 170 degrees SE. These sand
streaks of about 800 m wide are embedded in a low porous shale matrix
(phi < 5 %).
In layer 2 marine or lagoonal shales occur , in which distal mouthbar or
distal lagoonal delta occur. They translate into a low-porous (phi < 5%),
shaly sediment, with some irregular patches of somewhat higher porosity
(phi > 5%).
Layer 4 contains mouthbars or lagoonal deltas within lagoonal clays, so a
flow unit is expected which consists of an intermediate porosity region
(phi ~ 15%) with an approximate lobate shape embedded in a low-porosity
matrix (phi < 5%). The lobate shape is usually expressed as an ellipse (ratio
of the axes= 3:2) with the longest axis perpendicular to the paleocurrent
(which is between 110 and 170 degrees SE).
Opus Terra : Optimization & Uncertainty Solutions
PUNQ-S3 Geological Description
Expected sedimentary facies with estimates for width and
spacing for major flow units for each layer
Opus Terra : Optimization & Uncertainty Solutions
PUNQ-S3 Geological Modeling
The model contains 19x28x5 grid blocks, of which 1761 blocks
are active.
Layer 1, 3 and 5 : As describe in the geological description, we
consider two facies:
An high-porous sands (phi > 20 %);
A low porous shale matrix (phi < 5 %).
The geological modeling was performed by using an adaptive channel
modeling using the geological description and the hard observed data.
Opus Terra : Optimization & Uncertainty Solutions
10
PUNQ-S3 Geological Modeling
Layer 2 : we consider two facies:
A low porous shaly sediment (phi < 5%);
A high porous shaly sediment (phi > 5 %).
The geological modeling was performed by using ellipses as body
shape modeling using the geological description and the hard
observed data.
Layer 4 : we consider two facies:
An intermediate porosity region (phi ~ 15%);
A low-porosity matrix (phi < 5%).
The geological modeling was performed by using ellipses as body
shape using the geological description and the hard observed data.
Opus Terra : Optimization & Uncertainty Solutions
11
PUNQ-S3 Geological Modeling
The uncertain geological parameters of PUNQ-S3 are the
porosities, the vertical and horizontal permeabilities.
The parameterization of PUNQ-S3 model is based on the
geological description. We consider each facies, as describe in
the previous section, and we estimates the constante
properties for each facies: 18 parameters.
Opus Terra : Optimization & Uncertainty Solutions
12
PUNQ-S3 History Matching
The production scheduling has been inspired by the original
model
A first year of extended well testing, followed by a three year shut-in
period, before field production commences. The well testing year
consists of four three-monthly production periods, each having its own
production rate.
During field production, two weeks of each year are used for each well
to do a shut-in test to collect shut-in pressure data. The wells operate
under production constraint. After falling below a limiting bottom hole
pressure, they will switch to BHP-constraint.
13
PUNQ-S3 History Matching
Available data points used in history matching
Date
PRO-1
BHP
GOR
PRO-4
WC
BHP
GOR
PRO-5
WC
BHP
GOR
PRO-11
WC
BHP
GOR
PRO-12
WC
BHP
GOR
PRO-15
WC
BHP
GOR
1.01
224.0
225.2
228.7
219.3
231.0
217.1
91
211.7
210.6
222.9
202.7
218.7
193.5
182
215.6
216.6
223.4
208.4
220.8
209.0
274
219.6
224.4
230.0
219.1
224.7
216.8
366
226.3
229.5
230.7
228.4
229.8
228.1
1461 233.2
234.2
235.9
235.3
234.6
234.6
1826 201.0 135
190.4
215.7 63.6
203.8 67.2
209.1
1840 222.2
224.5
226.5
225.1
225.4
223.4
207.2 62.5
194.4 59.1
200.5 74.5
181.3 56.3
220.0
218.8
219.6
217.2
197.2 67.8
175.9 58.1
215.5
212.2
1642
72.6
1841
67
191.1 63.7
82.1
2008
191.7
165.3
2192 190.6
177.7
2206 215.4
2373
217.9
147.1
106
2557 184.4
170.7
202.6 62.5
186.0 62.4
2571 210.8
213.4
215.7
213.3
2572
2738
190.1
2922 178.5
74.8
0
2936 206.1
Sigma
WC
167.8
0.022
0
210.3
34
0.2
196.4 59.9 0.002 169.9 65.2 0.098 195.7 76.4
212.4
21
0.2
Opus Terra : Optimization & Uncertainty Solutions
208.9
6.2
0.2
212.4
6.3
0.2
170.0 49.8
208.1
0.2
5.7
0.2
14
PUNQ-S3 History Matching
Opus Terra : Optimization & Uncertainty Solutions
15
PUNQ-S3 History Matching
Sensitivity analysis by variable
Equal spacing sampler: 4 simulations by variable = 72 simulations.
Selected
Parameters
Opus Terra : Optimization & Uncertainty Solutions
16
PUNQ-S3 History Matching
Proxy model of the objective function with Sirenn
Training data
Experimental design : Fractionnal factorial sampler : 32 simulations +
central point.
Simulation performed for sensitivity analysis : 72 simulations
Minimization of the objective function with Glhis using the
Sirenn proxy
Opus Terra : Optimization & Uncertainty Solutions
17
PUNQ-S3 History Matching Results
18
PUNQ-S3 History Matching Results
19
PUNQ-S3 History Matching Results
20
PUNQ-S3 History Matching Results
21
PUNQ-S3 History Matching Results
22
PUNQ-S3 History Matching Results
23
PUNQ-S3 - Prediction
The next step consists to predict the ultimate recovery after
16.5 years.
Prediction obtains with 10 different solutions.
History Matching
Prediction
24
PUNQ-S3 - Conclusion
Opus Terra allows to:
Build a proxy of the objective function with a minimal number of
simulations
Perform a global optimization
The different solutions fit the observed data (Pressure, Watercut and Gas-oil ratio)
The predictions of differents solutions are very close to the
truth case.
Opus Terra : Optimization & Uncertainty Solutions
25
Useful Links
Opus Terra
Ocean Store: http://www.ocean.slb.com/Pages/Product.aspx?category=all&cat=Ocean&id=POTA-B1
Leaflet: http://www.terra3e.com/Docs/OpusTerra_Leaflet.pdf
Tutorial: http://terra3e.com/Docs/OpusTerra.avi
PUNQ-S3:
Imperial College: http://www3.imperial.ac.uk/earthscienceandengineering/research/perm/punq-s3model
CMA-ES :
Wikipedia: http://en.wikipedia.org/wiki/CMA-ES
Other products:
VolTerra: http://www.ocean.slb.com/Pages/Product.aspx?category=all&cat=Ocean&id=PVTE-B1
Scenarium: http://www.ocean.slb.com/Pages/Product.aspx?category=all&cat=Ocean&id=PSCN-B1
Sirenn: http://www.ocean.slb.com/Pages/Product.aspx?category=all&cat=Ocean&id=PSRN-B1
Glhis: http://www.ocean.slb.com/Pages/Product.aspx?category=all&cat=Ocean&id=PGLH-B1
Opus Terra : Optimization & Uncertainty Solutions
26
www.Terra3E.com
Dominique Gurillot : [email protected]
Jrmie Bruyelle : [email protected]
