Application of Weighted Early-Arrival Waveform Inversion to Shallow Land Data
|
|
- Geoffrey Rogers
- 5 years ago
- Views:
Transcription
1 Application of Weighted Early-Arrival Waveform Inversion to Shallow Land Data Han Yu, Dongliang Zhang and Xin Wang King Abdullah University of Science and Technology, Division of Physical Science and Engineering, Thuwal Saudi Arabia ABSTRACT Recent studies have shown that inverting traces weighted by the energy of the early-arrivals can improve the accuracy of estimating shallow velocities. This is explained by showing that the associated misfit gradient function tends to be sensitive to the kinematics of wave propagation and insensitive to the dynamics. A synthetic example verifies the theoretical predictions and show that the effects of noise and unpredicted amplitude variations in the inversion are reduced using this weighted early arrival waveform inversion (WEWI). We also apply this method to a D land data set for estimating the near-surface velocity distribution. The reverse time migration images suggest that, compared to the tomogram inverted directly from the early arrival waveforms, the WEWI tomogram provides a more convincing velocity model and more focused reflections in the deeper part of the image.. INTRODUCTION The near-surface velocity distribution is crucial for imaging the deeper parts of the Earth. Complex velocity variations at the near surface are often associated with undulating topography or irregular geology in the near-surface weathered layers (Amorium et al., 987; Taner et al., 998). If the near-surface velocity distribution is not accurately estimated, the coherency of the deeper migrated reflections can be strongly degraded (White, 989; Marsden, 993). To partly remedy this problem, the near-surface velocity model with smooth variations can be estimated by traveltime tomography (Zhu and McMechan, 989; Pratt and Goulty, 99; Aki and Richards, 00) that inverts the first-arrival traveltimes. However, in geologically complex areas, a more highly resolved velocity model is needed for imaging deeper reflectors. In this regard, waveform inversion (Tarantola, 984; Mora, 987; Zhou et al., 995) was developed to invert for more accurate tomograms by finite-frequency 6 seismic wave propagation. To reduce the computational time and local minima problems (Sirgue and Pratt, 004), early-arrival waveform inversion (EWI) was proposed by Sheng et al. (006) in the space-time domain and later applied to marine data (Boonyasiriwat et al., 00). In this work, we carry out the inversion on land data by following the conventional EWI method but using a recently developed objective misfit function (Shen, 00), whose gradient is more robust and focuses more on matching the phase rather than the amplitude in the data. However, the associated gradient does not have an important energy normalization term which is important for optimal imaging. In this work, the gradient associated with this weighted early arrival waveform inversion (WEWI) is properly normalized and shown to significantly improve the accuracy of the final tomogram. Instead of replacing the amplitude spectrum of a calculated trace with that of the corresponding observed trace (Sun and Schuster, 993), we implement WEWI in the time domain by normalizing both the observed and calculated early arrivals using the L norm of the trace, where this approach avoids the phase wrapping problem in the frequency domain (Shin and Min, 006). Our synthetic results demonstrate that compared to EWI, WEWI can mitigate the effects of noise and unpredicted amplitude variations in the data and robustly invert for highly resolved near-surface tomogram. Moreover, a land data test illustrates that WEWI produces a more accurate shallow subsurface tomogram where the energy is focused in the deeper part. This paper is organized into four sections. The first part is the introduction, and the second part analyzes the misfit function associated with its gradient in our approach. In section 3, numerical results are shown for inverting data associated with the Marmousi model and a field experiment in Saudi Arabia. The last section presents the conclusions.. THEORY In many field, particularly land data sets, there are strong elastic arrivals such as surface waves that cannot be modeled by the
2 6 Yu et al. acoustic wave equation. In addition, the amplitudes of some traces are distorted due to unexplained environmental sources and not explained by geometric spreading. Some elastic effects in the data can be reduced by applying an early arrival window to mute the later arrivals. Therefore, the conventional waveform inversion misfit function is modified by Shen (00) and expressed in the time domain as E= = p(x r, t x s ), s, r s, r p calc (x r, t x s ) p obs(x r, t x s ) p calc (x r, t x s ) p obs (x r, t x s ), () which indicates that p pt p calc p p calc calc eliminates the phase information of p calc in p obs. Therefore the common phases in p calc and p obs will not be selected to match again in the next iteration. Figure shows fact the fact this new virtual source is perpendicular to the recorded data. It thus weakens the effects of p calc and strengthens the phase difference between p obs and p calc in the back propagated wave fields, therefore making WEWI robust and significant. where p(x r, t x s ) denotes the pressure field trace recorded at the receiver position x r, with listening time t and a source at x s ; p denotes the L norm of the N vector p, namely p T p, where N is the number of time samples in the trace. Here, p obs represents the recorded trace with windowed early arrivals and p calc represents the synthetic early arrivals. The synthetic data are calculated by solving the constant-density acoustic wave equation, p(x, t x s ) c (x) t p(x, t x s )= s(x, t x s ), () where c(x) represents the velocity model at position x. The solution to equation is calculated by a second order in time and eighth order in space staggered-grid method (Levander, 988). Equation normalizes the observed and the synthetic early arrivals so that their energy can be compared at the same scale, and the waveform inversion in this case is more sensitive to phase differences in the misfit function. To unveil this fact, the Fréchet derivative [grad(x)= E c(x) ] of the functional E with respect to c(x) is calculated by grad(x)= p T (p calc/ p calc ) c = ( p pt p calc p calc p calc p calc ) T p calc c. (3) Here, if the scalar / p calc is ignored, the first term in equation 3 exactly corresponds to the gradient of conventional waveform misfit function. In this paper, we invert the field data using this gradient but ignore / p calc because it is included in calculating the step length when updating the velocity model. Previous synthetic tests (Shen, 00) missed the superscript term in the denominator p calc of equation 3, thus making the phase match less accurate. This misfit function becomes more significant if the term pt p calc is not small, which can be caused p calc by complicated geological conditions, because it attaches more importance to accurately predicting the phases rather than the amplitudes. Note, that the gradient term p pt p calc p p calc calc is orthogonal to p calc since ( p ( p)t p calc p calc p calc ) T p calc = p T p calc p T p calc = 0, (4) Figure : Construction of data residuals as back propagating sources for EWI and WEWI. The velocity model is estimated by an iterative conjugate gradient method where c k+ (x)=c k (x)+λ k d k (x), (5) and the conjugate directions are defined by d k = P k g k +β k d k, (6) for iterations k=,,..., k max, g=[grad(x)], and P is the conventional geometrical-spreading preconditioner (Causse et al., 999). The scalarλ k is the step length which can be determined by a quadratic line-search method (Nocedal and Wright, 999), and d k (x) is the component of the direction vector d k (x) indexed by x. For the first iteration, we set d 0 = g 0. The parameterβ k is calculated by the Polak-Ribiére formula (Nocedal and Wright, 999) β k = gt k (P kg k P k g k ) g T k P k g k. (7) To compute the gradient direction at each iteration reduces to computing the reverse time migration operation. Additional forward modelings are required for the line search. The initial velocity model c 0 (x) is the traveltime tomogram inverted by picked first arrivals (Nemeth et al., 997), and equation 5 is iteratively applied until the objective functional E satisfies a stopping criterion.
3 Weighted EWI to Shallow Land Data NUMERICAL TESTS OF WEWI 3. Synthetic Data Example The Marmousi model is used to test the robustness and quality of WEWI before it is applied to a land data set in the next section. First, a synthetic data set is generated based on the Marmousi model (Figure (a)) with a 576 (horizontal) by 84 (vertical) gridded mesh with a 6.0 m grid interval. There are 60 shots with a 54 m shot spacing, and for each shot, the number of receivers is 90 with a 6 m receiver spacing. The recording length is.5 s with a sampling rate 0.5 ms. White noise is added to each trace of every shot gather. The nonzero meanvalued noise consists of two parts: random receiver noise is added to every CSG and a random static amplitude shift is also applied to each trace. Figures 3(a) and 3(b) show one common shot gather before and after adding the noise, and they are displayed with the same amplitude scale. Figure 3(c) shows the noise mask added to Figure 3(a) that simulates dead traces. Finally, WEWI and conventional EWI are used to invert the noisy data up to 0.5 s after the first arrival using the same initial velocity model (see Figure (b)). The resulting tomograms after the 30 th iteration are presented in Figures 4(a) and 4(b), which proves that WEWI is less sensitive to the noise; the matched phase information in the predicted data is excluded in calculating the gradient using reverse time migration (RTM). Their misfit gradient Figures 5(a) and 5(b) for the first iteration further validate the advantage of the WEWI method, followed by two obviously different convergence rates of the data residual shown in Figure 5(c). The comparison of the two tomograms apparently indicates WEWI can invert for a more accurate velocity model with unpredicted noise in the data, and imply the robustness of the WEWI compared to EWI (a) Marmousi Velocity Model (b) Initial Velocity Model Figure : The Marmousi models with (a) the true velocity distribution and (b) the its smoothed version as the initial velocity model for waveform inversion. km/s (a) CSG #0 before Adding the Noise (b) CSG #0 after Adding the Noise (c) The Noise Mask for CSG # Figure 3: The CSG #0 generated by the Marmousi model (a) before, (b) after adding the noise, and (c) its noise mask. 3. Land Data Example 3.. Acquisition and Processing A D seismic survey is was carried out near KAUST with the acquisition geometry illustrated in Figure 6. The D acquisition line consists of 79 shots and 40 vertical component geophones per shot, with a uniform spacing of 30 m for both shots and receivers. For each channel, the record length is s with a 4 ms sampling rate. For a common shot gather (CSG), the shot position is in the middle of the 40 receivers, so the largest offset is 0 m. WEWI is applied to the first CSGs, and the horizontal distance for our inversion is restricted between shot # and # (red crosses in Figure 6). The topography of these shots is shown in Figure 7, which is almost flat if the horizontal distance X and the depth Z are of the same scale. In this case, we choose Z= 5 m as the surface for the inversion and ignore the elevation variations of the geophones and shotpoints. Prior to data processing, it is useful to estimate the nearsurface velocity distribution from the picked first arrivals. Figure 8 shows CSG #7 with picked first arrivals marked by red crosses, and the direct wave and the refractions are respectively marked by the white and green dashed lines. The slopes of the
4 64 Yu et al. (a) WEWI Tomogram (a) WEWI Gradient for the First Iteration km/s (b) Conventional EWI Tomogram (b) EWI Gradient for the First Iteration Figure 4: The inverted tomograms using (a) WEWI and (b) conventional EWI. two lines suggest that the velocity corresponding to the first layer is approximately v = 750 m/ s=3750 m/s, whereas the velocity for the second layer is about v = 800 m/0.5 s= 5600 m/s. We first apply the F-K filter to remove apparent surface waves in all the CSGs. As an example, the raw CSG # is shown in Figure 9(a), with most of the surface waves eliminated in Figure 9(b) by F-K filtering. The field data are then transformed from 3D to D format by applying the filter i/ω to all the traces in the frequency domain and scaled by t for 3D geometrical spreading in the time domain. In the meantime, the data are also filtered with a pass band from 0 to 0 Hz. The results after filtering of CSG # are shown in Figure 9(c). The spectra of a trace in CSG # before and after the processing steps are presented in Figure 0. Normalized L Data Residual (c) Data Residual Convergence Rate Comparison for Synthetic Data Iteration Number WEWI Convetional EWI Figure 5: (a) WEWI and (b) conventional EWI gradients for the first interation, and (c) the comparison of their convergence rate. Acquisition Geometry with Equal Scales in Distance 3.. Applying WEWI to The Land Data Set WEWI is implemented using a starting model calculated by first-arrival traveltime tomgoraphy with a 757 (horizontal) by 97 (vertical) gridded mesh and a 7.5 m grid interval. First, a ray-based traveltime tomography method is used to invert the first arrival traveltimes for a smooth velocity model. The traveltime tomogram is shown in Figure, where the high-velocity layer is about 80 m deep, as illustrated by the refraction data in Figure 8. Second, WEWI is used to invert the data restricted by dynamically increasing time windows. The time window stretches from the beginning of the record to 4 s after the first arrivals for the first 0 iterations and its size quadratically increases up to 0.5 s in the next 0 iterations. All the traces bounded by their time windows in the observed and calculated data are normalized according to equation during each iteration. The / p calc term in equation 3 is dropped and the rest of the terms are exactly followed when calculating the gradient. Y [km].6 0 Red Crosses: the first CSGs for inversion X [km] Figure 6: The acquisition geometry with equal scales in the horizontal and the vertical directions. The inverted WEWI tomogram at the 30 th iteration is shown in Figure (a). At the 0 th iteration, the observed data and cal-
5 Weighted EWI to Shallow Land Data 65 0 The Topography of Shot to Figure 7: Topography of the first shots.. (a) Raw CSG # CSG # 7 with Picked First Arrivals Figure 8: A common shot gather with its picked first arrivals. Red crosses: picked first arrivals; green dashed line: the estimated direct wave; white dash line: the estimated refraction. culated data recorded 4 s after the first arrivals mostly agree with one another in Figures 3(a) and 3(b). Note that nearoffset traces within 50 m from the source and far-offset traces more than 3300 m from the sources are not inverted due to the poor quality of the recorded data. After the dynamic windows are applied, many later arrival events also match according to Figures 4(a) and 4(b) at the 30 th iteration without harming the matched previous arrivals. Figure 5 depicts the residual vs. iteration number plot when applying WEWI for 30 iterations. Note that the residual gradually includes more data for comparison while the window size grows very slowly. To verify the inverted velocity model, we migrate the high frequency portion ( 45 Hz) of the observed data using the traveltime and the WEWI tomograms as migration velocity models. The reverse time migration (RTM) images above m in depth are exhibited in Figures 6(a) and 6(b) and their associated common image gathers are also shown in Figures 7(a) and 7(b). The RTM image based on the WEWI tomogram shows a more continuous structure in the shallow part and more focused energy in the deep part compared to the RTM image computed with the traveltime tomogram. The CIGs based on the WEWI tomogram also are flatter for both the near-surface and the deeper reflectors. The EWI is tomogram shown in Figure (b). Some shallow reflectors can still be detected according to this tomogram but they are much stronger compared to the tomogram in Figure (a). Figure (b) could be less realistic because the estimated near-surface velocity at Z = 80m is around 5600 m/s and not higher than 6000 m/s according to the picked first arrivals in Figure 8. The calculated CSG # for the degraded tomogram presented in Figure 3(c) does not match well with (b) F K Filtered CSG # (c) CSG # after Bandpass Filtering and 3D to D Transformation Figure 9: A raw (a) CSG #, (b) with its surface waves removed by the dip filter, and (c) the CSG # after bandpass filtering and 3D to D transformation. Figure 3(a) even for the early arrivals at the intermediate 0 th iteration. When the time window increases to 0.5 s in the 30 th iteration, the calculated CSG # with conventional EWI (4(c)) shows greater difference from the observed CSG # than its counterpart with WEWI (Figure 4(b)). The RTM image (Figure 6(c)) using this tomogram is also inferior compared to Figure 6(b), and so are the CIGs (Figures 7(c) and 8(c)) associated with it. Moreover, the normalized EWI data residual only decreases to 5 after 30 iterations as shown in Figure CONCLUSIONS A modified misfit function is applied and its associated Fréchet derivative is exactly followed to calculate the gradient for updating the velocity model. This modified function using traceby-trace normalization attaches more importance to accurately
6 66 Yu et al. (a) Spectrum of a Trace of Raw CSG # (a) Observed CSG # Normalized Amplitude Frequency (Hz) (b) Spectrum of the Trace in (a) after Filtering Normalized Amplitude Frequency (Hz) (b) Calculated CSG # with WEWI Figure 0: The spectrum of a trace in CSG # (a) before and (b) after processing. Traveltime Tomogram km/s (c) Calculated CSG # with Conventional EWI Figure : The traveltime tomogram for the picked first arrivals from CSGs # #. (a) WEWI Tomogram (b) Conventional EWI Tomogram Figure : (a) The WEWI, and (b) the conventional EWI tomograms inverted from the early arrivals of CSGs # #. predicting the phase than the amplitude of the recorded data. The gradient associated with this objective function mitigates km/s Figure 3: (a) The processed CSG #, (b) the calculated CSG # with the WEWI gradient, and (c) the calculated CSG # computed with conventional EWI. the effects of noise and unpredicted amplitude variations, and it strengthens the frequency difference between the observed and calculated data in the back propagated wave fields. The robustness and the phase matching property of WEWI are then validated by inverting for the Marmousi model using a polluted data set. This method is also tested on a real case by inverting CSGs of a land data set after regular processing. A dynamic time window is also used to invert the early arrivals. So that both refractions and some early arriving reflections are included in the WEWI approach. Our field data results suggest that WEWI can generate a more accurate and highly resolved velocity model compared to the conventional EWI tomogram. Although the data with peak frequency around 5 Hz are first used in the inversion, the drawbacks of WEWI for land data largely come from the lack of lower frequency data from -5 Hz. Part of this problem can come from the filters that remove low frequency information polluted by surface waves or other noise, and this might be remedied by better processing
7 Weighted EWI to Shallow Land Data 67 (a) Observed CSG # Muted by A Dynamic Time Window (a) RTM Image Using Traveltime Tomogram (b) RTM Image using WEWI Tomogram (b) Calculated CSG # Muted by A Dynamic Time Window (c) RTM Image Using Conventional EWI Tomogram (c) Calculated CSG # with Conventional EWI Figure 6: RTM images based on (a) the traveltime tomogram, (b) the WEWI tomogram, and (c) the conventional EWI tomogram. to continue widening the time window or inverting the events from deeper reflectors Figure 4: (a) The processed CSG #, (b) the calculated CSG # with WEWI, and (c) the calculated CSG # with EWI muted by a dynamic time window. Normalized L Data Residual 0.9 Data Residual Convergence Rate Comparison for Land Data WEWI Conventional EWI Iteration Nunber Figure 5: The convergence rate of WEWI for 30 iterations. techniques. Also, WEWI may still encounter the cycle skipping problem although it better utilizes the phase information in the data. Cycle skipping problems can be partly overcome by only inverting the near offset traces, and then with later iterations invert the larger offset traces. For reflections events from deep reservoir geology in oil industry, it is also necessary 5. ACKNOWLEDGMENTS We would like to thank the 03 sponsors of the CSIM Consortium ( for their financial support. The computation resource Shaheen ( for inversion provided by the high performance computing (HPC) center of King Abdullah University of Science and Technology (KAUST) is greatly appreciated. We also thank Prof. Schuster and anonymous CSIM members for their professional comments in the development of this work. REFERENCES Aki, K. and P. G. Richards, 00, Quantitative seismology: nd edition. University Science Books. Amorium, W. N. D., P. Hubral, and M. Tygel, 987, Computing field statics with the help of seismic tomography: Geophysical Prospecting, 35, Boonyasiriwat, C., G. T. Schuster, P. Valasek, and W. Cao, 00, Applications of multiscale waveform inversion to marine data using a flooding technique and dynamic earlyarrival windows: Geophysics, 75, R9 R36. Causse, E., R. Mittet, and B. Ursin, 999, Preconditioning for full-waveform inversion in viscoacoustic media: Geophysics, 64, Levander, A. R., 988, Fourth-order finite-difference p-sv seismograms: Geophysics, 53,
8 68 Yu et al. (a) CIGs based on Traveltime Tomogram (a) CIGs based on Traveltime Tomogram (b) CIGs based on WEWI Tomogram (b) CIGs based on WEWI Tomogram (c) CIGs based on Conventional EWI Tomogram (c) CIGs based on Conventional EWI Tomogram Figure 7: One zoomed view of common image gathers based on (a) the traveltime tomogram, (b) the WEWI tomogram, and (c) the conventional EWI tomogram. Marsden, D., 993, Statics corrections-a review: The Leading Edge,, Mora, P. R., 987, Nonlinear two-dimensional elastic inversion of multioff-set seismic data: Geophysics, 5, 8. Nemeth, T., E. Normark, and F. Qin, 997, Dynamic smoothing in crosswell traveltime tomography: Geophysics, 6, Nocedal, J. and S. J. Wright, 999, Numerical optimization: Springer series in operations research and financial engineering. Pratt, R. G. and N. R. Goulty, 99, Combining wave-equation imaging with traveltime tomography to form high-resolution images from crosshole data: Geophysics, 56, Shen, X., 00, Near-surface velocity estimation by weighted early-arrival waveform inversion: 80th SEG Technical Program Expanded Abstracts, Sheng, J., A. Leeds, M. Buddensiek, and G. T. Schuster, 006, Early arrival waveform tomography on near-surface refraction data: Geophysics, 7, U47 U57. Shin, C. and D. J. Min, 006, Waveform inversion using a logrithmic wavefield: Geophysics: Geophysics, 7, R3 R4. Figure 8: Another zoomed view of common image gathers based on (a) the traveltime tomogram, (b) the WEWI tomogram, and (c) the conventional EWI tomogram. Sirgue, L. and R. G. Pratt, 004, Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies: Geophysics, 69, Sun, Y. and G. T. Schuster, 993, Time-domain phase inversion: SEG Technical Program Expanded Abstracts 993, Taner, M. T., D. E. Wagner, E. Baysal, and L. Lu, 998, A unified method for -d and 3-d refraction statics: Geophysics, 63, Tarantola, A., 984, Inversion of seismic reflection data in the acoustic approximation: Geophysics, 49, White, D. J., 989, Two-dimensional seismic refraction tomography: Geophysics,, Zhou, C., W. Cai, Y. Luo, G. T. Schuster, and S. Hassanzadeh, 995, Acoustic wave-equation traveltime and waveform inversion of crosshole seismic data: Geophysics, 60, Zhu, X. and G. A. McMechan, 989, Estimation of a twodimensional seismic compressional-wave velocity distribution by iterative tomographic imaging: International Journal of Imaging System and Technology,, 3 7.
Multi-source Least-squares Migration of Gulf of Mexico Data
Multi-source Least-squares Migration of Gulf of Mexico Data Xin Wang. King Abdullah University of Science and Technology, Thuwal 955-69, Kingdom of Saudi Arabia Corresponding author is Xin Wang. E-mail
More informationDownloaded 09/16/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Time-domain incomplete Gauss-Newton full-waveform inversion of Gulf of Mexico data Abdullah AlTheyab*, Xin Wang, Gerard T. Schuster, King Abdullah University of Science and Technology Downloaded 9// to
More informationFull waveform inversion of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC)
of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC) Summary (FWI) is a promising technology for next generation model building. However, it faces with many
More informationF031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data
F031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data D Yingst (ION Geophysical), C. Wang* (ION Geophysical), J. Park (ION Geophysical), R. Bloor (ION Geophysical), J. Leveille
More informationImprovements in time domain FWI and its applications Kwangjin Yoon*, Sang Suh, James Cai and Bin Wang, TGS
Downloaded 0/7/13 to 05.196.179.38. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Improvements in time domain FWI and its applications Kwangjin Yoon*,
More informationDownloaded 10/26/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
Reflection Full-waveform Inversion for Inaccurate Starting Models Abdullah AlTheyab* and G. T. Schuster, King Abdullah University of Science and Technology (KAUST) SUMMARY We propose a seismic-reflection
More informationU043 3D Prestack Time Domain Full Waveform Inversion
U043 3D Prestack Time Domain Full Waveform Inversion D.V. Vigh* (WesternGeco), W.E.S. Starr (WesternGeco) & K.D. Kenneth Dingwall (WesternGeco) SUMMARY Despite the relatively high computational demand,
More informationElastic full waveform Inversion for land walkaway VSP data
Elastic full waveform Inversion for land walkaway VSP data Olga Podgornova, Scott Leaney, arwan Charara, Schlumberger; Eric von Lunen, Nexen Energy ULC Summary Full waveform inversion (FWI) is capable
More informationCrosswell Imaging by 2-D Prestack Wavepath Migration
Crosswell Imaging by 2-D Prestack Wavepath Migration Hongchuan Sun ABSTRACT Prestack wavepath migration (WM) is applied to 2-D synthetic crosswell data, and the migrated images are compared to those from
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationChapter 5. 3D data examples REALISTICALLY COMPLEX SYNTHETIC INVERSION. Modeling generation and survey design
Chapter 5 3D data examples In this chapter I will demonstrate the e ectiveness of the methodologies developed in the previous chapters using 3D data examples. I will first show inversion results of a realistically
More informationBuilding starting model for full waveform inversion from wide-aperture data by stereotomography
Building starting model for full waveform inversion from wide-aperture data by stereotomography Vincent Prieux 1, G. Lambaré 2, S. Operto 1 and Jean Virieux 3 1 Géosciences Azur - CNRS - UNSA, France;
More informationRobustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
More informationH003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography
H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography M. Cavalca* (Schlumberger - Westerngeco) & R.P. Fletcher (Schlumberger - Westerngeco) SUMMARY Estimation of the
More informationHigh definition tomography brings velocities to light Summary Introduction Figure 1:
Saverio Sioni, Patrice Guillaume*, Gilles Lambaré, Anthony Prescott, Xiaoming Zhang, Gregory Culianez, and Jean- Philippe Montel (CGGVeritas) Summary Velocity model building remains a crucial step in seismic
More informationStrategies for elastic full waveform inversion Espen Birger Raknes and Børge Arntsen, Norwegian University of Science and Technology
Strategies for elastic full waveform inversion Espen Birger Raknes and Børge Arntsen, Norwegian University of Science and Technology SUMMARY Ocean-bottom cables (OBC) have become common in reservoir monitoring
More informationSUMMARY LEAST-SQUARES MIGRATION THEORY
Making the most out of the least (squares migration) Gaurav Dutta, Yunsong Huang, Wei Dai, Xin Wang, and G.T. Schuster King Abdullah University of Science and Technology SUMMARY Standard migration images
More informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
More informationSUMMARY INTRODUCTION THEORY
Illumination compensation for subsalt image-domain wavefield tomography Tongning Yang, Center for Wave Phenomena, Colorado School of Mines, Jeffrey Shragge, School of Earth and Environment, University
More information3D image-domain wavefield tomography using time-lag extended images
CWP-748 3D image-domain wavefield tomography using time-lag extended images Tongning Yang and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT Image-domain wavefield tomography is
More informationRefraction Full-waveform Inversion in a Shallow Water Environment
Refraction Full-waveform Inversion in a Shallow Water Environment Z. Zou* (PGS), J. Ramos-Martínez (PGS), S. Kelly (PGS), G. Ronholt (PGS), L.T. Langlo (PGS), A. Valenciano Mavilio (PGS), N. Chemingui
More informationSUMMARY ELASTIC SCALAR IMAGING CONDITION
Robust 3D scalar imaging condition for elastic RTM Yuting Duan, presently at Shell International Exploration and Production Inc., formerly at Center for Wave Phenomena, Colorado School of Mines Paul Sava,
More informationProgress Report on: Interferometric Interpolation of 3D SSP Data
Progress Report on: Interferometric Interpolation of 3D SSP Data Sherif M. Hanafy ABSTRACT We present the theory and numerical results for interferometrically interpolating and extrapolating 3D marine
More informationSUMMARY INTRODUCTION THEORY. Objective function of ERWI
Extended Reflection Waveform Inversion via Differential Semblance Optimization Yujin Liu, William W. Symes and Zhenchun Li, China University of Petroleum (Huadong), Rice University SUMMARY Reflection-based
More informationDownloaded 09/03/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Full-waveform inversion in a shallow water environment: A North Sea 3D towed-streamer data example Kathy Zou*, Lars Tore Langlo, Grunde Rønholt, Jaime Ramos-Martinez and Steve Kelly, PGS Summary We apply
More informationWave-equation migration from topography: Imaging Husky
Stanford Exploration Project, Report 123, October 31, 2005, pages 49 56 Short Note Wave-equation migration from topography: Imaging Husky Jeff Shragge 1 INTRODUCTION Imaging land seismic data is wrought
More informationSummary. Introduction
Dmitry Alexandrov, Saint Petersburg State University; Andrey Bakulin, EXPEC Advanced Research Center, Saudi Aramco; Pierre Leger, Saudi Aramco; Boris Kashtan, Saint Petersburg State University Summary
More informationSUMMARY. earth is characterized by strong (poro)elasticity.
with illumination compensation Tongning Yang, Center for Wave Phenomena, Colorado School of Mines, Jeffrey Shragge, School of Earth and Environment, University of Western Australia, and Paul Sava, Center
More informationTh C 02 Model-Based Surface Wave Analysis and Attenuation
Th C 02 Model-Based Surface Wave Analysis and Attenuation J. Bai* (Paradigm Geophysical), O. Yilmaz (Paradigm Geophysical) Summary Surface waves can significantly degrade overall data quality in seismic
More informationMulticomponent land data pre-processing for FWI: a benchmark dataset
Multicomponent land data pre-processing for FWI: a benchmark dataset Raul Cova, Bernie K. Law and Kris Innanen CRWES/University of Calgary Summary Successful full-waveform inversion (FWI) studies using
More informationAPPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA.
APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA. BY: ISAAC KUMA YEBOAH. Department of Engineering, Regent University
More informationTh LHR5 08 Multi-modal Surface Wave Inversion and Application to North Sea OBN Data
Th LHR5 08 Multi-modal Surface Wave Inversion and pplication to North Sea ON Data S. Hou (CGG), D. Zheng (CGG), X.G. Miao* (CGG) & R.R. Haacke (CGG) SUMMRY Surface-wave inversion (SWI) for S-wave velocity
More informationSuccessful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment.
Successful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment. Sergey Birdus 1, Dean Criddle 2, Alexey Artyomov 1, Li Li 1, Qingbing Tang 1,
More informationTarget-oriented wavefield tomography: A field data example
Target-oriented wavefield tomography: A field data example Yaxun Tang and Biondo Biondi ABSTRACT We present a strategy for efficient migration velocity analysis in complex geological settings. The proposed
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationLeast squares Kirchhoff depth migration: important details
Least squares Kirchhoff depth migration: important details Daniel Trad CREWES-University of Calgary Summary Least squares migration has been an important research topic in the academia for about two decades,
More informationAccelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method
Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,
More informationLeast squares Kirchhoff depth migration with. preconditioning
Least squares Kirchhoff depth migration with preconditioning Aaron Stanton University of Alberta, Department of Physics, 4-83 CCIS, Edmonton AB T6G E (April, 3) Running head: Least Squares Migration ABSTRACT
More informationG004 Robust Frequency-domain Multi-parameter Elastic Full Waveform Inversion
G004 Robust Frequency-domain Multi-parameter Elastic Full Waveform Inversion R. Brossier* (Joseph Fourier University), S. Operto (CNRS/University of Nice-Sophia-Antipolis/Géoazur) & J. Virieux (Joseph
More informationGeometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT. August Lau and Chuan Yin.
Geometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT August Lau and Chuan Yin January 6, 2017 Abstract The goal of seismic processing is to convert input data
More informationGG450 4/5/2010. Today s material comes from p and in the text book. Please read and understand all of this material!
GG450 April 6, 2010 Seismic Reflection I Today s material comes from p. 32-33 and 81-116 in the text book. Please read and understand all of this material! Back to seismic waves Last week we talked about
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Full Waveform Inversion Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 18, 2014 Bee Bednar (Panorama Technologies)
More information1D internal multiple prediction in a multidimensional world: errors and recommendations
1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate
More informationSUMMARY NOISE REDUCTION BY SPATIAL SLOPE ESTIMATION
: steerable and non-linear filter approach Abdulaziz M. Almuhaidib and M. Nafi Toksöz, Earth Resources Laboratory, MIT SUMMARY We present an approach based on local-slope estimation for the separation
More informationDownward propagation: using true near-surface velocity model, one-way wave- equation extrapolator:
Class 9: Beyond Seismic Statics Corrections Wed, Oct 7, 2009 Dynamic corrections: wave-equation and wavefield datuming methods Skipping the near-surface by interferometry (no model needed) Common-Focus
More informationLab 3: Depth imaging using Reverse Time Migration
Due Wednesday, May 1, 2013 TA: Yunyue (Elita) Li Lab 3: Depth imaging using Reverse Time Migration Your Name: Anne of Cleves ABSTRACT In this exercise you will familiarize yourself with full wave-equation
More informationWe B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study
We B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study E. Zabihi Naeini* (Ikon Science), N. Kamath (Colorado School of Mines), I. Tsvankin (Colorado School of Mines), T. Alkhalifah
More informationDeconvolution in the radial trace domain
R-T domain deconvolution Deconvolution in the radial trace domain David C. Henley ABSTRACT The radial trace (R-T) domain has been shown to be useful for coherent noise attenuation and other seismic wavefield
More informationWe N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions
We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)
More informationMigration from a non-flat datum via reverse-time extrapolation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 50 Migration from a non-flat datum via reverse-time extrapolation Gopal Palacharla 1 ABSTRACT Land surveys usually have elevation changes,
More informationSEG/New Orleans 2006 Annual Meeting
3-D tomographic updating with automatic volume-based picking Dimitri Bevc*, Moritz Fliedner, Joel VanderKwaak, 3DGeo Development Inc. Summary Whether refining seismic images to evaluate opportunities in
More informationFull-waveform inversion for reservoir characterization: A synthetic study
CWP-889 Full-waveform inversion for reservoir characterization: A synthetic study Nishant Kamath, Ilya Tsvankin & Ehsan Zabihi Naeini ABSTRACT Most current reservoir-characterization workflows are based
More informationSUMMARY THEORY INTRODUCTION
Double-Difference Waveform Inversion of 4D Ocean Bottom Cable Data: Application to Valhall, North Sea Di Yang, Michael Fehler and Alison Malcolm, MIT, Faqi Liu and Scott Morton, Hess Corporation SUMMARY
More informationDownloaded 01/20/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
D Multi-source Least-squares Reverse Time Migration Wei Dai, C. Boonyasiriwat and Gerard T. Schuster, King Abdullah University of Science and Technology Downloaded 1//1 to 19.171.17.1. Redistribution subject
More informationHeadwave Stacking in Terms of Partial Derivative Wavefield
Geosystem Engineering, 7(1), 21-26 (March 2004) Headwave Stacking in Terms of Partial Derivative Wavefield Changsoo Shin School of Civil, Urban and Geosystem Engineering, Seoul National University, San
More informationSimultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco
Simultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco Summary In this paper, we review the near-surface challenges
More informationGeogiga Seismic Pro 8.3 Release Notes
Geogiga Seismic Pro 8.3 Release Notes Copyright 2017, All rights reserved. Table of Contents Introduction...1 Part 1 Utility Modules...2 Part 2 Reflection Modules...4 Updates in SF Imager...5 Updates in
More informationProcessing converted-wave data in the tau-p domain: rotation toward the source and moveout correction
τ-p domain converted-wave processing Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction Raul Cova and Kris Innanen ABSTRACT The asymmetry of the converted-wave
More informationA least-squares shot-profile application of time-lapse inverse scattering theory
A least-squares shot-profile application of time-lapse inverse scattering theory Mostafa Naghizadeh and Kris Innanen ABSTRACT The time-lapse imaging problem is addressed using least-squares shot-profile
More informationSUMMARY. min. m f (m) s.t. m C 1
Regularizing waveform inversion by projections onto convex sets application to the D Chevron synthetic blind-test dataset Bas Peters *, Zhilong Fang *, Brendan Smithyman #, Felix J. Herrmann * * Seismic
More informationSUMMARY THEORY INTRODUCTION
Acoustic 3D least-squares reverse time migration using the energy norm Daniel Rocha, Paul Sava & Antoine Guitton Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose a least-squares reverse
More informationLeast Squares Kirchhoff Depth Migration: potentials, challenges and its relation to interpolation
Least Squares Kirchhoff Depth Migration: potentials, challenges and its relation to interpolation Daniel Trad CGG Summary Least Squares Migration (LSM), like interpolation, has the potential to address
More informationE044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media
E044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media M. Cavalca* (WesternGeco), I. Moore (WesternGeco), L. Zhang (WesternGeco), S.L. Ng (WesternGeco), R.P. Fletcher (WesternGeco)
More informationStochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston
Stochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston Summary With the development of computational power, there has been an increased
More informationOn the Scattering Effect of Lateral Discontinuities on AVA Migration
On the Scattering Effect of Lateral Discontinuities on AVA Migration Juefu Wang* and Mauricio D. Sacchi Department of Physics, University of Alberta, 4 Avadh Bhatia Physics Laboratory, Edmonton, AB, T6G
More informationReverse-time migration imaging with/without multiples
Reverse-time migration imaging with/without multiples Zaiming Jiang, John C. Bancroft, and Laurence R. Lines Imaging with/without multiples ABSTRACT One of the challenges with reverse-time migration based
More informationSource deghosting for synchronized multi-level source streamer data Zhan Fu*, Nan Du, Hao Shen, Ping Wang, and Nicolas Chazalnoel (CGG)
Source deghosting for synchronized multi-level source streamer data Zhan Fu*, Nan Du, Hao Shen, Ping Wang, and Nicolas Chazalnoel (CGG) Summary Ghost wavefield elimination is pivotal for improving the
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 5: Refraction Migration Methods Wed, April 13, 2011
GEOPHYS 242: Near Surface Geophysical Imaging Class 5: Refraction Migration Methods Wed, April 13, 2011 Migration versus tomography Refraction traveltime and wavefield migration The theory of interferometry
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationG017 Beyond WAZ - A Modeling-based Evaluation of Extensions to Current Wide Azimuth Streamer Acquisition Geometries
G017 Beyond WAZ - A Modeling-based Evaluation of Extensions to Current Wide Azimuth Streamer Acquisition Geometries M. Cvetkovic* (ION Geophysical), Z. Zhou (ION Geophysical / GXT Imaging Solutions) &
More informationA projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines
A projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines SUMMARY A Hessian matrix in full waveform inversion (FWI) is difficult
More informationM. Warner* (S-Cube), T. Nangoo (S-Cube), A. Umpleby (S-Cube), N. Shah (S-Cube), G. Yao (S-Cube)
Tu A12 15 High-Resolution Reflection FWI M. Warner* (S-Cube), T. Nangoo (S-Cube), A. Umpleby (S-Cube), N. Shah (S-Cube), G. Yao (S-Cube) Summary We demonstrate reflection FWI on a less-than-ideal 3D narrow-azimuth
More informationMitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic
Mitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic A.S. Long* (PGS), S. Lu (PGS), D. Whitmore (PGS), H. LeGleut (PGS), R. Jones (Lundin
More informationIterative resolution estimation in Kirchhoff imaging
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 369?? Iterative resolution estimation in Kirchhoff imaging Robert G. Clapp, Sergey Fomel, and Marie Prucha 1 ABSTRACT We apply iterative
More informationA comparison between time domain and depth domain inversion to acoustic impedance Laurence Letki*, Kevin Darke, and Yan Araujo Borges, Schlumberger
Laurence Letki*, Kevin Darke, and Yan Araujo Borges, Schlumberger Summary Geophysical reservoir characterization in a complex geologic environment remains a challenge. Conventional amplitude inversion
More informationDownloaded 10/23/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
ACQUISITION APERTURE CORRECTION IN ANGLE-DOMAIN TOWARDS THE TRUE- REFLECTION RTM Rui Yan 1*, Huimin Guan 2, Xiao-Bi Xie 1, Ru-Shan Wu 1, 1 IGPP, Earth and Planetary Sciences Department, University of California,
More informationSUMMARY. amounts to solving the projected differential equation in model space by a marching method.
Subsurface Domain Image Warping by Horizontal Contraction and its Application to Wave-Equation Migration Velocity Analysis Peng Shen, Shell International E&P, William W. Symes, Rice University SUMMARY
More informationG012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry
G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry D.F. Halliday* (Schlumberger Cambridge Research), P.J. Bilsby (WesternGeco), J. Quigley (WesternGeco) & E. Kragh (Schlumberger
More informationRay based tomography using residual Stolt migration
Stanford Exploration Project, Report 11, November 11, 00, pages 1 15 Ray based tomography using residual Stolt migration Robert G. Clapp 1 ABSTRACT In complex areas, residual vertical movement is not an
More informationModel-space vs. data-space normalization for recursive depth migration
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Model-space vs. data-space normalization for recursive depth migration James Rickett 1 ABSTRACT Illumination problems caused by finite-recording
More informationAnisotropy-preserving 5D interpolation by hybrid Fourier transform
Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,
More informationThe dynamics of statics
ACQUISITION/PROCESSING Coordinated by Jeff Deere The dynamics of statics M. TURHAN TANER, Rock Solid Images, Houston, USA A. J. BERKHOUT, Technical University of Delft, The Netherlands SVEN TREITEL, TriDekon,
More informationGeogiga Seismic Pro 8.0 Release Notes
Geogiga Seismic Pro 8.0 Release Notes Copyright 2015, All rights reserved. Table of Contents Introduction...1 Part 1 General Enhancements...3 Trace Display...4 Color Section Display...6 Curve Plotting...8
More informationSeismic Reflection Method
Seismic Reflection Method 1/GPH221L9 I. Introduction and General considerations Seismic reflection is the most widely used geophysical technique. It can be used to derive important details about the geometry
More informationEfficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers
Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou
More informationWe G High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation
We G103 05 High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation G. Hilburn* (TGS), Y. He (TGS), Z. Yan (TGS) & F. Sherrill (TGS) SUMMARY An approach
More informationY015 Complementary Data-driven Methods for Interbed Demultiple of Land Data
Y015 Complementary Data-driven Methods for Interbed Demultiple of Land Data S. Sonika* (WesternGeco), A. Zarkhidze (WesternGeco), J. Heim (WesternGeco) & B. Dragoset (WesternGeco) SUMMARY Interbed multiples
More informationComplex-beam migration: non-recursive and recursive implementations
Comple-beam migration: non-recursive and recursive implementations Tianfei Zhu*, CGGVeritas, Calgar AB, Canada Tianfei.Zhu@cggveritas.com Summar We have developed a comple-beam method for shot-domain prestack
More informationPractical implementation of SRME for land multiple attenuation
Practical implementation of SRME for land multiple attenuation Juefu Wang* and Shaowu Wang, CGGVeritas, Calgary, Canada juefu.wang@cggveritas.com Summary We present a practical implementation of Surface
More informationBuilding accurate initial models using gain functions for. waveform inversion in the Laplace domain
Building accurate initial models using gain functions for waveform inversion in the Laplace domain Wansoo Ha 1 and Changsoo Shin 2 1 Department of Energy Resources Engineering, Pukyong National University
More informationDownloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
High fidelity imaging using reflections, refractions and multiples: North Sea example, Johan Sverdrup Field Grunde Rønholt*, Øystein Korsmo, Barbara Danielsen, Samuel Brown, Sverre Brandsberg-Dahl, Alejandro
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Aliasing Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 16, 2014 Bee Bednar (Panorama Technologies) Seismic
More informationTh G Surface-offset RTM Gathers - What Happens When the Velocity is Wrong?
Th G103 01 Surface-offset RTM Gathers - What Happens When the Velocity is Wrong? J.P. Montel* (CGG) SUMMARY Surface offset migrated common image gathers (SOCIGs) built by wave equation migration (WEM)
More informationSUMMARY. u s(e, x,t)u r(e, x +,t), (1)
Image-domain and data-domain waveform tomography: a case study Esteban Díaz 1, Yuting Duan 1, Gerhard Pratt 2, and Paul Sava 1 1 Center for Wave Phenomena, Colorado School of Mines, 2 University of Western
More informationu = v is the Laplacian and represents the sum of the second order derivatives of the wavefield spatially.
Seismic Elastic Modelling Satinder Chopra, Arcis Corporation, Calgary, Canada Summary Seismic modeling experiments were performed to understand the seismic responses associated with a variety of subsurface
More informationA case study for salt model building using CFP full azimuth data W. Gao*, Z. Guo, M. Guo, Q. Zhang, S. Hightower, G. Cloudy Jr. and Z.
case study for salt model building using CFP full azimuth data W. Gao*, Z. Guo, M. Guo, Q. Zhang, S. Hightower, G. Cloudy Jr. and Z. Li, TGS Summary We present a case study of the salt model building for
More informationResidual Moveout Analysis in a 3D dataset from Amplitude Variations with Offfset
P-270 Residual Moveout Analysis in a 3D dataset from Amplitude Variations with Offfset Summary Subhendu Dutta, ONGC The amplitude variations with offset (AVO) is a very useful tool which compliments the
More informationReverse time migration of multiples: Applications and challenges
Reverse time migration of multiples: Applications and challenges Zhiping Yang 1, Jeshurun Hembd 1, Hui Chen 1, and Jing Yang 1 Abstract Marine seismic acquisitions record both primary and multiple wavefields.
More informationMaximizing the value of the existing seismic data in Awali field Bahrain, by utilizing advanced 3D processing technology.
Maximizing the value of the existing seismic data in Awali field Bahrain, by utilizing advanced 3D processing technology. Eduard Maili* (OXY - Tatweer), Scott Burns (OXY), Neil Jones (Consultant, OXY)
More informationThe determination of the correct
SPECIAL High-performance SECTION: H i gh-performance computing computing MARK NOBLE, Mines ParisTech PHILIPPE THIERRY, Intel CEDRIC TAILLANDIER, CGGVeritas (formerly Mines ParisTech) HENRI CALANDRA, Total
More information