Oh, Juwon; Alkhalifah, Tariq Ali(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09-06)[Conference Paper]
For the purpose of extracting higher resolution information from a 3D field data set, we apply a 3D elastic orthorhombic (ORT) anisotropic full waveform inversion (FWI) to hopefully better represent the physics of the Earth. We utilize what we consider as the optimal parameterization for surface acquired seismic data over a potentially orthorhombic media. This parameterization admits the possibility of incorporating a hierarchical implementation moving from higher anisotropy symmetry to lower ones. From the analysis of the radiation pattern of this new parameterization, we focus the inversion of the 3D data on the parameters that may have imprint on the data with minimal tradeoff, and as a result we invert for the horizontal P-wave velocity model, an ε1 model, its orthorhombic deviation, and the shear wave velocity. The inverted higher resolution models provide reasonable insights of the medium.
Dutta, Gaurav; Schuster, Gerard T.(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
A wave-equation gradient optimization method is presented that inverts for the subsurface Q distribution by minimizing a skeletonized misfit function ε. Here, ε is the sum of the squared differences between the observed and the predicted peak/centroid frequency shifts of the early-arrivals. The gradient is computed by migrating the observed traces weighted by the frequency-shift residuals. The background Q model is perturbed until the predicted and the observed traces have the same peak frequencies or the same centroid frequencies. Numerical tests show that an improved accuracy of the inverted Q model by wave-equation Q tomography (WQ) leads to a noticeable improvement in the migration image quality.
Feng, Zongcai; Schuster, Gerard T.(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
Elastic least-squares reverse time migration (LSRTM) is used to invert synthetic particle-velocity data and crosswell pressure field data. The migration images consist of both the P- and Svelocity perturbation images. Numerical tests on synthetic and field data illustrate the advantages of elastic LSRTM over elastic reverse time migration (RTM). In addition, elastic LSRTM images are better focused and have better reflector continuity than do the acoustic LSRTM images.
Hanafy, Sherif; Schuster, Gerard T.(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
We present parsimonious refraction interferometry where a densely populated refraction data set can be obtained from just two shot gathers. The assumptions are that the first arrivals are comprised of head waves and direct waves, and a pair of reciprocal shot gathers is recorded over the line of interest. The refraction traveltimes from these reciprocal shot gathers can be picked and decomposed into O(N2) refraction traveltimes generated by N virtual sources, where N is the number of geophones in the 2D survey. This enormous increase in the number of virtual traveltime picks and associated rays, compared to the 2N traveltimes from the two reciprocal shot gathers, allows for increased model resolution and better condition numbers in the normal equations. Also, a reciprocal survey is far less time consuming than a standard refraction survey with a dense distribution of sources.
Liu, Zhaolun; AlTheyab, Abdullah; Hanafy, Sherif M.; Schuster, Gerard T.(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
We demonstrate that near-surface heterogeneities can be imaged by natural migration of backscattered surface waves in common shot gathers. No velocity model is required because the data are migrated onto surface points with the virtual Green's functions computed from the shot gathers. Migrating shot gathers recorded by 2D and 3D land surveys validates the effectiveness of detecting nearsurface heterogeneities by natural migration. The implication is that more accurate hazard maps can be created by migrating surface waves in land surveys.
Li, Vladimir; Wang, Hui; Tsvankin, Ilya; Diaz, Esteban; Alkhalifah, Tariq Ali(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
Wavefield tomography can handle complex subsurface geology better than ray-based techniques and, ultimately, provide a higher resolution. Here, we implement forward and adjoint wavefield extrapolation for VTI (transversely isotropic with a vertical symmetry axis) media using a pseudospectral operator that employes a separable approximation of the P-wave dispersion relation. This operator is employed to derive the gradients of the differential semblance optimization (DSO) and modified stack-power objective functions. We also obtain the gradient expressions for the data-domain objective function, which can incorporate borehole information necessary for stable VTI velocity analysis. These gradients are compared to the ones obtained with a space-time finite-difference (FD) scheme for a system of coupled wave equations. Whereas the kernels computed with the two wave-equation operators are similar, the pseudospectral method is not hampered by the imprint of the shear-wave artifact. Numerical examples also show that the modified stack-power objective function produces cleaner gradients than the more conventional DSO operator.
Alkhalifah, Tariq Ali; Wu, Zedong(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
Addressing anisotropy in full wavenumber inversion (FWI) is crucial to obtaining credible models, and it is extremely challenging considering the multi parameter nature of the inversion. A successful FWI in anisotropic media takes into account the sensitivity of the data (or the wave) to the long and short wavelength components of the anisotropic parameters. Considering the low sensitivity of FWI to the anellipticity parameter ? when parametrizing the acoustic transversely isotropic model with the horizontal velocity, η and ε, we develop a combined FWI and reflection waveform inversion (RWI) to invert for the anisotropic parameters that influence surface seismic data. This practical waveform inversion (PWI) separates the parameters to their resolvable scales, with information accessed from the data fitting (FWI) and the image focusing (RWI) objectives. With this parametrization, the RWI role is to obtain a smooth ηmodel, as well as velocity, while FWI focusses on the scattering potential of the horizontal velocity. The parameter η is used to produce the Born scattered wavefield for the RWI part and eventually fit the amplitude for the imperfect physics in the FWI part.
Wang, Hanchen(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
Using full waveform inversion (FWI) to locate microseismic and image microseismic events allows for an automatic process (free of picking) that utilizes the full wavefield. However, waveform inversion of microseismic events faces incredible nonlinearity due to the unknown source location (space) and function (time). We develop a source independent FWI of microseismic events to invert for the source image, source function and the velocity model. It is based on convolving reference traces with the observed and modeled data to mitigate the effect of an unknown source ignition time. The adjoint-state method is used to derive the gradient for the source image, source function and velocity updates. The extended image for source wavelet in z axis is extracted to check the accuracy of the inverted source image and velocity model. Also the angle gather is calculated to see if the velocity model is correct. By inverting for all the source image, source wavelet and the velocity model, the proposed method produces good estimates of the source location, ignition time and the background velocity for part of the SEG overthrust model.
Feng, Shihang; Schuster, Gerard T.(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
The wave-equation traveltime and waveform inversion (WTW) methodology is developed to invert for anisotropic parameters in a vertical transverse isotropic (VTI) meidum. The simultaneous inversion of anisotropic parameters v0, ε and δ is initially performed using the wave-equation traveltime inversion (WT) method. The WT tomograms are then used as starting background models for VTI full waveform inversion. Preliminary numerical tests on synthetic data demonstrate the feasibility of this method for multi-parameter inversion.
Alkhalifah, Tariq Ali; Masmoudi, Nabil; Oh, Juwon(SEG Technical Program Expanded Abstracts 2016, Society of Exploration Geophysicists, 2016-09)[Conference Paper]
Multi parameter full waveform inversion (FWI) usually suffers from the inherent tradeoffin the multi parameter nature of the model space. In orthorhombic anisotropy, such tradeoffis magnified by the large number of parameters involved in representing the elastic or even the acoustic approximation of such a medium. However, using a new parameterization with distinctive scattering features, we can condition FWI to invert for the parameters the data are sensitive to at different stages, scales, and locations in the model. Specifically, with a combination made up of a velocity and particular dimensionless ratios of the elastic coefficients, the scattering potential of the anisotropic parameters have stationary scattering radiation patterns as a function of the type of anisotropy. With our new parametrization, the data is mainly sensitive to the scattering potential of 4 parameters: the horizontal velocity in the x direction, x, which provides scattering mainly near zero offset in the x vertical plane, εd, which is the ratio of the horizontal velocity squared in the x and x direction, and δ3 describing the anellipticity in the horizontal plane. Since, with this parametrization, the radiation pattern for the horizontal velocity and ε is azimuth independent, we can perform an initial VTI inversion for these two parameters, and then use the other two parameters to fit the azimuth variation in the data. This can be done at the reservoir level or any region of the model. Including the transmission from reflections, the migration velocity analysis (MVA) component, into the picture, the multi azimuth surface seismic data are mainly sensitive to the long wavelength components of uh, δ3, and εd through the diving waves, and η1, ηd, and δ3, in the transmission to or from reflectors (especially, in the presence of large offsets). They are also sensitive to the short wavelength component of uh and ε.
The export option will allow you to export the current search results of the entered query to a file. Different
formats are available for download. To export the items, click on the button corresponding with the preferred download format.
By default, clicking on the export buttons will result in a download of the allowed maximum amount of items.
For anonymous users the allowed maximum amount is 50 search results.
To select a subset of the search results, click "Selective Export" button and make a selection of the items you want to export.
The amount of items that can be exported at once is similarly restricted as the full export.
After making a selection, click one of the export format buttons. The amount of items that will be exported is indicated in the bubble next to export format.