Linearized inversion of 3-D multi-offset data: background reconstruction and AVO inversion

Abstract

We present a linearized inversion algorithm for 3-D multi-offset seismic reflection data. The fundamental assumptions are that the background medium is only depth-dependent and that the 3-D multi-offset seismic data are regularly sampled as a series of 2-D multi-offset profiles. Our algorithm can be performed in the following three steps.(1) The pressure field P(x, y, h, t ) is 4-D Fourier transformed to P(k,, k,, k,,, o) with respect to the midpoint coordinates (x, y), the half-offset coordinates (h, 0), and the time t.(2) A 2-D inversion operator is applied to P(k,, k,, k,, o) for each given midpoint wavenumber k,. This computation can be carried out in parallel for all the k, wavenumbers so that this step becomes numerically equivalent to a 2-D multi-offset inversion in the Fourier domain.(3) A deconvolution by the source and for spatial effects is applied to correct for geometrical spreading and to separate the elastic parameters.About 98 per cent of the computation time of this algorithm takes place in the second step. Therefore it is worthwhile to note that this step is suitable for a data-parallel computing model.Most of the central memory and/or storage device required in this algorithm is used to store the 4-D FFT of the pressure field. By taking advantage of the limited frequency content of the seismic data, we can reduce its storage to about 4 Gbyte for a section of 256 lines, each containing up to 256 CMP.The input for this algorithm is a smooth, depth-dependent, background velocity model (the low spatial frequencies), and the results are the contrasts of P-wave impedance and acoustic velocity (the high spatial frequencies). The inversion for high spatial frequencies of P-wave impedance is essentially based on small-offset data, while that of acoustic velocity is essentially based on large-offset data. By using the equivalence between acoustic scattering and P-P scattering, the algorithm can also output the classical AVO parameters: P-wave intercept and gradient coefficients.This algorithm can alternatively be used to estimate the background velocity as follows. First, the pressure field is 4-D Fourier transformed, and inversions with various background velocity models are performed. For isotropic media, the good models for the background are then extracted by analysing outputs of P-wave impedance or acoustic velocity, as both yield the same background velocity. The anisotropic behaviour of the background model is depicted by differences between the background model obtained from the analysis of P-wave impedance outputs and that obtained from the analysis of acoustic velocity outputs. The computation time, in this process, can be substantially reduced either by wavenumber filtering because we are seeking a coarse grid model or by using analytic background velocity models to simplify the inversion solution further. 507 508 L. T. Ikelle 1 0 15 Travel times & Geom. Spread.Figure 1.In the linearized inverse problem theory, the model is decomposed into two parts: the low-spatial-frequency component and the...

Prestack amplitude versus angle inversion for Young's modulus and Poisson's ratio based on the exact Zoeppritz equations

“…Stewart, Zhang and Guthoff () showed that the Lamé parameters (λ and μ) of a medium better differentiate rock properties. The Zoeppritz equations are the basis of the prestack AVO/AVA inversion, but scholars have used numerous approximate equations for inversion because the Zoeppritz equations are strongly nonlinear (Ikelle ; Buland and Omre ; Li et al . ; Yin, Yang and Zhang ).…”

Formulation of the linearized forward problems for multicomponent OBS data in a water/solid configuration using the reciprocity theorem

“…Furthermore, its closed form makes Born approximation very attractive for constructing inversion algorithms (e.g. Snieder 1990; Ikelle 1995de). However, the adaptation of Born‐based solutions to the forward problem of multicomponent OBS data requires (1) the explicit introduction of the boundary conditions at a water/solid interface, (2) the use of the reciprocity theorem for a fluid/solid configuration and (3) the use of P ‐ and S ‐wave potentials.…”

Linearized inversion of data of multi-offset data for vertically inhomogeneous background

“…In this paper, the system of integral equations for small perturbations of elastic parameters was deduced and the uniqueness of the inverse problems was proved. Further, the similar approach was applied in [15] for the linearized inverse problem for elastic wave equation in the frequency domain. In [16] it was generalized for anisotropic elastic medium.…”

“…It should be noted that in these papers, maybe except [15,16], the main attention was paid to formal derivation of some mathematical relations which later were used in order to deduce some numerical algorithm of linearized inversion. Unfortunately, influence of such input parameters of the method as size of apertures, noise level, proximity of a priori known background to the exact structure of the medium under consideration and so on onto stability and quality of the gained results happened to be here out of interests.…”

“…In particular, papers [14][15][16] develop techniques when the initial multidimensional inverse problem is reduced to the set of one-dimensional integral equations of the first kind. As the correspondent integral operators have continuous kernels, they are compact operators; therefore, they have no bounded inverses.…”