Multi-scale anisotropic diffusion filtering method based on pre-stack CRP trace sets

Abstract

The invention relates to a multi-scale anisotropic diffusion filtering method based on pre-stack CRP trace sets. The method includes the steps that firstly, the input pre-stack CRP trace sets are regularized; secondly, multi-scale decomposition is conducted through a two-dimensional Mallat algorithm, and each decomposed sub section is initialized; thirdly, a diffusion threshold value is obtained through a diffusion coefficient method, and the sixth step is executed; after diffusion tensor parameters are obtained based on anisotropic diffusion of a diffusion tensor method, the parameters are substituted into a nonlinear anisotropic equation so as to conduct anisotropic diffusion filtering, and the fourth step is executed; fourthly, the SNR, the MSE and the PSNR of the sub sections obtained after each iteration are calculated, and then the optimal earthquake sub section is optimized; fifthly, processing of the fourth step is executed on all the sub sections obtained after anisotropic diffusion filtering, and the sixth step is executed; sixthly, pre-stack CRP trace set nondestructive reconstruction is conducted on information of all the iterated sub sections through a two-dimensional Mallat reconstruction algorithm, and then the optimal pre-stack CRP trace set is output after reconstruction. The multi-scale anisotropic diffusion filtering method can be widely applied to the processing process of various oil exploration earthquake data.

Description

technical field [0001] The invention relates to a seismic data processing method in the oil exploration process, in particular to a multi-scale anisotropic diffusion filtering method based on pre-stack CRP (common reflection point) gathers. Background technique [0002] In actual production, with the improvement of seismic exploration precision, people have higher and higher requirements for the quality of seismic data. The "three highs" (high resolution, high signal-to-noise ratio, and high fidelity) processing technology of seismic data has always been a hot spot and difficulty in the field of geophysical exploration. Among them, the signal-to-noise ratio restricts the improvement of resolution, and the fidelity It directly affects the reliability of data interpretation results. As the exploration target gradually changes from conventional oil and gas reservoirs to lithological and subtle oil and gas reservoirs, pre-stack inversion has become an important method. The main...

AI Technical Summary

Problems solved by technology

Weickert introduced the diffusion tensor instead of the diffusion coefficient to extract the directional information of the image structure, which improved the P-M model's defect of losing effective information on the edge details, and later proposed a method of constructing the diffusion tensor with rotation invariance, which overcomes the time step Long limit, but the construction template is more complicated

Images