Method for velocity analysis using waveform inversion in Laplace domain for geophysical imaging
원문보기
IPC분류정보
국가/구분
United States(US) Patent
등록
국제특허분류(IPC7판)
G01V-001/30
G01V-001/28
G01V-001/34
G06F-019/00
출원번호
UP-0942352
(2007-11-19)
등록번호
US-7797110
(2010-10-04)
발명자
/ 주소
Shin, Changsoo
출원인 / 주소
Shin's Geophysics
대리인 / 주소
Conley Rose, P.C.
인용정보
피인용 횟수 :
9인용 특허 :
7
초록▼
A Laplace transform system comprising a processor, a measured time domain wavefield, a velocity model, and Laplace damping constants, wherein the processor is programmed to calculate a substantially about zero frequency component of a Fourier transform of a time domain damped wavefield, wherein the
A Laplace transform system comprising a processor, a measured time domain wavefield, a velocity model, and Laplace damping constants, wherein the processor is programmed to calculate a substantially about zero frequency component of a Fourier transform of a time domain damped wavefield, wherein the time domain damped wavefield is damped by the Laplace damping constants to obtain long wavelength velocity information for deeper subsurface regions.
대표청구항▼
What is claimed is: 1. A processor configured to execute a Laplace transform system comprising: a measured time domain wavefield; a velocity model; and Laplace damping constants, wherein the processor is programmed to calculate a substantially about zero frequency component of a Fourier transform o
What is claimed is: 1. A processor configured to execute a Laplace transform system comprising: a measured time domain wavefield; a velocity model; and Laplace damping constants, wherein the processor is programmed to calculate a substantially about zero frequency component of a Fourier transform of a time domain damped wavefield, wherein the time domain damped wavefield is damped by the Laplace damping constants to obtain long wavelength velocity information for deeper subsurface regions. 2. The processor of claim 1 wherein a Laplace transform of the measured time domain wavefield is equivalent to a Poisson's equation when the Laplace damping constant is equal to zero. 3. The processor of claim 1 wherein the measured time domain wavefield is measured using seismic sources and receivers on a surface of a geological site of interest. 4. The processor of claim 1 wherein the measured time domain wavefield is prestack reflection data from subsurface structure. 5. The processor of claim 1 wherein a Laplace transform of the measured time domain wavefield is numerically calculated by integrating the damped time domain wavefield in a limit of a maximum recording time. 6. The processor of claim 1 wherein optimum values from 0 to 100 are chosen for the Laplace damping constant of the measured time domain wavefield to minimize integration error. 7. The processor of claim 1 wherein a small Laplace damping constant is used in a Laplace transform of the measured time domain wavefield to obtain information of deep parts of the waveform velocity model. 8. The processor of claim 1 wherein a large Laplace damping constant is used in a Laplace transform of the measured time domain wavefield to resolve information of shallow parts of the waveform velocity model. 9. A processor configured to execute a Laplace domain inversion system comprising: a damped recorded wavefield in a Laplace domain; a wave equation; an objective function; and a velocity model, wherein the processor is programmed to solve the wave equation in the Laplace domain, minimize the objective function, and calculate the velocity model that corresponds to the damped recorded wavefield to analyze subsurface regions. 10. The processor of claim 9 further comprising a source wavelet, wherein the processor is programmed to calculate the source wavelet and the velocity model. 11. The processor of claim 9 wherein the objective function is one of a logarithmic, an integral, and a power function. 12. The processor of claim 9 wherein the damped recorded wavefield in a Laplace domain is converted from data in a time domain. 13. The processor of claim 9 wherein the wave equation in a time domain comprises a wavefield component, a second-order time derivative of the wavefield component, a mass matrix component, and a stiffness matrix component. 14. The processor of claim 9 wherein the wave equation with proper boundary conditions, and using one of a real Laplace damping constant, and a complex Laplace damping constant, is solved in the Laplace domain by a factorization of some components of the wave equation and using forward or backward mathematical substitutions to obtain a damped calculated wavefield. 15. The processor of claim 9 wherein the objective function is equal to a summed squares of all differences between one of logarithmic, integral or power values of the damped recorded wavefield and a calculated damped wavefield. 16. The processor of claim 9 wherein the objective function is smooth and comprises no local minima or a small number of local minima. 17. The processor of claim 9 wherein one of a steepest-descent method, and a Gauss-Newton method is used to minimize the objective function by solving for a velocity model value that reduces a derivative value of the logarithmic objective to zero. 18. The processor of claim 9 wherein the processor is programmed to iterate starting with an initially calculated damped wavefield. 19. The processor of claim 9 wherein the processor is programmed to converge to the velocity model that corresponds to the damped recorded wavefield even if a different damped wavefield is initially calculated. 20. A method, implemented by one or more processors, for analyzing subsurface regions, comprising: receiving, at the one or more processors, collected prestack reflection data in a time domain; transforming, by the one or more processors, the time domain prestack reflection data to Laplace domain reflection data; initializing, by the one or more processors, a preliminary velocity model that represents a subsurface structure; calculating, by the one or more processors, one of a logarithmic, an integral, and a power objective function; verifying, by the one or more processors, whether the objective function satisfies a convergence criterion; updating, by the one or more processors, the velocity model if the convergence criterion is not met; and generating, by the one or more processors, a subsurface image from the velocity model if the convergence criterion is met to analyze subsurface regions. 21. The method of claim 20 further comprising: using, by the one or more processors, both the transformed prestack reflection data in the Laplace domain and the updated velocity model to minimize one of the logarithmic, the integral, and the power objective function in an iterative manner until the convergence criterion is met. 22. The method of claim 21 further comprising: initializing, by the one or more processors, a preliminary source wavelet; and updating, by the one or more processors, the source wavelet. 23. The method of claim 20 further comprising: confirming, by the one or more processors, that the convergence criterion is met when the value of one of the logarithmic, the integral, and the power objective function at a final iteration is found below a certain minimum and with a negligible difference from the values of the preceding iterations. 24. The method of claim 20 further comprising: generating, by the one or more processors, a subsurface image from the final updated velocity model using prestack depth migration models.
연구과제 타임라인
LOADING...
LOADING...
LOADING...
LOADING...
LOADING...
이 특허에 인용된 특허 (7)
Wang Shein S. (Ponca City OK) Cavanaugh Thomas D. (Ponca City OK) Hanson Douglas W. (Ponca City OK), Common offset depth migration with residual moveout correction.
Handley George E. (Golden CO) Moore A. M. (Houston TX) Dickman Robert E. (Evergreen CO), Method and apparatus for obtaining seismic vibrator reflection data.
Carrazzone James J. (Houston TX) Chang David (Houston TX) Lewis Catherine (Bellaire TX) Shah Pravin M. (Houston TX) Wang David Y. (Houston TX), Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.