Angle domain generalized Radon transform for elastic multi-parameter inverse scattering inversion

Geophysics ◽  
2021 ◽  
pp. 1-71
Author(s):  
Xuelei Li ◽  
Yanjie Wei ◽  
Wei Ouyang

Linearized algorithms based on the Born approximation are well-known and popular techniques for quantitative seismic imaging and inversion. However, linearization methods usually suffer from some significant problems, such as computational cost for the required number of iterations, requirement for background models, and uncertain and unstable multi-parameter extraction, which make the methods difficult to implement in practical applications. To avoid these problems, we propose an angle-domain generalized Radon transform (AD-GRT) inversion in 2D elastic isotropic media. This AD-GRT is an approximate transform between the seismic data and an angle-domain model, which acts as a scattering function, and the seismic data can be reconstructed accurately, even when the background models are incorrect. The density and Lam短oduli perturbation parameters can be extracted stably from the inverted angle-domain scattering function. Deconvolution of the source wavelet is taken into account to remove the effect of the wavelet and improve the resolution and accuracy of the inversion results. The derived AD-GRT inversion is non-iterative and is as efficient as the traditional elastic GRT method. The additional dimension of the angle domain has little effect on the computational cost of the AD-GRT, as opposed to other extended-domain inversion/migration methods. Our method also can be used to solve non-linear Born inversion problems using iteration, which can significantly improve their convergence rate. Three numerical examples illustrate that the angle-domain scattering function inversion, data reconstruction, and multi-parameter extraction using the presented AD-GRT inversion are effective.

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. V317-V328
Author(s):  
Jitao Ma ◽  
Guoyang Xu ◽  
Xiaohong Chen ◽  
Xiaoliu Wang ◽  
Zhenjiang Hao

The parabolic Radon transform is one of the most commonly used multiple attenuation methods in seismic data processing. The 2D Radon transform cannot consider the azimuth effect on seismic data when processing 3D common-depth point gathers; hence, the result of applying this transform is unreliable. Therefore, the 3D Radon transform should be applied. The theory of the 3D Radon transform is first introduced. To address sparse sampling in the crossline direction, a lower frequency constraint is introduced to reduce spatial aliasing and improve the resolution of the Radon transform. An orthogonal polynomial transform, which can fit the amplitude variations in different parabolic directions, is combined with the dealiased 3D high-resolution Radon transform to account for the amplitude variations with offset of seismic data. A multiple model can be estimated with superior accuracy, and improved results can be achieved. Synthetic and real data examples indicate that even though our method comes at a higher computational cost than existing techniques, the developed approach provides better attenuation of multiples for 3D seismic data with amplitude variations.


Author(s):  
Jimmy Ming-Tai Wu ◽  
Qian Teng ◽  
Shahab Tayeb ◽  
Jerry Chun-Wei Lin

AbstractThe high average-utility itemset mining (HAUIM) was established to provide a fair measure instead of genetic high-utility itemset mining (HUIM) for revealing the satisfied and interesting patterns. In practical applications, the database is dynamically changed when insertion/deletion operations are performed on databases. Several works were designed to handle the insertion process but fewer studies focused on processing the deletion process for knowledge maintenance. In this paper, we then develop a PRE-HAUI-DEL algorithm that utilizes the pre-large concept on HAUIM for handling transaction deletion in the dynamic databases. The pre-large concept is served as the buffer on HAUIM that reduces the number of database scans while the database is updated particularly in transaction deletion. Two upper-bound values are also established here to reduce the unpromising candidates early which can speed up the computational cost. From the experimental results, the designed PRE-HAUI-DEL algorithm is well performed compared to the Apriori-like model in terms of runtime, memory, and scalability in dynamic databases.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. A17-A21 ◽  
Author(s):  
Juan I. Sabbione ◽  
Mauricio D. Sacchi

The coefficients that synthesize seismic data via the hyperbolic Radon transform (HRT) are estimated by solving a linear-inverse problem. In the classical HRT, the computational cost of the inverse problem is proportional to the size of the data and the number of Radon coefficients. We have developed a strategy that significantly speeds up the implementation of time-domain HRTs. For this purpose, we have defined a restricted model space of coefficients applying hard thresholding to an initial low-resolution Radon gather. Then, an iterative solver that operated on the restricted model space was used to estimate the group of coefficients that synthesized the data. The method is illustrated with synthetic data and tested with a marine data example.


Author(s):  
Francisco González ◽  
Pierangelo Masarati ◽  
Javier Cuadrado ◽  
Miguel A. Naya

Formulating the dynamics equations of a mechanical system following a multibody dynamics approach often leads to a set of highly nonlinear differential-algebraic equations (DAEs). While this form of the equations of motion is suitable for a wide range of practical applications, in some cases it is necessary to have access to the linearized system dynamics. This is the case when stability and modal analyses are to be carried out; the definition of plant and system models for certain control algorithms and state estimators also requires a linear expression of the dynamics. A number of methods for the linearization of multibody dynamics can be found in the literature. They differ in both the approach that they follow to handle the equations of motion and the way in which they deliver their results, which in turn are determined by the selection of the generalized coordinates used to describe the mechanical system. This selection is closely related to the way in which the kinematic constraints of the system are treated. Three major approaches can be distinguished and used to categorize most of the linearization methods published so far. In this work, we demonstrate the properties of each approach in the linearization of systems in static equilibrium, illustrating them with the study of two representative examples.


