water-saturation estimation from seismic and rock-property
TRANSCRIPT
![Page 1: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/1.jpg)
Water-saturation estimation from seismic and rock-property trends
Zhengyun Zhou Fred J. HiltermanHaitao Ren Mritunjay Kumar
Center for Applied Geosciences and EnergyUniversity of Houston
![Page 2: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/2.jpg)
Objective
The objective of this research is to differentiate Fizz saturation from Gassaturation based on seismic and rock property trends.
“Fizz” ⇒ Low gas saturation“Gas” ⇒ Economic gas saturation
![Page 3: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/3.jpg)
MotivationWill you drill the well at prospect location?
Gas well
Dry well
?
Discovery Prospect
Amplitude map
![Page 4: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/4.jpg)
Outline
• Seismic Field Calibration• Rock-Property Transforms• Application of Transforms• Conclusions
![Page 5: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/5.jpg)
Seismic Field Calibration
GasWet
Well B
Well B is a known gas reservoir. At down-dip equivalent, we assume reservoir is wet.
Fairfield data
![Page 6: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/6.jpg)
Well B – Resistivity and Sonic in Pay Zone
densityTotal Pay = 65 ft Gas Sand Interval
sonic = 100µsec10,000 ft/sec
![Page 7: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/7.jpg)
Well B – Rock Properties
Gas Sand 65 ft 10000 ft/sec
2.00 gm/cc
Shale 9090 ft/sec2.35 gm/cc
Shale 9090 ft/sec2.35 gm/cc
Depth10000 ft
Model created from well-log curves
![Page 8: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/8.jpg)
Thin-Bed Synthetic Match With Migrated CDP Data
0 5 10 15 20 25 30 35Incident angle (degrees)
0 5 10 15 20 25 30 35Incident angle (degrees)
Does the prospect have the same AVO response?
Field CDP data
Thin-Bed Synthetics
Fairfield data
![Page 9: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/9.jpg)
Discovery versus Prospect AVO SignaturesDiscovery Prospect
Prospect has same AVO response as Discovery.
![Page 10: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/10.jpg)
Drilling Results: Hard Shale Over ProspectDiscovery Prospect
Fizz Fizz
Fizz and gas are differentiated by down-dip wet response.
Can we quantify seismic AVO to water saturation ?
![Page 11: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/11.jpg)
Outline
• Seismic Field Calibration• Rock-Property Transforms• Application of Transforms• Conclusions
![Page 12: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/12.jpg)
Assumption
For one particular hydrocarbon reservoir, the rock matrix is assumed to be the same as its down-dip equivalent. But the prospect and down-dip equivalent can have different thicknesses.
![Page 13: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/13.jpg)
Local Rock Properties – GOM - Louisiana
ShipShoal
Main Pass
S. Pass
New OrleansField A
600 km
80 wells Above onset of abnormal pressure
Locally, what are the velocity and density variations for sand and shale?
GDC data
![Page 14: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/14.jpg)
Sand Velocity (ft/sec)
Average VelocityAfter rejecting outliers
Standard Deviation
GDC data
Rock-Property Variations
![Page 15: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/15.jpg)
Rock-Properties: ± 1 Standard Deviation
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Shaleα+1std, ρ+1std-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Shaleα-1std, ρ-1std
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Base Model
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Shale and Sandα+1std, ρ+1std
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Shale and Sandα-1std, ρ-1std
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Base Model
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Sandα-1std, ρ-1std
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Base Model
WetFizzGas
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Ampl
itude
Sandα+1std, ρ+1std
![Page 16: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/16.jpg)
Rock-Property Transforms
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Shale Properties: +1 Std. Deviation
-0.2
-0.1
0
0.1
0 10 20 30 40Incident Angle
Am
plitu
de
Shale Properties: -1 Std. Deviation
1. NI(wet) – NI(gas) ≈ Constant
2. Slope is proportional to NIMore positive NI ⇒ Larger SlopeMore negative NI ⇒ Smaller slope
Let’s quantify these two observations
Two Observations:
![Page 17: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/17.jpg)
ShipShoal
Main Pass
S. Pass
New OrleansField A
Miocene Rock Properties – GOM - Louisiana
600 km
50 shelf wells, 239 Miocene sand packages• 149 hydrocarbon charged• 90 brine filled• 1600-6500m depths
• Measure sand and encasing shale rock properties• Fluid substitution:
Wet-, gas-, and fizz-saturated rock properties GDC data
![Page 18: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/18.jpg)
NI for 239 Miocene sand packages
-0.3
-0.2
-0.1
0.0
0.1
0.2
-0.1 0 0.1 0.2NI with Brine
NI w
ith G
as a
nd B
rine
07 0
Wet
Wet
NI(Gas) = -.08 + 1.25 NI(Wet)R2 = 0.8 GDC data
![Page 19: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/19.jpg)
Quantifying Local Reflectivity TransformsObservation 1
-0.08
-0.06
-0.04
-0.02
0
0.02
0.04
0.06
0 0.02 0.04 0.06 0.08 0.1NI (Wet)
NI (
Hyd
roca
rbon
s)
NIGAS=-.08+1.12 NIWET
NIFIZZ = -.05 +1.06 NIWET
Pore Fluid Transforms
Observation 2
Slope Transforms
y = 0.4219x - 0.0397
y = 0.3834x - 0.0936
y = 0.3847x - 0.1189
-0.16
-0.14
-0.12
-0.1
-0.08
-0.06
-0.04
-0.02
0
-0.1 -0.05 0 0.05 0.1Near
Far
RC(30°)WET=-.04+.42 RC(0°)WET
RC(30°)FIZZ = -.09 + .38 RC(0°)FIZZ
RC(30°)GAS = -.12 + .38 RC(0°)GAS
How are these transforms applied?
![Page 20: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/20.jpg)
Outline
• Seismic Field Calibration• Rock-Property Transforms• Application of Transforms• Conclusions
![Page 21: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/21.jpg)
Seismic Amplitude Model
θ
θTRb
A(θ) ≈ k RCZP_TOP(θ)λ
cos(θTR)4πb
Seismic amplitude
Reflection coefficient
Wavelengthb< λ/4
![Page 22: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/22.jpg)
Field Measurements – Amplitude Notation
A(0°)WET A(30°)WET A(0°)ANOMALYA(30°)ANOMALY
Prospect AVODown Dip AVO
How are seismic amplitudes converted into Reflection Coefficients?
![Page 23: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/23.jpg)
Transforming Seismic Amplitude into NISlope Transform
RC(30°) = b0 + b1 NI + Thin-bed amplitude response
NI =A(0°) * b0
A(30°)/cos(30°) – b1 A(0°)
RC(30°)WET = -.04 + .42 NIWETRC(30°)FIZZ = -.09 + .38 NIFIZZRC(30°)GAS = -.12 + .38 NIGAS
Rock propertymeasurements of b0 and b1
=NIGAS1.15 A(30°)ANOMALY – 0.38 A(0°)ANOMALY
-0.12 A(0°)ANOMALY
NI Test Statistic for Gas
![Page 24: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/24.jpg)
WET
HC HC
Near-Offset
WET
Far-Offset
HC HC
NI (Wet) = -.04 Near/ (1.15 Far - .42 Near)
NI of Wet
NI (Gas) = -.08 + 1.12 NI(Wet)
Standard for NI of gas
NI of gas
NI (Fizz) = -.05 + 1.06 NI (Wet)
Standard for NI of fizz
NI of Fizz
15.0
-131.2
-180.0
-82.5
-33.75
Amplitude
Estimation in Wet Area
Amplitude Maps
NI of Wet
Slope transformfor wet
NI of Gas and Fizz
Pore fluidtransforms:Wet GasWet Fizz
NI
![Page 25: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/25.jpg)
Test NI of Gas and Fizz
Slope transformfor gas and fizz
WET
HC HC
Near-Offset
WET
Far-Offset
HC HC
NI (gas) = f ( Near, Far )
HC HC
Gas test of NI
NI (fizz) = f ( Near, Far )
HC HC
Fizz test of NI
NI of Gas and Fizz
Compare tests with
standards
Standard for NI of gas
NI of gas
Standard for NI of fizz
NI of Fizz
15.0
-131.2
-180.0
-82.5
-33.75
Amplitude
NI
Estimation in Hydrocarbon Area
Amplitude Maps
![Page 26: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/26.jpg)
Test NI of Gas and Fizz
Slope transformfor gas and fizz
WET
HC HC
Near-Offset
WET
Far-Offset
HC HC
NI (gas) = f ( Near, Far )
Gas Fizz
Gas test of NI
NI (fizz) = f ( Near, Far )
Gas Fizz
Fizz test of NI
NI of Gas and Fizz
Compare tests with
standards
Standard for NI of gas
NI of gas
Standard for NI of fizz
NI of Fizz
15.0
-131.2
-180.0
-82.5
-33.75
Amplitude
NI
Estimation in Hydrocarbon Area
Amplitude MapsHC
![Page 27: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/27.jpg)
Outline
• Seismic Field Calibration• Rock-Property Transforms• Application of Transforms• Conclusions
![Page 28: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/28.jpg)
Conclusions• Fizz and gas can have the same AVO responses, but
down-dip water-saturated AVO response discriminates fizz from gas reservoirs.
• NI of wet, gas and fizz saturations vary with changing rock properties. However, the values of (NIWET –NIGAS) and of (NIWET –NIFIZZ) remain fairly stable. In the Pore-Fluid Transforms, linear relationships are used to predict NIGAS and NIFIZZ from NIWET .
• Near and far amplitude maps combined with Slope Transforms estimate the Reflection Coefficients for various pore fluids. Water saturation can be determined by comparing the NI values predicted in wet area and in prospect area.
Needs field verification !
![Page 29: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/29.jpg)
AcknowledgementsRQL sponsors
Permission to show data:Geophysical Development Corporation – PetrophysicalFairfield - Seismic
PetroSeismic Inc – Development of TIPS software
Rocky Roden – Assistance in preparing seismic data throughSMT AVOPAK.
![Page 30: Water-saturation estimation from seismic and rock-property](https://reader031.vdocument.in/reader031/viewer/2022012103/616a0b2111a7b741a34e2fb2/html5/thumbnails/30.jpg)
Thank you for your attention!Thank you for your attention!Thank you for your attention!