Numerical simulation method and device for earthquake vector wave field

Abstract

The invention provides a numerical simulation method and device for an earthquake vector wave field. The numerical simulation method comprises the following steps: performing L-times auto-convolution operation on an original window function, thereby acquiring an auto-convolution window function; performing weighted array operation on the auto-convolution window function and the original window function, thereby acquiring an auto-convolution combined window function; cutting off a space convolution sequence of a pseudo-spectral method by the auto-convolution combined window function, thereby acquiring a finite difference operator; performing numerical simulation on the earthquake vector wave field according to the finite difference operator. According to the numerical simulation method provided by the invention, the technical problems of the prior art that the value dispersion of a numerical simulation result of the finite difference earthquake vector wave field is larger and the numerical simulation precision is low because the cutting of the imminent finite difference operator through the window function cannot consider the spectrum coverage area and precision error stability are solved, and the technical effect of effectively increasing the simulation precision of the numerical simulation for the earthquake vector wave field is achieved.

Description

technical field [0001] The invention relates to the technical field of numerical simulation, in particular to a method and device for numerical simulation of seismic vector wave field. Background technique [0002] With the deepening of oil and gas exploration, the demand for higher-precision imaging and inversion is becoming more and more urgent, and the numerical simulation technology of seismic wave field, which is the basis of imaging and inversion methods, has become particularly important. Numerical simulation of seismic wave field has always been a hot field of seismology research. It has important theoretical and practical significance for people to understand the propagation law of seismic wave, interpret actual seismic data, characterize subsurface medium structure domain lithology, and exploit earth resources. [0003] There are three main types of commonly used numerical simulation methods: geometric ray method, integral equation method, wave equation method, and...

AI Technical Summary

Problems solved by technology

[0012] At present, commonly used window functions include generalized weighted Hanning (Hamming) window, Gaussian window, improved binomial window, etc. In this case, it is difficult to guarantee the performance of the side lobe, which makes the numerical simulation of the seismic wave field not very effective.

For example, for the generalized weighted Hanning window, its essence is a trigonometric window. The amplitude response of the trigonometric window has a relatively narrow main lobe, but the attenuation of the side lobe is not high enough. The truncation approximation of the trigonometric window function is limited. For the difference operator, the precision error fluctuation of the truncation approximation is relatively large. Even if there is a large spectral coverage, the fluctuation of the precision error has a great influence on the effect of the approximation. The improved binomial window is an adjustable window function, but its main lobe is still too wide, resulting in too long transition zone, using the improved binomial window function to truncate and approximate the finite difference operator, its accuracy error spectrum coverage is small, especially for low-order finite difference operators (4, 8, 12th order), the accuracy error has large fluctuations

Images