Transcript Research Poster 36 x 48 - A

NS23A-1637
STATIC AND DYNAMIC RESERVOIR CHARACTERIZATION USING
HIGH RESOLUTION P-WAVE SEISMIC VELOCITY DATA IN DELHI FIELD, LA
Authors’ Emails:
[email protected]
[email protected]
Sidra (Shahid) Hussain and Dr. Thomas Davis
Reservoir Characterization Project, Geophysics Department, Colorado School of Mines
Abstract
2. Static Characterization
Static and dynamic reservoir characterization was done on high resolution P-wave
seismic data in Delhi Field, LA to study the complex stratigraphy of the Holt-Bryant
sands and to delineate the CO2 flow path. The interpretation was done on bandwidthextended seismic data. Acoustic impedance inversion done on monitor and base
surveys helped in delineating CO2 flowpaths and the channel geometry in the
reservoir.
To monitor the flow of CO2 within Paluxy and Tuscaloosa sandstone formations, dynamic characterization was done using merged
2008-10 survey as the base survey and RCP (June 2010) survey as the monitor survey. The surveys were cross-equalized to
increase the repeatability. Model-based acoustic impedance (AI) inversion was performed on both the surveys to quantify the
changes in acoustic impedance with the addition of CO2 and water in the field. Fluid substitution modeling was done on well data
to model the expected changes in AI with CO2 and water.
Overlap
RCP
June 2010
Mar
2010
Acknowledgements
3.1 Cross-Equalization
2008
1. Introduction to Delhi Field
Permanent
Patch
(Jan 2010)
Survey Re-gridding
1 mile
Location of the seismic surveys in the field. The merged survey of 2008 and March
2010 was used as the base survey for 3D interpretation and the RCP survey from
June 2010 was used as the monitor survey for time-lapse interpretation in this study.
(Sun Oil Co. and Denbury Resources Inc.)
Static Time Shift
Time-Variant Time
Shift
Bandwidth extension was applied by Geotrace on the 3D seismic datasets using the
method of continuous wavelet transform to recover the lower and the higher
frequencies in the seismic bandwidth that were lost from earth’s reflectivity during
transmission. The method increased the dominant frequency in the data and
decreased the tuning effect.
After
Bandwidth
Extension
T= 20ms
T= 12ms
F=1/T=83.3 Hz
F=1/T=50Hz
Areal extent is 6200 acres
15 miles long
Delhi OOIP 357 MMBbls
RCP Study Area 7 sq miles
The steps taken to cross-equalize the base and the monitor surveys
2
2
0
Before Cross-equalization
After Cross-equalization
Positive Seismic Amplitude
NRMS maps and their respective histograms within the reservoir zone. The repeatability has improved with cross-equalization.
However, some zones of low repeatability around the wells are observed where one expects to see changes with production and
injection. The middle of the reservoir shows high repeatability. This could mean that the middle of the reservoir is being bypassed.
3.2 Model-based Acoustic
Impedance Inversion
Statistical constant phase wavelets extracted from the merged (2008-2010) survey
within 800-1100 ms, Inline=1019 and Xline=1001-1379
0
Positive Seismic Amplitude
3.3 Fluid Substitution Modeling
0
AI % change
-5
GR log
-10
5% CO2
-15
10% CO2
-20
-25
1400
Post-BE data
traces
overlain on
pre-BE data
After BE
Before BE
QC for Phase-shift: For most peaks and troughs, there is no phase shift.
Hierarchy of a model-based acoustic impedance
inversion (Modified from Young, 2006).
Location of Delhi Field in relation to the surrounding structural features
(Modified after Mancini, et al. 2008a).
MCF/D
BPD
100,000
BPD
QC for Time Shift: The time
difference falls within the range of -1
and +1 (see histogram) which is
normal for a trace by trace
computational continuous wavelet
transform process.
20,000
WATER
INJ
100,000Prior
1970
Production data not
available
Prior
1970
Primary
Production
RF ~14%
Secondary
Production
RF ~ 40%
Production
data not
available
Primary
Productio
n RF
~14%
OIL
PROD
Secondary
Production
RF ~ 40%
50
0 45
65
60
55
WATE
WATER
R INJ
PROD
OIL
Time in years
PROD
55
60
65
0
N
Injector Wells
Producer Wells
1 mile
2005
2400
Percentage change in acoustic
impedance with the addition of
CO2 under different pressure
conditions.
10% CO2
1600
1800
2000
Pressure (psi)
2200
2400
Percentage change in acoustic
impedance with the increase in
effective pressure for CO2
replacing brine.
4. Results
35
1 mile
GAS
PROD
1800
2000
2200
Pressure (psi)
5% CO2
N
GAS WATER
PROD PROD
Production history data of Delhi Field before the tertiary recovery started
0
(Modified from
Time in years
2005
45 Denbury Resources Inc.).
50
QC of well-to-seismic ties: Synthetic
seismograms (red) created using the
wavelets shown below the ties; The
correlation coefficient of pre-BE data is
94% while that of the post-BE data is
73% which is reasonable for a high
frequency data.
1600
8
7
6
5
4
3
2
1
0
1400
2.2 Structural Interpretation
MCF/
D20,00
0
AI% change with effective pressure for CO2
replacing brine
AI % change at different effective pressures
for CO2 replacing brine
2.1.1 QC of Bandwidth Extension
A generalized stratigraphic section of the Cretaceous formations in Delhi Field
(Denbury Resources Inc.).
Cross-Normalization
2.1 Bandwidth Extension
Before
Bandwidth
Extension
Location of Delhi Field (Evolution Petroleum Corporation).
Spectral Shaping
AI % change
I would like to thank Dr. Tom Davis and my thesis committee for helping me
with this research, Denbury Resources Inc. for sponsoring the research and
providing us with the data, Geotrace for doing the bandwidth extension and
Colorado School of Mines, Geophysics Department’s Reservoir
Characterization Project (RCP) for funding my research.
•
•
•
•
3. Dynamic Characterization
-25
Injector Wells
Producer Wells
Amplitude map of Paluxy from the
Amplitude map of the top of TUSC 7 from
the merged BE survey showing bright sandstonemerged BE survey showing SW-NE
trending meandering channel-like
bodies in red.
features in hot colors.
-4
An arbitrary dip line showing amplitude
difference through four of the wells in the
phase-1 injection pattern. The effect of CO2
can be seen around the injectors.
AI % difference map with a 2 ms window centered at
with the production data from June 2010 overlain; the
black dashed line is the oil-water contact; notice the
negative impedance change below the OWC. The
lighter yellow color shows area where Paluxy is not
being swept completely.
AI % difference map at TUSC 7 top; the white
triangles are TUSC 7 injectors and the white
circles are TUSC 7 producers; the black
polygons show the flow paths of CO2
illuminating channel-like features in the
sandstone.