Method of simulating the transmission rules of seismic waves in a TI medium

Abstract

The invention provides a method of simulating the transmission rules of seismic waves in a TI medium, and belongs to the seismic prospecting field. The method includes the steps: firstly, determining a six-order matrix coefficient on the basis of a to-be-solved problem; secondly, solving the first-order speed, i.e., a stress equation; thirdly, solving a convolution differential operator in the second step; fourthly, calculating a window function coefficient, and conducting windowing cross cutting on the convolution differential operator; fifthly, setting the parameter, the epicenter parameter and the epicenter frequency of a model, and carrying out discretization on the model according to geological bodies; and sixthly, conducting time iteration, obtaining the speed and the stress value within a whole area, and ending a simulation process and outputting a wave field snapshot and an earthquake record if the setting time is up.

Description

technical field [0001] The invention belongs to the field of seismic exploration, and specifically relates to a method for simulating the propagation law of seismic waves in a TI medium, which can simulate the propagation law of seismic waves in a transversely isotropic (TI) medium, and is used for forward simulation of wave equations to simulate seismic waves Propagation laws in the medium. Background technique [0002] Studying the propagation law of seismic waves in complex media is an important content in geophysics and an effective simulation method to help people understand the earth's media. [0003] Seismic wave numerical simulation is based on the known underground structure and physical parameters to simulate and study the propagation law of seismic waves in various underground media. It is widely used in the field of seismic exploration and natural earthquakes, but the anisotropy of actual underground media is ubiquitous. The wave equation can be used to describe...

AI Technical Summary

Problems solved by technology

[0003] Seismic wave numerical simulation is based on the known underground structure and physical parameters to simulate and study the propagation law of seismic waves in various underground media. It is widely used in the field of seismic exploration and natural earthquakes, but the anisotropy of actual underground media is ubiquitous. The wave equation can be used to describe the wave propagation process in the subsurface medium, so the numerical solution methods based on the wave equation, such as finite difference method, finite element method, pseudospectral method, spectral element method, etc., are widely used in the simulation of anisotropic media. These methods have their own advantages and disadvantages: the theoretical basis of the finite difference method is based on Taylor expansion, the calculation speed is fast and the efficiency is high, but it is difficult to deal with problems in irregular regions; the finite element method is based on the variational theory, the simulation accuracy is high, and it can handle Geological bodies with complex shapes, but the calculation process is complicated and the calculation amount is large; the pseudo-spectral method is also one of the numerical solutions to solve partial differential equations. Fast Fourier transform is introduced in the calculation process, and the accuracy can reach the accuracy of the spectrum. However, due to the Fourier transform Global information is required, and the pseudospectral method is incapable of dealing with wave field propagation in complex media; the spectral element method that combines the advantages of finite element and pseudospectral method also has problems such as the boundary and the new method are difficult to match.

Images