2014 ◽  
Vol 672-674 ◽  
pp. 1964-1967
Author(s):  
Jun Qiu Wang ◽  
Jun Lin ◽  
Xiang Bo Gong

Vibroseis obtained the seismic record by cross-correlation detection calculation. compared with dynamite source, cross-correlation detection can suppress random noise, but produce more correlation noise. This paper studies Radon transform to remove correlation noise produced by electromagnetic drive vibroseis and impact rammer. From the results of processing field seismic records, we can see that Radon transform can remove correlation noise by vibroseis, the SNR of vibroseis seismic data is effectively improved.


2016 ◽  
Vol 12 (3) ◽  
pp. 145
Author(s):  
Subarsyah Subarsyah ◽  
Tumpal Benhard Nainggolan

Interferensi water-bottom multipel terhadap reflektor primer menimbulkan efek bersifat destruktif yang menyebabkan penampang seismik menjadi tidak tepat akibat kehadiran reflektor semu. Teknik demultiple perlu diaplikasikan untuk mengatenuasi multipel. Transformasi parabolic radon merupakan teknik atenuasi multipel dengan metode pemisahan dalam domain radon. Multipel sering teridentifikasi pada penampang seismik. Untuk memperbaiki penampang seismik akan dilakukan dengan metode transformasi parabolic radon. Penerapan metode ini mengakibatkan reflektor multipel melemah dan tereduksi setelah dilakukan muting dalam domain radon terhadap zona multipel. Beberapa reflektor primer juga ikut melemah akibat pemisahan dalam domain radon yang kurang optimal, pemisahan akan optimal membutuhkan distribusi offset yang lebar. Kata kunci: Parabolic radon, multipel, atenuasi Water-bottom mutiple interference often destructively interfere with primary reflection that led to incorrect seismic section due to presence apparent reflector. Demultiple techniques need to be applied to attenuate the multiple. Parabolic Radon transform is demultiple attenuation technique that separate multiple and primary in radon domain. Water-bottom mutiple ussualy appear and easly identified on seismic data, parabolic radon transform applied to improve the seismic section. Application of this method to data showing multiple reflectors weakened and reduced after muting multiple zones in the radon domain. Some of the primary reflector also weakened due to bad separation in radon domain, optimal separation will require a wide distribution of offsets. Keywords: Parabolic radon, multiple, attenuation


Recent applications of conventional iterative coordinate descent (ICD) algorithms to multislice helical CT reconstructions have shown that conventional ICD can greatly improve image quality by increasing resolution as well as reducing noise and some artifacts. However, high computational cost and long reconstruction times remain as a barrier to the use of conventional algorithm in the practical applications. Among the various iterative methods that have been studied for conventional, ICD has been found to have relatively low overall computational requirements due to its fast convergence. This paper presents a fast model-based iterative reconstruction algorithm using spatially nonhomogeneous ICD (NH-ICD) optimization. The NH-ICD algorithm speeds up convergence by focusing computation where it is most needed. The NH-ICD algorithm has a mechanism that adaptively selects voxels for update. First, a voxel selection criterion VSC determines the voxels in greatest need of update. Then a voxel selection algorithm VSA selects the order of successive voxel updates based upon the need for repeated updates of some locations, while retaining characteristics for global convergence. In order to speed up each voxel update, we also propose a fast 3-D optimization algorithm that uses a quadratic substitute function to upper bound the local 3-D objective function, so that a closed form solution can be obtained rather than using a computationally expensive line search algorithm. The experimental results show that the proposed method accelerates the reconstructions by roughly a factor of three on average for typical 3-D multislice geometries.


2019 ◽  
Vol 34 (1) ◽  
Author(s):  
Tumpal Bernhard Nainggolan ◽  
Said Muhammad Rasidin ◽  
Imam Setiadi

Multiple often and always appear in marine seismic data due to very high acoustic impedance contrasts. These events have undergone more than one reflection. This causes the signal to arrive back at the receiver at an erroneous time, which, in turn, causes false results and can result in data misinterpretation. Several types of multiple suppression have been studied in literature. Methods that attenuate multiples can be classified into three broad categories: deconvolution methods; filtering methods and wavefield prediction subtraction methods. The study area is situated on Seram Sea in between 131°15’E – 132°45’E and 3°0’S – 4°0’S, Seram Trough which is located beneath Seram Sea at northern part of the Banda-Arc – Australian collision zone and currently the site of contraction between Bird’s Head and Seram. This research uses predictive deconvolution and FK-filter to attenuate short period multiple from their move out, then continued by SRME method to predict multiple that cannot be attenuated from previous method, then followed by Radon transform to attenuate multiple that still left and cannot be attenuated by SRME method. The result of each method then compared to each other to see how well multiple attenuated. Predictive deconvolution and F-K filter could not give satisfactory result especially complex area where multiple in dipping event is not periodic, SRME method successfully attenuate multiple especially in near offset multiple without need subsurface information, while SRME method fails to attenuate long offset multiple, combination of SRME method and Radon transform can give satisfactory result with careful selection of the Radon transform parameters because it can obscure some primary reflectors. Based on geological interpretation, Seram Trough is built by dominant structural style of deposited fold and thrust belt. The deposited fold and thrust belt has a complexly fault geometry from western zone until eastern of seismic line.


Sign in / Sign up

Export Citation Format

Share Document