Denoising methods and systems for preserving amplitude variation with offset (AVO) features of pre-stack seismic data

Patent metadata

Field	Value
Publication number	US-12487373-B2
Application number	US-202318329564-A
Country	US
Kind code	B2
Filing date	Jun 5, 2023
Priority date	Aug 2, 2022
Publication date	Dec 2, 2025
Grant date	Dec 2, 2025

Abstract

The present disclosure discloses a denoising method and a system for preserving amplitude variation with offset (AVO) features of pre-stack seismic data. The method includes: acquiring seismic wave data at at least one receiving point through a seismic signal acquisition device, and storing the seismic wave data in a memory; in response to ending of an acquisition operation of the seismic signal with the seismic signal acquisition device, determining, based on the seismic wave data in the memory, pre-denoising angle gather data through a processing device; and obtaining and outputting a denoised angle gather signal through inputting the pre-denoising angle gather data to a denoising device and performing denoising processing on the pre-denoising angle gather data based on the denoising device. In the present disclosure, an abnormal signal is identified in the memory, then the angle gather signal is denoised, and a user also use a user terminal to view signals before and after denoising and compare a denoising effect, which can improve denoising efficiency and accuracy.

First claim

What is claimed is:

1 . A denoising method for preserving amplitude variation with offset (AVO) features of pre-stack seismic data based on extended Radon transform, wherein the method comprises:
acquiring seismic wave data at at least one receiving point through a seismic signal acquisition device, and storing the seismic wave data in a memory, wherein the seismic signal acquisition device includes a seismic exploration instrument and a seismic geophone; in response to ending of an acquisition operation for a seismic signal with the seismic signal acquisition device, determining, based on the seismic wave data in the memory, pre-denoising angle gather data through a processing device, including:
identifying abnormal signal in the memory;
eliminating the abnormal signal; and
designating entire remaining seismic wave data as the pre-denoising angle gather data; and
obtaining and outputting a denoised angle gather signal through inputting the pre-denoising angle gather data to a denoising device and performing denoising processing on the pre-denoising angle gather data based on the denoising device, including: determining a range of t in an angle gather through performing division processing for overlapping time windows on an input angle gather; converting the input angle gather from an angle domain to a feature domain by a compressive sensing algorithm based on the range of t in the angle gather; obtaining a model parameter m(t,R 0 ,b) through inverting a formula (1) by a formula (2) based on the range of t in the angle gather, and converting the angle gather in the pre-stack seismic data from the angle domain to the feature domain using the model parameter m(t, R 0 , b);
D ( t ,θ)=Σ R 0 Σ b m ( t,R 0 ,b )*( R 0 +b sin 2 θ)  (1);

∥ Rm−d∥ 2 2 +λ∥m∥ 1   (2);

wherein D(t, θ) denotes angle gather data, R 0 +b sin 2 θ characterizes an AVO effect template that a seismic reflection amplitude changes with an angle, and the model parameter m(t, R 0 , b) in a transform space comprises AVO features of data; t denotes time, θ denotes an angle; d denotes the data; and R 0 , b is used to characterize the AVO features; obtaining an effective signal maintaining AVO features of original data by performing filtering processing on the angle gather in the feature domain, including:
selecting at least one set of candidate filtering parameters through vector matching based on a signal frequency distribution feature and a data noise feature in the pre-denoising angle gather data;
generating an evaluation result by evaluating a filtering effect of the at least one set of candidate filtering parameters, including:
determining the evaluation result by using an evaluation model based on the at least one set of candidate filtering parameter, an angle gather feature, and an overlapping parameter, wherein the evaluation model is a machine learning model, the evaluation model includes a filtering feature layer and a result evaluation layer, wherein
an input of the filtering feature layer includes the at least one set of candidate filtering parameter, the angle gather feature in the feature domain before filtering, and the overlapping parameter, and an output of the filtering feature layer includes the filtering feature;
an input of the result evaluation layer includes the filtering feature, the pre-denoising angle gather data, and smoothness of the overlapping parameter, and an output of the result evaluation layer includes a noise elimination rate and a signal loss rate corresponding to the candidate filtering parameter;
the output of the filtering feature layer is the input of the result evaluation layer, the result evaluation layer and the filtering feature layer are obtained through joint training; and
the evaluation result includes the noise elimination rate and the signal loss rate;

