Discretization method and system of angle domain generalized Radon transformation, terminal and medium

Abstract

The invention relates to a discretization method and system of angle domain generalized Radon transformation, a terminal and a storage medium. The method comprises the following steps: reading seismic trace data, setting imaging point coordinates of the seismic trace data, and reading travel time, amplitude and direction table of each imaging point coordinate; calculating a scattering angle and a discrete interval; calculating a sampling point theta0 closest to the scattering angle and a partition boundary of the sampling point theta0, calculating a distance l from the scattering angle to the partition boundary thetabound, judging whether the partition boundary thetabound splits a scattering angle discrete surface or not according to the distance l, if so, splitting a discrete unit near the partition boundary thetabound into two units by using a discrete unit splitting algorithm, and if not, splitting the discrete unit near the partition boundary thetabound into two units; calculating the areas of the two split units according to the discrete interval; and respectively using the two split units to calculate imaging accumulation at the sampling point theta0 and at the position of theta0 + / -delta theta. According to the invention, a continuous and smooth angle domain inversion result can be obtained.

Description

technical field [0001] The present application belongs to the technical field of information extraction in angle domain, and in particular relates to a discrete method, system, terminal and storage medium of generalized Radon transform in angle domain. Background technique [0002] The inventor proposed and established the angle-domain generalized Radon transform (AD-GRT) inversion method (Li Xuelei et al., 2020), which can realize multi-parameter inversion of acoustic waves and elastic media without iteration. Among them, the theoretical framework of the generalized Radon transform in the angle domain is: [0003] Based on the Born approximation, the wave field of the acoustic wave medium can be expressed as a generalized Radon transform (GRT) form (frequency domain): [0004] p'(r,s,ω)=ω 2 ∫dx 1 dx 3 A s A r exp[iω(φ s +φ r )]κ 0 F(ω)f(x,θ) (1) [0005] Among them, p′(r,s,ω) is the acoustic pressure scattering wave field, A s and φ s represent the amplitude and...

AI Technical Summary

Problems solved by technology

However, due to the discrete integer features, this conventional algorithm leads to different θ 0 f(y,θ 0 ) calculation is very unbalanced, and it is easy to appear with θ 0 jagged oscillating ups and downs

In particular, the distribution of source and receiver points is relatively sparse, or θ 0 When the value sampling is relatively dense, this kind of shock fluctuation is very serious, making it difficult to extract accurate f(y,θ 0 ) distribution relationship, which makes the angle domain inversion results invalid

Even if some smoothing factors can slow down the shock, there are still obvious fluctuations, and the smoothing factor will bring additional information loss

Images