We apply a rock physics analysis to well log data from the North-East Gulf of Mexico

Size: px
Start display at page:

Download "We apply a rock physics analysis to well log data from the North-East Gulf of Mexico"

Transcription

1 Rock Physics for Fluid and Porosity Mapping in NE GoM JACK DVORKIN, Stanford University and Rock Solid Images TIM FASNACHT, Anadarko Petroleum Corporation RICHARD UDEN, MAGGIE SMITH, NAUM DERZHI, AND JOEL WALLS, Rock Solid Images May 07, 2003 We apply a rock physics analysis to well log data from the North-East Gulf of Mexico to establish an effective-medium transform between the acoustic and elastic impedance on the one hand and lithology, pore fluid, and porosity on the other hand. These transforms are upscaled and applied to acoustic and elastic impedance inversion volumes to map lithology and porosity. INTRODUCTION. The ultimate goal of the geoscientist is to determine the reservoir properties (lithology and porosity) and conditions (pore fluid and pressure) from seismic data. This goal can be achieved by applying rock physics transforms to a volume of seismically-derived elastic properties. The basis for establishing rock physics transforms are controlled experiments where the reservoir and elastic sediment properties are measured on the same sample of earth at the same conditions. Such experimental data come from well logs and cores. The rock physics diagnostic is aimed at finding a deterministic effective-medium model appropriate for the rock under examination, calibrating it to the data, and then using it as the ultimate link (transform) between various rock properties. We select this rational deterministic approach, as opposed to pure statistics, in order to (a) understand and generalize trends seen in the data and use them outside of ranges available in the data; and (b) determine the domains of applicability of the trends. Below, the rock physics diagnostic is applied to well log data from the NE GoM and then used to map fluid and quantify porosity from acoustic and elastic impedance inversion volumes. WELL LOG DATA. The interval under examination includes shales and several high- 1

2 porosity sand packages (Figure 1). The upper sands are water-saturated while the lower sand contains hydrocarbons. The hydrocarbon-filled sand layer includes high-impedance calcite streaks. The S-wave data are not available in the well. One goal of this study is to quantify the porosity of the reservoir from seismicallyderived elastic parameters. This is why in Figure 1 we display an impedance-porosity cross-plot together with the depth plots of gamma-ray (GR), saturation, porosity, and impedance. The shale and wet sand form two separate branches in the impedanceporosity frame. In the same range of the total porosity, the impedance in the shale is smaller than in the wet sand. In both lithologies, evident are distinctive impedanceporosity trends. The impedance in the high-porosity hydrocarbon sand is smaller than in the wet sand of comparable porosity mostly due to the effect of the pore fluid on the bulk modulus. This is why the impedance-porosity trend present in the hydrocarbon reservoir is different from that observed in the wet sand located above. ROCK PHYSICS DIAGNOSTIC. In order to find an effective-medium model that will provide a porosity-lithology-impedance transform, we first theoretically saturate, using the V p -only fluid substitution, the entire interval under examination with the same pore fluid, namely, the formation brine. This pore-fluid equalization (we also call this procedure common fluid denominator ) allows us to concentrate purely on the effects of lithology and porosity on the elastic properties. The impedance in the brine-substituted interval is plotted versus porosity in Figure 2, top. It is color-coded by GR and the in-situ water saturation. The model-based impedance-porosity curves from the uncemented sand/shale model are superimposed upon the data points in the same figure. The sand data points are bound by the curves drawn for pure quartz sand and sand with 25% clay content. Model-based impedance-porosity curves are also drawn for rock with 75% and 100% clay content. These curves serve at lower bounds for the shale data. Most of the shale data lie above the 75% clay content curve which simply means that the shale is not pure clay. 2

3 Figure 1. Well log curves (from left to right, GR, water saturation, the total porosity, and impedance) and the impedance-porosity cross-plot. In the second, third, and fourth rows, we highlight (in red) the shale, wet sand, and hydrocarbon sand, respectively. Depth is in fictitious units. The original (not fluid-substituted) impedance is plotted versus porosity in Figure 2, 3

4 bottom. The uncemented sand/shale curve drawn for sand with 10% clay and for the insitu reservoir fluid is superimposed on the data. It can serve as an average impedanceporosity transform in the reservoir sand. The corresponding best-fit equations are: I p = f f 2, f = I p I p 2, (1) where the impedance ( I p ) is in km/s g/cc and porosity (f) is in fraction. Figure 2. Impedance versus total porosity color-coded by GR (left) and water saturation (right). Top: Impedance calculated for the brine-saturated interval, including the reservoir sand. The model curves come from the uncemented sand/shale model and are calculated for brine-saturated conditions. The curves are (top to bottom) for zero, 25%, 75%, and 100% clay content, respectively. Bottom: The insitu impedance. The yellow curve is for 10% clay content and the in-situ reservoir fluid. A comparison between the in-situ velocity and impedance data and the velocity and 4

5 impedance calculated using the uncemented sand/shale model with the in-situ clay content, porosity, and saturation as input (Figure 3) indicates that the selected rock physics transform is appropriate for the interval under examination. Figure 3. Comparison between sonic data and modeled data. Well log curves (from left to right, GR, water saturation, P-wave velocity, and impedance). In the velocity and impedance frames, the black curves represent the well log data while the red curves are computed from the uncemented sand/shale model using the clay content, porosity, and saturation curves as input. Depth is in fictitious units. In Figure 4 we compare the porosity measured in the pay zone with that predicted from the measured impedance using the transform from the second Equation (1). The prediction is very close to the data. Figure 4. Comparison between porosity data and prediction in the pay zone. In the porosity frame, the black curve is measured porosity while the red curve is porosity calculated from the measured impedance using the transform from Equation (1). Based on these results, we adopt the uncemented sand/shale model as appropriate for 5

