"uuid","repository link","title","author","contributor","publication year","abstract","subject topic","language","publication type","publisher","isbn","issn","patent","patent status","bibliographic note","access restriction","embargo date","faculty","department","research group","programme","project","coordinates"
"uuid:88951b0d-8b6a-421d-a232-914a620a847e","http://resolver.tudelft.nl/uuid:88951b0d-8b6a-421d-a232-914a620a847e","Surface-wave inversion for a P-velocity profile with a constant depth gradient of the squared slowness","Ponomarenko, A. V. (St. Petersburg State University); Kashtan, B. M. (St. Petersburg State University); Troyan, V. N. (St. Petersburg State University); Mulder, W.A. (TU Delft Applied Geophysics and Petrophysics; Shell Global Solutions International B.V.)","","2017","Surface waves are often used to estimate a near-surface shear-velocity profile. The inverse problem is solved for the locally one-dimensional problem of a set of homogeneous horizontal elastic layers. The result is a set of shear velocities, one for each layer. To obtain a P-wave velocity profile, the P-guided waves should be included in the inversion scheme. As an alternative to a multi-layered model, we consider a simple smooth acoustic constant-density velocity model, which has a negative constant vertical depth gradient of the squared P-wave slowness and is bounded by a free surface at the top and a homogeneous half-space at the bottom. The exact solution involves Airy functions and provides an analytical expression for the dispersion equation. If the v_{s}/v_{p} ratio is sufficiently small, the dispersion curves can be picked from the seismic data and inverted for the continuous P-wave velocity profile. The potential advantages of our model are its low computational cost and the fact that the result can serve as a smooth starting model for full-waveform inversion. For the latter, a smooth initial model is often preferred over a rough one. We test the inversion approach on synthetic elastic data computed for a single-layer P-wave model and on field data, both with a small v_{s}/v_{p} ratio. We find that a single-layer model can recover either the shallow or deeper part of the profile but not both, when compared with the result of a multi-layer inversion that we use as a reference. An extension of our analytic model to two layers above a homogeneous half-space, each with a constant vertical gradient of the squared P-wave slowness and connected in a continuous manner, improves the fit of the picked dispersion curves. The resulting profile resembles a smooth approximation of the multi-layered one but contains, of course, less detail. As it turns out, our method does not degrade as gracefully as, for instance, diving-wave tomography, and we can only hope to fit a subset of the dispersion curves. Therefore, the applicability of the method is limited to cases where the v_{s}/v_{p} ratio is small and the profile is sufficiently simple. A further extension of the two-layer model to more layers, each with a constant depth gradient of the squared slowness, might improve the fit of the modal structure but at an increased cost.","P-wave velocity profile; Squared-slowness gradient; Surface-wave inversion","en","journal article","","","","","","","","2018-07-31","","","Applied Geophysics and Petrophysics","","",""
"uuid:e5221293-3aa4-4b5f-9ba7-00f6ae7e4251","http://resolver.tudelft.nl/uuid:e5221293-3aa4-4b5f-9ba7-00f6ae7e4251","Estimating a continuous p-wave velocity profile with constant squared-slowness gradient models from seismic field data","Ponomarenko, A.V.; Kashtan, B.M.; Troyan, V.N.; Mulder, W.A.","","2015","We inverted seismic field data for a continuous, laterally invariant P-wave velocity profile. Instead of the usual approach that involves horizontal layers with piecewise constant densities and velocities, we consider models of one or two layers with a constant gradient of the squared slowness above a homogeneous halfspace. With a single layer above a halfspace, there are three parameters. With two layers, there are five. We solve the inverse problem by a direct grid search over a wide range of parameters. The results were compared to that of a piecewise-constant multi-layer inversion result. In the single-layer case, either the shallow or the deeper part of the model would match the multi-layer case, depending on which modes of the surface waves were selected. With two layers, a considerably better agreement is obtained over a larger depth range. Our method is limited to cases with a small Vs/Vp-ratio but has only 5 parameters. It could be a useful alternative to piecewise-constant multi-layer inversion, in particular if continuous P-velocity profiles are sought. These are sometimes better suited as a starting model for full waveform inversion than models with many discontinuities.","","en","conference paper","EAGE","","","","","","","","Civil Engineering and Geosciences","Geoscience and Engineering","","","",""
"uuid:74b69402-43b9-4530-bfad-1332c26bee0f","http://resolver.tudelft.nl/uuid:74b69402-43b9-4530-bfad-1332c26bee0f","Free-surface Multiples and full-waveform inversion spectral resolution","Kazei, V.V.; Kashtan, B.M.; Troyan, V.N.; Mulder, W.A.","","2015","Low frequencies play a crucial role in the convergence of full-waveform inversion to the correct model in most of its current implementations. However, the lower the frequencies, the bigger are the amplitudes of the surface waves, causing the inversion to be driven by the latter. If they are not blanked out or removed, this may lead to convergence problems. To analyze this situation, we consider the simplest case where surface waves are present: an acoustic layer over a halfspace. We earlier analyzed the contributions of various wave types to the wavenumber spectrum of a velocity perturbation above a reflecting halfspace, without a free surface. Here, we extend this spectral sensitivity analysis to the case with a free surface, which generates multiples and ghosts. In this setting, the surface guided P-waves can be considered as a superposition of free-surface multiples. Our analysis shows that the conditioning of the linearized inverse problem, which is solved at each iteration of full-waveform inversion, becomes worse when multiples are taken into account. At the same time the inclusion of multiples increases the sensitivity to some low wavenumbers in the model spectrum, which should be beneficial for full-waveform inversion once a suitable preconditioner has been found.","","en","conference paper","EAGE","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:97b8d0b7-72c5-45bf-850f-0909393ce5e9","http://resolver.tudelft.nl/uuid:97b8d0b7-72c5-45bf-850f-0909393ce5e9","Coupling of Elastic Isotropic Medium Parameters in Iterative Linearized Inversion","Anikiev, D.V.; Kashtan, B.M.; Mulder, W.A.; Troyan, V.N.","","2014","An elastic isotropic medium is described with three parameters. In seismic migration the perturbation of one elastic parameter affects the images of all the three, which means that these parameters are coupled. For an effective quantitative reconstruction of the true elastic medium reflectivity one can apply an iterative linearized migration/inversion, where minimization of the misfit functional is done by the conjugate gradient method. The final result of the iterative approach can be obtained directly by Newton’s method, using the pseudo-inverse Hessian matrix. Calculation of this matrix for a realistic model is an extremely resource-intensive problem, but for a model of a scatterer in a homogeneous elastic medium it is quite feasible. This paper presents the numerical results of elastic linearized inversion for this simple model, calculated both with iterative approach and Newton’s method. Experiments show that in the both cases the elastic parameters have coupled weaker than in the case of migration. The iterative approach allows achieving acceptable quality of the inversion, but requires a large number of iterations. For faster convergence it is necessary to use the preconditioned conjugate gradient method. The optimal preconditioning will improve the convergence of the method as well as the quality of inversion.","linearized inversion; seismic migration; elastic isotropic medium; coupling of elastic parameters","en","journal article","","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:f869b091-a713-4c3c-84c7-4a7372db92c0","http://resolver.tudelft.nl/uuid:f869b091-a713-4c3c-84c7-4a7372db92c0","Estimation of the p-wave velocity profile of elastic real data based on surface wave inversion","Ponomarenko, A.V.; Kashtan, B.M.; Troyan, V.N.; Mulder, W.A.","","2013","Recently, we proposed an analytical approach to invert for a smoothly varying near-surface P-wave velocity profile that has a squared slowness linearly decreasing with depth. The exact solution for such a velocity profile in the acoustic approximation can be expressed in terms of Airy functions and leads to a dispersion equation. The method was successfully applied to synthetic elastic data with small Vs/Vp-ratio. Here, we apply the method to land data. The result agrees with that of multi-layered inversion, confirming its potential to provide an initial P-wave velocity model for acoustic full waveform inversion. Compared to multi-layered inversion, the method is simpler to use and produces a smooth model characterized by three parameters. In some cases, having a smooth rather than a blocky initial model for full waveform inversion is more appropriate.","","en","conference paper","EAGE","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:a33951a6-b588-470d-887e-e2924871c0e7","http://resolver.tudelft.nl/uuid:a33951a6-b588-470d-887e-e2924871c0e7","Surface wave inversion for a p-wave velocity profile: Estimation of the squared slowness gradient","Ponomarenko, A.V.; Kashtan, B.M.; Troyan, V.N.; Mulder, W.A.","","2013","Surface waves can be used to obtain a near-surface shear wave profile. The inverse problem is usually solved for the locally 1-D problem of a set of homogeneous horizontal elastic layers. The output is a set of shear velocity values for each layer in the profile. P-wave velocity profile can be estimated if higher modes and P-guided waves are used in the inversion scheme. Here, we use an exact acoustic solution to invert for the P-velocity profile in an elastic model with a decreasing constant vertical gradient of the squared P-wave slowness, bounded by a free surface on the top and a homogeneous halfspace at the bottom. The exact acoustic solution can be expressed in Airy functions and leads to a dispersion equation. We can invert several modes of the dispersion equation for the single gradient parameter of the squared P-wave slowness from elastic data. As a first test case, we invert for the P-wave velocity profile of synthetic 2-D isotropic elastic data with a small Vs/Vp-ratio, using the first two dispersive P-wave modes. The method does not require any picking and should be able to provide an initial model for full waveform inversion when applied to real data.","","en","conference paper","EAGE","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:e6c2357e-e8a5-41ca-9759-ca41638d438a","http://resolver.tudelft.nl/uuid:e6c2357e-e8a5-41ca-9759-ca41638d438a","Spectral sensitivity analysis of FWI in a constant-gradient background velocity model","Kazei, V.; Kashtan, B.M.; Troyan, V.N.; Mulder, W.A.","","2013","Full waveform inversion suffers from local minima, due to a lack of low frequencies in the data. A reflector below the zone of interest may help in recovering the long-wavelength components of a velocity perturbation, as demonstrated in a paper by Mora. Because smooth models are more popular as initial guesses for FWI, we consider the Born approximation for a perturbation in a reference model with a constant velocity gradient. Analytic expressions are found that describe the spatial wavenumber spectrum of the recorded seismic signal as a function of the spatial spectrum of the inhomogeneity. We study this spectrum in more detail in terms of sensitivities. Since the velocity model is inhomogeneous near the perturbation, we need to specify its depth. We compare these sensitivities and find that low frequencies are extremely useful for the first stages of inversion – a well-known fact. However, also the high-frequency data contain some information about the low spatial wavenumbers in the perturbation, which offers opportunities for inversion in the absence of low frequencies in the data. We observe that the longer wavenumbers are better resolved in the deeper parts of the model if large enough offsets are available.","","en","conference paper","EAGE","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:f945fd0c-180b-4b4a-995c-7b84d37c1bc2","http://resolver.tudelft.nl/uuid:f945fd0c-180b-4b4a-995c-7b84d37c1bc2","FWI sensitivity analysis in the presence of free-surface multiples","Kazei, V.V.; Kashtan, B.M.; Troyan, V.N.; Mulder, W.A.","","2013","It is generally believed that full waveform inversion needs very low frequencies in the data to avoid convergence to a local minimum, which would lead to an incorrect velocity model. Often, low frequencies are not present in the data. Then, a kinematically accurate initial model is required that contains the long-wavelength structures that cannot be reconstructed from the data. Mora (1989), however, showed that reflector below the target area may help in recovering some of the longwavelength information. Here, we extend his analysis to the case with a free surface, generating multiples. At first sight, multiples should improve the sensitivity and resolution. Indeed, our study shows that sensitivity becomes higher in some parts of the model spectrum when multiples are included. If we consider thin layers and only first-order multiples, then the sensitivity is improved. However, when the layer is thick enough to allow for normal modes of higher orders, areas of high sensitivity appear and this leads to a poorer conditioning of the inverse problem in the wavenumber domain.","","en","conference paper","","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:0f0b7743-c611-4754-ac69-505651b052d8","http://resolver.tudelft.nl/uuid:0f0b7743-c611-4754-ac69-505651b052d8","On the Contribution of Head Waves to Full Waveform Inversion","Kazei, V.V.; Ponomarenko, A.V.; Troyan, V.N.; Kashtan, B.M.; Mulder, W.A.","","2012","Full waveform inversion suffers from local minima, due to a lack of low frequencies in the data. A reflector below the zone of interest may, however, help in recovering the long-wavelength components of a velocity perturbation, as demonstrated in a paper by Mora. With the Born approximation for the perturbation in a reference model consisting of two homogeneous isotropic acoustic halfspaces, analytic expressions can be found that describe the spatial spectrum of the recorded seismic signal as a function of the spatial spectrum of the inhomogeneity. We study this spectrum in more detail by separately considering direct, reflected and head waves. Taking the reflection coefficient of the deeper reflector into account, we obtain sensitivity estimates for each of these types of waves. Although the head waves have a relatively small contribution to the reconstruction of the velocity perturbation, compared to the other waves, they contain reliable long-wavelength information that can be beneficial for full waveform inversion.","","en","conference paper","","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""