AVO for one- and two-fracture set models
|
|
- Linette Hodge
- 5 years ago
- Views:
Transcription
1 GEOPHYSICS, VOL. 70, NO. 2 (MARCH-APRIL 2005); P. C1 C5, 7 FIGS., 3 TABLES / AVO for one- and two-fracture set models He Chen 1,RaymonL.Brown 2, and John P. Castagna 1 ABSTRACT A theoretical comparison is made of PP and PS angle dependent reflection coefficients at the top of two fractured-reservoir models using exact, general, anisotropic reflection coefficients. The two verticalfracture models are taken to have the same total crack density. The primary issue investigated is determination of the fracture orientation using azimuthal AVO analysis. The first model represents a single-fracture set and the second model has an additional fracture set oblique to the first set at an angle of 60. As expected, the PP-wave anisotropy is reduced when multiple fracture sets are present, making the determination of orientation more difficult than for the case of a singlefracture set. Long offsets are required for identification of dominant fracture orientations using PP-wave AVO. PS-wave AVO, however, is quite sensitive to fracture orientations, even at short offsets. For multiple-fracture sets, PS signals can potentially be used to determine orientations of the individual sets. INTRODUCTION The use of reflected P-wave and P-SV-wave amplitudeversus-offset analysis (AVO) can potentially aid in the characterization of fractured reservoirs (Thomsen, 1988). Cores and well logs offer methods of describing fractures at a very fine scale. The seismic method potentially offers a method of characterizing the reservoir at a larger scale that can be used for predicting reservoir properties. Typically, the estimation of fracture properties, especially those for multiple-fracture sets, requires multiple types of data (e.g., Shen et al., 2002; Grechka et al., 2000). A recent study by DeVault et al. (2002) indicates that S-waves have great potential application for fractured reservoirs. The utility of PS-waves is examined here for the particular problem of fracture-set orientation when multiplefracture sets are present. In this paper, the offset-dependent reflectivity for the top of a fractured reservoir is studied with a view toward identifying those approaches best suited for determining fracture-set orientations within the reservoir. We examine the azimuthal AVO variations of two fractured models using exact anisotropic reflection coefficients (Rokhlin et al., 1986). Model 1 is designed containing a single-fracture set, and model 2 consists of two-fracture sets. The singlefracture-set model exhibits a stronger azimuthal AVO variation. The results of the modeling study illustrate the potential problems and approaches that can be used in the interpretation of multiple-fracture sets. SYNTHETIC FRACTURED-RESERVOIR MODELS Two models of fractured reservoirs are set up for synthetic azimuthal AVO studies in this paper. The models studied are simple and limited to a constant fracture density so that the effects of fracture orientation could be studied. See the paper by Liu et al. (2000) for a discussion of multiple-fracture sets. The matrix material of the reservoirs is assumed to be limestone with 10% porosity, and the overlying seal layer is assumed to be shale. Both the incident shale and background (matrix) for the fractured reservoir layers are assumed to be isotropic. Only the presence of fractures contributes to the anisotropy of the models examined. Table 1 lists the petrophysical parameters for the isotropic portions of the models. All of these parameters are set within the range of the laboratory and empirical relationships between V p, V s,andρ described by Castagna et al. (1993) for the rocks considered. Model 1 is assumed to have one set of vertical fractures whose planes have an azimuth measured from north of 0.The X 2 -axis points north, and the X 1 -axis points east. A positive azimuth angle is measured toward the east from the X 2 -axis. This means that the normals to the first set of fractures are pointed parallel to the X 1 -axis (Figure 1b). Model 2 has two sets of fractures whose planes are at azimuths of 0 and 60 (Figure 1c). The crack density of the fracture set in model 1 is set equal to 0.1. For the purpose of comparison between models, the total crack density in model 2 is also set as 0.1, with the Manuscript received by the Editor August 25, 2002; revised manuscript received August 8, 2004; published online March 22, University of Oklahoma, Institute of Exploration and Development Geosciences, Norman, Oklahoma che@gcn.ou.edu; jcastagna@ou.edu. 2 Oklahoma Geological Survey, 100 E. Boyd Street, Norman, Oklahoma raybrown@ou.edu. c 2005 Society of Exploration Geophysicists. All rights reserved. C1
2 C2 Chen et al. crack density of each fracture set equal to In this way, the total crack density is the same for models 1 and 2, while their anisotropy differs. Table 2 lists the fracture parameters described above. The stiffness matrices (Table 3) are calculated assuming that the excess compliance caused by each fracture set can Table 1. Model matrix information. Reservoir matrix Overlying shale V p (m/s) V s (m/s) Density (g/cm 3 ) Porosity 10% simply be added to find the change in compliance caused by the presence of both fracture sets. This assumption neglects any interaction between the fractures. As a first-order approximation, we use the scalar cracks described by Schoenberg and Sayers (1995) to estimate the compliance for each fracture set. Using scalar cracks simplifies the stiffness tensor for the models, but it is not expected to change the major conclusions for the paper. The reflection-coefficient calculations are made using the method described by Chen (2000) and Rokhlin et al. (1986). First, the horizontal slowness of the incident wave is determined. Next, this horizontal slowness is used to set up two sixth-order polynomials that can be solved to find the reflected and refracted slownesses for the two media across the reflecting boundary. The eigenvectors for each of the reflected and refracted waves are determined next. Finally, the boundary conditions are set up, and the reflection and refraction coefficients are determined. According to Hudson (1980, 1981), one vertical set of fractures can introduce transverse isotropy with a horizontal symmetry axis (HTI), while multiple sets of fractures can cause arbitrary anisotropy. An example of what happens with multiple fracture sets can be examined in two simple models discussed below. The reader is referred to the inversion work by Bakulin et al. (2002) for a study of two orthogonal-fracture sets. An incident P-wave is assumed for each model. The vertical planes of incidence for the P-waves will be described by the azimuths of the planes. These planes will all contain the X 3 -axis, and the azimuth is measured positively in a clockwise movement from north (Figure 1a) (the same direction as the X 2 -axis). The three expected reflections are a quasi-p (PP)- wave, a quasi-in-plane S (PSI)-wave and a quasi-out-of-plane S-wave. The incident phase and/or ray (group) angle (both are the same for incident medium, which is isotropic) range for the study varies from 0 to 45. This angle is measured from the X 3 -axis in the plane of incidence. The variation of reflection amplitudes with angle of incidence is observed at four different azimuth directions: 0, 30, 60,and 90, respectively. In other words, the plane of incidence is varied, although the fracture geometry is the same; i.e., the stiffness tensor remains the same while the azimuthal angle changes. AZIMUTHAL INTERPRETATION OF AVO FROM TOP OF RESERVOIR The first wave mode to be considered is a downgoing P-wave reflected from the reservoir as a P-wave (PP). Figures 2 and 3 show the exact azimuthal PP reflection amplitudes from models 1 (single-fracture set) and 2 (two-fracture Table 2. Model fracture information. Figure 1. (a) Azimuthal direction of incident plane, described by angle φ. (b)x 1 X 2 -plane sketch of fracture orientation in model 1. (c) X 1 X 2 -plane sketch of fracture orientations in model 2. Crack density Azimuthal angle Dip angle Model 1 Fracture set no Model 2 Fracture set no Fracture set no
3 Fractured AVO C3 sets), respectively. First, there are only slight differences on the normal-incidence PP reflection coefficients, showing that the total crack density may be the primary factor related to fractures that influence the zero-offset PP reflections. Secondly, the azimuthal variations of the PP reflection can only be clearly observed within each curve family when the incidence angles exceed 25. This is why long offsets are expected to be required to detect vertical-fracture orientations with P-waves. Note that model 1 (Figure 2), with a single fracture set, shows more azimuthal AVO than model 2 (Figure 3). Note also that, in Figure 3, the azimuth 90 curve is identical to the azimuth 30 curve. For these two situations, the incidence plane happens to be perpendicular to one of the two fracture sets. As a result, the P-wave sees the same reflection response at these azimuths, because the two fracture sets have the same crack density. The next reflected phase to be considered is the incident P-wave, which reflects from the reservoir as an S-wave (PS). The polarization of S-waves for anisotropic media can be expressed in many ways. To simplify the description, we describe Table 3. Stiffness matrix information. C ij (GPa) Isotropic background: Single-fracture set (model 1): Two-fracture sets (model 2): Figure 2. Azimuthal P-wave reflection coefficients for model 1. Figure 3. Azimuthal P-wave reflection coefficients for model 2.
4 C4 Chen et al. the S-wave polarization as in-plane (near-sv) and out-ofplane (near-sh). In this way, the reader gets an intuitive feeling for the polarization without actually plotting the polarization or having to keep track of fracture alignments. Figures 4 and 5 show the out-of-plane shear reflection (or near-sh reflection) amplitudes from models 1 and 2, respectively. This mode is unique to the anisotropic circumstances and is caused by coupling between the in-plane (near-sv) and out-of-plane (near-sh) components of motion at the interface. It can be seen that: 1) The azimuthal variation of the converted PS reflections is much more sensitive to fractures than PP reflections (Figures 2 and 3). Even at small incidence angles, Figures 4 and 5 show a distinct azimuthal AVO response for both models 1 and 2. 2) The normal incidence reflections are always zero for all quasi-ps-waves, because all the fracture sets are vertical. There is no converted shear wave when the P-wave is normally incident upon the reflection surface for the model assumed. 3) In Figure 4, for model 1, the azimuth 90 and 0 AVO curves are zero. With only one set of vertical fractures at azimuth 90, the reservoir model is an HTI medium with a single vertical-symmetry plane perpendicular to the fractures (Tsvankin, 1995, 1996). One of the symmetry planes is the isotropic plane at azimuth 0 (X 2 X 3 plane) and the other is an effective VTI plane at azimuth 90 direction (X 1 X 3 plane) (Tsvankin, 1995, 1996). These models show that the converted shear-wave reflections can be used to find the orientation of the fracture set if only one fracture set exists. 4) In Figure 5, for model 2, a zero-reflection coefficient curve is observed at an azimuth of 60, which equally splits the two fracture systems. The two reflection-coefficient curves for azimuths 90 and 30 show up as mirror images of each other reflected about the azimuth 60 curve. These observations indicate that the azimuth 60 planeactsasa symmetry plane, but that it is not parallel or perpendicular to any fracture sets. In this two-fracture-set model, the 60 azimuth plane can be misinterpreted as an apparent effective orientation of one fracture set. Figure 4. Azimuthal P-SH reflection coefficients for model 1. Figures 6 and 7 show the quasi-in-plane shear reflections (or near-sv) for models 1 and 2, respectively. Once again, as in the out-of-plane S-wave case, the in-plane S-wave reflection curves show strong azimuthal variations for fracture models under small incidence-angle conditions. Even more importantly, the azimuthal AVO variations are larger than those of out-of-plane S-wave curves (compare to Figures 4 and 5). Another observation is that the azimuthal AVO from the singlefracture-set model (model 1) is larger than those of the twofracture-set model (model 2). Figure 5. Azimuthal P-SH reflection coefficients for model 2. Figure 6. Azimuthal P-SV coefficients for model 1.
5 Fractured AVO C5 If the ideal fracture sets contribute equally to fracture permeability, then one of the symmetry planes identified by PS reflections is in the direction of maximum horizontal permeability, while the other represents a minimum direction of horizontal permeability. If there is no PS reflection at normal incidence, the existing fracture sets can be assumed to be normal to the reflection boundary. ACKNOWLEDGMENTS H. Chen received financial support from the OU Integrated Reservoir Characterization Consortium sponsors. Appreciation is expressed to Bill Lamb for his discussions of the topic. H. Chen and R. Brown acknowledge the support of this work by the U.S. Department of Energy through project DE-AC26-99BC Figure 7. Azimuthal P-SV coefficients for model 2. CONCLUSIONS Exact reflection-coefficient versus angle-of-incidence curves for P- and PS-waves were calculated for models with one or two sets of fractures. We find that: 1) The azimuthal variation of converted PS reflections is more sensitive to fractures than PP reflections. Even at small angles of incidence, the PS reflection coefficients show distinct azimuthal AVO variations. The azimuthal variations of the PP reflection caused by fractures can be clearly observed only when the incidence angle is large (greater than 25 for the models in this synthetic study). 2) For the same total crack density, a model containing a single-fracture set exhibits the largest variation of azimuthal AVO. Multiple-fracture sets tend to weaken the azimuthal anisotropic effects (especially for PP-waves). 3) With large-offset data, the azimuthal variation in P-wave reflections may indicate the existence of fracture zones. For larger relative-amplitude variations with azimuth, there is a greater possibility for a single-fracture set or a dominant fracture orientation. 4) The out-of-plane/in-plane PS-wave azimuthal AVOs can be used, when they can be observed, to identify symmetry planes as well as the orientation of individual fracture sets. This offers a powerful tool for understanding the fracture makeup of a reservoir. REFERENCES Bakulin, A., V. Grechka, and I. Tsvankin, 2000, Estimation of fracture parameters from reflection seismic data, Part III: Fractured models with monoclinic symmetry: Geophysics, 65, , 2002, Seismic inversion for the parameters of two orthogonal fracture sets in a VTI background medium: Geophysics, 67, Castagna, J. P., M. L. Batzle, J. E. Daiser, K. M. Tubman, and M. L. Burnett, 1993, Shearwave velocity control, in J. P. Castagna and G. E. Backus, eds., Offset-dependent reflectivity: Theory and practice of AVO analysis: SEG. Chen, H., 2000, Anisotropic effects upon amplitude-vs-offset response in realistic earth models: Ph.D. dissertation, University of Oklahoma. DeVault, B., T. L. Davis, I. Tsvankin, R. Verm, and F. Hilterman, 2002, Multicomponent AVO analysis, Vacuum Field, New Mexico: Geophysics, 67, Grechka, V., P. Contreras, and I. Tsvankin, 2000, Inversion of normal moveout for monoclinic media: Geophysical Prospecting, 48, Hudson, J. A., 1980, Overall properties of a cracked solid: Mathematical Proceedings of the Cambridge Philosophical Society, 88, , 1981, Wave speeds and attenuation of elastic waves in material containing cracks: Geophysical Journal of the Royal Astronomical Society, 64, Liu, E., J. A. Hudson, and T. Pointer, 2000, Equivalent medium representation of fractured rock: Journal of Geophysical Research, 105, Rokhlin, S. I., T. K. Bolland, and L. Adler, 1986, Reflection and refraction of elastic waves on a plane interface between two generally anisotropic media: Journal of the Acoustic Society of America, 79, no. 4, Schoenberg, M., and C. Sayers, 1995, Seismic anisotropy of fractured rock: Geophysics, 60, Shen, F., X. Zhu, and M. N. Toksoz, 2002, Effects of fractures on NMO velocities and P-wave azimuthal AVO response: Geophysics, 67, Thomsen, L., 1988, Reflection seismology over azimuthally anisotropic media: Geophysics, 53, Tsvankin, I., 1995, Body-wave radiation patterns and AVO in transversely isotropic media: Geophysics, 60, , 1996, P-wave signatures and notation for transversely isotropic media: An overview: Geophysics, 61,
HTI anisotropy in heterogeneous elastic model and homogeneous equivalent model
HTI anisotropy in heterogeneous elastic model and homogeneous equivalent model Sitamai, W. Ajiduah*, Gary, F. Margrave and Pat, F. Daley CREWES/University of Calgary. Summary In this study, we compared
More informationAzimuthal Anisotropy Investigations for P and S Waves: A Physical Modelling Experiment
zimuthal nisotropy Investigations for P and S Waves: Physical Modelling Experiment Khaled l Dulaijan, Gary F. Margrave, and Joe Wong CREWES Summary two-layer model was built using vertically laminated
More informationJoint Anisotropic Parameter Inversion Based on Nonhyperbolic Moveout and Azimuthal Amplitude Variance
URTe: 2455484 Joint Anisotropic Parameter Inversion Based on Nonhyperbolic Moveout and Azimuthal Amplitude Variance Fangyu Li*, Jie Qi, Bin Lyu, Kurt J. Marfurt, University of Oklahoma. opyright 2 16,
More informationAnisotropic 3D Amplitude Variation with Azimuth (AVAZ) Methods to Detect Fracture-Prone Zones in Tight Gas Resource Plays*
Anisotropic 3D Amplitude Variation with Azimuth (AVAZ) Methods to Detect Fracture-Prone Zones in Tight Gas Resource Plays* Bill Goodway 1, John Varsek 1, and Christian Abaco 1 Search and Discovery Article
More informationUNDERSTANDING FRACTURE ORIENTATION BY REMOVING POLARIZATION DISTORTION FOR DIRECT SHEAR WAVES. Terence Campbell
UNDERSTANDING FRACTURE ORIENTATION BY REMOVING POLARIZATION DISTORTION FOR DIRECT SHEAR WAVES Terence Campbell Department of Geological Sciences The University of Texas at Austin Abstract The progressive
More informationAnisotropic 3D Amplitude Variation with Azimuth (AVAZ) Methods to Detect Fracture Prone Zones in Tight Gas Resource Plays
Antropic 3D Amplitude Variation with Azimuth (AVAZ) Methods to Detect Fracture Prone Zones in Tight Gas Resource Plays Bill Goodway* EnCana Corp., Calgary, AB, Canada william.goodway@encana.com and John
More informationSummary. Introduction
Chris Davison*, Andrew Ratcliffe, Sergio Grion (CGGeritas), Rodney Johnston, Carlos Duque, Jeremy Neep, Musa Maharramov (BP). Summary Azimuthal velocity models for HTI (Horizontal Transverse Isotropy)
More informationWe B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study
We B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study E. Zabihi Naeini* (Ikon Science), N. Kamath (Colorado School of Mines), I. Tsvankin (Colorado School of Mines), T. Alkhalifah
More informationMC 1.2. SEG/Houston 2005 Annual Meeting 877
PS-wave moveout inversion for tilted TI media: A physical-modeling study Pawan Dewangan and Ilya Tsvankin, Center for Wave Phenomena, Colorado School of Mines (CSM), Mike Batzle, Center for Rock Abuse,
More informationFermat s principle and ray tracing in anisotropic layered media
Ray tracing in anisotropic layers Fermat s principle and ray tracing in anisotropic layered media Joe Wong ABSTRACT I consider the path followed by a seismic signal travelling through velocity models consisting
More informationComparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media
Raytracing in VTI media Comparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media Jianli Yang and Don C. Lawton ABSTRACT Determination of the conversion
More informationAzimuthal Fourier Coefficient Elastic Inversion
Azimuthal Fourier Coefficient Elastic Inversion Benjamin Roure*, Hampson-Russell Software & Services, a CGGVeritas Company, Calgary, Canada Benjamin.Roure@cggveritas.com and Jon Downton, Hampson-Russell
More informationElastic physical seismic modeling Don C. Lawton,ScottP. Cheadle, EricV. Gallantand MalcolmB. Bertram
Elastic physical seismic modeling Don C. Lawton,ScottP. Cheadle, EricV. Gallantand MalcolmB. Bertram ABSTRACT A physieal seismic modeling system for studying elastic waves in solid models has been established.
More informationEstimating interval shear-wave splitting from multicomponent virtual shear check shots
GEOPHYSICS, VOL. 73, NO. 5 SEPTEMBER-OCTOBER 2008 ; P. A39 A43, 5 FIGS. 10.1190/1.2952356 Estimating interval shear-wave splitting from multicomponent virtual shear check shots Andrey Bakulin 1 and Albena
More informationSUMMARY INTRODUCTION METHOD: AVOZ WITH SVD FOR FRACTURE INVER- SION
Fracture density inversion from a physical geological model using azimuthal AVO with optimal basis functions Isabel Varela, University of Edinburgh, Sonja Maultzsch, TOTAL E&P, Mark Chapman and Xiang-Yang
More informationMinimizing Fracture Characterization Uncertainties Using Full Azimuth Imaging in Local Angle Domain
P-237 Minimizing Fracture Characterization Uncertainties Using Full Azimuth Imaging in Local Angle Domain Shiv Pujan Singh*, Duane Dopkin, Paradigm Geophysical Summary Shale plays are naturally heterogeneous
More informationNumerical modeling of shear-wave splitting and azimuthal velocity analysis in fractured media
Shear-wave splitting analysis Numerical modeling of shear-wave splitting and azimuthal velocity analysis in fractured media Zimin Zhang, Don C. Lawton, and Robert R. Stewart ABSTRACT This report presents
More informationFoolproof AvO. Abstract
Foolproof AvO Dr. Ron Masters, Geoscience Advisor, Headwave, Inc Copyright 2013, The European Association of Geoscientists and Engineers This paper was prepared for presentation during the 75 th EAGE Conference
More informationComparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media
Comparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media Jianli Yang* and Don C. Lawton CREWES, University of Calgary, 5 UniversityDR.NW, Calgary, AB, TN
More informationFull-waveform inversion for reservoir characterization: A synthetic study
CWP-889 Full-waveform inversion for reservoir characterization: A synthetic study Nishant Kamath, Ilya Tsvankin & Ehsan Zabihi Naeini ABSTRACT Most current reservoir-characterization workflows are based
More informationFractured Reservoir Characterization using Azimuthal AVO
Fractured Reservoir Characterization using Azimuthal AVO Burke J. Minsley, Daniel R. Burns, and Mark E. Willis Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary Sciences Massachusetts
More informationTesseral 2.5D-3C Modeling
Tesseral 2.5D-3C Modeling Outline Basic theory and validation tests for different models Basic theory of 2.5D-3C modeling 2-Layer model with the lower HTI layer 2-layer model with both HTI layers Multi-layered
More informationEfficient Double-Beam Characterization for Fractured Reservoir. Yingcai Zheng, Xinding Fang, Michael C. Fehler and Daniel R. Burns
Efficient Double-Beam Characterization for Fractured Reservoir Yingcai Zheng, Xinding Fang, Michael C. Fehler and Daniel R. Burns Abstract We proposed an efficient target-oriented method to characterize
More informationFracture Quality from Integrating Time-Lapse VSP and Microseismic Data
Fracture Quality from Integrating Time-Lapse VSP and Microseismic Data Mark E. Willis, Daniel R. Burns, Rongrong Lu, M. Nafi Toksöz, Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary
More informationm=[a,b,c,d] T together with the a posteriori covariance
zimuthal VO analysis: stabilizing the model parameters Chris Davison*, ndrew Ratcliffe, Sergio Grion (CGGVeritas), Rodney Johnston, Carlos Duque, Musa Maharramov (BP). solved using linear least squares
More informationThe azimuth-dependent offset-midpoint traveltime pyramid in 3D HTI media
The azimuth-dependent offset-midpoint traeltime pyramid in 3D HTI media Item Type Conference Paper Authors Hao Qi; Stoas Alexey; Alkhalifah Tariq Ali Eprint ersion Publisher's Version/PDF DOI.9/segam3-58.
More informationMulticomponent wide-azimuth seismic data for defining fractured reservoirs
Multicomponent wide-azimuth seismic data for defining fractured reservoirs Evaluating and exploiting azimuthal anisotropy Data Processing Figure 1 A typical surface outcrop showing aligned fractures Figure
More informationRobustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
More informationCOMPUTATIONALLY EFFICIENT RAY TRACING ALGORITHM FOR SIMULATION OF TRANSDUCER FIELDS IN ANISOTROPIC MATERIALS
Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation NDE 2011, December 8-10, 2011 COMPUTATIONALLY EFFICIENT RAY TRACING ALGORITHM FOR SIMULATION OF TRANSDUCER FIELDS IN ANISOTROPIC
More informationAnisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS
Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS Summary Anisotropic depth model building using surface seismic data alone is non-unique and
More informationPrestack Kirchhoff time migration for complex media
Stanford Exploration Project, Report 97, July 8, 998, pages 45 6 Prestack Kirchhoff time migration for complex media Tariq Alkhalifah keywords: time migration, anisotropy ABSTRACT Constructing the seismic
More informationObstacles in the analysis of azimuth information from prestack seismic data Anat Canning* and Alex Malkin, Paradigm.
Obstacles in the analysis of azimuth information from prestack seismic data Anat Canning* and Alex Malkin, Paradigm. Summary The azimuth information derived from prestack seismic data at target layers
More informationInversion of Fracture Parameters by Using the Artificial Neural Network
IC -AI 99 International Conference 57 Inversion of Fracture Parameters by Using the Artificial Neural Network Hengchang Dai, Xiang-Yang Li, and Colin MacBeth British Geological Survey, Murchison House,
More informationNon-hyperbolic moveout inversion of wide-azimuth P-wave data for orthorhombic media
Geophysical Prospecting, 2006, 54, 535 552 Non-hyperbolic moveout inversion of wide-azimuth P-wave data for orthorhombic media Ivan Vasconcelos and Ilya Tsvankin Center for Wave Phenomena, Department of
More informationSUMMARY INTRODUCTION THEORY. The anisotropic elastic wave equation is
Passive wavefield imaging using the energy norm Daniel Rocha and Paul Sava, Center for Wave Phenomena, Colorado School of Mines, Jeffrey Shragge, Centre for Energy Geoscience, The University of Western
More informationSeparation of diffracted waves in TI media
CWP-829 Separation of diffracted waves in TI media Yogesh Arora & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines Key words: Diffractions, Kirchhoff, Specularity, Anisotropy ABSTRACT
More informationOn the Scattering Effect of Lateral Discontinuities on AVA Migration
On the Scattering Effect of Lateral Discontinuities on AVA Migration Juefu Wang* and Mauricio D. Sacchi Department of Physics, University of Alberta, 4 Avadh Bhatia Physics Laboratory, Edmonton, AB, T6G
More informationRotation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J. Fowler, WesternGeco.
otation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J Fowler, WesternGeco Summary Symmetry groups commonly used to describe seismic anisotropy
More informationFracture Detection using Amplitude versus Offset and Azimuth Analysis of a 3D P-wave Seismic Dataset and Synthetic Examples
Fracture Detection using Amplitude versus Offset and Azimuth Analysis of a 3D P-wave Seismic Dataset and Synthetic Examples Burke J. Minsley, Mark E. Willis, Mary Krasovec, Daniel R. Burns, and M. Nafi
More informationSeismic physical modeling II: VVAZ and AVAZ effects observed on reflections from isolated HTI targets
Seismic physical modeling II: VVAZ and AVAZ effects observed on reflections from isolated HTI targets Joe Wong ABSTRACT Seismic physical modeling was used to investigate VVAZ and AVAZ phenomena associated
More informationPS wave AVO aspects on processing, inversion, and interpretation
PS wave AVO aspects on processing, inversion, and interpretation Yong Xu, Paradigm Geophysical Canada Summary Using PS wave AVO inversion, density contrasts can be obtained from converted wave data. The
More informationCross-validation of 3C and 9C prestack shear wave splitting analysis at Big Sky
Cross-validation of 3C and 9C prestack shear wave splitting analysis at Big Sky Jeff P. Grossman* and Gulia Popov, Sensor Geophysical, a Global Geophysical Company Summary We present a case study (reproduced
More informationEarthStudy 360. Full-Azimuth Angle Domain Imaging and Analysis
EarthStudy 360 Full-Azimuth Angle Domain Imaging and Analysis 1 EarthStudy 360 A New World of Information for Geoscientists Expanding the Frontiers of Subsurface Exploration Paradigm EarthStudy 360 is
More informationSUMMARY. predicted by Rüger s equation (Rüger, 2001) using the known elastic parameters of the model.
Azimuthal AVO over a simulated fractured medium: A physical modeling experiment Faranak Mahmoudian, Joe Wong, and Gary F. Margrave, University of Calgary, CREWES SUMMARY A seismic physical model experiment
More informationSeismic modelling with the reflectivity method
Yongwang Ma, Luiz Loures, and Gary F. Margrave ABSTRACT Seismic modelling with the reflectivity method Numerical seismic modelling is a powerful tool in seismic imaging, interpretation and inversion. Wave
More informationAVO Analysis with Multi-Offset VSP Data
AVO Analysis with Multi-Offset VSP Data Z. Li*, W. Qian, B. Milkereit and E. Adam University of Toronto, Dept. of Physics, Toronto, ON, M5S 2J8 zli@physics.utoronto.ca T. Bohlen Kiel University, Geosciences,
More information( ) characterize fractured reservoirs include shear wave splitting (Vetri et al., 2003) and the amplitude-versusangle-and-azimuths.
Seismic characterization of fractured reservoirs using 3D double beams Yingcai Zheng*, Xinding Fang, Michael C. Fehler and Daniel R. Burns, Department of Earth, Atmospheric and Planetary Sciences, Massachusetts
More informationSeismic Data Preconditioning for Improved Reservoir Characterization (Inversion and Fracture Analysis)
Seismic Data Preconditioning for Improved Reservoir Characterization (Inversion and Fracture Analysis) Darren Schmidt *, Hampson-Russell, Calgary, AB, Canada Darren.Schmidt@cggveritas.com Alicia Veronesi,
More informationFLUID SUBSTITUTION AND AVAZ ANALYSIS OF FRACTURED DOMAINS: AN ULTRASONIC EXPERIMENTAL STUDY
FLUID SUBSTITUTION AND AVAZ ANALYSIS OF FRACTURED DOMAINS: AN ULTRASONIC EXPERIMENTAL STUDY A Thesis Presented to the Faculty of the Department of Earth and Atmospheric Sciences University of Houston In
More informationMulti-Azimuth PSDM Processing in the Presence of Orthorhombic Anisotropy -- A Case History Offshore North West Australia Summary
Multi-Azimuth PSDM Processing in the Presence of Orthorhombic Anisotropy -- A Case History Offshore North West Australia Sergey Birdus, Jason Sun*, Wenting Sun and Yi Xie, CGGVeritas, Marco Gazzoli and
More informationwhere f (ω) is the transfer function for fracture scattering. where r 2 (ω) is the reflectivity of the lower reflectors.
A robust method for fracture orientation and density detection from seismic scattered energy Xinding Fang, Mike Fehler, Zhenya Zhu, Tianrun Chen, and Daniel Burns, MIT Earth Resources Laboratory SUMMARY
More informationModeling in Seismic. Tesseral Geo Modeling. Capability to compute 2D-2C and 3D- 3C gathers Enables Geo-scientists to make:
www.tesseral-geo.com Multiparameter Numerical Medium for Seismic Modeling, Planning, Imaging & Interpretation Worldwide Tesseral Geo Modeling Modeling in Seismic Ray-tracing, and especially, finite-difference
More informationcv R z design. In this paper, we discuss three of these new methods developed in the last five years.
Nick Moldoveanu, Robin Fletcher, Anthony Lichnewsky, Darrell Coles, WesternGeco Hugues Djikpesse, Schlumberger Doll Research Summary In recent years new methods and tools were developed in seismic survey
More information2012 Tesseral Technologies
In the past, full-wavefield seismic modeling was not widely utilized in geophysical studies due to the lack of PC computational power. Less computationally intensive and simplified methods of seismic modeling
More informationProcessing of anisotropic data in the τ -p domain: II Common-conversion-point sorting
GEOPHYSICS, VOL. 70, NO. 4 (JULY-AUGUST 2005); P. D29 D36, 6 FIGS., 1 TABLE. 10.1190/1.1993708 Processing of anisotropic data in the τ -p domain: II Common-conversion-point sorting Mirko van der Baan 1
More informationAzimuthal binning for improved fracture delineation Gabriel Perez*, Kurt J. Marfurt and Susan Nissen
Azimuthal binning for improved fracture delineation Gabriel Perez*, Kurt J. Marfurt and Susan issen Abstract We propose an alternate way to define azimuth binning in Kirchhoff prestack migration. This
More informationFourth-order NMO velocity for P-waves in layered orthorhombic media vs. offset-azimuth
Fourt-order NMO velocity for P-waves in layered orrombic media vs. set-azimut Zvi Koren* and Igor Ravve Paradigm Geopysical Summary We derive te fourt-order NMO velocity of compressional waves for a multi-layer
More information3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines
3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose 3D angle decomposition methods from elastic reverse
More informationSEG/New Orleans 2006 Annual Meeting
Accuracy improvement for super-wide angle one-way waves by wavefront reconstruction Ru-Shan Wu* and Xiaofeng Jia, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary To
More informationAVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz
AVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz Summary We propose an alternative approach of AVA analysis based on RTM angle-domain
More informationReflectivity modeling for stratified anelastic media
Reflectivity modeling for stratified anelastic media Peng Cheng and Gary F. Margrave Reflectivity modeling ABSTRACT The reflectivity method is widely used for the computation of synthetic seismograms for
More informationSUMMARY ELASTIC SCALAR IMAGING CONDITION
Robust 3D scalar imaging condition for elastic RTM Yuting Duan, presently at Shell International Exploration and Production Inc., formerly at Center for Wave Phenomena, Colorado School of Mines Paul Sava,
More informationInternational Symposium Non-Destructive Testing in Civil Engineering (NDT-CE) September 15-17, 2015, Berlin, Germany
More Info at Open Access Database www.ndt.net/?id=18355 Effect of Surface Unevenness on In Situ Measurements and Theoretical Simulation in Non-Contact Surface Wave Measurements Using a Rolling Microphone
More informationProcessing converted-wave data in the tau-p domain: rotation toward the source and moveout correction
τ-p domain converted-wave processing Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction Raul Cova and Kris Innanen ABSTRACT The asymmetry of the converted-wave
More informationWe N Fracture and Stress from in Situ, Full-azimuth Seismic Data - Application to a Kuwait Oil Field Dataset
http://www.earthdoc.org/publication/publicationdetails/?publication=80685 1-4 June 2015 IFEMA Madrid We N105 13 Fracture and Stress from in Situ, Full-azimuth Seismic Data - Application to a Kuwait Oil
More informationTrue Amplitude Extraction Technology and Its Application in Pre-Stack Seismic Data
Journal of Earth Science, Vol. 26, No. 4, p. 522 529, August 25 ISSN 674-487X Printed in China DOI:.7/s2583-5-567-2 True Amplitude Extraction Technology and Its Application in Pre-Stack Seismic Data Xiaohui
More information1D internal multiple prediction in a multidimensional world: errors and recommendations
1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate
More informationStrategies for elastic full waveform inversion Espen Birger Raknes and Børge Arntsen, Norwegian University of Science and Technology
Strategies for elastic full waveform inversion Espen Birger Raknes and Børge Arntsen, Norwegian University of Science and Technology SUMMARY Ocean-bottom cables (OBC) have become common in reservoir monitoring
More informationElastic wave mode separation for TTI media
CWP-628 Elastic wave mode separation for TTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes for isotropic elastic wavefields is typically
More information2010 SEG SEG Denver 2010 Annual Meeting
Localized anisotropic tomography with checkshot : Gulf of Mexico case study Andrey Bakulin*, Yangjun (Kevin) Liu, Olga Zdraveva, WesternGeco/Schlumberger Summary Borehole information must be used to build
More information3D S wave statics from direct VSP arrivals Michael J. O'Brien, Allied Geophysics, and Paritosh Singh, Colorado School of Mines
3D S wave statics from direct VSP arrivals Michael J. O'Brien, Allied Geophysics, and Paritosh Singh, Colorado School of Mines Summary We observe that direct S wave arrivals in 3D VSP data contain the
More informationA Short Narrative on the Scope of Work Involved in Data Conditioning and Seismic Reservoir Characterization
A Short Narrative on the Scope of Work Involved in Data Conditioning and Seismic Reservoir Characterization March 18, 1999 M. Turhan (Tury) Taner, Ph.D. Chief Geophysicist Rock Solid Images 2600 South
More informationSeismic Reflection Method
Seismic Reflection Method 1/GPH221L9 I. Introduction and General considerations Seismic reflection is the most widely used geophysical technique. It can be used to derive important details about the geometry
More informationAVO modeling in TESSERAL package
AVO modeling in TESSERAL package AVO modeling can be performed for complex media: isotropic, anisotropic, absorbing, gradient, thin-layered, for plain and curved reflecting boundaries. Package estimates
More informationElastic full waveform Inversion for land walkaway VSP data
Elastic full waveform Inversion for land walkaway VSP data Olga Podgornova, Scott Leaney, arwan Charara, Schlumberger; Eric von Lunen, Nexen Energy ULC Summary Full waveform inversion (FWI) is capable
More informationDownloaded 05/15/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Downloaded 5/15/1 to 184.7... Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ VVAZ vs. AVAZ: practical implementation and comparison of two fracture detection
More informationRaypath and traveltime computations Main for 2D transversely isotropic media with Heading dipping symmetry axes Bing Zhou Authors Stewart Greenhalgh
Exploration Geophysics (2006) 37, 150 159 Raypath and traveltime computations Main for 2D transversely isotropic media with Heading dipping symmetry axes Bing Zhou Authors Stewart Greenhalgh Key Words:
More informationHIGH RESOLUTION STACKING OF SEISMIC DATA. Marcos Ricardo Covre, Tiago Barros and Renato da Rocha Lopes
HIGH RESOLUTION STACKING OF SEISMIC DATA Marcos Ricardo Covre, Tiago Barros and Renato da Rocha Lopes School of Electrical and Computer Engineering, University of Campinas DSPCom Laboratory, Department
More information3-D vertical cable processing using EOM
Carlos Rodriguez-Suarez, John C. Bancroft, Yong Xu and Robert R. Stewart ABSTRACT Three-dimensional seismic data using vertical cables was modeled and processed using equivalent offset migration (EOM),
More informationu = v is the Laplacian and represents the sum of the second order derivatives of the wavefield spatially.
Seismic Elastic Modelling Satinder Chopra, Arcis Corporation, Calgary, Canada Summary Seismic modeling experiments were performed to understand the seismic responses associated with a variety of subsurface
More informationDispersion Curve Extraction (1-D Profile)
PS User Guide Series 2015 Dispersion Curve Extraction (1-D Profile) Working with Sample Data ("Vs1DMASW.dat") Prepared By Choon B. Park, Ph.D. January 2015 Table of Contents Page 1. Overview 2 2. Bounds
More informationWe G High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation
We G103 05 High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation G. Hilburn* (TGS), Y. He (TGS), Z. Yan (TGS) & F. Sherrill (TGS) SUMMARY An approach
More informationSEG/New Orleans 2006 Annual Meeting
3-D tomographic updating with automatic volume-based picking Dimitri Bevc*, Moritz Fliedner, Joel VanderKwaak, 3DGeo Development Inc. Summary Whether refining seismic images to evaluate opportunities in
More informationMulticomponent f-x seismic random noise attenuation via vector autoregressive operators
Multicomponent f-x seismic random noise attenuation via vector autoregressive operators Mostafa Naghizadeh and Mauricio Sacchi ABSTRACT We propose an extension of the traditional frequency-space (f-x)
More informationTransactions, SMiRT-22 San Francisco, California, August 18-23, 2013 Division 5
Transactions, SMiRT-22 San Francisco, California, August 18-23, 2013 LACK OF CORRELATION (INCOHERENCE) MODELING AND EFFECTS FROM REALISTIC 3D, INCLINED, BODY AND SURFACE SEISMIC MOTIONS N. Tafazzoli 1,
More informationShear Wave Seismic Study Comparing 9C3D SV and SH Images with 3C3D C-Wave Images FINAL REPORT
Shear Wave Seismic Study Comparing 9C3D SV and SH Images with 3C3D C-Wave Images FINAL REPORT Reporting Period Start Date: January 1, 2003 Reporting Period End Date: June 30, 2004 John Beecherl, Principal
More informationMulti-component seismic-resolution analysis using finite-difference acquisition modelling*
CSIRO PUBLISHING www.publish.csiro.au/journals/eg Exploration Geophysics, 28, 39, 189 197 Multi-component seismic-resolution analysis using finite-difference acquisition modelling* Shaun Strong 1,2 Steve
More informationA A14 A24 A34 A44 A45 A46
Computation of qp- and qs-wave Rays raveltimes Slowness Vector and Polarization in General Anisotropic Media Junwei Huang 1 * Juan M. Reyes-Montes 1 and R. Paul Young 2 1 Applied Seismology Consultants
More informationSEG/San Antonio 2007 Annual Meeting
Yaofeng He* and Ru-Shan Wu Modeling and Imaging Laboratory, Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, CA, 95064, USA Summary A one-way and one-return boundary
More information10 th Biennial International Conference & Exposition P-080. Wojciech Kobusinski, Geofizyka Torun S.A. Summary. Introduction
10 th Biennial International Conference & Exposition P-080 Improving knowledge about earth subsurface based on anisotropic depth imaging. Case of land, full azimuth seismic survey for unconventional hydrocarbons
More informationAnisotropy-preserving 5D interpolation by hybrid Fourier transform
Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,
More informationWe LHR5 06 Multi-dimensional Seismic Data Decomposition for Improved Diffraction Imaging and High Resolution Interpretation
We LHR5 06 Multi-dimensional Seismic Data Decomposition for Improved Diffraction Imaging and High Resolution Interpretation G. Yelin (Paradigm), B. de Ribet* (Paradigm), Y. Serfaty (Paradigm) & D. Chase
More informationMain components of full-waveform inversion for reservoir characterization
Main components of full-waveform inversion for reservoir characterization Ehsan Zabihi Naeini 1*, Tariq Alkhalifah 2, Ilya Tsvankin 3, Nishant Kamath 3 and Jiubing Cheng 4 Abstract Building a 3D reservoir
More informationPressure and shear wave separation of ocean bottom seismic data
Pressure and shear wave separation of ocean bottom seismic data Ohad Barak ABSTRACT PZ summation is a common method for separating the upgoing wavefield from the downgoing wavefield in data acquired by
More informationEstimation of interval anisotropy parameters using velocity-independent layer stripping
GEOPHYSICS, VOL. 74, NO. 5 SEPTEMBER-OCTOBER 009; P. WB117 WB17, 10 FIGS., 8 TABLES. 10.1190/1.315746 Estimation of interval anisotropy parameters using velocity-independent layer stripping Xiaoxiang Wang
More informationComputational Ocean Acoustics
Computational Ocean Acoustics Ray Tracing Wavenumber Integration Normal Modes Parabolic Equation Wavenumber Integration Range-independent Integral Transform solution Exact depth-dependent solution Global
More informationReflection seismic Method - 2D
Reflection seismic Method - 2D Acoustic Impedance Seismic events Wavelets Convolutional model Resolution Stacking and Signal/Noise Data orders Reading: Sheriff and Geldart, Chapters 6, 8 Acoustic Impedance
More informationConsiderations in 3D depth-specific P-S survey design
Considerations in 3D depth-specific P-S survey design Don C. Lawton and Peter W. Cary 3D P-S survey design ABSTRACT A new sparse-shot design for 3D P-S surveys is introduced. In the sparse shot design
More informationNSE 1.7. SEG/Houston 2005 Annual Meeting 1057
Finite-difference Modeling of High-frequency Rayleigh waves Yixian Xu, China University of Geosciences, Wuhan, 430074, China; Jianghai Xia, and Richard D. Miller, Kansas Geological Survey, The University
More informationEARS2232 Exploration Seismics. Exploration. 2006/2007 Semester 1. Teaching Team. Dr. Sebastian Rost (Module Leader)
EARS3 Exploration Seismics Teaching Team Dr. Sebastian Rost (Module Leader) Dr. Graham Stuart (Seismic Interpretation) 006/007 Semester 1 Ben Dando (Demonstrator) Objectives On completion of this module
More information