performing a plurality of rounds of iterative updating on the at least one set of the candidate filtering parameters based on the evaluation result to obtain an optimal filtering parameter;
performing a two-dimensional filtering processing on the pre-denoising angle gather data in the feature domain with the optimal filtering parameter to generate an effective signal that meets a preset condition, wherein the optimal filtering parameter includes at least one of a filter type and filtering strength;
obtaining optimized angle gather data in a corresponding time window through performing inverse conversion processing on the effective signal; obtaining a complete angle gather signal through performing splicing processing for overlapping time windows on the optimized angle gather data in a plurality of time windows, to realize denoising processing for preserving the AVO features of the pre-stack seismic data; and sending the denoised angle gather signal and the pre-denoising angle gather data to a user terminal for display and confirmation, wherein the user terminal displays the denoised angle gather signal and the pre-denoising angle gather data for users by displaying images and videos and according to information on the user terminal, the user visually views changes in angle gather signals before and after denoising to confirm a denoising effect to determine whether a denoising result is available.

2 . The denoising method according to claim 1 , wherein the model parameter m(t,R 0 ,b) is obtained through inverting the formula (1) based on an inversion algorithm of a compressive sensing theory.

3 . The denoising method according to claim 1 , wherein the effective signal maintaining the AVO features of the original data is obtained through performing the two-dimensional filtering processing on the angle gather in the feature domain.

4 . The denoising method according to claim 1 , wherein the optimized angle gather data in the corresponding time window is obtained through performing the inverse conversion processing on the effective signal using the formula (1).

5 . A denoising system for preserving AVO features of pre-stack seismic data based on extended Radon transform, wherein the system is applied to the denoising method of claim 1 , and the system comprises:
a seismic signal acquisition device configured to acquire seismic wave data at at least one receiving point, wherein the seismic signal acquisition device includes a seismic exploration instrument and a seismic geophone; a memory configured to store the seismic wave data; a processing device configured to determine pre-denoising angle gather data based on the seismic wave data, wherein to determine pre-denoising angle gather data based on the seismic wave data, the processing device is further configured to:
identify abnormal signal in the memory;
eliminate the abnormal signal; and
designate entire remaining seismic wave data as the pre-denoising angle gather data;
a denoising device configured to perform denoising processing on the pre-denoising angle gather data to obtain and output a denoised angle gather signal, including: a first processing module configured to determine a range of t in an angle gather through performing division processing for overlapping time windows on an input angle gather; a conversion module configured to convert the input angle gather from an angle domain to a feature domain by a compressive sensing algorithm based on the range of t in the angle gather; a filtering module configured to obtain an effective signal maintaining AVO features of original data by performing filtering processing on the angle gather in the feature domain, wherein to obtain the effective signal maintaining AVO features of the original data, the filtering module is further configured to:
select at least one set of candidate filtering parameters through vector matching based on a signal frequency distribution feature and a data

Assignees

• Univ Southwest Petroleum

Inventors

• Cao Weiping
• Yang Ran
• Huang Xuri
• Cui Xiaoqing
• Li Mengcheng
• Li Luo
• Li Haoyuan
• Ren Mengyu
• Yang Sheng
• Li Moyan
• Li Xinwang

Classifications

• G01V1/307 Primary

  for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity · CPC title

• G01V1/364 Primary

  Seismic filtering (G01V1/37 takes precedence) · CPC title

• G01V1/362 Primary

  Effecting static or dynamic corrections; Stacking · CPC title

• G01V2210/512

  Pre-stack · CPC title

• G01V2210/324

  Filtering · CPC title