...

Advanced Subsurface Characterization for CO Geologic Sequestration and Induced Seismicity Evaluations

by user

on
Category: Documents
25

views

Report

Comments

Transcript

Advanced Subsurface Characterization for CO Geologic Sequestration and Induced Seismicity Evaluations
Advanced Subsurface Characterization for CO2 Geologic
Sequestration and Induced Seismicity Evaluations
Tiraz Birdie, Lynn Watney, Aimee Scheffer, Jason Rush, Eugene Holubnyak, Mina Fazelalavi,
John Doveton, Jennifer Raney, Saugata Datta, Dennis Hedke, and Jennifer Roberts
Carbon Management Technology Conference
November 19th 2015
Sugarland, TX
Sedimentary Basins in the US
Suitability of Sedimentary Basins for Geologic Sequestration
Evolution of Sedimentary Basins
Carbonates
Evaporites
Shales
(Caprock)
Present Day Evolution of Evaporative Cap - Arabian Gulf
Analog to portions of the Arbuckle CO2 injection zone at Wellington Field, KS
Tidal Effects Reflected in Geologic Logs
Time
GR/CGR/
SP/Cal
Microresistivity
Neutron-Den-Pe
Sonic
Impedance
Reflection
Coefficient
Synthetic
100 Hz
Top Cherokee Gp.
Top Mississippian
Secondary caprock
Pay
CO2-EOR pilot
Depth
Equiv.
Pierson Fm.
~ 50 million years
Caprock
Carbonate
Evaporite
ØN
Primary caprock
Interval
Chattanooga Sh.
Simpson Group
Top Arbuckle
Jefferson CityCotter
Baffle/barrier
-Tight, dense
- High
impedance
Roubidoux Fm.
Gasconade Dol.
Gunter Ss.
CO2 Injection
zone
Importance of Characterization on CO2 Plume and Pressure
Projections
Coarse Vertical Characterization
Meter-scale Vertical Characterization
Co2 Plume
Induced Pressure
Co2 Plume
Induced Pressure
Wellington CO2Sequestration and EOR Site
26,000 tons
Mississippian
26,000 tons
(Saline Aquifer)
Nuclear Magnetic Resonance Image Logging
Effect of Pore Space on Time to Return to Equilibrium
NMR Variables and Their Influence on
Petropohphsical Properties
T1 Time: Time to align the protons with the magnetic field
T2 Time: Time for protons to recover, through bulk, diffusion, and surface relaxivity
Mississippian MRI
Caprock MRI
#1-28, lower Miss to top Arbuckle
Magnetic resonance imaging analysis
low high
small large
T2
(pore size)
Permeability
Lower Mississippian
argillaceous dolosiltone,
small pores
50 ft
high
low
Porosity
Chattanooga Shale
Smallest pores
Simpson shales,
Smallest pores
Top Arbuckle
Caprock evidence of lower Miss. :
• Micro-nano darcy perm
• Quiet fracture wise
• Organic matter ~2% TOC
Flow Units in the Lower Arbuckle Injection Zone
KGS #1-32
0
2
Porosity%
4
6
8
10
12
Porosity%
14
4900
0
50 ft
Connected vugs
 Solution &
