2 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.