--> An 'Elastic Impedance' Approach, by Satinder Chopra and Ritesh Kumar Sharma; #41082 (2012).
[First Hit]

Datapages, Inc.Print this page

Click to view article in PDF format.

 

GCAn 'Previous HitElasticNext Hit Previous HitImpedanceNext Hit' Approach*

 

Satinder Chopra1and Ritesh Kumar Sharma1

 

Search and Discovery Article #41082 (2012)

Posted November 12, 2012

 

*Adapted from the Geophysical Corner column, prepared by the author, in AAPG Explorer, October, 2012. Editor of Geophysical Corner is Satinder Chopra ([email protected]). Managing Editor of AAPG Explorer is Vern Stefanic; Larry Nation is Communications Director. AAPG©2012

 

1 Arcis Corp., Calgary, Canada ([email protected])

 

General Statement

A detailed investigation of seismic amplitudes can yield information pertaining to lithological variation in subsurface sedimentary rock formations and the existence and extent of some hydrocarbon zones. This objective can be facilitated in a process called seismic inversion, which transforms seismic amplitudes into acoustic Previous HitimpedanceNext Hit values. In doing so, the seismic reflection response gets transformed into layered Previous HitimpedanceNext Hit response, which makes the interpretation of the lithological and fluid information more convenient – each transformed Previous HitimpedanceNext Hit trace can now be considered as an Previous HitimpedanceNext Hit log curve and the seismic volume as logs recorded in wells drilled at every seismic trace location. Just as the changes in the character of Previous HitimpedanceNext Hit log curves are indicative of changes in lithology, porosity and fluid content, similar changes seen on inverted Previous HitimpedanceNext Hit traces are interpretable of these properties in a lateral sense over an area and so over a volume.

Acoustic Previous HitimpedanceNext Hit inversion has now become an integral part of most interpretation projects today. While this is a beneficial tool for the seismic interpreter, acoustic Previous HitimpedanceNext Hit inversion is usually run on stacked seismic traces – that is, the individual prestack time migrated offset gathers are stacked and then transformed into Previous HitimpedanceNext Hit. To better exploit the fluid effects that manifest on prestack gathers as variation of amplitudes with offset or angle, prestack Previous HitimpedanceNext Hit inversion also can be carried out. Of course, it would take longer - and so the trade-off is usually between the cost, time and the method to be used.

General statement

Figures

Examples

Conclusion

Acknowledgments

 

 

 

 

 

 

 

 

 

General statement

Figures

Examples

Conclusion

Acknowledgments

 

 

 

 

 

 

 

 

 

General statement

Figures

Examples

Conclusion

Acknowledgments

 

 

 

 

 

 

 

 

 

General statement

Figures

Examples

Conclusion

Acknowledgments

Figure Captions

Figure 1. Display of P-wave, S-wave, Density and Gamma-Ray log curves from a well in the Magdelena Valley, Colombia. The low values for all the curve displays are to the left and high values to the right. To the far right is shown a comparison of an Acoustic Previous HitImpedanceNext Hit (AI) curve with the computed Previous HitElasticNext Hit Previous HitImpedanceNext Hit (EI) (30°) curve. Notice the decrease in Previous HitimpedanceNext Hit (deviation in the blue and red curves) at the gas-producing zone as indicated with the purple arrow. (Data courtesy: PetroNorte, Colombia).

Figure 2. Synthetic seismogram tie for the low-angle (~10°) near-stack. The synthetic seismogram was generated by using the Previous HitimpedanceNext Hit log curve. The correlation seems to be reasonably good. (Data courtesy: PetroNorte, Colombia).

Figure 3. Synthetic seismogram tie for the high-angle (~30°) far-stack. The synthetic seismogram was generated by using the computed EI (30°) log curve. The correlation seems to be reasonably good. Notice the weakening of the amplitudes at the location of the orange arrows. (Data courtesy: PetroNorte, Colombia).

Figure 4. (Left) Segment of the acoustic Previous HitimpedanceNext Hit section (with the overlaid AI log) shows anomalously low values of Previous HitimpedanceNext Hit at the gas-producing zone (in yellow highlighted zone). (Right) Equivalent segment from EI (30°) section showing the anomaly as much more pronounced. (Data courtesy: PetroNorte, Colombia).

Examples

A simple way to examine the variations of amplitude as a function of offset is to generate the offset-limited seismic volumes, such as the near-, mid- and far-offset (or angle) volumes. Variations seen on these volumes in desired zones could then be indicative of the fluid information.