fracture
Nonconnected
vugs
5020
5030
5040
5050
5070
5080
5100
5110
5130
5010
5090
5100
5120
5110
Interparticle/matrix
5120
5140
5130
5150
5140
5160
Flow unit boundaries
5060
5070
5090
4950
Wellington #1-28
5060
5080
4940
5000
Wellington #1-32
5050
14
4990
5000
5040
12
4980
4990
5030
10
4970
4980
5020
8
4930
4960
4970
5010
6
4920
4930
4960
4
4910
4920
4950
2
Ø
4900
4910
4940
Wells 3500 ft apart
KGS #1-28
Ø
Doveton and Fazelalavi, 2012
Utilize whole core
analysis, NMR,
spectral sonic, and
resistivity logs
Porosity and Hydraulic Conductivity Comparsion with
Core Data
Good match with core data
MRI and PHND Estimates of Porosity
MRI effective porosity
and neutron-density
crossplot porosity in
the Arbuckle of
Wellington #1-32
Conclusion: there is a good match
between MRI porosity and lithologycorrected neutron-density porosity
which is a useful cross-validation of
these logs
Permeability Profile of Arbuckle
(NMR)
Redox reactive ions
reflect changes in
biogeochemistry
(microbial) occurring
between upper and
lower Arbuckle, in
turn attributed to lack
of hydraulic
communication
Scheffer, KGS
Upper Arbuckle
(High Permeability)
Mid Arbuckle
(Low Permeability)
Lower Arbuckle
(High Permeability)
Ion Based Verification of Baffle Zone and Caprock
Scheffer, KGS
• Brine of Lower Arbuckle vary substantially from Upper Arbuckle
• Lower Arbuckle brines cluster together
• Upper Arbuckle values more spaced out, suggests smaller baffles
Isotopic Verification of Baffle Zone and Caprock
Oxygen & Hydrogen
isotopes of brines from
DST and perf & swabbing
Mississippian
Brines
(distinct from
Arbuckle)
Upper Arbuckle
-- distinct
Lower Arbuckle
injection interval
-Waters distinct from upper Arbuckle and Miss
- Lower intervals are also geochemically
homogeneous  infer fracture connectivity
Scheffer, 2012
Microbial Ecology and Validation of Baffle Zone
Baffle top
•
Lowest biomass coincides with low perm zone (Lower JCC) and low DOC
•
Highest biomass coincides with high perm and high concentrations of sulfate
•
Same 9 genera were found in brine from Upper Arbuckle depths
•
Brine from tight zone had 7 genera; 3 less and 1 unique
•
Supports mixing of Upper Arbuckle and some degree of separation below
Seismic Profile Confirms Permeability Stratification in Arbuckle
South
Impedance = ρ x Ø
East
KGS #1-32
KGS #1-28
Top Oread
Thick
Lansing Group
Shales
Top Kansas City Ls.
Top Mississippian
Lower Pierson
Top Arbuckle
Baffle or potential barrier to vertical flow
(high impedance)
Low impedance injection interval
Bot Arbuckle
Hedke, 2012
Seismic Structure Mapping Confirms Regional
Presence of Caprock
Entry Pressure Analysis
Rhomaa-Umma Analysis
Depth-Constrained Clustering
using Potassium, Uranium, Thorium
condensed
condensed
condensed
condensed
Oil show
condensed
Significant flooding surface
from core description
K- RhomaaU- Umma
Th Lith Ø
Rt
Injection zone
Simulated Plume and Pressure Projections
Holubnyack, Rush, Fazelalavi (KGS)
•
Baffle Zones keep CO2 plume and pressures
confined deep in the Arbuckle injection zone
Earthquake Trends in Southern Kansas
Disposed Volume
140
0
120
20
100
40
80
60
60
80
40
Number of Quakes (>M2.0)
Disposal Volumes (MMBL)
Annual Earthquake Count
100
20
0
120
1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 2012 2014
Induced Seismicity - Physical Mechanisms
Subsurface stress field
Faults slips when principal
stresses exceed threshold
Drilling Induced Fractures and Leak-off Test
Used to Estimate Principal Stresses
XRMI Log
Drilling Induced
Fractures
𝑆𝐻𝑚𝑎𝑥 ~ 3𝑆𝐻𝑚𝑖𝑛 - 2𝑃𝑝 + 0.1(𝑆ℎ𝑚𝑖𝑛 - 𝑃𝑝 )
Fault Identification
• Extensive data acquisition required to
identify faults and assess seismic risk
regionally
3D Stress Analysis Used to Estimate Fault Slip Tendency
Slip Tendency Plot
𝝈𝒏
ST = 𝝉/ 𝝈𝒏
•
ST= 0.3 (lower than of 0.5 for
fault slippage)
𝝉
Fly UP