Residual migration in VTI media using anisotropy continuation
|
|
- Joel Craig
- 6 years ago
- Views:
Transcription
1 Stanford Exploration Project, Report SERGEY, Noember 9, 2000, pages 671?? Residual migration in VTI media using anisotropy continuation Tariq Alkhalifah Sergey Fomel 1 ABSTRACT We introduce anisotropy continuation as a process which relates changes in seismic images to perturbations in the anisotropic medium parameters. This process is constrained by two kinematic equations, one for perturbations in the normal-moeout NMO) elocity the other for perturbations in the dimensionless anisotropy parameter η. We consider separately the case of post-stack migration show that the kinematic equations in this case can be soled explicitly by conerting them to ordinary differential equations by the method of characteristics. Comparing the results of kinematic computations with synthetic numerical experiments confirms the theoretical accuracy of the method. INTRODUCTION A well-known paradox in seismic imaging is that the detailed information about the subsurface elocity is required before a reliable image can be obtained. In practice, this paradox leads to an iteratie approach to building the image. It looks attractie to relate small changes in elocity parameters to inexpensie operators perturbing the image. This approach has been long known as residual migration. A classic result is the theory of residual post-stack migration Rothman et al., 1985), extended to the prestack case by Etgen 1990). In a recent paper, Fomel 1996) introduced the concept of elocity continuation as the continuous model of the residual migration process. All these results were based on the assumption of the isotropic elocity model. Recently, emphasis has been put on the importance of considering anisotropy its influence on data. Alkhalifah Tsankin 1995) demonstrated that, for TI media with ertical symmetry axis VTI media), just two parameters are sufficient for performing all time-related processing, such as normal moeout NMO) correction including non-hyperbolic moeout correction, if necessary), dip-moeout DMO) correction, prestack poststack time migration. One of these two parameters, the short-spread NMO elocity for a horizontal reflector, is gien by nmo 0) = 1 + 2δ, 1) where is the ertical P-wae elocity, δ is one of Thomsen s anisotropy parameters 1 tariq@sep.stanford.edu,sergey@sep.stanford.edu 671
2 672 Alkhalifah Fomel SEP 94 Thomsen, 1986). Taking h to be the P-wae elocity in the horizontal direction, the other parameter, η, is gien by η h ɛ δ nmo 2 1) = 0) 1 + 2δ, 2) where ɛ is another of Thomsen s parameters. In addition, Alkhalifah 1997) has showed that the dependency on just two parameters becomes exact when the ertical shear wae elocity V S0 ) is set to zero. Setting V S0 = 0 leads to remarkably accurate kinematic representations. It also results in much simpler equations that describe P-wae propagation in VTI media. Throughout this paper, we use these simplified, yet accurate, equations, based on setting V S0 = 0, to derie the continuation equations. Because we are only considering time sections, for the sake of simplicity, we denote nmo by. Thus, time processing in VTI media, depends on two parameters η), whereas in isotropic media only counts. In this paper, we generalize the elocity continuation concept to hle VTI media. We define anisotropy continuation as the process of seismic image perturbation when either or η change as migration parameters. This approach is especially attractie, when the initial image is obtained with isotropic migration that is with η = 0). In this case, anisotropy continuation is equialent to introducing anisotropy in the model without the need of repeating the migration step. For the sake of simplicity, we start from the post-stack case purely kinematic description. We define howeer the guidelines for moing to the more complicated interesting cases of prestack migration dynamic equations. The results are preliminary, but they open promising opportunities for seismic data processing in presence of anisotropy. THE GENERAL THEORY In the case of zero-offset reflection, the ray trael distance, l, from the source to the reflection point is related to the two-way zero-offset time, t, by the simple equation l = g t, 3) where g is the half of the group elocity, best expressed in terms of its components, as follows: g = 2 gx gτ. Here gx denotes the horizontal component of group elocity, is the ertical P-wae elocity, gτ is the -normalized ertical component of the group elocity. Under the assumption of zero shear-wae elocity in VTI media, these components hae the following analytic expressions: gx = 2 p x 1 2η + 2η pτ 2 ) η) p x 2 + p τ 2, 4)
3 SEP 94 Anisotropy Continuation 673 gτ = η p x 2 ) p τ η) p x 2 + p τ 2, 5) where p x is the horizontal component of slowness, p τ is the normalized again by the ertical P-wae elocity ) ertical component of slowness. The two components of the slowness ector are related by the following eikonal-type equation Alkhalifah, 1997): p τ = 1 2 p x η p x 2. 6) Equation 6) corresponds to a normalized ersion of the dispersion relation in VTI media. If we consider η as imaging parameters migration elocity migration anisotropy coefficient), the ray length l can be taken as an imaging inariant. This implies that the partial deriaties of l with respect to the imaging parameters are zero. Therefore, l = g t + g l η = g η t + g t = 0, 7) t = 0. 8) η Applying the simple chain rule to equations 7) 8), we obtain t = t, t η = t η, 9) where t = p τ, the two-way ertical traeltime is gien by τ = gτ t. Combining equations 7-9) eliminates the two-way zero-offset time t, which leads to the equations = g τ p τ gτ g, 10) η = g η τ p τ gτ g. 11) After some tedious algebraic manipulation, we can transform equations 5) 6) to the general form = τ F p x,,η), 12)
4 674 Alkhalifah Fomel SEP 94 η = τ F η p x,,η). 13) Since the residual migration is applied to migrated data, with the time axis gien by τ the reflection slope gien by x, instead of t p x, respectiely, we need to eliminate p x from equations 12) 13). This task can be achieed with the help of the following explicit relation, deried in Appendix A, p 2 x = 2τ x η) τ x 2 + S, 14) where τ x = x, S = 8 2 η τ x η) τ x 2 ) 2. Inserting equation 14) into equations 12) 13) yields exact, yet complicated equations, describing the continuation process for η. In summary, these equations hae the form ) = τ f x,,η 15) ) η = τ f η x,,η. 16) Equations of the form 15) 16) contain all the necessary information about the kinematic laws of anisotropy continuation in the domain of zero-offset migration. Linearization A useful approximation of equations 15) 16) can be obtained by simply setting η equal to zero in the right-h side of the equations. Under this approximation, equation 15) leads to the kinematic elocity-continuation equation for elliptically anisotropic media, which has the following relatiely simple form: = τ ) τ x τ x 2 ) τ x 2. 17) It is interesting to note that setting =, yields Fomel s expression for isotropic media Fomel, 1996) gien by = τ τ x 2. 18)
5 SEP 94 Anisotropy Continuation 675 Alkhalifah Tsankin 1995) hae shown that time-domain processing algorithms for elliptically anisotropic media should be the same as those for isotropic media. Howeer, in anisotropic continuation, elliptical anisotropy isotropy differ by a ertical scaling factor that is related to the difference between the ertical NMO elocities. In isotropic media, when elocity is continued, both the ertical NMO elocities which are the same) are continued together, whereas in anisotropic media including elliptically anisotropic) the NMO-elocity continuation is separated from the ertical-elocity one, equation 19) corresponds to continuation only in the NMO elocity. This also implies that equation 19) is more flexible than equation 18), in that we can isolate the ertical-elocity continuation a parameter that is usually ambiguous in surface processing) from the rest of the continuation process. Using τ = z, where z is depth, we immediately obtain the equation = τ, which represents the ertical-elocity continuation. Setting η = 0 = in equation 16) leads to the following kinematic equation for η-continuation: η = τ4 4 τ x τ 2. 19) x We include more discussion about different aspects of linearization in Appendix B. The next section presents the analytic solution of equation 19). Later in this paper, we compare the analytic solution with a numerical synthetic example. ORDINARY DIFFERENTIAL EQUATION REPRESENTATION: ANISOTROPY RAYS According to the classic rules of mathematical physics, the solution of the kinematic equations 15) 16) can be obtained by soling the following system of ordinary differential equations: dx dm = τ f m x, dτ m dm = τ f m m + τ m f m, dτ dm = ττ f m x + τ m, x dτ x dm = τ x f m. 20) Here m sts for either or η, τ x = x, f = m. To trace the η rays, we must first identify the initial alues x 0, τ 0, τ x0, τ m0 from the boundary conditions. The ariables x 0 τ 0 describe the initial position of a reflector in a time-migrated section, τ x0 describes its migrated slope, τ m0 is simply f m 0,τ 0,τ x0 ). Using the exact kinematic expressions for f results in rather complicated representations of the ordinary differential equations. The linearized expressions, on the other h, are simple allow for a straightforward analytical solution.
6 676 Alkhalifah Fomel SEP 94 From kinematics to dynamics The kinematic η-continuation equation 19) corresponds to the following linear fourth-order dynamic equation 4 P t 3 η + 4 P 2 x 2 t η + t4 4 P x 4 = 0, 21) where the t coordinate refers to the ertical traeltime τ, Pt, x,η) is the migrated image, parameterized in the anisotropy parameter η. To find the correspondence between equations 19) 21), it is sufficient to apply a ray-theoretical model of the image Pt, x,η) = Ax,η) f t τx,η)) 22) as a trial solution to 21). Here the surface t = τx,η) is the anisotropy continuation waefront - the image of a reflector for the corresponding alue of η, the function A is the amplitude. Substituting the trial solution into the partial differential equation 21) considering only the terms with the highest asymptotic order those containing the fourth-order deriatie of the waelet f ), we arrie at the kinematic equation 19). The next asymptotic order the third-order deriaties of f ) gies us the linear partial differential equation of the amplitude transport, as follows: τx 2 ) A η + 22 τ x τη 2 2 ττx 2 ) A x + 2 A 2τ x τ xη + τ η τ xx 6 2 ττx 2 τ ) xx = 0. 23) We can see that when the reflector is flat τ x = 0 τ xx = 0), equation 23) reduces to the equality A η = 0, the amplitude remains unchanged for different η. This is of course a reasonable behaior in the case of a flat reflector. It doesn t guarantee though that the amplitudes, defined by 23), behae equally well for dipping cured reflectors. The amplitude behaior may be altered by adding low-order terms to equation 21). According to the ray theory, such terms can influence the amplitude behaior, but do not change the kinematics of the wae propagation. An appropriate initial-alue condition for equation 21) is the result of isotropic migration that corresponds to the η = 0 section in the t, x,η) domain. In practice, the initial-alue problem can be soled by a finite-difference technique. SYNTHETIC TEST Residual post-stack migration operators can be obtained by generating synthetic data for a model consisting of diffractors for gien medium parameters then migrating the same data with different medium parameters. For example, we can generate diffractions for isotropic media migrate those diffractions using an anisotropic migration. The resultant operator
7 SEP 94 Anisotropy Continuation 677 Distance km) Distance km) Time s) 1.0 Time s) Figure 1: Residual post-stack migration operators calculated by soling equation 19), oerlaid aboe synthetic operators. The synthetic operators are obtained by applying TI poststack migration with η = 0.1 left) η = 0.2 right) to three diffractions generated considering isotropic media. The NMO elocity for the modeling migration is 2.0 km/s. anico-impres [NR]
8 678 Alkhalifah Fomel SEP 94 describes the correction needed to transform an isotropically migrated section to an anisotropic one, that is the anisotropic residual migration operator. Figure 1 shows such synthetic operators oerlaid by kinematically calculated operators that were computed with the help of equation 19) the continuation equations for the case of small η). Despite the inherent accuracy of the synthetic operators, they suffer from the lack of aperture in modeling the diffractions, therefore, beyond a certain angle the operators anish. The agreement between the synthetic calculated operators for small angles, especially for the η = 0.1 case, promises reasonable results in future dynamic implementations. CONCLUSIONS We hae extended the concept of elocity continuation in isotropic media to continuations in both the NMO elocity the anisotropy parameter η for VTI media. Despite the fact that we hae considered the simple case of post-stack migration separately, the exact kinematic equations describing the continuation process are anything but simple. Howeer, useful insights into this problem are deduced from linearized approximations of the continuation equations. These insights include the following conclusions: The leading order behaior of the elocity continuation is proportional to τ 2 x, which corresponds to small or moderate dips. The leading order behaior of the η continuation is proportional to τx 4, which corresponds to moderate or steep dips. Both leading terms are independent of the strength of anisotropy η). In practical applications, the initial migrated section is obtained by isotropic migration,, therefore, the residual process is used to correct for anisotropy. Setting η = 0 in the continuation equations for this type of an application is a reasonable approximation, gien that η=0 is the starting point we consider only weak to moderate degrees of anisotropy η 0.1). Numerical experiments with synthetically generated operators confirm this conclusion. Future directions of research will focus on a practical implementation of anisotropy continuation as well as on extending the theory to the case of residual prestack migration. REFERENCES Alkhalifah, T., Tsankin, I., 1995, Velocity analysis for transersely isotropic media: Geophysics, 60, Alkhalifah, T., 1997, Acoustic approximations for processing in transersely isotropic media: submitted to Geophysics. Etgen, J., 1990, Residual prestack migration interal elocity estimation: Ph.D. thesis, Stanford Uniersity.
9 SEP 94 Anisotropy Continuation 679 Fomel, S., 1996, Migration elocity analysis by elocity continuation: SEP 92, Press, W. H., Teukolsky, S. A., Vetterling, W. T., Flannery, B. P., 1992, Numerical recipes, the art of scientific computing: Cambridge Uniersity Press, New York. Rothman, D. H., Lein, S. A., Rocca, F., 1985, Residual migration applications limitations: Geophysics, 50, no. 1, Thomsen, L., 1986, Weak elastic anisotropy: Geophysics, 51, no. 10, APPENDIX A RELATING THE ZERO-OFFSET AND MIGRATION SLOPES The chain rule of differentiation leads to the equality p x = t x = p τ x, where p τ = t. It is conenient to transform equality A-1) to the form x = p x p τ. A-1) A-2) Using the expression for p τ from the main text, we can write equation A-2) as a quadratic polynomial in px 2 as follows apx 4 + bp2 x + c = 0, A-3) where a = 2 2 η, b = x ) η) + 1, c = x )2. Since η can be small as small as zero for isotropic media), we use the following form of solution to the quadratic equation px 2 = 2c b ± b 2 4ac Press et al., 1992). This form does not go to infinity as η approaches 0. We choose the solution with the negatie sign in front of the square root, because this solution complies with the isotropic result when η is equal to zero.
10 680 Alkhalifah Fomel SEP 94 APPENDIX B LINEARIZED APPROXIMATIONS Although the exact expressions might be sufficiently constructie for actual residual migration applications, linearized forms are still useful, because they gie us aluable insights into the problem. The degree of parameter dependency for different reflector dips is one of the most obious insights in the anisotropy continuation problem. Perturbation of a small parameter proides a general mechanism to simplify functions by recasting them into power-series expansion oer a parameter that has small alues. Two ariables can satisfy the small perturbation criterion in this problem: The anisotropy parameter η η << 1) the reflection dip τ x τ x << 1 or p x << 1). Setting η = 0 yields equation 19) for the elocity continuation in elliptical anisotropic media η = 4 4 τ τ x τ 2 x τ )) x τ ) 2 x τ ). B-1) 2 x which represents the case when we initially introduce anisotropy into our model. Because p x the zero-offset slope) is typically lower than τ x the migrated slope), we perform initial expansions in terms of y = p x. Applying the Taylor series expansion of equations 12) 13) in terms of y dropping all terms beyond the fourth power in y, we obtain = τ p 2 x ) 3 4 τ p x ) η) ), B-2) 2 4 η = 4 4 τ p x ). B-3) 2 Although both equations are equal to zero for p x =0, the leading term in the elocity continuation is proportional to px 2, whereas the the leading term in the η continuation is proportional to px 4. As a result the elocity continuation has greater influence at lower angles than the η continuation. It is also interesting to note that both leading terms are independent of the size of anisotropy η). Despite the typically lower alues of p x, expansions in terms of τ x are more important, but less accurate. For small τ x, p x τ x,, therefore, the leading-term behaior of τ x expansions is the same as that of p x As a result, we arrie at the equation = τ ) 2 τ x + 4 τ ) + τ ) + 12η τ ) ) τ x, B-4)
11 SEP 94 Anisotropy Continuation 681 η = 4 τ ) 4 τ x. B-5) 2 Most of the terms in equations B-4) B-5) are functions of the difference between the ertical NMO elocities. Therefore, for simplicity without a loss of generality, we set = keep only the terms up to the eighth power in τ x. The resultant expressions take the form = τ τ x η τ τ x η 4 25η) τ τ x η 5 83η + 144η 2) τ τ x 8 η = 4 τ τ x η) τ τ x η + 156η 2) τ τ x 8. B-6) B-7) Curiously enough, the second term of the η continuation heaily depends on the size of anisotropy 20η). The first term of equation B-6) τx 2 ) is the isotropic term; all other terms in equations B-6) B-7) are induced by the anisotropy.
12 682 SEP 94
An acoustic wave equation for orthorhombic anisotropy
Stanford Exploration Project, Report 98, August 1, 1998, pages 6?? An acoustic wave equation for orthorhombic anisotropy Tariq Alkhalifah 1 keywords: Anisotropy, finite difference, modeling ABSTRACT Using
More informationv v Downloaded 01/11/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
The pseudo-analytical method: application of pseudo-laplacians to acoustic and acoustic anisotropic wae propagation John T. Etgen* and Serre Brandsberg-Dahl Summary We generalize the pseudo-spectral method
More informationThe Importance of Anisotropy for Prestack Imaging
The Importance of Anisotropy for Prestack Imaging Xiaogui Miao*, Daid Wilkinson and Scott Cheadle Veritas DGC Inc., 715-5th Ae. S.W., Calgary, AB, T2P 5A2 Calgary, AB Xiaogui_Miao@eritasdgc.com ABSTRACT
More informationVelocity continuation and the anatomy of residual prestack time migration
GEOPHYSICS, VOL. 68, NO. 5 SEPTEMBER-OCTOBER 003; P. 1650 1661, 11 FIGS. 10.1190/1.160639 Velocity continuation and the anatomy of residual prestack time migration Sergey Fomel ABSTRACT Velocity continuation
More informationTraveltime approximations for inhomogeneous transversely isotropic media with a horizontal symmetry axis
Geophysical Prospecting, 2013, 61, 495 503 doi: 10.1111/j.1365-2478.2012.01067.x Traveltime approximations for inhomogeneous transversely isotropic media with a horizontal symmetry axis Tariq Alkhalifah
More informationG022 Multi-azimuth Seismic Data Imaging in the Presence of Orthorhombic Anisotropy
G0 Multi-azimuth Seismic Data Imaging in the Presence of Orthorhomic Anisotropy Y. Xie* (CGGVeritas), S. Birdus (CGGVeritas), J. Sun (CGGVeritas) & C. Notfors (CGGVeritas) SUMMARY The presence of orthorhomic
More informationStanford Exploration Project, Report 97, July 8, 1998, pages
Stanford Exploration Project, Report 97, July 8, 998, pages 7 6 Stanford Exploration Project, Report 97, July 8, 998, pages 7 The offset-midpoint traveltime pyramid in transversely isotropic media Tariq
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (04) 98 653 66 GJI Seismology doi: 0.093/gji/ggu9 Effective orthorhombic anisotropic models for wavefield extrapolation Wilson Ibanez-Jacome Tariq Alkhalifah
More informationNonhyperbolic Reflection Moveout for Orthorhombic Media
Nonhyperbolic Reflection Moveout for Orthorhombic Media AbdulFattah Al-Dajani and M. Nafi Toksöz Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary Sciences Massachusetts Institute of
More informationAcoustic anisotropic wavefields through perturbation theory
GEOPHYSICS, VOL. 78, NO. 5 (SEPTEMBER-OCTOBER 2013); P. WC41 WC50, 22 FIGS. 10.1190/GEO2012-0391.1 Acoustic anisotropic wavefields through perturbation theory Tariq Alkhalifah 1 ABSTRACT Solving the anisotropic
More informationPrestack exploding reflector modelling and migration for anisotropic media
Geophysical Prospecting, 015, 63, 10 doi: 10.1111/1365-478.1148 Prestack exploding reflector modelling and migration for anisotropic media T. Alkhalifah Physical Sciences and Engineering, King Abdullah
More informationDepth-domain velocity analysis in VTI media using surface P-wave data: Is it feasible?
GEOPHYSICS, VOL. 66, NO. 3 (MAY-JUNE 2001); P. 897 903, 4 FIGS., 1 TABLE. Depth-domain velocity analysis in VTI media using surface P-wave data: Is it feasible? Yves Le Stunff, Vladimir Grechka, and Ilya
More informationStanford Exploration Project, Report 115, May 22, 2004, pages
Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 248 Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 First-order lateral interval velocity estimates without picking
More informationReflection moveout and parameter estimation for horizontal transverse isotropy
GEOPHYSICS, VOL. 62, NO. 2 (MARCH-APRIL 1997); P. 614 629, 7 FIGS. Reflection moveout and parameter estimation for horizontal transverse isotropy Ilya Tsvankin ABSTRACT Transverse isotropy with a horizontal
More informationAn efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models
An efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models Item Type Conference Paper Authors Waheed Umair bin; Alkhalifah Tariq Ali Eprint version Pre-print
More informationTTI Anisotropic Depth Migration: Which Tilt Estimate Should We Use?
TTI Anisotropic Depth Migration: Which Tilt Estimate Should We Use? Francois Audebert* CGG Americas, Calgary, Alberta, Canada faudebert@cgg.com and Volker Dirks CGG Americas, Calgary, Alberta, Canada Abstract
More informationP185 TTI Anisotropic Depth Migration - Which Tilt Estimate Should We Use?
P18 TTI Anisotropic Depth Migration - Which Tilt Estimate Should We Use? F.S. Audebert* (CGG Americas Inc.), A. Pettenati (Ecole Supérieure d' Electricité) & V. Dirks (CGG Americas Inc) SUMMARY We perform
More informationChapter 1. Introduction EARTH MODEL BUILDING
Chapter 1 Introduction Seismic anisotropy in complex earth subsurface has become increasingly important in seismic imaging due to the increasing offset and azimuth in modern seismic data. To account for
More informationA synthetic example of anisotropic P -wave. processing for a model from the Gulf of Mexico
A synthetic example of anisotropic P -wave processing for a model from the Gulf of Mexico Baoniu Han, Tagir Galikeev, Vladimir Grechka, Jérôme Le Rousseau and Ilya Tsvankin Center for Wave Phenomena, Department
More informationAnalysis of image gathers in factorized VTI media
GEOPHYSICS, VOL. 68, NO. 6 NOVEMBER-DECEMBER 2003); P. 206 2025, 3 FIGS. 0.90/.635055 Analysis of image gathers in factorized VTI media Debashish Sarkar and Ilya Tsvankin ABSTRACT Because events in image
More informationE : Ground-penetrating radar (GPR)
Geophysics 3 March 009 E : Ground-penetrating radar (GPR) The EM methods in section D use low frequency signals that trael in the Earth by diffusion. These methods can image resistiity of the Earth on
More informationUNDERSTAND MOTION IN ONE AND TWO DIMENSIONS
SUBAREA I. COMPETENCY 1.0 UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS MECHANICS Skill 1.1 Calculating displacement, aerage elocity, instantaneous elocity, and acceleration in a gien frame of reference
More informationWave-equation tomography for anisotropic parameters
Wave-equation tomography for anisotropic parameters Yunyue (Elita) Li and Biondo Biondi ABSTRACT Anisotropic models are recognized as more realistic representations of the subsurface where complex geological
More informationNMO residual reduction on medium anisotropy in field X using Fomel and Stovas method
Journal of Physics: Conference Series PAPER OPEN ACCESS NMO residual reduction on medium anisotropy in field X using Fomel and Stovas method To cite this article: M S Rosid et al 018 J. Phys.: Conf. Ser.
More information(a) Taking the derivative of the position vector with respect to time, we have, in SI units (m/s),
Chapter 4 Student Solutions Manual. We apply Eq. 4- and Eq. 4-6. (a) Taking the deriatie of the position ector with respect to time, we hae, in SI units (m/s), d ˆ = (i + 4t ˆj + tk) ˆ = 8tˆj + k ˆ. dt
More informationMigration velocity analysis in factorized VTI media
GEOPHYSICS, VOL. 69, NO. 3 (MAY-JUNE 2004); P. 708 718, 17 FIGS. 10.1190/1.1759457 Migration velocity analysis in factorized VTI media Debashish Sarkar and Ilya Tsvankin ABSTRACT One of the main challenges
More informationStatus: Unit 2, Chapter 3
1 Status: Unit, Chapter 3 Vectors and Scalars Addition of Vectors Graphical Methods Subtraction of Vectors, and Multiplication by a Scalar Adding Vectors by Components Unit Vectors Vector Kinematics Projectile
More informationImprovement of stacking image by anisotropic velocity analysis using P-wave seismic data
P-46 Improvement of stacking image by anisotropic velocity analysis using P-wave seismic data Laxmidhar Behera*, and Prakash Khare, National Geophysical Research Institute (NGRI), Hyderabad Summary Anisotropy
More informationStanford Exploration Project, Report SERGEY, November 9, 2000, pages 683??
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 683?? 682 Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 683?? Velocity continuation by spectral methods Sergey
More informationPlane-wave migration in tilted coordinates
Stanford Exploration Project, Report 124, April 4, 2006, pages 1 16 Plane-wave migration in tilted coordinates Guojian Shan and Biondo Biondi ABSTRACT Plane-wave migration in tilted coordinates is powerful
More information3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie Zhang, Wei Zhang, University of Science and Technology of China (USTC)
Downloaded 01/03/14 to 16.01.198.34. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ 3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie
More informationSUPPLEMENTARY MATERIAL. Authors: Alan A. Stocker (1) and Eero P. Simoncelli (2)
SUPPLEMENTARY MATERIAL Authors: Alan A. Stocker () and Eero P. Simoncelli () Affiliations: () Dept. of Psychology, Uniersity of Pennsylania 34 Walnut Street 33C Philadelphia, PA 94-68 U.S.A. () Howard
More informationImplicit 3-D depth migration by wavefield extrapolation with helical boundary conditions
Stanford Exploration Project, Report 97, July 8, 1998, pages 1 13 Implicit 3-D depth migration by wavefield extrapolation with helical boundary conditions James Rickett, Jon Claerbout, and Sergey Fomel
More informationFowler DMO and time migration for transversely isotropic media
GEOPHYSICS, VOL. 61, NO. 3 (MAY-JUNE 1996), P. 835-845, 6 FIGS. Fowler DMO and time migration for transversely isotropic media John Anderson*, Tariq Alkhalifaht, and Ilya Tsvankint ABSTRACT The main advantage
More informationDoppler shifts in astronomy
7.4 Doppler shift 253 Diide the transformation (3.4) by as follows: = g 1 bck. (Lorentz transformation) (7.43) Eliminate in the right-hand term with (41) and then inoke (42) to yield = g (1 b cos u). (7.44)
More informationMin Chen Department of Mathematics Purdue University 150 N. University Street
EXISTENCE OF TRAVELING WAVE SOLUTIONS OF A HIGH-ORDER NONLINEAR ACOUSTIC WAVE EQUATION Min Chen Department of Mathematics Purdue Uniersity 150 N. Uniersity Street 47907-067 chen@math.purdue.edu Monica
More information3-D description of normal moveout in anisotropic inhomogeneous media
GEOPHYSICS, VOL. 63, NO. 3 (MAY-JUNE 1998); P. 1079-1092,10 FIGS. 3-D description of normal moveout in anisotropic inhomogeneous media Vladimir Grechka* and Ilya Tsvankin* ABSTRACT We present a new equation
More informationVISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION
VISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION Predict Obsere Explain Exercise 1 Take an A4 sheet of paper and a heay object (cricket ball, basketball, brick, book, etc). Predict what will
More informationMapping the conversion point in vertical transversely isotropic (VTI) media
Mapping the conversion point in vertical transversely isotropic (VTI) media Jianli Yang and Don. C. Lawton Conversion-point mapping ABSTRACT The important aspect of converted-wave (P-S) seismology is that
More informationCases of integrability corresponding to the motion of a pendulum on the two-dimensional plane
Cases of integrability corresponding to the motion of a pendulum on the two-dimensional plane MAXIM V. SHAMOLIN Lomonoso Moscow State Uniersity Institute of Mechanics Michurinskii Ae.,, 99 Moscow RUSSIAN
More informationChapter 1. The Postulates of the Special Theory of Relativity
Chapter 1 The Postulates of the Special Theory of Relatiity Imagine a railroad station with six tracks (Fig. 1.1): On track 1a a train has stopped, the train on track 1b is going to the east at a elocity
More informationLecture #8-6 Waves and Sound 1. Mechanical Waves We have already considered simple harmonic motion, which is an example of periodic motion in time.
Lecture #8-6 Waes and Sound 1. Mechanical Waes We hae already considered simple harmonic motion, which is an example of periodic motion in time. The position of the body is changing with time as a sinusoidal
More informationFourier finite-difference wave propagation a
Fourier finite-difference wave propagation a a Published in Geophysics, 76, T123-T129, (2011) Xiaolei Song and Sergey Fomel ABSTRACT We introduce a novel technique for seismic wave extrapolation in time.
More informationAVAZ inversion for fracture orientation and intensity: a physical modeling study
AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian*, Gary F. Margrave, and Joe Wong, University of Calgary. CREWES fmahmoud@ucalgary.ca Summary We present
More informationСollisionless damping of electron waves in non-maxwellian plasma 1
http:/arxi.org/physics/78.748 Сollisionless damping of electron waes in non-mawellian plasma V. N. Soshnio Plasma Physics Dept., All-Russian Institute of Scientific and Technical Information of the Russian
More informationChapter 2: 1D Kinematics Tuesday January 13th
Chapter : D Kinematics Tuesday January 3th Motion in a straight line (D Kinematics) Aerage elocity and aerage speed Instantaneous elocity and speed Acceleration Short summary Constant acceleration a special
More informationSingle soliton solution to the extended KdV equation over uneven depth
Eur. Phys. J. E 7) : DOI./epje/i7-59-7 Regular Article THE EUROPEAN PHYSICAL JOURNAL E Single soliton solution to the etended KdV equation oer uneen depth George Rowlands, Piotr Rozmej,a, Eryk Infeld 3,
More informationCompensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS
Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy
More informationInterval anisotropic parameters estimation in a least squares sense Case histories from West Africa
P-263 Summary Interval anisotropic parameters estimation in a least squares sense Patrizia Cibin*, Maurizio Ferla Eni E&P Division (Milano, Italy), Emmanuel Spadavecchia - Politecnico di Milano (Milano,
More information4-vectors. Chapter Definition of 4-vectors
Chapter 12 4-ectors Copyright 2004 by Daid Morin, morin@physics.harard.edu We now come to a ery powerful concept in relatiity, namely that of 4-ectors. Although it is possible to derie eerything in special
More informationGroup and Renormgroup Symmetry of Quasi-Chaplygin Media
Nonlinear Mathematical Physics 1996, V.3, N 3 4, 351 356. Group and Renormgroup Symmetry of Quasi-Chaplygin Media Vladimir F. KOVALEV Institute for Mathematical Modelling, Miusskaya 4-A, Moscow, 125047,
More informationDirect nonlinear traveltime inversion in layered VTI media Paul J. Fowler*, Alexander Jackson, Joseph Gaffney, and David Boreham, WesternGeco
Paul J. Fowler*, Alexander Jackson, Joseph Gaffney, and David Boreham, WesternGeco Summary We present a scheme for direct nonlinear inversion of picked moveout traveltimes in block-layered isotropic or
More informationP S-wave polarity reversal in angle domain common-image gathers
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? P S-wave polarity reversal in angle domain common-image gathers Daniel Rosales and James Rickett 1 ABSTRACT The change in the reflection
More informationSemi-implicit Treatment of the Hall Effect in NIMROD Simulations
Semi-implicit Treatment of the Hall Effect in NIMROD Simulations H. Tian and C. R. Soinec Department of Engineering Physics, Uniersity of Wisconsin-Madison Madison, WI 5376 Presented at the 45th Annual
More information10. Yes. Any function of (x - vt) will represent wave motion because it will satisfy the wave equation, Eq
CHAPER 5: Wae Motion Responses to Questions 5. he speed of sound in air obeys the equation B. If the bulk modulus is approximately constant and the density of air decreases with temperature, then the speed
More informationMagnetic Fields Part 3: Electromagnetic Induction
Magnetic Fields Part 3: Electromagnetic Induction Last modified: 15/12/2017 Contents Links Electromagnetic Induction Induced EMF Induced Current Induction & Magnetic Flux Magnetic Flux Change in Flux Faraday
More informationSection 6: PRISMATIC BEAMS. Beam Theory
Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam
More informationReversal in time order of interactive events: Collision of inclined rods
Reersal in time order of interactie eents: Collision of inclined rods Published in The European Journal of Physics Eur. J. Phys. 27 819-824 http://www.iop.org/ej/abstract/0143-0807/27/4/013 Chandru Iyer
More informationPhysics 2A Chapter 3 - Motion in Two Dimensions Fall 2017
These notes are seen pages. A quick summary: Projectile motion is simply horizontal motion at constant elocity with ertical motion at constant acceleration. An object moing in a circular path experiences
More informationChapter 1: Kinematics of Particles
Chapter 1: Kinematics of Particles 1.1 INTRODUCTION Mechanics the state of rest of motion of bodies subjected to the action of forces Static equilibrium of a body that is either at rest or moes with constant
More informationSUMMARY ANGLE DECOMPOSITION INTRODUCTION. A conventional cross-correlation imaging condition for wave-equation migration is (Claerbout, 1985)
Comparison of angle decomposition methods for wave-equation migration Natalya Patrikeeva and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Angle domain common image gathers offer
More informationOn anelliptic approximations for qp velocities in VTI media
Geophysical Prospecting, 2004, 52, 247 259 On anelliptic approximations for qp velocities in VTI media Sergey Fomel Bureau of Economic Geology, The University of Texas, Austin, University Station, Box
More information(a) During the first part of the motion, the displacement is x 1 = 40 km and the time interval is t 1 (30 km / h) (80 km) 40 km/h. t. (2.
Chapter 3. Since the trip consists of two parts, let the displacements during first and second parts of the motion be x and x, and the corresponding time interals be t and t, respectiely. Now, because
More informationMoveout approximation for P waves in a homogeneous VTI medium
Moveout approximation for P waves in a homogeneous VTI medium Véronique Farra 1 and Ivan Pšenčík 2 1 Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, UMR 7154 CNRS,
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
FINITE DIFFERENCE TIME DOMAIN MODELING OF INFRASOUND PROPAGATION: APPLICATION TO SHADOW ZONE ARRIVALS AND REGIONAL PROPAGATION Catherine de Groot-Hedlin Uniersity of California at San Diego Sponsored by
More informationAVAZ and VVAZ practical analysis to estimate anisotropic properties
AVAZ and VVAZ practical analysis to estimate anisotropic properties Yexin Liu*, SoftMirrors Ltd., Calgary, Alberta, Canada yexinliu@softmirrors.com Summary Seismic anisotropic properties, such as orientation
More informationAzimuthal AVO and Curvature. Jesse Kolb* David Cho Kris Innanen
Azimuthal AVO and Curvature Jesse Kolb* David Cho Kris Innanen Azimuthal AVO What is azimuthal AVO?: Analysis of incidence angle and azimuthal amplitude variations of reflection coefficients; Measures
More informationy (m)
4 Spring 99 Problem Set Optional Problems Physics February, 999 Handout Sinusoidal Waes. sinusoidal waes traeling on a string are described by wae Two Waelength is waelength of wae?ofwae? In terms of amplitude
More informationMotion in Two and Three Dimensions
PH 1-A Fall 014 Motion in Two and Three Dimensions Lectures 4,5 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More informationLinear Momentum and Collisions Conservation of linear momentum
Unit 4 Linear omentum and Collisions 4.. Conseration of linear momentum 4. Collisions 4.3 Impulse 4.4 Coefficient of restitution (e) 4.. Conseration of linear momentum m m u u m = u = u m Before Collision
More informationElastic wavefield separation for VTI media
CWP-598 Elastic wavefield separation for VTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes from isotropic elastic wavefields is typically
More informationA possible mechanism to explain wave-particle duality L D HOWE No current affiliation PACS Numbers: r, w, k
A possible mechanism to explain wae-particle duality L D HOWE No current affiliation PACS Numbers: 0.50.-r, 03.65.-w, 05.60.-k Abstract The relationship between light speed energy and the kinetic energy
More informationVelocity, Acceleration and Equations of Motion in the Elliptical Coordinate System
Aailable online at www.scholarsresearchlibrary.com Archies of Physics Research, 018, 9 (): 10-16 (http://scholarsresearchlibrary.com/archie.html) ISSN 0976-0970 CODEN (USA): APRRC7 Velocity, Acceleration
More informationVelocity analysis and imaging in transversely isotropic media:
Velocity analysis and imaging in transversely isotropic media: Methodology and a case study TARIQ ALKHALIFAH, ILYA TSVANKIN, KEN LARNER Colorodo School of Mines Golden, Colorado JOHN TOLDI Chevron Overseas
More informationGeneral Lorentz Boost Transformations, Acting on Some Important Physical Quantities
General Lorentz Boost Transformations, Acting on Some Important Physical Quantities We are interested in transforming measurements made in a reference frame O into measurements of the same quantities as
More informationAmplitude preserving AMO from true amplitude DMO and inverse DMO
Stanford Exploration Project, Report 84, May 9, 001, pages 1 167 Amplitude preserving AMO from true amplitude DMO and inverse DMO Nizar Chemingui and Biondo Biondi 1 ABSTRACT Starting from the definition
More informationCorrespondence between the low- and high-frequency limits for anisotropic parameters in a layered medium
GEOPHYSICS VOL. 74 NO. MARCH-APRIL 009 ; P. WA5 WA33 7 FIGS. 10.1190/1.3075143 Correspondence between the low- and high-frequency limits for anisotropic parameters in a layered medium Mercia Betania Costa
More informationPhysics 4A Solutions to Chapter 4 Homework
Physics 4A Solutions to Chapter 4 Homework Chapter 4 Questions: 4, 1, 1 Exercises & Problems: 5, 11, 3, 7, 8, 58, 67, 77, 87, 11 Answers to Questions: Q 4-4 (a) all tie (b) 1 and tie (the rocket is shot
More informationA Petroleum Geologist's Guide to Seismic Reflection
A Petroleum Geologist's Guide to Seismic Reflection William Ashcroft WILEY-BLACKWELL A John Wiley & Sons, Ltd., Publication Contents Preface Acknowledgements xi xiii Part I Basic topics and 2D interpretation
More informationDO PHYSICS ONLINE. WEB activity: Use the web to find out more about: Aristotle, Copernicus, Kepler, Galileo and Newton.
DO PHYSICS ONLINE DISPLACEMENT VELOCITY ACCELERATION The objects that make up space are in motion, we moe, soccer balls moe, the Earth moes, electrons moe, - - -. Motion implies change. The study of the
More informationNMO-velocity surfaces and Dix-type formulas in anisotropic heterogeneous media
GEOPHYSICS, VOL. 67, NO. 3 (MAY-JUNE 2002); P. 939 951, 6 FIGS., 1 TABLE. 10.1190/1.1484536 NMO-velocity surfaces and Dix-type formulas in anisotropic heterogeneous media Vladimir Grechka and Ilya Tsvankin
More informationAVAZ inversion for fracture orientation and intensity: a physical modeling study
AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian and Gary F Margrave ABSTRACT AVAZ inversion We present a pre-stack amplitude inversion of P-wave data
More informationAn Improvement of Velocity Variation with Offset (VVO) Method in Estimating Anisotropic Parameters
DOI 10.7603/s40875-014-0003-0 GSTF Journal of Geological Sciences (JGS) Vol. No.1, July 015 An Improvement of Velocity Variation with Offset (VVO) Method in Estimating Anisotropic Parameters Ida Herawati,
More informationTransmission lines using a distributed equivalent circuit
Cambridge Uniersity Press 978-1-107-02600-1 - Transmission Lines Equialent Circuits, Electromagnetic Theory, and Photons Part 1 Transmission lines using a distributed equialent circuit in this web serice
More informationSPACE-TIME HOLOMORPHIC TIME-PERIODIC SOLUTIONS OF NAVIER-STOKES EQUATIONS. 1. Introduction We study Navier-Stokes equations in Lagrangean coordinates
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 218, pp. 1 5. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SPACE-TIME HOLOMORPHIC
More informationA-level Mathematics. MM03 Mark scheme June Version 1.0: Final
-leel Mathematics MM0 Mark scheme 660 June 0 Version.0: Final Mark schemes are prepared by the Lead ssessment Writer and considered, together with the releant questions, by a panel of subject teachers.
More informationKinetic plasma description
Kinetic plasma description Distribution function Boltzmann and Vlaso equations Soling the Vlaso equation Examples of distribution functions plasma element t 1 r t 2 r Different leels of plasma description
More informationSHEAR-WAVE REFLECTION MOVEOUT FOR AZIMUTHALLY ANISOTROPIC MEDIA
SHEAR-WAVE REFLECTION MOVEOUT FOR AZIMUTHALLY ANISOTROPIC MEDIA AbdulFattah AI-Dajani and M. Nafi Toksoz Earth Resources Laboratory Department of Earth, Atmospheric, and Planetary Sciences Massachusetts
More informationMCAT Physics - Problem Drill 06: Translational Motion
MCAT Physics - Problem Drill 06: Translational Motion Question No. 1 of 10 Instructions: (1) Read the problem and answer choices carefully () Work the problems on paper as 1. An object falls from rest
More informationSupplementary Information Microfluidic quadrupole and floating concentration gradient Mohammad A. Qasaimeh, Thomas Gervais, and David Juncker
Mohammad A. Qasaimeh, Thomas Gerais, and Daid Juncker Supplementary Figure S1 The microfluidic quadrupole (MQ is modeled as two source (Q inj and two drain (Q asp points arranged in the classical quardupolar
More informationP137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients
P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients M. Ayzenberg (StatoilHydro), A. Aizenberg (Institute of Petroleum Geology and Geophysics),
More informationMotion in Two and Three Dimensions
PH 1-1D Spring 013 Motion in Two and Three Dimensions Lectures 5,6,7 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More informationChapter 7. Seismic imaging. 7.1 Assumptions and vocabulary
Chapter 7 Seismic imaging Much of the imaging procedure was already described in the previous chapters. An image, or a gradient update, is formed from the imaging condition by means of the incident and
More informationPhysics Kinematics: Projectile Motion. Science and Mathematics Education Research Group
F FA ACULTY C U L T Y OF O F EDUCATION E D U C A T I O N Department of Curriculum and Pedagogy Physics Kinematics: Projectile Motion Science and Mathematics Education Research Group Supported by UBC Teaching
More informationSeismic inversion for the parameters of two orthogonal fracture sets in a VTI background medium
GEOPHYSICS, VOL. 67, NO. 1 (JANUARY-FEBRUARY 2002); P. 292 299, 3 FIGS. 10.1190/1.1451801 Seismic inversion for the parameters of two orthogonal fracture sets in a VTI background medium Andrey Bakulin,
More informationLECTURE 2: CROSS PRODUCTS, MULTILINEARITY, AND AREAS OF PARALLELOGRAMS
LECTURE : CROSS PRODUCTS, MULTILINEARITY, AND AREAS OF PARALLELOGRAMS MA1111: LINEAR ALGEBRA I, MICHAELMAS 016 1. Finishing up dot products Last time we stated the following theorem, for which I owe you
More informationLowrank RTM for converted wave imaging
Lowrank RM for converted wave imaging Lorenzo Casasanta * (CGG), Zhiguang Xue (U Austin), Sam Gray (CGG) Summary his paper is an attempt to fill the technology gap existing between pure P- and PS-wave
More informationS 1 S 2 A B C. 7/25/2006 Superposition ( F.Robilliard) 1
P S S S 0 x S A B C 7/5/006 Superposition ( F.Robilliard) Superposition of Waes: As we hae seen preiously, the defining property of a wae is that it can be described by a wae function of the form - y F(x
More informationUNIVERSITY OF CALGARY. A comparison of different methods for estimating Thomsen's anisotropy parameters. Chunyan Xiao A THESIS
Important Notice This copy may be used only for the purposes of research and private study, and any use of the copy for a purpose other than research or private study may require the authorization of the
More informationInversion of normal moveout for monoclinic media
Inversion of normal moveout for monoclinic media Vladimir Grechka, Pedro Contreras, and Ilya Tsvankin Center for Wave Phenomena, Department of Geophysics, Colorado School of Mines, Golden, CO 80401-1887
More information