For example, a low-Previous HitimpedanceNext Hit gas-sand sandwiched between shale would yield an increase of amplitude with offset. Such a variation can be detected on comparing the near-offset seismic volume with the equivalent far-offset volume, and noticing high amplitude anomalies on the latter corresponding to the gas samples.

As amplitudes of the near-offset traces are related to the changes in acoustic Previous HitimpedanceNext Hit they can be calibrated with well log curves or synthetic seismograms. However, if a far-offset or a far-angle stack has to be calibrated with the log data or synthetic seismograms, there is no analogous set of log curves that could be used for the purpose.

Back in 1999, Patrick Connolly from BP pointed this out and suggested the generalization of acoustic Previous HitimpedanceNext Hit for variable incidence angle using a linearized version of Zoeppritz equations. He called this "Previous HitelasticNext Hit Previous HitimpedanceNext Hit," and it provides the framework to calibrate and invert non-zero offset data.

The Previous HitelasticNext Hit Previous HitimpedanceNext Hit approach is strongly dependent on the medium parameters (VP, VS, and density and angle of incidence), and so is often regarded as the rock attribute analog of acoustic Previous HitimpedanceNext Hit for varying angles of incidence.

In actual practice, the CMP gather at the position of the well is picked up, different angle ranges are selected and angle stacks generated. Given the VP, VS and density log curves, the Previous HitelasticNext Hit Previous HitimpedanceNext Hit is calculated for different angles of incidence. The angle stack traces from the gather and those derived from the log curves (Previous HitelasticNext Hit Previous HitimpedanceNext Hit, or "EI") are compared for a visual assessment and interpretation.

Another useful and meaningful display is the comparison of the acoustic Previous HitimpedanceNext Hit log curve with the Previous HitelasticNext Hit Previous HitimpedanceNext Hit curve at the far-angle that is admissible for the given data. In Figure 1, the acoustic Previous HitimpedanceNext Hit log is compared with the EI (30 degree) log curve for a discovery well from Colombia. The target is related to Eocene fluvial deposits, mainly composed of interbedded medium- to fine-grained quartz sandstones and claystones.

The gas was detected during mud logging and on the electric log curve; however, the density and neutron curve crossover is not as high as expected, probably due to low saturation as well as its position. The saturation is expected to increase in the up-dip direction. Notice that there is a decrease of Previous HitimpedanceNext Hit at the gas-sand interface, and so it will show up as higher amplitudes on the seismic data.

It may be mentioned that the Previous HitelasticNext Hit Previous HitimpedanceNext Hit values vary significantly with the incidence angle – and because of this, when Previous HitelasticNext Hit Previous HitimpedanceNext Hit logs have to be displayed alongside with or overlaid on acoustic Previous HitimpedanceNext Hit, they have to be scaled in such a way that the EI values for all angles fall in the range of the normal acoustic Previous HitimpedanceNext Hit values. In Figure 1, such a scaling has been applied to EI (30 degrees).

As stated, Previous HitelasticNext Hit Previous HitimpedanceNext Hit also provides a convenient way of producing synthetic seismograms for variable angles of incidence. The computed EI (30 degree) log curve can be used for producing synthetic seismograms, which could now be used for correlation with far-offset/angle stack.

In Figure 2 we show the correlation of a synthetic seismogram (generated with the acoustic Previous HitimpedanceNext Hit log curve) with the near-stack from the Magdelena Valley, Colombia. The tie seems to be reasonably good. Figure 3 shows a similar synthetic seismogram (generated by using the EI log curve) tie with the far-offset stack. This tie again seems to be good – but, as expected, the amplitudes at the indicated locations on the far-stack are weaker and seem to tie accordingly with the synthetic.

In Figure 4 we show a comparison of a segment of an acoustic Previous HitimpedanceNext Hit section derived from post-stack AI inversion and the equivalent EI (30 degree) section. Notice the differences in the yellow highlighted zones that enclose the gas-producing reservoir. While the hydrocarbon-bearing zone is indicated on the AI section, it appears more pronounced and convincing – and its correlation with the overlaid Previous HitimpedanceNext Hit curve also is much better than the correlation with the AI section.

Conclusion

Thus Previous HitelasticNext Hit Previous HitimpedanceTop attribute serves to combine the benefits of working with inverted data with far-angle data where the fluid information resides.

Acknowledgments

We thank PetroNorte, Colombia, for giving us permission for presentation of the results shown in this study. We also thank Arcis Seismic Solutions for permission to present this work.

Return to top.