Youzuo Lin and Lianjie Huang
|
|
- Michael Gardner
- 5 years ago
- Views:
Transcription
1 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 and Steeply-Dipping Fault Zones Using Elastic-Waveform Inversion with Edge-Guided Regularization Youzuo Lin and Lianjie Huang Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA s: (YL); (LH) Keywords: Edge information, Elastic-waveform inversion, Modified total-variation regularization, Regularization techniques, Velocity estimation. ABSTRACT Imaging correct locations of fault zones is important for geothermal exploration and optimizing enhanced geothermal systems. Elastic reverse-time migration has the potential to directly image steeply-dipping fault zones. However, it requires both accurate compressional and shear wave velocity models. Elastic-waveform inversion is a promising tool for velocity estimation. Because of the ill-posedness and complexity of elastic-waveform inversion, it becomes much more challenging to accurately obtain velocity estimation compared to acoustic-waveform inversion, particularly in the deep regions and fault zones. To improve velocity estimation, we develop an elastic-waveform inversion method with an edge-guided modified total-variation regularization scheme to improve the inversion accuracy and robustness, particularly for steeply-dipping fault zones with widths of much smaller than the seismic wavelength. The new regularization scheme incorporates the edge information into waveform inversion. The edge information of subsurface interfaces is obtained iteratively using migration imaging during elastic-waveform inversion. Our new elastic-waveform inversion takes advantage of the robustness of migration imaging to improve velocity estimation. We validate the improved capability of our new elastic-waveform inversion method using synthetic seismic data for a complex model containing several steeply-dipping fault zones. Our inversion results are much better than those produced without using edge-guided regularization. Elastic-waveform inversion with an edge-guided modified total-variation regularization scheme has the potential to accurately estimate both the compressional- and shear-wave velocities within steeply-dipping fault zones, which would significantly improve direct imaging of these faults. 1. INTRODUCTION Direct imaging steeply-dipping fault zones is important because these fault zones may provide paths for hydrothermal flow or confine the boundaries of geothermal reservoirs. Huang et al. (2011) demonstrated that reverse-time migration has the potential to image steeply-dipping fault zones. The method requires an accurate velocity model for direct imaging of steeply-dipping fault zones. Full-waveform inversion (FWI) is a promising tool for velocity estimation of subsurface structures (Tarantola 1984, Mora 1987, Pratt et al. 1998, Sirgue and Pratt 2004). However, one of the difficulties with FWI is the convergence to the local minima that makes the technique very sensitive to the starting velocity model. Besides the local minima issue, the computational complexity in solving an elastic-waveform inversion problem becomes much more expensive than the acoustic-waveform inversion. Regularization techniques are often used to alleviate the ill-posedness problem of waveform inversion caused by the limited data coverage (Hu et al. 2009, Burstedde and Ghattas 2009, Ramirez and Lewis 2010, Guitton 2011). Tomographic inversion usually gives smoothed image reconstructions. Dewaraja et al. (2010) used the boundary information to improve the reconstruction quality for SPECT imaging. Guo and Yin (2012) encoded the edge information in MRI reconstructions, which also yields better imaging results. Baritaux and Unser (2010) applied the edge information to the Fluorescence Diffuse Optical Tomography, and obtained enhanced reconstructions compared to results yielded without using the edge information. We have recently employed an edge-guided regularization scheme and reverse-time migration (RTM) in acoustic-waveform (or full-waveform) inversion to obtain accurate and high-resolution velocity reconstructions of compressional-wave velocities in complex geophysical models containing sharp interfaces and steeply-dipping fault zones (Lin and Huang 2013). Elastic-waveform inversion is needed for estimating both compressional- and shear-wave velocities using multi-component seismic data. We develop an elastic-waveform inversion (EWI) method with edge-guided regularization to improve the accuracy of reconstructions of both compressional- and shear-wave velocities using multi-component seismic data. The edge-guided regularization can be combined with any regularization techniques. In this paper, we study the performance of elastic-waveform inversion with the edge-guided regularization and a modified total-variation regularization scheme. The method uses the interface information from elastic-wave reverse-time migration to regularize elastic-waveform inversion. Combining with the modified totalvariation regularization scheme, our new elastic-waveform inversion is able to reconstruct velocities for complex geophysical models containing sharp interfaces and steeply-dipping fault zones. We use synthetic multi-component seismic reflection data to validate the improved capability of our new elastic-waveform inversion for estimating both the compressional-wave and shear-wave velocities of a complex geophysical model from the Soda Lake geothermal field. The geophysical model is constructed using a prestack depth migration image at the Soda Lake geothermal field. The model contains several steeply-dipping fault zones with widths of much smaller than the seismic wavelength. Our results demonstrate that our new elastic-waveform inversion method produces much more accurate estimation of both compressional- and shear-wave velocities in the regions with sharp interfaces and steeply-dipping fault zones compared to those obtained without using the edge-guided regularization technique. 1
2 2. ELASTIC-WAVEFORM INVERSION The elastic-wave equation in the time domain is given by 2 U T ρ ( K U ) ( μ ( U U )) s, (1) 2 t where is the density, K and μ are the elastic moduli, s is the source term and U is the displacement wavefield. The solution to the above elastic wave equation can be written as U f ( K, μ, ρ, s ) (2) where the function of f is the wave-propagation operator. Numerical techniques, such as finite-difference and spectral element methods, can be used to solve for the forward problem. Let m be the model parameter, Eq. (2) becomes U f ( m ). (3) To invert for the model parameter m, we need to solve the minimization problem 2 E ( m ) m in p f ( m ) 2, (4) m where E ( m ) is the misfit function, 2 stands for the L 2 norm, and p represents recorded multi-component waveforms. The resulting model m minimizes the square difference between observed and synthetic waveforms. 3. ELASTIC-WAVEFORM INVERSION WITH EDGE-GUIDED MODIFIED TV REGULARIZATION The cost function with a modified total-variation (TV) regularization scheme is given by Huang et al. (2008) E ( m, u ) m in{ p f ( m ) 2 1 m u u }, (5) T V m, u where 1 and 2 are both positive regularization parameters, u is an auxiliary vector with a dimension of m, and the TV term u for a 2D model is defined as the L 1 norm given by TV u T V u 1 ( x u ) i, j ( z u ) i, j, (6) i, j with ( x u ) i, j u i 1, j u and i, j ( z u ) i, j u i, j 1 u i, j. To incorporate the edge information, we reformulate the TV term given by Eq. (6) as u E T V w u 1 wi, j ( ( x u ) i, j ( z u ) i, j ), (7) i, j where the weighting parameter w controls the amount of regularization among adjacent spatial grid points. We set the weighting value as the following: 0 if p o in t (i, j) is o n th e e d g e w i, j. 1 if p o in t (i, j) is n o t o n th e e d g e (8) We then obtain elastic-waveform inversion with the edge-guided modified TV regularization scheme given by E ( m, u ) m in{ p f ( m ) 2 1 m u w u }. (9) 1 m, u We rewrite elastic-waveform inversion with the edge-guided modified TV regularization in Eq. (9) as E ( m, u ) m in { m in { p f ( m ) 2 1 m u } u }. (10) E T V u m 4. EDGE DETECTION During each iteration step of elastic-waveform inversion, we compute forward-propagation wavefields from sources and backwardpropagation wavefields from receivers. Therefore, we exploit these wavefields to obtain the edges of heterogeneities or interfaces of subsurface structures using elastic reverse-time migration. Consequently, we gain the edge information during elastic-waveform inversion. After the edges are determined, we apply the weighting coefficients according to Eq. (8). 2
3 5. NUMERICAL RESULTS We use synthetic surface seismic data for a model with the compressional and shear velocities in Figs. 1a and 1b to demonstrate the improvement of our new elastic-waveform inversion method with the edge-guided modified TV regularization scheme for velocity estimation. The compressional and shear models are constructed using geologic features found at the Soda Lake geothermal field. It contains several steeply-dipping fault zones. There are three basalt regions in Fig. 1 with a high velocity value. A total of 165 common-shot gathers of synthetic seismic data with 720 receivers at the top surface of the model are used to invert for velocity values of the model. The shot interval is 20 m and the receiver interval is 5 m. A Ricker wavelet with a center frequency of 25 Hz is used as the source time function. We smooth the original velocity models in Figs. 1a and 1b by averaging the slowness within two wavelengths, resulting in a model in Figs. 1c and 1d. We use these two smoothed models as the starting models for EWI reconstructions. Figure 2 shows EWI reconstructions with no regularization (Figs. 2a and 2b) and with edge-guided modified TV regularization (Figs. 2c and 2d). The EWI results obtained with the edge-guided modified TV regularization scheme in Figs. 2c and 2d shows accurate velocity reconstruction, particularly in the deep regions of the models that are poorly reconstructed without using any regularization in Figs. 2a and 2b. In addition, our new method greatly improves reconstructions of the steeply-dipping fault zones. For comparison of the velocity values estimated using the two EWI methods, we display two profiles: a horizontal profile and a vertical profile. The horizontal profile cuts through two basalt regions and all five fault zones. The vertical profile is along the center of the largest basalt region. Figure 3 displays the horizontal profiles of the two EWI reconstructions along the depth of 0.7 km. Both methods reconstruct the locations of all five fault zones. EWI with edge-guided modified TV regularization significantly improves the reconstructed velocity values of the fault zones compared to EWI without regularization. The true values of the fault zones are 2125 m/s for compressional velocity and 1062 m/s for shear velocity. The estimated velocity values of EWI with edge-guided modified TV regularization are approximately 2180 m/s for the compressional velocity and 1068 m/s for the shear velocity. However the approximated reconstructed compressional and shear velocity values of EWI without regularization are far away from the true values. Figure 4 shows a comparison of the vertical profiles of EWI reconstructions in Fig. 2. EWI without regularization produces oscillated profiles in both compressional and shear velocity results. By contrast, EWI with edge-guided modified TV regularization accurately reconstructs both compressional and shear velocity values and the interfaces. 6. CONCLUSIONS We have developed a novel elastic-waveform inversion method with edge-guided modified total-variation regularization. The method employs the edge information in combination with a modified total-variation regularization scheme. We have validated the capability of our new elastic-waveform inversion method for accurate reconstructions of velocities for steeply-dipping fault zones and the deep regions of the model. Our elastic-waveform inversion results of synthetic seismic data for a Soda Lake geothermal velocity model demonstrate that our new method can accurately reconstruct not only velocity values but also sharp interfaces. Our novel elastic-waveform inversion method with edge-guided modified total-variation regularization could be a powerful tool for accurate velocity estimation, which could lead to significantly improved elastic reverse-time migration for direct imaging of steeply-dipping fault zones. ACKNOWLEDGEMENTS This work was supported by the Geothermal Technologies Program of the U.S. Department of Energy through contract DE-AC52-06NA25396 to Los Alamos National Laboratory (LANL). The computation was performed using super-computers of LANL's Institutional Computing Program. We thank Magma Energy (U.S.) Corp. for providing us with migration images and James Echols for his help in constructing the velocity model for the Soda Lake geothermal field. 3
4 (a) True compressional-wave velocity model (b) True shear-wave velocity model (c) Initial compressional-wave velocity model (d) Initial shear-wave velocity model Figure 1: Velocity models used for EWI: (a) compressional velocity model and (b) shear velocity model from the Soda Lake geothermal field for generating synthetic seismic data; (c) smoothed compressional velocity model and (d) smoothed shear velocity model used as the starting models for EWI. 4
5 (a) Reconstruction of compressional velocity (b) Reconstruction of shear velocity (c) Reconstruction of compressional velocity (d) Reconstruction of shear velocity Figure 2: Reconstructions of the velocity model in Fig. 1 produced using EWI without regularization (a, b) and with edgeguided modified TV regularization (c, d). EWI without using regularization yields poor velocity reconstructions. EWI with edge-guided modified TV regularization produces significantly improved velocity reconstructions, particularly in the deep regions of the model and within the steeply-dipping fault zones. 5
6 (a) Reconstruction of compressional velocity model (b) Reconstruction of shear velocity model (c) Reconstruction of compressional velocity model (d) Reconstruction of shear velocity model Figure 3: Comparison of horizontal profiles of EWI reconstructions in Figs. 2 along the depth of 0.7 km. In each panel, the red line shows the true velocity value, the green line is the initial guess, and the blue line is the EWI reconstruction result. EWI without regularization in (a) and (b) yields oscillated profiles. EWI with edge-guided modified TV regularization in (c) and (d) yields both accurate compressional and shear velocity values. The true compressional velocity value within the fault zones is 2125 m/s and the true shear velocity value is 1062 m/s, while the velocity values of EWI with edge-guided modified TV regularization in (c) and (d) are approximately 2180 m/s and 1068 m/s. By contrast, the reconstructed velocity values in (a) and (b) are far away from the true values. 6
7 (a) Reconstruction of compressional velocity model (b) Reconstruction of shear velocity model (c) Reconstruction of compressional velocity model (d) Reconstruction of shear velocity model Figure 4: Comparison of vertical profiles of EWI reconstructions in Figs. 2 along the horizontal position of 1.5 km. In each panel, the red line shows the true velocity value, the green line is the initial guess, and the blue line is the EWI reconstruction result. EWI without regularization in (a) and (b) yields oscillated profiles. EWI with edge-guided modified TV regularization displayed in (c) and (d) accurately reconstructs both velocity values and the interfaces. 7
8 REFERENCES Baritaux, J.-C., and Unser, M.: A priori guided reconstruction for FDOT using mixed norms, Proceedings, IEEE International Symposium on Biomedical Imaging: From Nano to Macro, (2010). Burstedde, C., and Ghattas, O.: Algorithmic strategies for full waveform inversion: 1d experiments, Geophysics, 74, (2009), Dewaraja, Y. K., Koral, K. F. and Fessler, J. A.: Regularized reconstruction in quantitative SPECT using CT side information from hybrid imaging, Physics in Medicine and Biology, (2010), Guitton, A.: A blocky regularization scheme for full waveform inversion, Expanded Abstracts, 81th SEG Annual Meeting, (2011), Guo,W., and Yin, W. : Edge guided reconstruction for compressive imaging, SIAM Journal on Imaging Sciences, 5, (2012), Hu, W., Abubakar, A. and Habashy, T.: Simultaneous multi-frequency inversion of full-waveform seismic data, Geophysics, 74, (2009), Huang, L., Kelley, S., Zhang, Z., Rehfeldt, K., Albrecht, M., and Kaufman, G.: Imaging faults with reverse-time migration for geothermal exploration at Jemez Pueblo in New Mexico, Transactions, Geothermal Resources Council, 35, (2011), Huang, Y., NG, M. and Wen, Y. : A fast total variation minimization method for image restoration, SIAM Journal on Multiscale Modeling and Simulation, 7, (2008), Lin, Y., and Huang, L.: Accurate velocity estimation for steep fault zones using full-waveform inversion with edge-guided regularization, Transactions, Geothermal Resources Council, (2013), Mora, P.: Nonlinear two-dimensional elastic inversion of multi-offset seismic data, Geophysics, 54, (1987), Pratt, R. G., Shin, C. and Hicks, G. J.: Gauss-newton and full newton methods in frequency space seismic waveform inversion, Geophysical Journal International, 13, (1998), Ramirez, A., and Lewis, W.: Regularization and full-waveform inversion: A two-step approach, Expanded Abstracts, 80th SEG Annual Meeting, (2010), Sirgue, L., and Pratt, R. G.: Efficient waveform inversion and imaging: A strategy for selecting temporal frequency, Geophysics, 69, (2004), Tarantola, A.: Inversion of seismic reflection data in the acoustic approximation, Geophysics, 49, (1984),
Unravel Faults on Seismic Migration Images Using Structure-Oriented, Fault-Preserving and Nonlinear Anisotropic Diffusion Filtering
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Unravel Faults on Seismic Migration Images Using Structure-Oriented,
More informationAnisotropic Seismic Imaging and Inversion for Subsurface Characterization at the Blue Mountain Geothermal Field in Nevada
PROCEEDINGS, 43rd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 12-14, 2018 SGP-TR-213 Anisotropic Seismic Imaging and Inversion for Subsurface Characterization
More informationW011 Full Waveform Inversion for Detailed Velocity Model Building
W011 Full Waveform Inversion for Detailed Velocity Model Building S. Kapoor* (WesternGeco, LLC), D. Vigh (WesternGeco), H. Li (WesternGeco) & D. Derharoutian (WesternGeco) SUMMARY An accurate earth model
More informationGeothermal Reservoir Imaging Using 2016 Walkaway VSP Data from the Raft River Geothermal Field
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Geothermal Reservoir Imaging Using 2016 Walkaway VSP Data from
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 informationFull-waveform inversion application in different geological settings Denes Vigh*, Jerry Kapoor and Hongyan Li, WesternGeco
Full-waveform inversion application in different geological settings Denes Vigh*, Jerry Kapoor and Hongyan Li, WesternGeco Summary After the synthetic data inversion examples, real 3D data sets have been
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 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 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 information2D Laplace-Domain Waveform Inversion of Field Data Using a Power Objective Function
Pure Appl. Geophys. Ó 213 Springer Basel DOI 1.17/s24-13-651-4 Pure and Applied Geophysics 2D Laplace-Domain Waveform Inversion of Field Data Using a Power Objective Function EUNJIN PARK, 1 WANSOO HA,
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 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 informationWe Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter
We-09-04 Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter J. Giraud* (WesternGeco Geosolutions), M. De Stefano (WesternGeco Geosolutions)
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 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 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 informationMigration with Implicit Solvers for the Time-harmonic Helmholtz
Migration with Implicit Solvers for the Time-harmonic Helmholtz Yogi A. Erlangga, Felix J. Herrmann Seismic Laboratory for Imaging and Modeling, The University of British Columbia {yerlangga,fherrmann}@eos.ubc.ca
More informationFWI with Compressive Updates Aleksandr Aravkin, Felix Herrmann, Tristan van Leeuwen, Xiang Li, James Burke
Consortium 2010 FWI with Compressive Updates Aleksandr Aravkin, Felix Herrmann, Tristan van Leeuwen, Xiang Li, James Burke SLIM University of British Columbia Full Waveform Inversion The Full Waveform
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 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 informationTOM 1.7. Sparse Norm Reflection Tomography for Handling Velocity Ambiguities
SEG/Houston 2005 Annual Meeting 2554 Yonadav Sudman, Paradigm and Dan Kosloff, Tel-Aviv University and Paradigm Summary Reflection seismology with the normal range of offsets encountered in seismic surveys
More informationTowards full waveform inversion: A torturous path
FWI Towards full waveform inversion: A torturous path J. Helen Isaac and Gary F. Margrave ABSTRACT Full waveform inversion (FWI) can be viewed as an iterative cycle involving forward modelling, pre-stack
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 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 informationA projected Hessian for full waveform inversion
CWP-679 A projected Hessian for full waveform inversion Yong Ma & Dave Hale Center for Wave Phenomena, Colorado School of Mines, Golden, CO 80401, USA (c) Figure 1. Update directions for one iteration
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 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 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. (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 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 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 informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 188, 1221 1242 doi: 10.1111/j.1365-246X.2011.05314.x Full waveform inversion strategy for density in the frequency domain Woodon Jeong, 1 Ho-Yong
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 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 informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Inverse Scattering Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 30, 2014 Bee Bednar (Panorama Technologies)
More informationSeismic Imaging of Supercritical Geothermal Reservoir Using Full-waveform Inversion Method
PROCEEDINGS, 44th Workshop on Geothermal Reservoir EngineeringPROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Seismic
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 informationElastic least-squares migration with two-way wave equation forward and adjoint operators
Elastic least-squares migration with two-way wave equation forward and adjoint operators Ke Chen and Mauricio D. Sacchi, Department of Physics, University of Alberta Summary Time domain elastic least-squares
More informationTIME-DEPENDENT SEISMIC TOMOGRAPHY OF GEOTHERMAL SYSTEMS
PROCEEDINGS, Thirty-Fourth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 9-11, 2009 SGP-TR-187 TIME-DEPENDENT SEISMIC TOMOGRAPHY OF GEOTHERMAL SYSTEMS
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 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 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 information3D VSP Experiment Design Study of the Raft River Geothermal Field
PROCEEDINGS, 41st Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 22-24, 2016 SGP-TR-209 3D VSP Experiment Design Study of the Raft River Geothermal Field
More informationA new second order absorbing boundary layer formulation for anisotropic-elastic wavefeld simulation
A new second order absorbing boundary layer formulation for anisotropic-elastic wavefeld simulation Junxiao Li, Kris Innanen and Bing Wang University of Calgary China University of Petroleum-Beijing Summary
More informationUncertainty quantification for inverse problems with a weak wave-equation constraint
Uncertainty quantification for inverse problems with a weak wave-equation constraint Zhilong Fang*, Curt Da Silva*, Rachel Kuske** and Felix J. Herrmann* *Seismic Laboratory for Imaging and Modeling (SLIM),
More informationMultiple Scenario Inversion of Reflection Seismic Prestack Data
Downloaded from orbit.dtu.dk on: Nov 28, 2018 Multiple Scenario Inversion of Reflection Seismic Prestack Data Hansen, Thomas Mejer; Cordua, Knud Skou; Mosegaard, Klaus Publication date: 2013 Document Version
More informationFull waveform inversion with wave equation migration and well control
FWI by WEM Full waveform inversion with wave equation migration and well control Gary F. Margrave, Robert J. Ferguson and Chad M. Hogan ABSTRACT We examine the key concepts in full waveform inversion (FWI)
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 informationUncertainty quantification for Wavefield Reconstruction Inversion
Uncertainty quantification for Wavefield Reconstruction Inversion Zhilong Fang *, Chia Ying Lee, Curt Da Silva *, Felix J. Herrmann *, and Rachel Kuske * Seismic Laboratory for Imaging and Modeling (SLIM),
More informationSparsity-promoting migration with multiples
Sparsity-promoting migration with multiples Tim Lin, Ning Tu and Felix Herrmann SLIM Seismic Laboratory for Imaging and Modeling the University of British Columbia Courtesy of Verschuur, 29 SLIM Motivation..
More informationInterval anisotropic parameters estimation in a least squares sense Case histories from West Africa
P-263 Summary Interval anisotropic parameters estimation in a least squares sense Patrizia Cibin*, Maurizio Ferla Eni E&P Division (Milano, Italy), Emmanuel Spadavecchia - Politecnico di Milano (Milano,
More 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 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 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 informationFull-Waveform Inversion with Gauss- Newton-Krylov Method
Full-Waveform Inversion with Gauss- Newton-Krylov Method Yogi A. Erlangga and Felix J. Herrmann {yerlangga,fherrmann}@eos.ubc.ca Seismic Laboratory for Imaging and Modeling The University of British Columbia
More informationP191 Bayesian Linearized AVAZ Inversion in HTI Fractured Media
P9 Bayesian Linearized AAZ Inversion in HI Fractured Media L. Zhao* (University of Houston), J. Geng (ongji University), D. Han (University of Houston) & M. Nasser (Maersk Oil Houston Inc.) SUMMARY A new
More informationNonlinear acoustic imaging via reduced order model backprojection
Nonlinear acoustic imaging via reduced order model backprojection Alexander V. Mamonov 1, Vladimir Druskin 2, Andrew Thaler 3 and Mikhail Zaslavsky 2 1 University of Houston, 2 Schlumberger-Doll Research
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 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 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 informationElastic versus acoustic inversion for marine surveys
Elastic versus acoustic inversion for marine surveys Item Type Article Authors Mora, Peter; Wu, Zedong Citation Mora P, Wu Z (2018) Elastic versus acoustic inversion for marine surveys. Geophysical Journal
More information11th Biennial International Conference & Exposition. Keywords Sub-basalt imaging, Staggered grid; Elastic finite-difference, Full-waveform modeling.
Sub-basalt imaging using full-wave elastic finite-difference modeling: A synthetic study in the Deccan basalt covered region of India. Karabi Talukdar* and Laxmidhar Behera, CSIR-National Geophysical Research
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 informationDownloaded 11/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
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
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 informationImaging complex structure with crosswell seismic in Jianghan oil field
INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Imaging complex structure with crosswell seismic in Jianghan oil field QICHENG DONG and BRUCE MARION, Z-Seis, Houston, Texas, U.S. JEFF MEYER, Fusion
More informationTime domain sparsity promoting LSRTM with source estimation
Time domain sparsity promoting LSRTM with source estimation Mengmeng Yang, Philipp Witte, Zhilong Fang & Felix J. Herrmann SLIM University of British Columbia Motivation Features of RTM: pros - no dip
More informationHigh Resolution Imaging of Fault Zone Properties
Annual Report on 1998-99 Studies, Southern California Earthquake Center High Resolution Imaging of Fault Zone Properties Yehuda Ben-Zion Department of Earth Sciences, University of Southern California
More informationOne-step extrapolation method for reverse time migration
GEOPHYSICS VOL. 74 NO. 4 JULY-AUGUST 2009; P. A29 A33 5 FIGS. 10.1190/1.3123476 One-step extrapolation method for reverse time migration Yu Zhang 1 and Guanquan Zhang 2 ABSTRACT We have proposed a new
More informationSummary. We present the results of the near-surface characterization for a 3D survey in thrust belt area in Sharjah, United Arab Emirates.
Near-surface characterization, challenges, and solutions for high-density, high-productivity, Alexander Zarkhidze*, Claudio Strobbia, Abdallah Ibrahim, WesternGeco; Luis Viertel Herrera, Abdulla Al Qadi,
More informationSeismic imaging and optimal transport
Seismic imaging and optimal transport Bjorn Engquist In collaboration with Brittany Froese, Sergey Fomel and Yunan Yang Brenier60, Calculus of Variations and Optimal Transportation, Paris, January 10-13,
More informationWe E Multiparameter Full-waveform Inversion for Acoustic VTI Medium with Surface Seismic Data
We E16 4 Multiarameter Full-waveform Inversion for Acoustic VI Medium with Surface Seismic Data X. Cheng* (Schlumberger) K. Jiao (Schlumberger) D. Sun (Schlumberger) & D. Vigh (Schlumberger) SUMMARY In
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 informationPART A: Short-answer questions (50%; each worth 2%)
PART A: Short-answer questions (50%; each worth 2%) Your answers should be brief (just a few words) and may be written on these pages if you wish. Remember to hand these pages in with your other exam pages!
More informationMain 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 informationTwo-Scale Wave Equation Modeling for Seismic Inversion
Two-Scale Wave Equation Modeling for Seismic Inversion Susan E. Minkoff Department of Mathematics and Statistics University of Maryland Baltimore County Baltimore, MD 21250, USA RICAM Workshop 3: Wave
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 informationNOTICE CONCERNING COPYRIGHT RESTRICTIONS
NOTICE CONCERNING COPYRIGHT RESTRICTIONS This document may contain copyrighted materials. These materials have been made available for use in research, teaching, and private study, but may not be used
More informationThe Coso Geothermal Area: A Laboratory for Advanced MEQ Studies for Geothermal Monitoring
The Coso Geothermal Area: A Laboratory for Advanced MEQ Studies for Geothermal Monitoring Bruce R. Julian U. S. Geological Survey, Menlo Park, CA 94025 USA julian@usgs.gov Gillian R. Foulger Dept. Earth
More informationAnisotropic anelastic seismic full waveform modeling and inversion: Application to North Sea offset VSP data
CONFERENCE PROCEEDINGS DOI: https://doi.org/10.18599/grs.2018.3.149-153 Anisotropic anelastic seismic full waveform modeling and inversion: Application to North Sea offset VSP data C. Barnes 1, M. Charara
More informationTh P7 02 A Method to Suppress Salt-related Converted Wave Using 3D Acoustic Modelling
Th P7 0 A Method to Suppress Salt-related Converted Wave Using 3D Acoustic Modelling J. Kumar* (Petroleum Geo-Services), M. Salem (ENI E&P), D.E. Cegani (ENI E&P) Summary Converted waves can be recorded
More informationData examples of logarithm Fourier domain bidirectional deconvolution
Data examples of logarithm Fourier domain bidirectional deconvolution Qiang Fu and Yi Shen and Jon Claerbout ABSTRACT Time domain bidirectional deconvolution methods show a very promising perspective on
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 informationTowards Modelling Elastic and Viscoelastic Seismic Wave Propagation in Boreholes
Towards Modelling Elastic and Viscoelastic Seismic Wave Propagation in Boreholes NA WANG, DONG SHI, BERND MILKEREIT Department of Physics, University of Toronto, Toronto, Canada M5S 1A7 Summary We are
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 informationF-K Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures
F-K Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures Yang Zhang 1, Xander Campman 1, Samantha Grandi 1, Shihong Chi 1, M. Nafi Toksöz 1, Mark E. Willis 1, Daniel R. Burns
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 informationTh P10 13 Alternative Misfit Functions for FWI Applied to Surface Waves
Th P0 3 Alternative Misfit Functions for FWI Alied to Surface Waves I. Masoni* Total E&P, Joseh Fourier University), R. Brossier Joseh Fourier University), J. Virieux Joseh Fourier University) & J.L. Boelle
More informationDownloaded 05/01/17 to Redistribution subject to SEG license or copyright; see Terms of Use at
Mapping Imbricate Structures in the Thrust Belt of Southeast Turkey by Large-Offset Seismic Survey Oz Yilmaz*, Anatolian Geophysical, Istanbul, Turkey; and Serdar Uygun, Ali Ölmez, and Emel Çalı, Turkish
More informationImproving subsurface imaging in geological complex area: Structure PReserving INTerpolation in 6D (SPRINT6D)
Improving subsurface imaging in geological complex area: Structure PReserving INTerpolation in 6D (SPRINT6D) Dan Negut* and Mark Ng, Divestco Inc. Copyright 2008, ACGGP This paper was selected for presentation
More informationAPPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION
METHODS AND APPLICATIONS OF ANALYSIS. c 2004 International Press Vol. 11, No. 3, pp. 345 352, September 2004 004 APPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION PING SHENG, GANG SUN,
More informationFinite difference elastic modeling of the topography and the weathering layer
Finite difference elastic modeling of the topography and the weathering layer Saul E. Guevara and Gary F. Margrave ABSTRACT Finite difference 2D elastic modeling is used to study characteristics of the
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 informationUncertainty quantification for Wavefield Reconstruction Inversion
using a PDE free semidefinite Hessian and randomize-then-optimize method Zhilong Fang *, Chia Ying Lee, Curt Da Silva, Tristan van Leeuwen and Felix J. Herrmann Seismic Laboratory for Imaging and Modeling
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 informationSeismic modeling evaluation of fault illumination in the Woodford Shale Sumit Verma*, Onur Mutlu, Kurt J. Marfurt, The University of Oklahoma
Seismic modeling evaluation of fault illumination in the Woodford Shale Sumit Verma*, Onur Mutlu, Kurt J. Marfurt, The University of Oklahoma Summary The Woodford Shale is one of the more important resource
More informationA nonlinear differential-semblance algorithm for waveform inversion
A nonlinear differential-semblance algorithm for waveform inversion Dong Sun The Rice Inversion Project Annual Review Meeting March 3, 212, 1/15 Agenda 1 Full Waveform Inversion 2 Extension and Differential
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 informationWavefield Reconstruction Inversion (WRI) a new take on wave-equation based inversion Felix J. Herrmann
Wavefield Reconstruction Inversion (WRI) a new take on wave-equation based inversion Felix J. Herrmann SLIM University of British Columbia van Leeuwen, T and Herrmann, F J (2013). Mitigating local minima
More informationWave-equation Hessian by phase encoding
Wave-equation Hessian by phase encoding Yaxun Tang ABSTRACT I present a method for computing wave-equation Hessian operators, also known as resolution functions or point-spread functions, under the Born
More information