Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform Inversion
In the presented study, multi-parameter inversion in the presence of attenuation is used for the reconstruction of the P- and the S- wave velocities and the density models of a synthetic shallow subsurface structure that contains a dipping high-velocity layer near the surface with varying thicknesse...
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
MDPI AG
2022-08-01
|
| Series: | Applied Sciences |
| Subjects: | |
| Online Access: | https://www.mdpi.com/2076-3417/12/15/7741 |
| _version_ | 1797442792337702912 |
|---|---|
| author | Niloofar Alaei Mehrdad Soleimani Monfared Amin Roshandel Kahoo Thomas Bohlen |
| author_facet | Niloofar Alaei Mehrdad Soleimani Monfared Amin Roshandel Kahoo Thomas Bohlen |
| author_sort | Niloofar Alaei |
| collection | DOAJ |
| description | In the presented study, multi-parameter inversion in the presence of attenuation is used for the reconstruction of the P- and the S- wave velocities and the density models of a synthetic shallow subsurface structure that contains a dipping high-velocity layer near the surface with varying thicknesses. The problem of high-velocity layers also complicates selection of an appropriate initial velocity model. The forward problem is solved with the finite difference, and the inverse problem is solved with the preconditioned conjugate gradient. We used also the adjoint wavefield approach for computing the gradient of the misfit function without explicitly build the sensitivity matrix. The proposed method is capable of either minimizing the least-squares norm of the data misfit or use the Born approximation for estimating partial derivative wavefields. It depends on which characteristics of the recorded data—such as amplitude, phase, logarithm of the complex-valued data, envelope in the misfit, or the linearization procedure of the inverse problem—are used. It showed that by a pseudo-viscoelastic time-domain full-waveform inversion, structures below the high-velocity layer can be imaged. However, by inverting attenuation of P- and S- waves simultaneously with the velocities and mass density, better results would be obtained. |
| first_indexed | 2024-03-09T12:47:13Z |
| format | Article |
| id | doaj.art-bb9e58c2ae294e5fb1210f6c10c85feb |
| institution | Directory Open Access Journal |
| issn | 2076-3417 |
| language | English |
| last_indexed | 2024-03-09T12:47:13Z |
| publishDate | 2022-08-01 |
| publisher | MDPI AG |
| record_format | Article |
| series | Applied Sciences |
| spelling | doaj.art-bb9e58c2ae294e5fb1210f6c10c85feb2023-11-30T22:11:00ZengMDPI AGApplied Sciences2076-34172022-08-011215774110.3390/app12157741Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform InversionNiloofar Alaei0Mehrdad Soleimani Monfared1Amin Roshandel Kahoo2Thomas Bohlen3Department of Mining Petroleum and Geophysics Engineering, Shahrood University of Technology, Shahrood 3619995161, IranDepartment of Mining Petroleum and Geophysics Engineering, Shahrood University of Technology, Shahrood 3619995161, IranDepartment of Mining Petroleum and Geophysics Engineering, Shahrood University of Technology, Shahrood 3619995161, IranGeophysical Institute (GPI), Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, GermanyIn the presented study, multi-parameter inversion in the presence of attenuation is used for the reconstruction of the P- and the S- wave velocities and the density models of a synthetic shallow subsurface structure that contains a dipping high-velocity layer near the surface with varying thicknesses. The problem of high-velocity layers also complicates selection of an appropriate initial velocity model. The forward problem is solved with the finite difference, and the inverse problem is solved with the preconditioned conjugate gradient. We used also the adjoint wavefield approach for computing the gradient of the misfit function without explicitly build the sensitivity matrix. The proposed method is capable of either minimizing the least-squares norm of the data misfit or use the Born approximation for estimating partial derivative wavefields. It depends on which characteristics of the recorded data—such as amplitude, phase, logarithm of the complex-valued data, envelope in the misfit, or the linearization procedure of the inverse problem—are used. It showed that by a pseudo-viscoelastic time-domain full-waveform inversion, structures below the high-velocity layer can be imaged. However, by inverting attenuation of P- and S- waves simultaneously with the velocities and mass density, better results would be obtained.https://www.mdpi.com/2076-3417/12/15/7741complex velocity modelfull waveform inversionwave attenuationpreconditioned conjugate gradientvibroseis sources |
| spellingShingle | Niloofar Alaei Mehrdad Soleimani Monfared Amin Roshandel Kahoo Thomas Bohlen Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform Inversion Applied Sciences complex velocity model full waveform inversion wave attenuation preconditioned conjugate gradient vibroseis sources |
| title | Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform Inversion |
| title_full | Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform Inversion |
| title_fullStr | Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform Inversion |
| title_full_unstemmed | Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform Inversion |
| title_short | Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform Inversion |
| title_sort | seismic imaging of complex velocity structures by 2d pseudo viscoelastic time domain full waveform inversion |
| topic | complex velocity model full waveform inversion wave attenuation preconditioned conjugate gradient vibroseis sources |
| url | https://www.mdpi.com/2076-3417/12/15/7741 |
| work_keys_str_mv | AT niloofaralaei seismicimagingofcomplexvelocitystructuresby2dpseudoviscoelastictimedomainfullwaveforminversion AT mehrdadsoleimanimonfared seismicimagingofcomplexvelocitystructuresby2dpseudoviscoelastictimedomainfullwaveforminversion AT aminroshandelkahoo seismicimagingofcomplexvelocitystructuresby2dpseudoviscoelastictimedomainfullwaveforminversion AT thomasbohlen seismicimagingofcomplexvelocitystructuresby2dpseudoviscoelastictimedomainfullwaveforminversion |