6 the well under examination. As we have already seen, this model can serve to predict porosity from impedance. It can also be used to calculate the S-wave velocity and Poisson s Ratio (PR) which are not available at the well. Our calculations indicate that, as expected, PR in the reservoir sand is low (about 0.13); PR in the wet sand is about 0.28; and PR in the shale is about UPSCALING. Elastic waves in the seismic frequency range usually sample large portions of the subsurface whose dimensions are on the order of the wavelength. As a result, the elastic structure of the subsurface extracted from seismic data can only represent the elastic properties averaged over relatively large intervals. The details that are apparent at the well log scale cannot be recovered from seismic data. To illustrate this effect, consider synthetic seismic traces modeled at the well under examination using a 50 Hz and 25 Hz Ricker wavelet and the velocity (V p and the modelpredicted V s ) and bulk density at the well (Figure 5). This synthetic seismic modeling is based on the convolution of the wavelet with the reflectivity series derived at the well. The reflectivity R pp (q) at an angle of incidence q is calculated according to Hilterman s approximation to the Zoeppritz equations: R pp (q) = R pp (0)cos 2 q Dn sin 2 q, R pp (0) = D(0.5ln I p ), (2) where n is PR and D indicates the difference between the lower and upper half-space at the reflecting interface. The pay zone is manifested by a positive reflection whose amplitude decreases with offset. The positive normal reflection is due to the gradual impedance increase from the shale above the pay zone, through the pay zone, and to the underlying hard shale. The AVO effect is due to the small Poisson s ratio characteristic of porous sand with hydrocarbon. Also in Figure 5 shown are the pseudo-impedance ( I pp ) and pseudo-pr (n p ) curves which were calculated from the synthetic traces by using Hilterman s AVO approximation the pseudo-impedance is simply the cumulative sum of the normalincidence trace while pseudo-pr is a linear combination of the cumulative sum of the far- 6

7 offset trace and the pseudo-impedance to produce Poisson s ratio reflectivity: I pp = t Ú T 0 (t)dt, n p = [cos -2 q Ú T q (t)dt - I pp ]/(2.25tan 2 q), (3) 0 0 t where T 0 is the normal-incidence trace; T q is the trace at angle q ; and t is the two-way travel time. Figure 5. Synthetic seismic at the well. From left to right: Synthetic gather (black) and stack (red); impedance and PR at the well; and pseudo-impedance and pseudo-pr calculated from the synthetic traces. Travel time is in fictitious units. The pseudo-impedance and pseudo-pr curves presented in Figure 5 capture the 7

8 character of the impedance and PR at the well but, at the same time, omit the details and, effectively, smooth the well log curves. Because of this smoothing effect, it is not clear whether the impedance-porosity transform expressed in Equation (1) can be useful at the seismic scale, how it has to be corrected (upscaled), and, finally, how to interpret porosity estimates that are based on the seismic-scale impedance values. To investigate the effects of upscaling on the elastic properties and also on rock physics transforms, we use the Backus average which estimates the seismic-scale response by harmonically averaging of the elastic moduli at the well. The bulk density is upscaled by means of arithmetic averaging. Then the upscaled impedance and PR are calculated from the Backus-upscaled elastic moduli and arithmetically-upscaled density according to the appropriate theory-of-elasticity equations. The upscaled elastic curves as well as the arithmetically-averaged porosity curve are shown in Figure 6. This display confirms what we already learned from the synthetic seismic traces: the log-scale details of the interval will be omitted in seismic data. Figure 6. Original (black) and upscaled (bold cyan) porosity, impedance, and PR curves at the well. Depth is in fictitious units. UPSCALED IMPEDANCE-POROSITY TRANSFORM. A question to address is whether the transform from the impedance to porosity derived from well log data -- 2 f = I p I p -- is valid at the seismic scale. To investigate this problem, let us plot the porosity versus impedance at the log scale and seismic scale (Figure 7). The upscaled trend generally follows the original porosity-impedance trend but slightly overestimates the porosity at a given impedance value. The best-linear-fit 8

9 transform between the upscaled impedance and averaged porosity is f u = (I pu - 6.5), (4) where the impedance is in km/s g/cc and porosity is in fraction, and the subscript u refers to the upscaled values. This equation is site-specific and should not be generalized because the results of upscaling may be strongly influenced by the size and shape of the reservoir and surrounding strata. In general, the validity of log-scale rock physics transforms at the seismic scale improves as the reservoir becomes thicker and the variations of the rock properties with depth around the reservoir become smoother. Figure 7. Log-scale porosity versus log-scale impedance (black); log-scale porosity versus upscaled impedance (blue); and upscaled porosity versus upscaled impedance (red). The yellow line is from Equation (4). INFLUENCE OF HIGH-IMPEDANCE LAYERS ON REFLECTION. Let us investigate is whether the high-impedance layers that bound the reservoir affect the seismic reflection and, eventually, porosity estimates from the seismic impedance. To address this question, we alter the elastic curves at the well by reducing the impedance in the three high-impedance layers to the surrounding impedance values (Figure 8). Then we calculate the Backus-upscaled impedance curves. 9

10 Figure 8. The impedance curve at the reservoir. Black is for the measured impedance; red is for the impedance artificially reduced within the high-impedance layers; cyan is for the upscaled original impedance curve; and yellow is for the upscaled reduced impedance curve. Depth is in fictitious units. The upscaling results indicate that the influence of the high-impedance layers on the upscaled impedance values is small. To further support this conclusion, we generate synthetic traces using the impedance curves at the well with and without the highimpedance layers (Figure 9). The reflections appear to be virtually unaffected by the high-impedance layers. This apparently counter-intuitive result is due to the fact that the averaging of the elastic properties in the subsurface is harmonic and as such gives larger weights to the softer parts of the interval. POROSITY FROM UPSCALED IMPEDANCE. Next we apply Equation (4) to the upscaled impedance curve to assess the accuracy of porosity estimates from the seismicscale impedance. The estimated porosity is compared to the log-scale porosity and also to the arithmetically-averaged porosity in Figure 10. The porosity curve derived from the upscaled impedance is essentially the same as the arithmetically averaged porosity curve. Both curves give a reasonably accurate average estimate to the log-scale porosity. 10

11 Figure 9. Synthetic traces with (top) and without (bottom) high-impedance layers in the reservoir zone. The display is the same as in Figure 5. 11

12 Figure 10. Left: Water saturation curve in the reservoir zone. Right: the log-scale porosity curve (black); arithmetically upscaled porosity curve (red); and porosity derived from the upscaled impedance (yellow). Depth is in fictitious units. POROSITY FROM SEISMIC IMPEDANCE INVERSION. Next we use acoustic and elastic impedance inversion sections obtained from real seismic data to identify the reservoir and estimate the total porosity within the reservoir. In the section displayed in Figure 11, the reservoir is identified by a low-value-pr strip, as is expected from the above rock physics analysis and pseudo-impedance and pseudo-pr inversion curves shown in Figures 5 and 9. The minimum-value strip in the PR inversion section does not precisely coincide in space with the maximum-value strip in the impedance-inversion section. This is expected from the impedance structure of the well log profile where the maximum of the impedance is located slightly below the reservoir while the minimum of the PR is located precisely at the reservoir. The reservoir in the inversion section is identified as the low-pr strip. The porosity in the reservoir is directly calculated from the impedance inversion using Equation (4). The calculated porosity section is displayed in Figure 11. The same approach as used in the vertical section is employed to create a seismic porosity cube from 3D acoustic and elastic impedance inversion (Figure 12). The dipping sand layer filled with hydrocarbon is apparent in the volume. The curved section of the sand by a constant-time horizon visible in the 3D display is an apparent indication of a three-way closure. 12

