Downloaded 11/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
|
|
- Sherman Haynes
- 6 years ago
- Views:
Transcription
1 Velocity tomography based on turning waves and finite-frequency sensitivity kernels Xiao-Bi Xie*, University of California at Santa Cruz, Jan Pajchel, Statoil, Norway Summary We present a velocity tomography method, which is based on finite-frequency sensitivity kernels and applied to turning waves in models with generally positive velocity gradient. The finite-frequency sensitivity kernels are calculated using the one-way propagator. They replace the ray tracing to create links between the velocity model errors and the observed travel time residuals. The new approach is a wave-equation based method and wave phenomena are naturally incorporated. To demonstrate the inversion process, numerical examples are presented. For 2D velocity models with 0% to 30% velocity errors, we generate the synthetic data set. The travel time differences between seismograms from the true velocity model and those from the initial (intermediate) model are used as input data. The inversion is solved in an iterative way and the velocity models can be successfully updated. Introduction The turning wave data are widely used to determine the subsurface velocities, both at depth and near surface. Traditionally, the turning wave tomography has been dominated by the turning-ray method which assumes an infinitely high frequency. When a broadband signal propagates through a complex region, the ray theory often poorly approximates the actual wave phenomenon. Given this situation, a robust velocity tomography method based on the wave equation theory is preferred. Under this category, the full waveform inversion (FWI) has recently become very attractive. However, if the initial model is drastically biased from the true velocity model, the cycle skipping may cause convergence problem in the FWI. The travel time data are formed by accumulated multiple forward scatterings. As a substitution, we can use wave equation based sensitivity kernel for travel time tomography, which can tolerant large velocity errors and deal with data composed of large time delays. The sensitivity of finite-frequency signal to velocity perturbation has been extensively investigated and successfully used in velocity tomography (e.g., Woodward, 992; Vasco et al., 995; Dahlen et al., 2000; Huang et al., 2000; Zhao et al., 2000; Spetzler and Snieder, 2004; Jocker, et al., 2006; and Liu et al., 2009) and migration velocity analysis (e.g., Shen et al., 2003; Sava and Biondi, 2004; de Hoop, 2006; Xie and Yang, 2007, 2008a, b; Sava and Vlad, 2008; Fliedner et al., 2007; and Xia and Jin 203). Sensitivity kernels are usually derived based on the scattering theory. Several methods can be used to calculate finite-frequency sensitivity kernels in heterogeneous models. The method based on the full-wave finitedifference calculation (e.g., Zhao et al. 2005) can deal with complex velocity models but is very time consuming. On the other hand, the one-way wave equation based propagator (e.g., Xie and Yang 2007, 2008a, b) or the Gaussian Beam method (Xie 20) can be used to generate sensitivity kernels as well. The one-way method is efficient and can properly handle multiple forward scatterings, consisting with the travel time data. Given the large data sizes involved, this method is particularly useful in exploration seismology. We propose a turning wave tomography method based on the finite-frequency sensitivity kernels calculated by one-way wave propagators. Numerical examples are calculated to validate the proposed method. Finite-Frequency Sensitivity Kernel and its Calculation The finite-frequency sensitivity kernel can be defined as how the changes in velocity model will affect the observed finite-frequency seismic signals, e.g., their amplitudes, phases or travel times, or inversely, how the observed signals can sense the changes in velocity model. For a monotonic signal, its phase change can be linked to the velocity model perturbation as rs, r m r K Fr, rs, r dv. () V For a broad band signal, its travel time can be linked to the velocity perturbation as t( rs, rg) mrkbr, rs, r Gdv, (2) V where rs, r and t( rs, r G) are the phase and travel time differences for waves propagating from r S to r G and due to changes in model parameter mr vr v 0 r, is the frequency, v r is the velocity perturbation, v0 r is the background velocity. The spatial integral dv includes all regions with velocity perturbations. V The frequency domain phase sensitivity kernel K can be F expressed as (Woodward, 992; Spetzler and Snieder, 2004; Jocker et al., 2006) 2 Grr ; S, Grr ; KFrr, S, r imag 2k, (3) 0 GrG; rs, G r ; r, is the Green s function calculated in where 2 v0 r and from r to r 2, k0 v0 r, imag[ ] denotes taking the imaginary part. The broadband sensitivity kernel K for travel time can be obtained by stacking monotonic B kernels (Xie and Yang, 2008a) 756
2 W ( ) KBrr, S, rg KFrr, S, r d, (4) where W ( ) is a weighting function related to the source function (Xie and Yang, 2008a). Based on equations (3) and (4), calculating the sensitivity kernels is mostly calculating Green s functions in background velocity models. The oil and gas industry involves processing huge data sets and the efficiency is often a priority. The primary advantage of the one-way propagator is its high efficiency, often order of magnitude faster than the full-way finite difference method. The oneway propagator keeps most of the wave phenomena but neglects the back scattering and reverberations. It is suitable for simulating the travel time information generated by accumulated multiple forward scatterings. Figure. The stored kernel parameters FK. The enlarged parts show their details. The Inversion System Integral equation (2) forms the basis of velocity tomography, in which the sensitivity kernel K plays an B important role to associate the observed travel time residuals with the subsurface velocity errors between the true and trial velocity models. With observed travel time data t( rs, r G), we can invert velocity error vv 0 and use it in velocity model updating. Another important issue is the storage of sensitivity kernels. Unlike seismic rays, the finite-frequency sensitivity kernels are volumetric. Huge space is required to store thousands of kernels. This can cause problems in both storage and input/output. Several techniques can be used to reduce the storage space. Here, following Xie and Yang (2008b), we first partition the integral in equation (2) into the summation of integrals in small rectangular cells V r k t rs, rg m r K Br, rs, r G dv. (5) V r k k Within each cell, we use a hyperbolic function m r af r a2f2 r a3f3 r a4f4 r (6) to interpolate the unknown velocity perturbation, where f r, f r x, f r y, f r xy. (7) For the 4 corners of each cell, r 4, we have 4 equations, which link the coefficient a to the velocity values at the 4 i corners mr fr f2r f3r f4r a m 2 f 2 f2 2 f3 2 f 4 2 a r r r r r 2 (8) mr3 fr3 f2r3 f3r3 f4r3a 3 mr4 fr4 f2r4 f3r4 f4r4 a4 Substituting the values of hyperbolic functions at four corners, we have mr a m 2 x 0 0 a r 2, (9) mr3 0 y 0 a 3 m r4 x y xya4 where x and y are horizontal and vertical sizes of the cell. Using the inverse matrix x x 0 0 P, (0) y 0 y 0 xy xy xy xy equation (9) can express as ai Pm ij r j, () where repeating subscripts denote summation. Substituting (6) - () into the cell integral in (5), we have mrkbr, rs, rgdv V rk afk i i Br, rs, rgdv (2) V rk afk i ik Pm ij rj FKik where FK f r K r dv. (3) ik V r k i B Substituting equation (2) into equation (5) we have the linear system trs, rg Pm ij r jfk. (4) ik k j i 757
3 In this approach, we first determine the cell size according to the required accuracy for inversion. Within cell k, calculate integrals in equation (3) and store only 4 parameters FK for each cell. The accuracy of the kernel 4 storage is adaptive to the accuracy of the velocity model partitioning. Illustrated in Figure are parameters FK for about 000 kernels calculated for a velocity model, where selected kernels are enlarged to show their details. Equation (4) is the discretized version of the integral equation (2) and will be used in the following numerical tests. Figure 2. Snapshots of turning waves passing through velocity models with/without perturbations. Top: in the background model. Middle: in a model with -0% perturbations. Bottom: differential wavefield overlapped on differential velocity models. Numerical Examples () Turning waves in the perturbed velocity model. To illustrate how the velocity perturbations affect the turning wave propagation, shown in Figure 2 are snapshots of turning waves overlapped on the velocity models. These wavefields are calculated using the finite-difference method and a 5 Hz Ricker source. Shown in the top panel is the background model with a positive velocity gradient. In the middle panel, the model has a -0% perturbation which causes the wave front distortion and the amplitude focusing. Shown in the bottom panel is differential wavefield overlapped on differential velocity model. Note, even for such small perturbations, the wave front shift could be comparable to its wavelength. This may cause convergence difficulties for the conventional FWI. Figure 3. Top: synthetic seismograms in models with/without velocity perturbations. Middle: the velocity model with positive anomalies. Bottom: picked travel time differences from the synthetic data. Figure 4. Top panel: the background velocity which also serves as the initial model. Bottom panel: the true velocity model which has a positive velocity inclusion located at the center of the model. The middle panel is the inverted result after 4 iterations. 758
4 (2) Travel time data. To generate a synthetic data set, finite-difference method is used to calculate synthetic seismograms in the true velocity model and the initial (or intermediate) model. The travel time differences between two sets of seismograms are calculated using a cross correlation method. Illustrated in Figure 3 is such an example. The top panel is synthetic shot records in the perturbed and background models. Note the distance range where the travel times are affected by the velocity anomaly. Shown in the middle panel is the velocity model which has positive perturbations near the center part. The schematic wave path indicates the range the perturbations affecting the turning wave propagation. Shown in the bottom panel are picked travel time differences in this model. The data is composed of 35 shots each with 30 receivers. Each column is the result for a shot record. The patterns of the travel time delay provide the information for turning wave tomography. Figure 5. Travel time differences during iterations. (a) to (d) correspond to travel time differences before iterations to 4. Note different color bars are used for different panels. (3) Inversion example. In this example, the background is a -D model with a positive gradient and the true velocity model is the background plus a small patch with +0% velocity perturbations. The inversion is starting from the background model, which is shown on the top panel in Figure 4. The true velocity model is shown at the bottom. The inverted velocity model after 4 iterations is shown in the middle panel. Illustrated in Figure 5 are changes of travel time differences during the iteration process. The maximum travel time differences change from 0.5 s to 0.02 s. (4) Inversion example 2. In this example, we test a velocity model with large perturbations. In Figure 6, the top panel is the background velocity model which also serves as the initial model. The model in the bottom panel is the true velocity model. The maximum velocity perturbations are approximately 30%. Shown in the middle panel is the inverted velocity model after two iterations. The inverted result reproduces large-scale features of the original model but there are some small-scale artifacts due to the footprint of the cell size in the model partitioning. Certain smoothing during the iteration may solve this problem. Figure 6. Top panel: the background velocity model which has a vertical velocity gradient. This model also serves as the initial velocity model. Bottom panel: the true velocity model which is composed of the background velocity plus the velocity perturbation. Shown in the middle panel is the inverted velocity model. Conclusion We derived the theory for the turning-wave tomography based on finite-frequency sensitivity kernels. We developed the numerical method to calculate sensitivity kernels for the turning waves using the one-way propagator. This approach is highly efficient. A technique is designed for velocity model partitioning and efficiently store sensitivity kernels. We represent the unknown velocity model with interpolation functions and project the sensitivity kernels onto these functions. Instead of store the kernel itself, we store only its coefficients. With this technique, the storage is tremendously reduced, as well as the required I/O time. Because it is automatically adapted to the model partition, this method will not cause any accuracy loss. With the proposed inversion system, we successfully tested the turning wave tomography using synthetic data set. This method uses only the travel time information. Thus can handle initial models drastically biased from the true velocity model. The current results are preliminary. 759
5 References Dahlen, F., S. Huang, and G. Nolet, 2000, Fréchet kernels for finite-frequency traveltimes I, Theory: Geophysical Journal International, 4, de Hoop, M.V., R.D. van der Hilst, and P. Shen, 2006, Wave equation reflection tomography: annihilators and sensitivity kernels: Geophysical Journal International, 67, Fliedner, M.M., M.P. Brown, D. Bevc, and B. Biondi, 2007, Wave path tomography for subsalt velocity model building: 77th Annual International Meeting, SE Expanded Abstracts, Huang, S., F. Dahlen, and G. Nolet, 2000, Frechet kernels for finite-frequency travel times II, Examples: Geophys. J. Int. 4, Jocker, J., J. Spetzler, D. Smeulders, and J. Trampert, 2006, Validation of first-order diffraction theory for the traveltimes and amplitudes of propagating waves: Geophysics, 7, T67-T77. Liu, Y., L. Dong, Y. Wang, J. Zhu, and Z. Ma, 2009, Sensitivity kernels for seismic Fresnel volume tomography: Geophysics, 74, U35-U46. Sava, P.C., and B. Biondi, 2004, Wave-equation migration velocity analysis - I: Theory: Geophysical Prospecting, 52, Sava, P. and I. Vlad, 2008, Numeric implementation of wave-equation migration velocity analysis operators: Geophysics, 73, VE45-VE59. Shen, P., W. Symes, and C. Stolk, 2003, Differential semplance velocity analysis by wave-equation migration: 73rd Annual International Meeting, SE Expanded Abstracts, Spetzler, J. and R. Snieder, 2004, The Fresnel volume and transmitted waves: Geophysics, 69, Vasco, D.W., J.E. Peterson, Jr., and E.L. Majer, 995, Beyond ray tomography: Wavepaths and Fresnel volumes, Geophysics, 60, Woodward, M. J., 992, Wave-equation tomography: Geophysics, 57, Xia, F., and Jin S., 203, Angle-domain wave equation tomography using RTM image gathers, 83rd Annual International Meeting, SE Expanded Abstract, Xie, X.B., 20, Calculating finite-frequency sensitivity kernels using the Gaussian beam method, 8th Annual International Meeting, SE Expanded Abstracts 30, Xie, X.B., and H. Yang, 2007, A migration velocity updating method based on the shot index common image gather and finite-frequency sensitivity kernel: 77th Annual International Meeting, SE Expanded Abstracts, Xie, X.B., and H. Yang, 2008a, The finite-frequency sensitivity kernel for migration residual moveout and its applications in migration velocity analysis, Geophysics, 73, S24-S249. Xie, X.B., and H. Yang, 2008b, A wave-equation migration velocity analysis approach based on the finitefrequency sensitivity kernel, 78th Annual International Meeting, SE Expanded Abstracts, Zhao, L., T.H. Jordan, and C.H. Chapman, 2000, Threedimensional Fréchet differential kernels for seismic delay times: Geophysical Journal International, 4, Zhao, L., T.H. Jordan, K.B. Olsen, and P. Chen, 2005, Frechet kernels for imaging regional earth structure based on three-dimensional reference models: Bull. Seism. Soc. Am., 95, , doi: 0.785/
Main Menu. SEG Houston 2009 International Exposition and Annual Meeting. rp = rr dr, (1) () r
elocity analysis for plane-wave source migration using the finite-frequency sensitivity kernel Hui Yang, ormerly niversity of California at anta Cruz; presently LandOcean Energy ervices Co., Ltd., eijing,
More informationTraveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a
Traveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a a Published in SEP report, 103, 61-68 (2000) James Rickett 1 INTRODUCTION Estimating an accurate velocity function is one of the
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 informationTomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima
1 2 3 Auxiliary materials for Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima nuclear power plant area 4 5 6 7 8 9 Ping Tong 1,2, Dapeng Zhao 1 and Dinghui Yang 2 [1] {Department of Geophysics,
More informationWe G Elastic One-Return Boundary Element Method and Hybrid Elastic Thin-Slab Propagator
We G05 0 Elastic One-Return Boundary Element Method and Hybrid Elastic Thin-lab Propagator R.. Wu* (University of California) & Z. Ge (Beijing University, China) UMMARY We developed the theory and algorithm
More informationComplex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada
Page 1 of 10 Home Articles Interviews Print Editions Complex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada DECEMBER 2012 FOCUS ARTICLE Summary Gaussian-beam depth migration
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 information2012 SEG SEG Las Vegas 2012 Annual Meeting Page 1
Wei Huang *, Kun Jiao, Denes Vigh, Jerry Kapoor, David Watts, Hongyan Li, David Derharoutian, Xin Cheng WesternGeco Summary Since the 1990s, subsalt imaging in the Gulf of Mexico (GOM) has been a major
More informationSummary. Introduction
Detailed velocity model building in a carbonate karst zone and improving sub-karst images in the Gulf of Mexico Jun Cai*, Hao Xun, Li Li, Yang He, Zhiming Li, Shuqian Dong, Manhong Guo and Bin Wang, TGS
More informationSummary. Introduction
with constrained iterative tomography: Methodology and application Hans Kristian Helgesen* (Statoil), Jun Tang (WesternGeco), Jinjun Liu (Statoil), Antoun Salama, Randolph Pepper, Sam Madden, Marta Woodward,
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 informationFull Waveform Inversion via Matched Source Extension
Full Waveform Inversion via Matched Source Extension Guanghui Huang and William W. Symes TRIP, Department of Computational and Applied Mathematics May 1, 215, TRIP Annual Meeting G. Huang and W. W. Symes
More informationVELOCITY REDUCTION - A POINT IN VELOCITY MODEL BUILDING FOR PRE-STACK KIRCHOFF DEPTH MIGRATION
Physical Geography; Cartography; Geographic Information Systems & Spatial Planing VELOCITY REDUCTION - A POINT IN VELOCITY MODEL BUILDING FOR PRE-STACK KIRCHOFF DEPTH MIGRATION Dr. Maya Grigorova University
More informationValidity of the Rytov Approximation in the Form of Finite-Frequency Sensitivity Kernels
Pure Appl. Geophys. 172 (2015), 1409 1427 Ó 2014 Springer Basel DOI 10.1007/s00024-014-0991-8 Pure and Applied Geophysics Validity of the Rytov Approximation in the Form of Finite-Frequency Sensitivity
More informationFinite-frequency tomography
Finite-frequency tomography Lapo Boschi (lapo@erdw.ethz.ch) October 14, 2009 Application of Born theory (scattering theory) to the ray-theory solution (the forward problem) In earlier lectures we have
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2013) 194, 1042 1049 Advance Access publication 2013 May 3 doi: 10.1093/gji/ggt144 Finite-frequency sensitivity kernels for two-station surface wave
More informationImage-space wave-equation tomography in the generalized source domain
Image-space wave-equation tomography in the generalized source domain Yaxun Tang, Claudio Guerra, and Biondo Biondi ABSTRACT We extend the theory of image-space wave-equation tomography to the generalized
More informationVollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich
Seminar des FA Ultraschallprüfung Vortrag 1 More info about this article: http://www.ndt.net/?id=20944 Vollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich Thomas
More informationAttenuation compensation in viscoacoustic reserve-time migration Jianyong Bai*, Guoquan Chen, David Yingst, and Jacques Leveille, ION Geophysical
Attenuation compensation in viscoacoustic reserve-time migration Jianyong Bai*, Guoquan Chen, David Yingst, and Jacques Leveille, ION Geophysical Summary Seismic waves are attenuated during propagation.
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 informationSUMMARY. (Sun and Symes, 2012; Biondi and Almomin, 2012; Almomin and Biondi, 2012).
Inversion Velocity Analysis via Differential Semblance Optimization in the Depth-oriented Extension Yujin Liu, William W. Symes, Yin Huang,and Zhenchun Li, China University of Petroleum (Huadong), Rice
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 informationValidation of first-order diffraction theory for the traveltimes and amplitudes of propagating waves
GEOPHYSICS, VO. 71, NO. 6 NOVEMBER-DECEMBER 26; P. T167 T177, 13 FIGS. 1.119/1.2358412 Validation of first-order diffraction theory for the traveltimes and amplitudes of propagating waves Jeroen Jocker
More informationWaveform inversion for attenuation estimation in anisotropic media Tong Bai & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines
Waveform inversion for attenuation estimation in anisotropic media Tong Bai & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines SUMMARY Robust estimation of attenuation coefficients remains
More informationSUMMARY INTRODUCTION. f ad j (t) = 2 Es,r. The kernel
The failure mode of correlation focusing for model velocity estimation Hyoungsu Baek 1(*), Henri Calandra 2, and Laurent Demanet 1 1 Dept. of Mathematics and Earth Resources Lab, Massachusetts Institute
More informationImaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges
2017 SCEC Proposal Report #17133 Imaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges Principal Investigator Zhongwen Zhan Seismological Laboratory, California
More informationComparison between least-squares reverse time migration and full-waveform inversion
Comparison between least-squares reverse time migration and full-waveform inversion Lei Yang, Daniel O. Trad and Wenyong Pan Summary The inverse problem in exploration geophysics usually consists of two
More informationElastic least-squares reverse time migration
CWP-865 Elastic least-squares reverse time migration Yuting Duan, Paul Sava, and Antoine Guitton Center for Wave Phenomena, Colorado School of Mines ABSTRACT Least-squares migration (LSM) can produce images
More informationModeling, migration, and inversion in the generalized source and receiver domain
Modeling, migration, and inversion in the generalized source and receiver domain Yaxun Tang ABSTRACT I extend the theory of Born modeling, migration and inversion to the generalized source and receiver
More informationRegistration-guided least-squares waveform inversion
Registration-guided least-squares waveform inversion Hyoungsu Baek 1, Henri Calandra, Laurent Demanet 1 1 MIT Mathematics department, TOTAL S.A. January 15 013 Abstract Full waveform inversion with frequency
More information3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES
3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES Barbara Romanowicz Department of Earth and Planetary Science, U. C. Berkeley Dr. Barbara Romanowicz, UC Berkeley (KITP Colloquium 9/11/02) 1 Cartoon
More informationSUMMARY MOVEOUT DERIVATION FOR SICP-WEMVA
The sensitivity of converted-phase extended images to P- and S-wavespeed variations Andrey H. Shabelansky, Earth Resources Laboratory, Massachusetts Institute of Technology; Alison Malcolm, Memorial University
More informationElastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC)
Elastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC) Summary We develop an elastic full waveform inversion
More informationThe Improvement of 3D Traveltime Tomographic Inversion Method
Advances in Petroleum Exploration and Development Vol. 5, No., 013, pp. 36-40 DOI:10.3968/j.aped.1955438013050.1368 ISSN 195-54X [Print] ISSN 195-5438 [Online] www.cscanada.net www.cscanada.org The Improvement
More informationThe Estimation of PS Wave Statics By Using SCDCC + Stack Power Method
The Estimation of PS Wave Statics By Using SCDCC + Stack Power Method ABSTRACT Summary Aimin Xue * Petrol Sound Geoservices. 8, MacEwan Ridge Close, N.W., Calgary, Alberta T3K3A7, Canada. aiminx@telus.net
More informationVáclav Bucha. Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague. SW3D meeting June 6-7, 2016 C OM S TR 3 D
Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department
More informationYouzuo Lin and Lianjie Huang
PROCEEDINGS, Thirty-Ninth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 24-26, 2014 SGP-TR-202 Building Subsurface Velocity Models with Sharp Interfaces
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
FINITE-FREQUENCY SEISMIC TOMOGRAPHY OF BODY WAVES AND SURFACE WAVES FROM AMBIENT SEISMIC NOISE: CRUSTAL AND MANTLE STRUCTURE BENEATH EASTERN EURASIA Yong Ren 2, Wei Zhang 2, Ting Yang 3, Yang Shen 2,and
More informationJ.A. Haugen* (StatoilHydro ASA), J. Mispel (StatoilHydro ASA) & B. Arntsen (NTNU)
U008 Seismic Imaging Below "Dirty" Salt J.A. Haugen* (StatoilHydro ASA), J. Mispel (StatoilHydro ASA) & B. Arntsen (NTNU) SUMMARY Base and sub salt seismic imaging is still an unresolved issue. To solve
More informationApproximate- vs. full-hessian in FWI: 1D analytical and numerical experiments
Approximate- vs. full-hessian in FWI: 1D analytical and numerical experiments Raul Cova and Kris Innanen ABSTRACT Feasibility of using Full Waveform Inversion (FWI) to build velocity models has been increasing
More informationTu N Fault Shadow Removal over Timor Trough Using Broadband Seismic, FWI and Fault Constrained Tomography
Tu N118 05 Fault Shadow Removal over Timor Trough Using Broadband Seismic, FWI and Fault Constrained Tomography Y. Guo* (CGG), M. Fujimoto (INPEX), S. Wu (CGG) & Y. Sasaki (INPEX) SUMMARY Thrust-complex
More informationSeismic tomography with co-located soft data
Seismic tomography with co-located soft data Mohammad Maysami and Robert G. Clapp ABSTRACT There is a wide range of uncertainties present in seismic data. Limited subsurface illumination is also common,
More informationFull Waveform Inversion (FWI) with wave-equation migration. Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010
Full Waveform Inversion (FWI) with wave-equation migration (WEM) and well control Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010 Outline The FWI cycle The fundamental theorem of FWI Understanding
More informationSeismogram Interpretation. Seismogram Interpretation
Travel times in the Earth Ray paths, phases and their name Wavefields in the Earth: SH waves, P-SV waves Seismic Tomography Receiver Functions Seismogram Example Long-period transverse displacement for
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2) 85, 385 42 doi:./j.365-246x.2.4945.x Tomographic errors from wave front healing: more than just afastbias Alison E. Malcolm and Jeannot Trampert 2
More informationMarchenko redatuming: advantages and limitations in complex media Summary Incomplete time reversal and one-sided focusing Introduction
Marchenko redatuming: advantages and limitations in complex media Ivan Vasconcelos*, Schlumberger Gould Research, Dirk-Jan van Manen, ETH Zurich, Matteo Ravasi, University of Edinburgh, Kees Wapenaar and
More informationWenyong Pan and Lianjie Huang. Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 019 SGP-TR-14 Adaptive Viscoelastic-Waveform Inversion Using the Local Wavelet
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 informationNonlinear seismic imaging via reduced order model backprojection
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, Vladimir Druskin 2 and Mikhail Zaslavsky 2 University of Houston, 2 Schlumberger-Doll Research Center Mamonov, Druskin,
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 informationStructural Cause of Missed Eruption in the Lunayyir Basaltic
GSA DATA REPOSITORY 2015140 Supplementary information for the paper Structural Cause of Missed Eruption in the Lunayyir Basaltic Field (Saudi Arabia) in 2009 Koulakov, I., El Khrepy, S., Al-Arifi, N.,
More informationThe 2D/3D i-stats Workflow for Image-Based Near-Surface Modeling for Statics Corrections
The 2D/3D i-stats Workflow for Image-Based Near-Surface Modeling for Statics Corrections Öz Yilmaz CTO, GeoTomo LLC With the i-stats workflow, you no longer need first-break picking as for traveltime tomography,
More informationTh Guided Waves - Inversion and Attenuation
Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information
More informationVelocity analysis using AB semblance a
Velocity analysis using AB semblance a a Published in Geophysical Prospecting, v. 57, 311-321 (2009) Sergey Fomel ABSTRACT I derive and analyze an explicit formula for a generalized semblance attribute,
More information( ) ( ), 1, 1. Downloaded 10/14/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Chao Zhang*, Jie Zhang, University of Science and Technology of China (USTC); Zhikun Sun, GeoTomo LLC Summary Refraction traveltimes have long been applied for deriving long-wavelength statics solutions.
More informationResearch Project Report
Research Project Report Title: Prediction of pre-critical seismograms from post-critical traces Principal Investigator: Co-principal Investigators: Mrinal Sen Arthur Weglein and Paul Stoffa Final report
More informationExtended isochron rays in prestack depth (map) migration
Extended isochron rays in prestack depth (map) migration A.A. Duchkov and M.V. de Hoop Purdue University, 150 N.University st., West Lafayette, IN, 47907 e-mail: aduchkov@purdue.edu (December 15, 2008)
More informationFOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES
FOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES Vilma Castillejos Hernandez Supervisor: Tatsuhiko Hara MEE10508 ABSTRACT We performed time domain moment tensor
More informationOverview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples
Overview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples Guanghui Huang Education University of Chinese Academy of Sciences, Beijing, China Ph.D. in Computational
More informationVelocity Update Using High Resolution Tomography in Santos Basin, Brazil Lingli Hu and Jianhang Zhou, CGGVeritas
Lingli Hu and Jianhang Zhou, CGGVeritas Summary The exploration interest in the Santos Basin offshore Brazil has increased with the large deep water pre-salt discoveries, such as Tupi and Jupiter. As the
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 informationThe i-stats: An Image-Based Effective-Medium Modeling of Near-Surface Anomalies Oz Yilmaz*, GeoTomo LLC, Houston, TX
The i-stats: An Image-Based Effective-Medium Modeling of Near-Surface Anomalies Oz Yilmaz*, GeoTomo LLC, Houston, TX Summary Near-surface modeling for statics corrections is an integral part of a land
More informationMain Menu. Summary. Introduction. stands for general space-time coordinates. γ
Yu Geng*, ILab, IGPP, University of California, Santa Cruz, visiting from Institute of Wave and Information, Xi an Jiaotong University, Ru-Shan Wu, ILab, IGPP, University of California, Santa Cruz, Jinghuai
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 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 informationStrain Green s Tensors, Reciprocity, and Their Applications to Seismic Source and Structure Studies
Bulletin of the Seismological Society of America, Vol. 96, No. 5, pp. 1753 1763, October 2006, doi: 10.1785/0120050253 Strain Green s Tensors, Reciprocity, and Their Applications to Seismic Source and
More informationSeismic Waves in Complex 3 D Structures, 26 (2016), (ISSN , online at
Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationReflection Seismic Method
Reflection Seismic Method Data and Image sort orders; Seismic Impedance; -D field acquisition geometries; CMP binning and fold; Resolution, Stacking charts; Normal Moveout and correction for it; Stacking;
More informationA 27-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method
Geophysical Prospecting 04 6 58 77 doi: 0./365-478.090 A 7-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method Jing-Bo Chen Key Laboratory of Petroleum Resources
More informationEvidence of an axial magma chamber beneath the ultraslow spreading Southwest Indian Ridge
GSA Data Repository 176 1 5 6 7 9 1 11 1 SUPPLEMENTARY MATERIAL FOR: Evidence of an axial magma chamber beneath the ultraslow spreading Southwest Indian Ridge Hanchao Jian 1,, Satish C. Singh *, Yongshun
More informationStatic Corrections for Seismic Reflection Surveys
Static Corrections for Seismic Reflection Surveys MIKE COX Volume Editors: Series Editor: Eugene F. Scherrer Roland Chen Eugene F. Scherrer Society of Exploration Geophysicists Tulsa, Oklahoma Contents
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 informationFull waveform inversion in the Laplace and Laplace-Fourier domains
Full waveform inversion in the Laplace and Laplace-Fourier domains Changsoo Shin, Wansoo Ha, Wookeen Chung, and Ho Seuk Bae Summary We present a review of Laplace and Laplace-Fourier domain waveform inversion.
More informationSeismic Response and Wave Group Characteristics of Reef Carbonate Formation of Karloff-Oxford Group in Asser Block
Seismic Response and Wave Group Characteristics of Reef Zeng zhongyu Zheng xuyao Hong qiyu Zeng zhongyu Zheng xuyao Hong qiyu Institute of Geophysics, China Earthquake Administration, Beijing 100081, China,
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
TRANSITION ZONE WAVE PROPAGATION: CHARACTERIZING TRAVEL-TIME AND AMPLITUDE INFORMATION Peter M. Shearer and Jesse F. Lawrence University of California San Diego, Institute of Geophysics and Planetary Physics
More informationDirty salt velocity inversion: The road to a clearer subsalt image
GEOPHYSICS. VOL. 76, NO. 5 (SEPTEMBER-OCTOBER 2011); P. WB169 WB174, 8 FIGS. 10.1190/GEO2010-0392.1 Dirty salt velocity inversion: The road to a clearer subsalt image Shuo Ji 1, Tony Huang 1, Kang Fu 2,
More informationAttenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation
Attenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation Gaurav Dutta, Kai Lu, Xin Wang and Gerard T. Schuster, King Abdullah University of Science and Technology
More informationANGLE-DEPENDENT TOMOSTATICS. Abstract
ANGLE-DEPENDENT TOMOSTATICS Lindsay M. Mayer, Kansas Geological Survey, University of Kansas, Lawrence, KS Richard D. Miller, Kansas Geological Survey, University of Kansas, Lawrence, KS Julian Ivanov,
More informationMatrix formulation of adjoint Kirchhoff datuming
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 377 Matrix formulation of adjoint Kirchhoff datuming Dimitri Bevc 1 ABSTRACT I present the matrix formulation of Kirchhoff wave equation datuming
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
ANALYSIS AND SIMULATION OF THREE-DIMENSIONAL SCATTERING DUE TO HETEROGENEOUS CRUSTAL STRUCTURE AND SURFACE TOPOGRAPHY ON REGIONAL PHASES; MAGNITUDE AND DISCRIMINATION Arben Pitarka 1, Don V. Helmberger
More informationAnisotropic tomography for TTI and VTI media Yang He* and Jun Cai, TGS
media Yang He* and Jun Cai, TGS Summary A simultaneous anisotropic tomographic inversion algorithm is developed. Check shot constraints and appropriate algorithm preconditioning play an important role
More informationDeconvolution imaging condition for reverse-time migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 83 96 Deconvolution imaging condition for reverse-time migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT The reverse-time
More informationResidual Statics using CSP gathers
Residual Statics using CSP gathers Xinxiang Li and John C. Bancroft ABSTRACT All the conventional methods for residual statics analysis require normal moveout (NMO) correction applied on the seismic data.
More informationWE SRS2 11 ADAPTIVE LEAST-SQUARES RTM FOR SUBSALT IMAGING
Technical paper WE SRS2 11 ADAPTIVE LEAST-SQUARES RTM FOR SUBSALT IMAGING Authors C. Zeng (TGS), S. Dong (TGS), B. Wang (TGS) & Z. Zhang* (TGS) 2016 TGS-NOPEC Geophysical Company ASA. All rights reserved.
More informationStolt residual migration for converted waves
Stanford Exploration Project, Report 11, September 18, 1, pages 1 6 Stolt residual migration for converted waves Daniel Rosales, Paul Sava, and Biondo Biondi 1 ABSTRACT P S velocity analysis is a new frontier
More informationA Nonlinear Differential Semblance Strategy for Waveform Inversion: Experiments in Layered Media Dong Sun and William W Symes, Rice University
A Nonlinear Differential Semblance Strategy for Waveform Inversion: Experiments in Layered Media Dong Sun and William W Symes, Rice University SUMMARY This paper proposes an alternative approach to the
More informationEstimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
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 informationAPPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY
APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY Timur TEZEL Supervisor: Takuo SHIBUTANI MEE07169 ABSTRACT In this study I tried to determine the shear wave velocity structure in the crust
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2010) doi: 10.1111/j.1365-246X.2010.04643.x Empirically determined finite frequency sensitivity kernels for surface waves Fan-Chi Lin and Michael H.
More informationAngle-domain common-image gathers from anisotropic Gaussian beam migration and its application to anisotropy-induced imaging errors analysis
J. Earth Syst. Sci. (217) 126: 1 c Indian Academy of Sciences DOI 1.17/s124-16-783-x Angle-domain common-image gathers from anisotropic Gaussian beam migration and its application to anisotropy-induced
More informationRough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical
Rough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical Summary This paper discusses rough-sea de-ghosting for variabledepth streamer data. The
More informationRAYFRACT IN MARINE SURVEYS. The data provided by the contractor needed total re-picking and reinterpretation.
RAYFRACT IN MARINE SURVEYS As part of a geotechnical assessment and feasibility planning of channel improvement a shallow marine seismic refraction survey was undertaken. The data was initially processed
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. 2014) 196, 1827 1843 Advance Access publication 2014 January 2 doi: 10.1093/gji/ggt496 Non-linear partial derivative and its De Wolf approximation for
More informationSource estimation for frequency-domain FWI with robust penalties
Source estimation for frequency-domain FWI with robust penalties Aleksandr Y. Aravkin, Tristan van Leeuwen, Henri Calandra, and Felix J. Herrmann Dept. of Earth and Ocean sciences University of British
More informationIntroduction to Seismic Imaging
Introduction to Seismic Imaging Alison Malcolm Department of Earth, Atmospheric and Planetary Sciences MIT August 20, 2010 Outline Introduction Why we image the Earth How data are collected Imaging vs
More informationScattering-based decomposition of sensitivity kernels for full-waveform inversion
Scattering-based decomposition of sensitivity kernels for full-waveform inversion Daniel Leal Macedo (Universidade Estadual de Campinas) and Ivan asconcelos (Schlumberger Cambridge Research) Copyright
More informationSUMMARY REVIEW OF THE FREQUENCY DOMAIN L2 FWI-HESSIAN
Efficient stochastic Hessian estimation for full waveform inversion Lucas A. Willemsen, Alison E. Malcolm and Russell J. Hewett, Massachusetts Institute of Technology SUMMARY In this abstract we present
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
DEVELOPMENT OF REGIONAL PHASE TOMOGRAPHIC ATTENUATION MODELS FOR EURASIA Thorne Lay 1, Xiao-Bi Xie 1, and Xiaoning (David) Yang 2 University of California, Santa Cruz 1 and Los Alamos National Laboratory
More information