L-BFGS initial matrix solving method applied to least square reverse time migration

A technology of L-BFGS and least squares, which is applied in the field of seismic exploration migration imaging, can solve problems such as accelerated convergence speed, insufficient geometric diffusion compensation, and unbalanced illumination, and achieves improved imaging accuracy, uniform underground illumination, and accurate correction factors Effect

Active Publication Date: 2020-06-16
OCEAN UNIV OF CHINA
View PDF4 Cites 5 Cited by
  • Summary
  • Abstract
  • Description
  • Claims
  • Application Information

AI Technical Summary

Problems solved by technology

The model experiment results show that this method can effectively overcome the defects of insufficient geometric diffusion compensation and unbalanced illumination in the acoustic least squares reverse time migration method based on the conventional method, and significantly improve the imaging accuracy of the least squares reverse time migration. Accelerated the convergence speed, thus laying a good foundation for subsequent seismic data interpretation work

Method used

the structure of the environmentally friendly knitted fabric provided by the present invention; figure 2 Flow chart of the yarn wrapping machine for environmentally friendly knitted fabrics and storage devices; image 3 Is the parameter map of the yarn covering machine
View more

Image

Smart Image Click on the blue labels to locate them in the text.
Viewing Examples
Smart Image
  • L-BFGS initial matrix solving method applied to least square reverse time migration
  • L-BFGS initial matrix solving method applied to least square reverse time migration
  • L-BFGS initial matrix solving method applied to least square reverse time migration

Examples

Experimental program
Comparison scheme
Effect test

Embodiment 1

[0039] The present invention adopts the L-BFGS algorithm of optimizing the initial matrix to realize the least squares reverse time migration, and the specific process is as follows figure 2 shown. The present invention adopts the Marmousi model to describe the specific implementation. The Marmousi model is one of the commonly used standard models in the field of seismic exploration, and it is often used to test the feasibility and accuracy of algorithms under complex media conditions. In order to improve the calculation efficiency, under the premise of ensuring the complexity of the model, the original Marmousi model is thinned out. The Marmousi true velocity model after thinning is as follows image 3 As shown, Gaussian smoothing can be performed on it to obtain the thinned Marmousi background velocity model (such as Figure 4 shown), the theoretical reflection coefficient model can be obtained from the real velocity model and the background velocity model (the Marmousi ...

specific Embodiment approach

[0041] (1) Based on the thinned-out Marmousi velocity model, the seismic data collection is carried out in a full-array observation mode. There are 116 shots in total, 576 traces are received for each shot, the shot interval is 40m, and the trace interval is 8m. The shot points and receiver points are at Surface, with a depth of 0m. The least squares reverse time migration imaging is performed based on the velocity model and the acquired seismic data.

[0042] (2) In the kth iterative calculation of the least squares reverse time migration, k>=1, based on the background velocity model v obtained after velocity model smoothing 0 , taking the source wavelet w as the disturbance, according to formula (1), we can get the timing wave field P 0 .

[0043]

[0044] In formula (1), v 0 is the background velocity model, t is time, x and z are space coordinates respectively, P 0 is the simulated positive wave field, w is the source wavelet;

[0045] (3) Based on the positive wav...

the structure of the environmentally friendly knitted fabric provided by the present invention; figure 2 Flow chart of the yarn wrapping machine for environmentally friendly knitted fabrics and storage devices; image 3 Is the parameter map of the yarn covering machine
Login to view more

PUM

No PUM Login to view more

Abstract

The invention relates to an L-BFGS initial matrix solving method applied to least square reverse time migration, which belongs to the field of seismic exploration migration imaging. The method comprises the following steps of: in each iteration process of least square reverse time migration, based on a background speed model and a reflection coefficient model, taking shot point seismic source wavelets as disturbance, adopting a second-order scalar acoustic wave equation finite difference linearization forward modeling technology to obtain a timing wave field and a seismic record, and calculating a shot point one-way illumination compensation factor based on the timing wave field, performing reverse time continuation by taking the seismic record as disturbance to obtain a reverse time wavefield, obtaining a detection point one-way illumination compensation factor, and finally obtaining a shot-detection two-way illumination compensation factor. On the basis, the compensation factor is introduced into an L-BFGS algorithm as an initial approximate matrix of a Hessian inverse matrix, and an L-BFGS initial matrix solving method applied to least square reverse time migration is developed. According to the method, the least square reverse time migration imaging precision can be effectively improved, and the convergence speed is increased.

Description

technical field [0001] The invention belongs to the field of seismic exploration migration imaging, and in particular relates to a method for obtaining an L-BFGS initial matrix applied to least square reverse time migration. Background technique [0002] Seismic exploration is currently the most effective geophysical exploration method for proving underground geological structures and finding field oil and gas reservoirs. It mainly includes the following three links: field seismic data acquisition, seismic data processing, and seismic data interpretation. Seismic data processing is an important link connecting field seismic data acquisition and seismic data interpretation, and the migration imaging of underground structures is the core content of seismic data processing. Therefore, high-precision migration imaging technology has always been an important part of seismic exploration. research hotspots in the field. Compared with conventional reverse time migration, the least ...

Claims

the structure of the environmentally friendly knitted fabric provided by the present invention; figure 2 Flow chart of the yarn wrapping machine for environmentally friendly knitted fabrics and storage devices; image 3 Is the parameter map of the yarn covering machine
Login to view more

Application Information

Patent Timeline
no application Login to view more
Patent Type & Authority Applications(China)
IPC IPC(8): G01V1/28G01V1/34
CPCG01V1/282G01V1/345G01V2210/52G01V2210/74
Inventor 解闯宋鹏谭军李金山夏冬明姜秀萍赵波张超张锐埼张洪洋王绍文毛士博
Owner OCEAN UNIV OF CHINA
Who we serve
  • R&D Engineer
  • R&D Manager
  • IP Professional
Why Eureka
  • Industry Leading Data Capabilities
  • Powerful AI technology
  • Patent DNA Extraction
Social media
Try Eureka
PatSnap group products

Browse by: Latest US Patents, China's latest patents, Technical Efficacy Thesaurus, Application Domain, Technology Topic.

© 2024 PatSnap. All rights reserved.Legal|Privacy policy|Modern Slavery Act Transparency Statement|Sitemap