13 Figure 11. Left: PR inversion from real seismic. Middle: Acoustic impedance inversion from real seismic. Right: Seismic porosity derived from the impedance inversion in the reservoir zone. Depth and cross-line distance are in fictitious units. The vertical bars indicate the location of the well. The well log curves of saturation (left); GR (middle); and porosity (right) are superimposed upon the vertical sections. Figure 12. Mapping fluid and porosity from 3D acoustic and elastic impedance inversion. From left to right: PR inversion; acoustic impedance inversion; and total porosity produced from the acoustic impedance according to Equation (4). CUMULATIVE POROSITY. An important desired reservoir characterization parameter is the cumulative porosity ( C f ) also known as the total pore volume in the reservoir, or net-to-gross. This quantity is often expressed as the product of the porosity (f) and reservoir thickness ( h): C f = fh. Strictly speaking, it is the integral of the porosity with respect to depth, calculated within the reservoir: C f = z B Ú f(z)dz, (5) z T 13

14 where the integration is between the top ( z T ) and bottom ( z B ) depths of the reservoir. Moreover, only the accumulated porosity within the pay zone is of immediate interest. To calculate the cumulative pay-zone porosity ( C fp ), Equation (5) should be modified as: C fp = z B Ú f(z) Sw <1 dz, (6) z T where condition S w <1 simply means that during the integration, porosity outside of the pay zone is not counted. The cumulative pay-zone porosity is plotted versus depth in Figure 13 for the original (log-scale) curve as well as for the upscaled porosity. The two curves are very close to each other which means that the seismically derived C fp and the log-scale C fp are essentially the same. Figure 13. Cumulative porosity (left); cumulative impedance (middle); and cumulative inverse impedance (right) versus depth. Black curves are for the log scale values while yellow curves are for the upscaled values. The cumulative porosity is in fraction times meter; the cumulative impedance is in km/s g/cc times meter; and the cumulative inverse impedance is in s/km cc/g times meter. CUMULATIVE ATTRIBUTES. As we have already established, porosity can be related to the impedance. Next question is which seismic attribute can be used to estimate the cumulative porosity. One possibility is that such a cumulative attribute (CATT) is the cumulative impedance in the pay zone (C I p P ) which can be defined in the same manner as C fp : 14

15 C I p P = z B Ú I p (z) Sw <1 dz. (7) z T C I p P is plotted versus depth in Figure 13 (middle) for the original log-scale curve as well as for the upscaled curve. Once again, the two curves are very close to each other which means that the proposed integration operator will produce approximately the same results at the seismic scale as at the log scale. The next question is whether the proposed CATT can be uniquely translated into the cumulative porosity within a seismic volume. To investigate this possibility, let us crossplot C fp versus C I p P at the log and seismic scale. The result shown in Figure 14 (left) indicates that the cross-plot curves at both scales are close to each other which means that the proposed CATT (C I p P ) can be used to calculate the cumulative porosity from seismic. Figure 14. Cumulative porosity versus cumulative impedance (left); and versus cumulative inverse impedance (right). Black curves are for the log scale values while yellow curves are for the upscaled values. Another possible CATT may be the inverse cumulative impedance (C I p -1 P ): C I p -1 P = z B Ú I -1 p (z) dz. Sw <1 (8) z T This CATT is plotted versus depth in Figure 13 (right). The cumulative porosity is cross-plotted versus C I in Figure 14, right. In this cross-plot, the match between the -1 p P log-scale and seismic-scale values is even better than for the cumulative impedance on 15

16 the left. This is likely due to the fact that the elastic moduli upscale by means of harmonic averaging and, in this particular case, such averaging rule may be approximately valid for the impedance. The approach proposed here essentially introduces a new class of seismic attributes which are calculated from a seismic trace by repeated integration -- the seismic impedance is due to a one-time integration of a trace which is approximately treated as the reflectivity series and the cumulative impedance is due to one-time integration of the seismic impedance or repeated integration of the trace. CAVEAT -- DELINEATING PAY ZONE. An important condition for calculating C fp from a CATT was to integrate the latter within the pay zone. To meet this condition we will have to employ an additional attribute, such as the seismic PR, to locate the pay zone. Finding a threshold for this attribute to accurately delineate a pay is important for a successful use of CATTS. CONCLUSION. The rock physics diagnostic approach enables us to establish a sitespecific effective-medium model to transform the elastic rock properties into lithology, fluid, and porosity. An effective-medium model allows us to systematically and consistently perturb rock properties and conditions to predict the elastic response away from a well. A model established for log-scale data has to be upscaled to the seismic scale in order to be used with seismic inversion cubes. Such upscaling may call for corrections in the model, especially so in a heterogeneous environment with large elastic contrasts. ACKNOWLEDGEMENT. Uwe Strecker helped with stratigraphic and geological understanding of the example presented in this study. SUGGESTED READING Dvorkin, J., and Nur, A., 1996, Elasticity of High-Porosity Sandstones: Theory for Two North Sea Datasets, Geophysics, 61, Mavko, G., Chan, C., and Mukerji, T., 1995, Fluid substitution: Estimating changes in Vp without knowing Vs, Geophysics, 60, Backus, G.F., 1962, Long-wave elastic anisotropy produced by horizontal layering, JGR, 67,

RESEARCH PROPOSAL. Effects of scales and extracting methods on quantifying quality factor Q. Yi Shen

RESEARCH PROPOSAL. Effects of scales and extracting methods on quantifying quality factor Q. Yi Shen RESEARCH PROPOSAL Effects of scales and extracting methods on quantifying quality factor Q Yi Shen 2:30 P.M., Wednesday, November 28th, 2012 Shen 2 Ph.D. Proposal ABSTRACT The attenuation values obtained

More information

Rock physics of a gas hydrate reservoir. Gas hydrates are solids composed of a hydrogen-bonded ROUND TABLE

Rock physics of a gas hydrate reservoir. Gas hydrates are solids composed of a hydrogen-bonded ROUND TABLE ROUND TABLE Rock physics of a gas hydrate reservoir JACK DVORKIN and AMOS NUR, Stanford University, California, U.S. RICHARD UDEN and TURHAN TANER, Rock Solid Images, Houston, Texas, U.S. Gas hydrates

More information

Estimating the hydrocarbon volume from elastic and resistivity data: A concept

Estimating the hydrocarbon volume from elastic and resistivity data: A concept INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Estimating the hydrocarbon volume from elastic and resistivity data: A concept CARMEN T. GOMEZ, JACK DVORKIN, and GARY MAVKO, Stanford University,

