Full Waveform Inversion Method in Time Domain Based on Zero-Mean Normalized Cross-Correlation Objective Function

A technology of full waveform inversion and objective function, applied in underground medium imaging and medium velocity calculation, time domain full waveform inversion, time domain full waveform inversion based on zero-mean normalized cross-correlation objective function, can solve Waveform matching error, not easy to skip cycle, low-frequency information cannot be received by the geophone, etc., to achieve the effect of eliminating low-frequency non-average noise, reducing amplitude energy error, and improving inversion accuracy

Inactive Publication Date: 2019-11-22
JILIN UNIV
View PDF5 Cites 0 Cited by
  • Summary
  • Abstract
  • Description
  • Claims
  • Application Information

AI Technical Summary

Problems solved by technology

[0004] The main problem faced by the full waveform inversion method is the cycle skip problem, that is, the waveform phase difference between the simulated record and the observed record is more than half a cycle, resulting in a waveform matching error, which leads to an error in the inversion result
There are two main reasons for the cycle skipping. One is that the initial model of the inversion is too far from the real underground model, resulting in a large difference between the simulated record waveform and the observed record waveform, which leads to the cycle skipping; the other is the lack of low-frequency information, the seismic records The low-frequency signal in the medium carries the information of the large-scale structure and is not easy to skip cycles
However, low-frequency information is often missing in actual observation records, and there are many reasons for the lack of low-frequency information, such as missing low-frequency information from seismic wavelets excited by the source, low-frequency information not being received by the geophone, or low-frequency information being artificially removed due to noise pollution, etc.
For the traditional full waveform inversion based on the least squares objective function, the requirements for amplitude information are relatively high, and the waveform amplitude error caused by source excitation or seismic record acquisition has a great impact on the inversion process

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
  • Full Waveform Inversion Method in Time Domain Based on Zero-Mean Normalized Cross-Correlation Objective Function
  • Full Waveform Inversion Method in Time Domain Based on Zero-Mean Normalized Cross-Correlation Objective Function
  • Full Waveform Inversion Method in Time Domain Based on Zero-Mean Normalized Cross-Correlation Objective Function

Examples

Experimental program
Comparison scheme
Effect test

Embodiment 1

[0082] According to the exploration requirements, Parallel Computing Toolbox and MATLAB Distributed Computing Server (R2016b) were installed under the Windows 10 Professional Edition system to build the MATLAB parallel platform. Full waveform inversion test with missing low-frequency information using Marmousi, the real model ( Figure 4 ) and the initial model ( Figure 5 ).

[0083] The model parameters are as follows:

[0084] Table 1 Test parameters of full waveform inversion with missing low frequency information

[0085]

[0086]

[0087] The grid size of the model is 69×192, the grid spacing is 12.5, the lateral distance is 2400m, and the vertical depth is 862.5m. The seismic wave velocity ranges from 1500m / s to 4000m / s in the model. The grid points are all one geophone, and the interval between the geophones is 12.5m, and a mixed source contains 12 sources. The seismic source selects the Reckel wavelet with a main frequency of 20Hz ( Figure 3a ), the sampli...

Embodiment 2

[0090] According to the exploration requirements, Parallel Computing Toolbox and MATLAB Distributed Computing Server (R2016b) were installed under the Windows 10 Professional Edition system to build the MATLAB parallel platform. Full waveform inversion test using Marmousi when the amplitude is wrong, the real model ( Figure 4 ) and the initial model ( Figure 5 ).

[0091] The model parameters are as follows:

[0092] Table 2. Full waveform inversion test parameters in the low frequency band of observation records with amplitude errors

[0093] grid size grid distance Line length vertical depth Fluc Bandwidth 69×192 12.5m 2400m 862.5m 20Hz 0~20Hz

[0094] The grid size of the model is 69×192, the grid spacing is 12.5, the lateral distance is 2400m, and the vertical depth is 862.5m. The seismic wave velocity ranges from 1500m / s to 4000m / s in the model. The grid points are all one geophone, and the interval between the geophones is 12.5m, ...

Embodiment 3

[0097] According to the exploration requirements, Parallel Computing Toolbox and MATLAB Distributed Computing Server (R2016b) were installed under the Windows 10 Professional Edition system to build the MATLAB parallel platform. Use Marmousi to perform full waveform inversion test when the observation record contains strong low-frequency noise, and the noise-added observation record ( Figure 9 ), the real model ( Figure 4 ) and the initial model ( Figure 5 ). .

[0098] The model parameters are as follows:

[0099] Table 3 Full waveform inversion test parameters in the low frequency band with strong low frequency noise

[0100] grid size grid distance Line length vertical depth Fluc Bandwidth 69×192 12.5m 2400m 862.5m 20Hz 0~20Hz

[0101] The grid size of the model is 69×192, the grid spacing is 12.5, the lateral distance is 2400m, and the vertical depth is 862.5m. The seismic wave velocity ranges from 1500m / s to 4000m / s in the model....

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 a time domain full-waveform inversion method based on zero-mean normalized cross-correlation objective function. Calculation of zero-mean normalized cross-correlation betweensimulation and observation records is used for replacing calculation of a least square objective function in the traditional full-waveform inversion. On the basis of the normalized cross-correlation,absolute values of amplitudes of the simulated data and the observed data are kept to be between -1 and 1, so that the inversion result error caused by the amplitude error is reduced effectively. Because of the normalization of the amplitude energy, the energy weight of the offset in the seismic record is improved; and since the offset data carry information of the large scale structure in the model, the occurrence of the cycle skipping is reduced. Because the actual seismic record is easy to affect by the low frequency noises and the noises are non-mean noises, a de-mean item is introduced; and the interference of low frequency noises is reduced by de-mean value calculation of the simulation and observation records. With the method provided by the invention, the velocity parameter of theunderground medium is inverted accurately under the circumstance of being lack of the low frequency information or noise pollution on the low frequency information.

Description

technical field [0001] The invention belongs to the technical field of seismic exploration, and relates to an underground medium imaging and medium velocity calculation method in the field of seismic exploration, in particular to a time-domain full waveform inversion method, in particular to a zero-mean normalized cross-correlation objective function The time domain full waveform inversion method. Background technique [0002] Non-renewable resources such as petroleum and minerals play a vital role in all aspects of human life. Therefore, the corresponding resource exploration technology is constantly innovating and developing. In recent years, with the increasing scale of land and offshore oil and gas development, exploration targets have changed from shallow simple structural geological areas to deep complex structural areas and hidden oil and gas reservoirs, and the difficulty of exploration has increased. In order to further solve the problem of deep energy exploration...

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 Patents(China)
IPC IPC(8): G01V1/30G01V1/36
CPCG01V1/303G01V1/364G01V2210/324G01V2210/6222
Inventor 董士琦韩立国陈瑞鼎胡勇张盼尹语晨
Owner JILIN UNIV
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