3 results about "Finite difference coefficient" patented technology

The application discloses a fractional-order viscous wave equation numerical simulation method, system and device, relates to the technical field of seismic wave propagation numerical simulation, and comprises the following steps: substituting a wave equation plane wave solution into a fractional-order Laplacian operator finite difference formula to obtain a generating function corresponding to the finite difference formula; replacing regular grid finite difference coefficients with optimized regular grid finite difference coefficients, performing inverse Fourier transform on the generating function to obtain optimized difference coefficients; based on the optimized difference coefficients of the finite difference formula and the fractional-order Laplacian operator finite difference formula, performing numerical dispersion on a preset fractional-order Laplacian viscous wave equation to obtain a fractional-order Laplacian viscous wave equation finite difference wave field recursion formula of seismic wave field data, and then simulating seismic response characteristics in an attenuation medium. The application can efficiently simulate numerical values of seismic response characteristics in the attenuation medium.

The invention discloses an elastic wave field forward modeling acceleration method and system based on parallel one-dimensional transverse convolution. The method comprises the following steps: S1, constructing a one-dimensional transverse convolution kernel according to a finite difference coefficient; s2, on the basis of the one-dimensional transverse convolution kernel, executing one-dimensional transverse convolution operation on the input matrix, and calculating to obtain a partial derivative in the x direction; s3, through matrix transposition processing, data arrangement in the z direction is converted into continuous data in the row direction, one-dimensional transverse convolution operation is executed based on the transposed input matrix, and a partial derivative in the z direction is obtained through calculation; and S4, on the basis of the calculated partial derivatives in the x direction and the z direction, elastic wave field forward modeling is executed. According to the technical scheme, high arithmetic strength and high calculation efficiency of convolution in the GPU are fully utilized, and the problem that the GPU reads data slowly according to columns is solved by using a matrix transposition technology.

The application discloses a seismic vector wave field numerical simulation method, device and medium, relates to the technical field of seismic vector wave field numerical simulation, and comprises the following steps: according to the Nyquist sampling theorem, a difference dispersion relation formula containing staggered grid finite difference coefficients is determined; an error limiting condition is introduced, and the staggered grid finite difference dispersion relation formula is converted into a target function to be optimized; an improved zebra optimization algorithm is used to solve the target function, and an optimized staggered grid finite difference operator is obtained; the first-order velocity-stress elastic wave equation is subjected to difference dispersion processing by using the optimized staggered grid finite difference operator, and according to the medium parameter model of a test region, a high-precision seismic vector wave field in a calculation region can be calculated. The improved zebra optimization algorithm is used to optimize the staggered grid finite difference operator, the seismic vector wave field numerical simulation precision can be improved, and a foundation is laid for a subsequent seismic imaging inversion process.