More information

Downloaded 11/02/16 to Redistribution subject to SEG license or copyright; see Terms of Use at Summary.

Downloaded 11/02/16 to Redistribution subject to SEG license or copyright; see Terms of Use at   Summary. in thin sand reservoirs William Marin* and Paola Vera de Newton, Rock Solid Images, and Mario Di Luca, Pacific Exploración y Producción. Summary Rock Physics Templates (RPTs) are useful tools for well

More information

ROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES ABSTRACT

ROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES ABSTRACT ROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES PER AVSETH, JACK DVORKIN, AND GARY MAVKO Department of Geophysics, Stanford University, CA 94305-2215, USA

More information

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad*

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad* Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad* Carl Reine 1, Chris Szelewski 2, and Chaminda Sandanayake 3 Search and Discovery Article #41899 (2016)** Posted September

More information

Rock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa

Rock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa Rock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa Franklin Ruiz, Carlos Cobos, Marcelo Benabentos, Beatriz Chacon, and Roberto Varade, Luis Gairifo,

More information

SEG Houston 2009 International Exposition and Annual Meeting. that the project results can correctly interpreted.

SEG Houston 2009 International Exposition and Annual Meeting. that the project results can correctly interpreted. Calibration of Pre-Stack Simultaneous Impedance Inversion using Rock Physics Scott Singleton and Rob Keirstead, Rock Solid Images Log Conditioning and Rock Physics Modeling Summary Geophysical Well Log

More information

RC 1.3. SEG/Houston 2005 Annual Meeting 1307

RC 1.3. SEG/Houston 2005 Annual Meeting 1307 from seismic AVO Xin-Gong Li,University of Houston and IntSeis Inc, De-Hua Han, and Jiajin Liu, University of Houston Donn McGuire, Anadarko Petroleum Corp Summary A new inversion method is tested to directly

More information

BPM37 Linking Basin Modeling with Seismic Attributes through Rock Physics

BPM37 Linking Basin Modeling with Seismic Attributes through Rock Physics BPM37 Linking Basin Modeling with Seismic Attributes through Rock Physics W. AlKawai* (Stanford University), T. Mukerji (Stanford University) & S. Graham (Stanford University) SUMMARY In this study, we

More information

Effects of fluid changes on seismic reflections: Predicting amplitudes at gas reservoir directly from amplitudes at wet reservoir

Effects of fluid changes on seismic reflections: Predicting amplitudes at gas reservoir directly from amplitudes at wet reservoir GEOPHYSICS, VOL. 77, NO. 4 (JULY-AUGUST 2012); P. D129 D140, 15 FIGS., 2 TABLES. 10.1190/GEO2011-0331.1 Effects of fluid changes on seismic reflections: Predicting amplitudes at gas reservoir directly

More information

Seismic reservoir characterization of a U.S. Midcontinent fluvial system using rock physics, poststack seismic attributes, and neural networks

Seismic reservoir characterization of a U.S. Midcontinent fluvial system using rock physics, poststack seismic attributes, and neural networks CORNER INTERPRETER S Coordinated by Linda R. Sternbach Seismic reservoir characterization of a U.S. Midcontinent fluvial system using rock physics, poststack seismic attributes, and neural networks JOEL

More information

Rock physics and AVO applications in gas hydrate exploration

Rock physics and AVO applications in gas hydrate exploration Rock physics and AVO applications in gas hydrate exploration ABSTRACT Yong Xu*, Satinder Chopra Core Lab Reservoir Technologies Division, 301,400-3rd Ave SW, Calgary, AB, T2P 4H2 yxu@corelab.ca Summary

More information

IDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note. / K s. + K f., G Dry. = G / ρ, (2)

IDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note. / K s. + K f., G Dry. = G / ρ, (2) IDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note JACK DVORKIN, DAN MOOS, JAMES PACKWOOD, AND AMOS NUR DEPARTMENT OF GEOPHYSICS, STANFORD UNIVERSITY January 5, 2001 INTRODUCTION Gassmann's (1951)

More information

Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties

Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties Amaro C. 1 Abstract: The main goal of reservoir modeling and characterization is the inference of

More information

An overview of AVO and inversion

An overview of AVO and inversion P-486 An overview of AVO and inversion Brian Russell, Hampson-Russell, CGGVeritas Company Summary The Amplitude Variations with Offset (AVO) technique has grown to include a multitude of sub-techniques,

More information

Reservoir Characterization using AVO and Seismic Inversion Techniques

Reservoir Characterization using AVO and Seismic Inversion Techniques P-205 Reservoir Characterization using AVO and Summary *Abhinav Kumar Dubey, IIT Kharagpur Reservoir characterization is one of the most important components of seismic data interpretation. Conventional

More information

Pressure and Compaction in the Rock Physics Space. Jack Dvorkin

Pressure and Compaction in the Rock Physics Space. Jack Dvorkin Pressure and Compaction in the Rock Physics Space Jack Dvorkin June 2002 0 200 Compaction of Shales Freshly deposited shales and clays may have enormous porosity of ~ 80%. The speed of sound is close to

More information

The role of seismic modeling in Reservoir characterization: A case study from Crestal part of South Mumbai High field

The role of seismic modeling in Reservoir characterization: A case study from Crestal part of South Mumbai High field P-305 The role of seismic modeling in Reservoir characterization: A case study from Crestal part of South Mumbai High field Summary V B Singh*, Mahendra Pratap, ONGC The objective of the modeling was to

More information

An empirical method for estimation of anisotropic parameters in clastic rocks

An empirical method for estimation of anisotropic parameters in clastic rocks An empirical method for estimation of anisotropic parameters in clastic rocks YONGYI LI, Paradigm Geophysical, Calgary, Alberta, Canada Clastic sediments, particularly shale, exhibit transverse isotropic

More information

A New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset*

A New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset* A New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset* Changcheng Liu 1 and Prasad Ghosh 2 Search and Discovery Article #41764 (2016) Posted January 25, 2016 *Adapted

More information

Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics

Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics Peipei Li 1 and Tapan Mukerji 1,2 1 Department of Energy Resources Engineering 2 Department of

More information

The elastic properties such as velocity, density, impedance,

The elastic properties such as velocity, density, impedance, SPECIAL SECTION: Rr ock Physics physics Lithology and fluid differentiation using rock physics template XIN-GANG CHI AND DE-HUA HAN, University of Houston The elastic properties such as velocity, density,

More information

Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc.

Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc. Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc. Summary Rock physics establishes the link between reservoir properties,

More information

Th SBT1 14 Seismic Characters of Pore Pressure Due to Smectite-to-illite Transition

