Seismic imaging system for acoustic-elastic coupled media using accumulated Laplace gradient direction

Patent metadata

Field	Value
Publication number	US-9229121-B2
Application number	US-201313795072-A
Country	US
Kind code	B2
Filing date	Mar 12, 2013
Priority date	Mar 13, 2012
Publication date	Jan 5, 2016
Grant date	Jan 5, 2016

Abstract

The description relates to a seismic imaging technology technique for modeling a subsurface structure through waveform inversion in the Laplace domain. The seismic imaging system comprises a scaled gradient calculating unit calculating a scaled gradient, a modeling parameter updating unit updating the model parameters using the scaled gradient direction, and an iteration control unit controlling the scaled gradient calculating unit and the modeling parameter updating unit to repeat processing iteratively until a stopping criteria is met.

First claim

What is claimed is:

1. A seismic imaging system using Laplace-domain wave equation in acoustic-elastic coupled media, comprising:
a scaled gradient calculating hardware circuit calculating a scaled gradient direction of an objective function by multiplying a gradient of the objective function and an accumulated gradient, wherein the accumulated gradient is calculated using a vertically accumulated square sum of original gradient values by summing an accumulated gradient of a previous vertical layer and a square of the gradient of the objective function, and the objective function is defined by a residual of logarithmic wavefield in Laplace domain;
a modeling parameter updating hardware circuit updating model parameters for a discretized finite element equation for the acoustic-elastic coupled media from a Laplace-transformed acoustic wave equation, a Laplace-transformed elastic wave equation and continuity conditions at an interface between an acoustic media and an elastic media of the acoustic-elastic coupled media;
an iteration control hardware circuit controlling the scaled gradient calculating circuit and the modeling parameter updating hardware circuit to repeat processing iteratively until a stopping criteria is met; and
an imaging hardware circuit converting Lamè constants into P-wave velocity, S-wave velocity, and density and generating a submarine subsurface image showing the distribution of P-wave and S-wave velocities.

2. The seismic imaging system in claim 1 , wherein the system further comprises a source wavelet estimation hardware circuit estimating an unknown source wavelet from an initial source using a full Newton method and outputting the result to the modeling parameter updating hardware circuit.

3. The seismic imaging system of claim 1 , wherein the accumulated gradient is always positive.

4. A seismic imaging method using Laplace-domain wave equation in acoustic-elastic coupled media, comprising:
calculating, using a scaled gradient calculation hardware circuit, a scaled gradient direction of an objective function by multiplying a gradient of the objective function and an accumulated gradient, the accumulated gradient being calculated using a vertically accumulated square sum of original gradient values by summing an accumulated gradient of a previous vertical layer and a square of the gradient of the objective function, wherein the objective function is defined by a residual of logarithmic wavefield in Laplace domain;
updating model parameters for a discretized finite element equation for the acoustic-elastic coupled media from a Laplace-transformed acoustic wave equation, a Laplace-transformed elastic wave equation and continuity conditions at an interface between an acoustic media and an elastic media of the acoustic-elastic coupled media;
repeating calculating the scaled gradient direction and updating the model parameters iteratively until a stopping criteria is met; and
converting, using an imaging hardware circuit, Lamè constants into P-wave velocity, S-wave velocity, and density and generating a submarine subsurface image showing the distribution of P-wave and S-wave velocities.

5. The seismic imaging method in claim 4 , wherein updating the model parameters includes estimating an unknown source wavelet from an initial source using a full Newton method.

Assignees

• Shin Changsoo
• Univ Seoul Nat R & Db Found

Inventors

• Shin Changsoo

Classifications

• G01V1/282 Primary

  Application of seismic models, synthetic seismograms · CPC title