Multiscale Seismic Inversion in the Data and Image Domains

Handle URI:
http://hdl.handle.net/10754/621206
Title:
Multiscale Seismic Inversion in the Data and Image Domains
Authors:
Zhang, Sanzong ( 0000-0002-3972-073X )
Abstract:
I present a general methodology for inverting seismic data in either the data or image domains. It partially overcomes one of the most serious problems with current waveform inversion methods, which is the tendency to converge to models far from the actual one. The key idea is to develop a multiscale misfit function that is composed of both a simplified version of the data and one associated with the complex part of the data. Misfit functions based on simple data are characterized by many fewer local minima so that a gradient optimization method can make quick progress in getting to the general vicinity of the actual model. Once we are near the actual model, we then use the gradient based on the more complex data. Below, we describe two implementations of this multiscale strategy: wave equation traveltime inversion in the data domain and generalized differential semblance optimization in the image domain. • Wave Equation Traveltime Inversion in the Data Domain (WT): The main difficulty with iterative waveform inversion is that it tends to get stuck in local minima associated with the waveform misfit function. To mitigate this problem and avoid the need to fit amplitudes in the data, we present a waveequation method that inverts the traveltimes of reflection events, and so is less prone to the local minima problem. Instead of a waveform misfit function, the penalty function is a crosscorrelation of the downgoing direct wave and the upgoing reflection wave at the trial image point. The time lag which maximizes the crosscorrelation amplitude represents the reflection-traveltime residual that is back-projected along the reflection wavepath to update the velocity. Shot- and angle-domain crosscorrelation functions are introduced to estimate the reflection-traveltime residual by semblance analysis and scanning. In theory, only the traveltime information is inverted and there is no need to precisely fit the amplitudes or assume a high-frequency approximation. Results with both synthetic data and field records reveal both the benefits and limitations of WT. • Generalized Differental Semblance Optimization in the Image Domain (GDSO): We now extend the multiscale physics approach to differential semblance optimization (DSO) in the image domain. That is, we identify the space-lag offset H(x, z, h) in the subsurface-offset domain as an implicit function of velocity. It describes the smoothly varying moveout H(x, z, h) of the migration image m(x, z, h) in the subsurface-offset domain, which is analogous to the smoothly varying traveltime residual ∆τ(x) of a reflection event in a shot gather. The velocity model is found that minimizes the objective function ∑x,z,h H(x, z, h)2m(x, z, h)2, where coherent noise is eliminated everywhere except along the picked curve H(x, z, h). This method is denoted as generalized DSO (GDSO) and mitigates the coherent noise problem with DSO. Numerical examples are presented that empirically demonstrate its effectiveness in providing more accurate velocity models compared to conventional DSO.
Advisors:
Schuster, Gerard T. ( 0000-0001-7532-1587 )
Committee Member:
Hanafy, Sherif; Sun, Shuyu ( 0000-0002-3078-864X ) ; Wu, Ying ( 0000-0002-7919-1107 ) ; Luo, Yi
KAUST Department:
Physical Sciences and Engineering (PSE) Division
Program:
Earth Sciences and Engineering
Issue Date:
Dec-2015
Type:
Dissertation
Appears in Collections:
Dissertations

Full metadata record

DC FieldValue Language
dc.contributor.advisorSchuster, Gerard T.en
dc.contributor.authorZhang, Sanzongen
dc.date.accessioned2016-10-27T05:37:04Z-
dc.date.available2016-10-27T05:37:04Z-
dc.date.issued2015-12-
dc.identifier.urihttp://hdl.handle.net/10754/621206-
dc.description.abstractI present a general methodology for inverting seismic data in either the data or image domains. It partially overcomes one of the most serious problems with current waveform inversion methods, which is the tendency to converge to models far from the actual one. The key idea is to develop a multiscale misfit function that is composed of both a simplified version of the data and one associated with the complex part of the data. Misfit functions based on simple data are characterized by many fewer local minima so that a gradient optimization method can make quick progress in getting to the general vicinity of the actual model. Once we are near the actual model, we then use the gradient based on the more complex data. Below, we describe two implementations of this multiscale strategy: wave equation traveltime inversion in the data domain and generalized differential semblance optimization in the image domain. • Wave Equation Traveltime Inversion in the Data Domain (WT): The main difficulty with iterative waveform inversion is that it tends to get stuck in local minima associated with the waveform misfit function. To mitigate this problem and avoid the need to fit amplitudes in the data, we present a waveequation method that inverts the traveltimes of reflection events, and so is less prone to the local minima problem. Instead of a waveform misfit function, the penalty function is a crosscorrelation of the downgoing direct wave and the upgoing reflection wave at the trial image point. The time lag which maximizes the crosscorrelation amplitude represents the reflection-traveltime residual that is back-projected along the reflection wavepath to update the velocity. Shot- and angle-domain crosscorrelation functions are introduced to estimate the reflection-traveltime residual by semblance analysis and scanning. In theory, only the traveltime information is inverted and there is no need to precisely fit the amplitudes or assume a high-frequency approximation. Results with both synthetic data and field records reveal both the benefits and limitations of WT. • Generalized Differental Semblance Optimization in the Image Domain (GDSO): We now extend the multiscale physics approach to differential semblance optimization (DSO) in the image domain. That is, we identify the space-lag offset H(x, z, h) in the subsurface-offset domain as an implicit function of velocity. It describes the smoothly varying moveout H(x, z, h) of the migration image m(x, z, h) in the subsurface-offset domain, which is analogous to the smoothly varying traveltime residual ∆τ(x) of a reflection event in a shot gather. The velocity model is found that minimizes the objective function ∑x,z,h H(x, z, h)2m(x, z, h)2, where coherent noise is eliminated everywhere except along the picked curve H(x, z, h). This method is denoted as generalized DSO (GDSO) and mitigates the coherent noise problem with DSO. Numerical examples are presented that empirically demonstrate its effectiveness in providing more accurate velocity models compared to conventional DSO.en
dc.language.isoenen
dc.subjectSeismic Inversionen
dc.subjectWave-equation inversionen
dc.subjectTraveltimeen
dc.subjectDSOen
dc.subjectFull wave form inversionen
dc.subjectMulti-scaleen
dc.titleMultiscale Seismic Inversion in the Data and Image Domainsen
dc.typeDissertationen
dc.contributor.departmentPhysical Sciences and Engineering (PSE) Divisionen
thesis.degree.grantorKing Abdullah University of Science and Technologyen_GB
dc.contributor.committeememberHanafy, Sherifen
dc.contributor.committeememberSun, Shuyuen
dc.contributor.committeememberWu, Yingen
dc.contributor.committeememberLuo, Yien
thesis.degree.disciplineEarth Sciences and Engineeringen
thesis.degree.nameDoctor of Philosophyen
dc.person.id102018en
All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.