Th SBT1 14 Seismic Characters of Pore Pressure Due to Smectite-to-illite Transition Th SBT1 14 Seismic Characters of Pore Pressure Due to Smectite-to-illite Transition X. Qin* (University of Houston) & D. Han (University of Houston) SUMMARY In this study, we strive to understand unloading

More information

Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field

Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field Downloaded 09/12/14 to 84.215.159.82. Redistribution subject to SEG license or copyright; see Terms of Use

More information

Fluid-property discrimination with AVO: A Biot-Gassmann perspective

Fluid-property discrimination with AVO: A Biot-Gassmann perspective Fluid-property discrimination with AVO: A Biot-Gassmann perspective Brian H. Russell, Ken Hedlin 1, Fred J. Hilterman, and Laurence R. Lines ABSTRACT This paper draws together basic rock physics, AVO,

More information

Edinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains

Edinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains Frequency-dependent AVO attribute: theory and example Xiaoyang Wu, 1* Mark Chapman 1,2 and Xiang-Yang Li 1 1 Edinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains Road, Edinburgh

More information

Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston

Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston ain enu Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston Summary When S-wave velocity is absent, approximate Gassmann

More information

AVO analysis of 3-D seismic data at G-field, Norway

AVO analysis of 3-D seismic data at G-field, Norway Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 2014, 5(5):148-164 ISSN: 0976-8610 CODEN (USA): AASRFC AVO analysis of 3-D seismic data at G-field, Norway Fredrick

More information

Recent advances in application of AVO to carbonate reservoirs: case histories

Recent advances in application of AVO to carbonate reservoirs: case histories Recent advances in application of AVO to reservoirs: case histories Yongyi Li, Bill Goodway*, and Jonathan Downton Core Lab Reservoir Technologies Division *EnCana Corporation Summary The application of

More information

Reservoir properties inversion from AVO attributes

Reservoir properties inversion from AVO attributes Reservoir properties inversion from AVO attributes Xin-gang Chi* and De-hua Han, University of Houston Summary A new rock physics model based inversion method is put forward where the shaly-sand mixture

More information

The reason why acoustic and shear impedances inverted

The reason why acoustic and shear impedances inverted SPECIAL The Rocky SECTION: Mountain The Rocky region Mountain region Comparison of shear impedances inverted from stacked PS and SS data: Example from Rulison Field, Colorado ELDAR GULIYEV, Occidental

More information

Shaly Sand Rock Physics Analysis and Seismic Inversion Implication

Shaly Sand Rock Physics Analysis and Seismic Inversion Implication Shaly Sand Rock Physics Analysis and Seismic Inversion Implication Adi Widyantoro (IkonScience), Matthew Saul (IkonScience/UWA) Rock physics analysis of reservoir elastic properties often assumes homogeneity

More information

THE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOA-BOMBAL FIELDS

THE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOA-BOMBAL FIELDS THE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOA-BOMBAL FIELDS Cuesta, Julián* 1, Pérez, Richard 1 ; Hernández, Freddy 1 ; Carrasquel, Williams 1 ; Cabrera,

More information

Rock physics integration of CSEM and seismic data: a case study based on the Luva gas field.

Rock physics integration of CSEM and seismic data: a case study based on the Luva gas field. Rock physics integration of CSEM and seismic data: a case study based on the Luva gas field. Peter Harris*, Zhijun Du, Harald H. Soleng, Lucy M. MacGregor, Wiebke Olsen, OHM-Rock Solid Images Summary It

More information

Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field

Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field Amit Suman and Tapan Mukerji Department of Energy Resources Engineering Stanford University

More information

Pre-Stack Seismic Inversion and Amplitude Versus Angle Modeling Reduces the Risk in Hydrocarbon Prospect Evaluation

Pre-Stack Seismic Inversion and Amplitude Versus Angle Modeling Reduces the Risk in Hydrocarbon Prospect Evaluation Advances in Petroleum Exploration and Development Vol. 7, No. 2, 2014, pp. 30-39 DOI:10.3968/5170 ISSN 1925-542X [Print] ISSN 1925-5438 [Online] www.cscanada.net www.cscanada.org Pre-Stack Seismic Inversion

More information

Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait

Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait Osman Khaled, Yousef Al-Zuabi, Hameed Shereef Summary The zone under study is Zubair formation of Cretaceous

More information

GRAIN SORTING, POROSITY, AND ELASTICITY. Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University ABSTRACT

GRAIN SORTING, POROSITY, AND ELASTICITY. Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University ABSTRACT GRAIN SORTING, POROSITY, AND ELASTICITY Jack Dvorkin and Mario A. Gutierrez Geophysics Department, Stanford University July 24, 2001 ABSTRACT Grain size distribution (sorting) is determined by deposition.

More information

Integration of rock attributes to discriminate Miocene reservoirs for heterogeneity and fluid variability

Integration of rock attributes to discriminate Miocene reservoirs for heterogeneity and fluid variability P 397 Integration of rock attributes to discriminate Miocene reservoirs for heterogeneity and fluid variability Debajyoti Guha*, Ashok Yadav, Santan Kumar, Naina Gupta, Ashutosh Garg *Reliance Industries

More information

The Stanford VI reservoir

The Stanford VI reservoir The Stanford VI reservoir Scarlet A. Castro, Jef Caers and Tapan Mukerji Stanford Center for Reservoir Forecasting Stanford Rock Physics and Borehole Geophysics Project May 2005 Abstract An exhaustive

More information

AVO responses for varying Gas saturation sands A Possible Pathway in Reducing Exploration Risk

AVO responses for varying Gas saturation sands A Possible Pathway in Reducing Exploration Risk P 420 AVO responses for varying Gas saturation sands A Possible Pathway in Reducing Exploration Risk Naina Gupta*, Debajyoti Guha, Ashok Yadav, Santan Kumar, Ashutosh Garg *Reliance Industries Limited,

More information

Seismic reservoir characterization using the multiattribute rotation scheme: Case study in the South Falkland Basin

Seismic reservoir characterization using the multiattribute rotation scheme: Case study in the South Falkland Basin Seismic reservoir characterization using the multiattribute rotation scheme: Case study in the South Falkland Basin Pedro Alvarez 1, Bruce Farrer 2, and Doyin Oyetunji 1 Downloaded 08/15/16 to 50.207.209.154.

More information

Statistical Rock Physics

Statistical Rock Physics Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics

More information

2011 SEG SEG San Antonio 2011 Annual Meeting 771. Summary. Method

2011 SEG SEG San Antonio 2011 Annual Meeting 771. Summary. Method Geological Parameters Effecting Controlled-Source Electromagnetic Feasibility: A North Sea Sand Reservoir Example Michelle Ellis and Robert Keirstead, RSI Summary Seismic and electromagnetic data measure

More information

Summary. Seismic Field Example

Summary. Seismic Field Example Water-saturation estimation from seismic and rock-property trends Zhengyun Zhou*, Fred J. Hilterman, Haitao Ren, Center for Applied Geosciences and Energy, University of Houston, Mritunjay Kumar, Dept.

More information

Predicting Gas Hydrates Using Prestack Seismic Data in Deepwater Gulf of Mexico (JIP Projects)

Predicting Gas Hydrates Using Prestack Seismic Data in Deepwater Gulf of Mexico (JIP Projects) Predicting Gas Hydrates Using Prestack Seismic Data in Deepwater Gulf of Mexico (JIP Projects) Dianna Shelander 1, Jianchun Dai 2, George Bunge 1, Dan McConnell 3, Niranjan Banik 2 1 Schlumberger / DCS

More information

QUANTITATIVE INTERPRETATION

QUANTITATIVE INTERPRETATION QUANTITATIVE INTERPRETATION THE AIM OF QUANTITATIVE INTERPRETATION (QI) IS, THROUGH THE USE OF AMPLITUDE ANALYSIS, TO PREDICT LITHOLOGY AND FLUID CONTENT AWAY FROM THE WELL BORE This process should make

More information

Quick Look : Rule Of Thumb of Rock Properties in Deep Water, Krishna-Godavari Basin

Quick Look : Rule Of Thumb of Rock Properties in Deep Water, Krishna-Godavari Basin 5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 576-581 Quick Look : Rule Of Thumb of Rock Properties in Deep Water, Krishna-Godavari Basin Acharya M. N. 1, Biswal A. K. 1

More information

Quantitative interpretation using inverse rock-physics modeling on AVO data

Quantitative interpretation using inverse rock-physics modeling on AVO data Quantitative interpretation using inverse rock-physics modeling on AVO data Erling Hugo Jensen 1, Tor Arne Johansen 2, 3, 4, Per Avseth 5, 6, and Kenneth Bredesen 2, 7 Downloaded 08/16/16 to 129.177.32.62.

More information

Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel

Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel Heriot-Watt University Heriot-Watt University Research Gateway Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel DOI: 10.3997/2214-4609.20132136

More information

INT 4.5. SEG/Houston 2005 Annual Meeting 821

INT 4.5. SEG/Houston 2005 Annual Meeting 821 Kathleen Baker* and Mike Batzle, Colorado School of Mines, Richard Gibson, Texas A&M University Summary There have been many studies of how Direct Hydrocarbon Indicators (DHIs) can help or hurt us when

More information

AFI (AVO Fluid Inversion)

AFI (AVO Fluid Inversion) AFI (AVO Fluid Inversion) Uncertainty in AVO: How can we measure it? Dan Hampson, Brian Russell Hampson-Russell Software, Calgary Last Updated: April 2005 Authors: Dan Hampson, Brian Russell 1 Overview

More information

Th LHR2 08 Towards an Effective Petroelastic Model for Simulator to Seismic Studies

Th LHR2 08 Towards an Effective Petroelastic Model for Simulator to Seismic Studies Th LHR2 08 Towards an Effective Petroelastic Model for Simulator to Seismic Studies A. Briceno* (Heriot-Watt University), C. MacBeth (Heriot-Watt University) & M.D. Mangriotis (Heriot-Watt University)

More information

Downloaded 10/02/18 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 10/02/18 to Redistribution subject to SEG license or copyright; see Terms of Use at Multi-scenario, multi-realization seismic inversion for probabilistic seismic reservoir characterization Kester Waters* and Michael Kemper, Ikon Science Ltd. Summary We propose a two tiered inversion strategy

More information

Petrophysical Study of Shale Properties in Alaska North Slope

Petrophysical Study of Shale Properties in Alaska North Slope Petrophysical Study of Shale Properties in Alaska North Slope Minh Tran Tapan Mukerji Energy Resources Engineering Department Stanford University, CA, USA Region of Interest 1.5 miles 20 miles Stratigraphic

More information

LINK BETWEEN ATTENUATION AND VELOCITY DISPERSION

LINK BETWEEN ATTENUATION AND VELOCITY DISPERSION LINK BETWEEN ATTENUATION AND VELOCITY DISPERSION Jack Dvorkin Stanford University and Rock Solid Images April 25, 2005 SUMMARY In a viscoelastic sample, the causality principle links the attenuation of

More information

Downloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at AVO crossplot analysis in unconsolidated sediments containing gas hydrate and free gas: Green Canyon 955, Gulf of Mexico Zijian Zhang* 1, Daniel R. McConnell 1, De-hua Han 2 1 Fugro GeoConsulting, Inc.,

More information

Stochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell

Stochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell Stochastic vs Deterministic Pre-stack Inversion Methods Brian Russell Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological

More information

Geophysical model response in a shale gas

Geophysical model response in a shale gas Geophysical model response in a shale gas Dhananjay Kumar and G. Michael Hoversten Chevron USA Inc. Abstract Shale gas is an important asset now. The production from unconventional reservoir like shale

More information

Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function

Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function Jason M. McCrank, Gary F. Margrave, and Don C. Lawton ABSTRACT The Bayesian inference is used to estimate

More information

Delineating a sandstone reservoir at Pikes Peak, Saskatchewan using 3C seismic data and well logs

Delineating a sandstone reservoir at Pikes Peak, Saskatchewan using 3C seismic data and well logs Delineating a sandston reservoir at Pikes Peak Delineating a sandstone reservoir at Pikes Peak, Saskatchewan using 3C seismic data and well logs Natalia L. Soubotcheva and Robert R. Stewart ABSTRACT To

More information

AVO Attributes of a Deep Coal Seam

AVO Attributes of a Deep Coal Seam AVO Attributes of a Deep Coal Seam Jinfeng Ma* State Key Laboratory of Continental Dynamics, Northwest University, China jinfengma@sohu.com Igor Morozov University of Saskatchewan, Saskatoon, SK, Canada

More information

Practical aspects of AVO modeling

Practical aspects of AVO modeling Practical aspects of AVO modeling YONGYI LI, Paradigm Geophysical, Calgary, Canada JONATHAN DOWNTON, Veritas, Calgary, Canada, YONG XU, Arcis Corporation, Calgary, Canada AVO (amplitude variation with

More information

Quantitative Interpretation

Quantitative Interpretation Quantitative Interpretation The aim of quantitative interpretation (QI) is, through the use of amplitude analysis, to predict lithology and fluid content away from the well bore. This process should make

More information

Fred Mayer 1; Graham Cain 1; Carmen Dumitrescu 2; (1) Devon Canada; (2) Terra-IQ Ltd. Summary

Fred Mayer 1; Graham Cain 1; Carmen Dumitrescu 2; (1) Devon Canada; (2) Terra-IQ Ltd. Summary 2401377 Statistically Improved Resistivity and Density Estimation From Multicomponent Seismic Data: Case Study from the Lower Cretaceous McMurray Formation, Athabasca Oil Sands Fred Mayer 1; Graham Cain

More information

Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin

Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin P-276 Summary Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin K.L.Mandal*, R.K.Srivastava, S.Saha, Oil India Limited M.K.Sukla, Indian Institute of Technology, Kharagpur

More information

Rock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E.

Rock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E. Rock Physics and Quantitative Wavelet Estimation for Seismic Interpretation: Tertiary North Sea R.W.Simm 1, S.Xu 2 and R.E.White 2 1. Enterprise Oil plc, Grand Buildings, Trafalgar Square, London WC2N

More information

Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration

Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration 66 Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration Kenneth Bredesen 1, Erling Hugo Jensen 1, 2, Tor Arne Johansen 1, 2, and Per Avseth 3, 4

More information

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY Ida Ayu PURNAMASARI*, Hilfan KHAIRY, Abdelazis Lotfy ABDELDAYEM Geoscience and Petroleum Engineering Department Universiti Teknologi

More information

Quantifying Bypassed Pay Through 4-D Post-Stack Inversion*

Quantifying Bypassed Pay Through 4-D Post-Stack Inversion* Quantifying Bypassed Pay Through 4-D Post-Stack Inversion* Robert Woock 1, Sean Boerner 2 and James Gamble 1 Search and Discovery Article #40799 (2011) Posted August 12, 2011 *Adapted from oral presentation

More information

Crosswell tomography imaging of the permeability structure within a sandstone oil field.

Crosswell tomography imaging of the permeability structure within a sandstone oil field. Crosswell tomography imaging of the permeability structure within a sandstone oil field. Tokuo Yamamoto (1), and Junichi Sakakibara (2) (1) University of Miami and Yamamoto Engineering Corporation, (2)

More information

Keywords. Unconsolidated Sands, Seismic Amplitude, Oil API

Keywords. Unconsolidated Sands, Seismic Amplitude, Oil API Role of Seismic Amplitude in Assessment of Oil API Ravi Mishra*(Essar Oil Ltd.), Baban Jee (Essar Oil Ltd.), Ashish Kumar (Essar Oil Ltd.) Email: ravimishragp@gmail.com Keywords Unconsolidated Sands, Seismic

More information

3D petrophysical modeling - 1-3D Petrophysical Modeling Usning Complex Seismic Attributes and Limited Well Log Data

3D petrophysical modeling - 1-3D Petrophysical Modeling Usning Complex Seismic Attributes and Limited Well Log Data - 1-3D Petrophysical Modeling Usning Complex Seismic Attributes and Limited Well Log Data Mehdi Eftekhari Far, and De-Hua Han, Rock Physics Lab, University of Houston Summary A method for 3D modeling and

More information

Shear Wave Velocity Estimation Utilizing Wireline Logs for a Carbonate Reservoir, South-West Iran

Shear Wave Velocity Estimation Utilizing Wireline Logs for a Carbonate Reservoir, South-West Iran Iranian Int. J. Sci. 4(2), 2003, p. 209-221 Shear Wave Velocity Estimation Utilizing Wireline Logs for a Carbonate Reservoir, South-West Iran Eskandari, H. 1, Rezaee, M.R., 2 Javaherian, A., 3 and Mohammadnia,

More information

Rock Physics of Shales and Source Rocks. Gary Mavko Professor of Geophysics Director, Stanford Rock Physics Project

Rock Physics of Shales and Source Rocks. Gary Mavko Professor of Geophysics Director, Stanford Rock Physics Project Rock Physics of Shales and Source Rocks Gary Mavko Professor of Geophysics Director, Stanford Rock Physics Project 1 First Question: What is Shale? Shale -- a rock composed of mud-sized particles, such

More information

Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics

Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics Summary Madhumita Sengupta*, Ran Bachrach, Niranjan Banik, esterngeco. Net-to-gross (N/G ) is

More information

P125 AVO for Pre-Resonant and Resonant Frequency Ranges of a Periodical Thin-Layered Stack

P125 AVO for Pre-Resonant and Resonant Frequency Ranges of a Periodical Thin-Layered Stack P125 AVO for Pre-Resonant and Resonant Frequency Ranges of a Periodical Thin-Layered Stack N. Marmalyevskyy* (Ukrainian State Geological Prospecting Institute), Y. Roganov (Ukrainian State Geological Prospecting

More information

SENSITIVITY ANALYSIS OF AMPLITUDE VARIATION WITH OFFSET (AVO) IN FRACTURED MEDIA

SENSITIVITY ANALYSIS OF AMPLITUDE VARIATION WITH OFFSET (AVO) IN FRACTURED MEDIA SENSITIVITY ANALYSIS OF AMPLITUDE VARIATION WITH OFFSET AVO) IN FRACTURED MEDIA Mary L. Krasovec, William L. Rodi, and M. Nafi Toksoz Earth Resources Laboratory Department of Earth, Atmospheric, and Planetary

More information

Reservoir And Elastic Property Prediction Away From Well Control Kyle T. Spikes and Jack P. Dvorkin; Stanford Rock Physics Laboratory

Reservoir And Elastic Property Prediction Away From Well Control Kyle T. Spikes and Jack P. Dvorkin; Stanford Rock Physics Laboratory Reservoir And Elastic Property Prediction Away From Well Control Kyle T. Spikes and Jack P. Dvorkin; Stanford Rock Physics Laboratory Abstract The ultimate goal of this effort is to quantify lithology,

More information

Framework for AVO gradient and intercept interpretation

Framework for AVO gradient and intercept interpretation GEOPHYSICS, VOL. 63, NO. 3 (MAY-JUNE 1998); P. 948 956, 10 FIGS., 2 TABLES. Framework for AVO gradient and intercept interpretation John P. Castagna, Herbert W. Swan, and Douglas J. Foster ABSTRACT Amplitude

More information

Oil and Natural Gas Corporation Ltd., VRC(Panvel), WOB, ONGC, Mumbai. 1

Oil and Natural Gas Corporation Ltd., VRC(Panvel), WOB, ONGC, Mumbai. 1 P-259 Summary Data for identification of Porosity Behaviour in Oligocene Lime Stone of D18 Area Of Western Offshore, India V.K. Baid*, P.H. Rao, P.S. Basak, Ravi Kant, V. Vairavan 1, K.M. Sundaram 1, ONGC

More information

Integrating rock physics modeling, prestack inversion and Bayesian classification. Brian Russell

Integrating rock physics modeling, prestack inversion and Bayesian classification. Brian Russell Integrating rock physics modeling, prestack inversion and Bayesian classification Brian Russell Introduction Today, most geoscientists have an array of tools available to perform seismic reservoir characterization.

More information

Amplitude variation with offset AVO. and. Direct Hydrocarbon Indicators DHI. Reflection at vertical incidence. Reflection at oblique incidence

Amplitude variation with offset AVO. and. Direct Hydrocarbon Indicators DHI. Reflection at vertical incidence. Reflection at oblique incidence Amplitude variation with offset AVO and Direct Hydrocarbon Indicators DHI Reflection at vertical incidence Reflection coefficient R(p) c α 1 S wavespeed β 1 density ρ 1 α 2 S wavespeed β 2 density ρ 2

More information

Reservoir connectivity uncertainty from stochastic seismic inversion Rémi Moyen* and Philippe M. Doyen (CGGVeritas)

Reservoir connectivity uncertainty from stochastic seismic inversion Rémi Moyen* and Philippe M. Doyen (CGGVeritas) Rémi Moyen* and Philippe M. Doyen (CGGVeritas) Summary Static reservoir connectivity analysis is sometimes based on 3D facies or geobody models defined by combining well data and inverted seismic impedances.

More information

Pluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO*

Pluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO* Pluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO* *This paper has been submitted to the EAGE for presentation at the June 2001 EAGE meeting. SUMMARY

More information

Reservoir Rock Properties COPYRIGHT. Sources and Seals Porosity and Permeability. This section will cover the following learning objectives:

Reservoir Rock Properties COPYRIGHT. Sources and Seals Porosity and Permeability. This section will cover the following learning objectives: Learning Objectives Reservoir Rock Properties Core Sources and Seals Porosity and Permeability This section will cover the following learning objectives: Explain why petroleum fluids are found in underground

More information

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at Electrical anisotropy drivers in the Snøhvit region of the Barents Sea Michelle Ellis*, Lucy MacGregor, Rolf Ackermann, Paola Newton, Robert Keirstead, Alberto Rusic, Slim Bouchrara, Amanda Geck Alvarez,

More information

AVO inversion and Lateral prediction of reservoir properties of Amangi hydrocarbon field of the Niger Delta area of Nigeria

AVO inversion and Lateral prediction of reservoir properties of Amangi hydrocarbon field of the Niger Delta area of Nigeria IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 2321 0990, p-issn: 2321 0982.Volume 2, Issue 2 Ver. II. (Mar-Apr. 2014), PP 08-17 AVO inversion and Lateral prediction of reservoir properties

More information

Post-stack inversion of the Hussar low frequency seismic data

Post-stack inversion of the Hussar low frequency seismic data Inversion of the Hussar low frequency seismic data Post-stack inversion of the Hussar low frequency seismic data Patricia E. Gavotti, Don C. Lawton, Gary F. Margrave and J. Helen Isaac ABSTRACT The Hussar

More information

Competing Effect of Pore Fluid and Texture -- Case Study

Competing Effect of Pore Fluid and Texture -- Case Study Competing Effect of Pore Fluid and Texture -- Case Study Depth (m) Sw Sxo. m Poisson's Ratio.. GOC.1 5 7 8 9 P-Impedance OWC 15 GR.. RHOB.5 1 Saturation...5. 1. 1.5 Vs (km/s).. Poisson's Ratio 5 7 P-Impedance

More information

Lithology prediction and fluid discrimination in Block A6 offshore Myanmar

Lithology prediction and fluid discrimination in Block A6 offshore Myanmar 10 th Biennial International Conference & Exposition P 141 Lithology prediction and fluid discrimination in Block A6 offshore Myanmar Hanumantha Rao. Y *, Loic Michel, Hampson-Russell, Kyaw Myint, Ko Ko,

More information

Neural Inversion Technology for reservoir property prediction from seismic data

Neural Inversion Technology for reservoir property prediction from seismic data Original article published in Russian in Nefteservice, March 2009 Neural Inversion Technology for reservoir property prediction from seismic data Malyarova Tatyana, Kopenkin Roman, Paradigm At the software

More information

Use of Seismic Inversion Attributes In Field Development Planning

Use of Seismic Inversion Attributes In Field Development Planning IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 2321 0990, p-issn: 2321 0982.Volume 6, Issue 2 Ver. II (Mar. Apr. 2018), PP 86-92 www.iosrjournals.org Use of Seismic Inversion Attributes

More information

AVO Crossplotting II: Examining Vp/Vs Behavior

AVO Crossplotting II: Examining Vp/Vs Behavior AVO Crossplotting II: Examining Vp/Vs Behavior Heath Pelletier* Talisman Energy, Calgary, AB hpelletier@talisman-energy.com Introduction The development of AVO crossplot analysis has been the subject of

More information

Estimation of density from seismic data without long offsets a novel approach.

Estimation of density from seismic data without long offsets a novel approach. Estimation of density from seismic data without long offsets a novel approach. Ritesh Kumar Sharma* and Satinder Chopra Arcis seismic solutions, TGS, Calgary Summary Estimation of density plays an important

More information

Pre-stack (AVO) and post-stack inversion of the Hussar low frequency seismic data

Pre-stack (AVO) and post-stack inversion of the Hussar low frequency seismic data Pre-stack (AVO) and post-stack inversion of the Hussar low frequency seismic data A.Nassir Saeed, Gary F. Margrave and Laurence R. Lines ABSTRACT Post-stack and pre-stack (AVO) inversion were performed

More information

We Density/Porosity Versus Velocity of Overconsolidated Sands Derived from Experimental Compaction SUMMARY

We Density/Porosity Versus Velocity of Overconsolidated Sands Derived from Experimental Compaction SUMMARY We 6 Density/Porosity Versus Velocity of Overconsolidated Sands Derived from Experimental Compaction S. Narongsirikul* (University of Oslo), N.H. Mondol (University of Oslo and Norwegian Geotechnical Inst)

More information

We P2 04 Rock Property Volume Estimation Using the Multiattribute Rotation Scheme (MARS) - Case Study in the South Falkland Basin

We P2 04 Rock Property Volume Estimation Using the Multiattribute Rotation Scheme (MARS) - Case Study in the South Falkland Basin We P2 04 Rock Property Volume Estimation Using the Multiattribute Rotation Scheme (MARS) - Case Study in the South Falkland Basin P.K. Alvarez* (Rock Solid Images), B. Farrer (Borders & Southern Petroleum),

More information