Search

CN-121995492-A - Automatic dynamic correction processing method, device and medium

CN121995492ACN 121995492 ACN121995492 ACN 121995492ACN-121995492-A

Abstract

The invention provides an automatic dynamic correction processing method, an automatic dynamic correction processing device and a medium, and belongs to the field of oil-gas geophysical exploration. The method comprises the steps of inputting a CMP gather of seismic data and arranging the CMP gather according to the sequence of offset distances, calculating the dynamic matching relation between the seismic channel with the minimum offset distance and the seismic channels of other channels, calculating the initial moveout static correction value of each seismic channel, constructing an objective function of the moveout correction value, and searching the optimal moveout correction value, wherein the step is 100. The invention shortens the time for picking up the velocity spectrum, improves the precision of picking up the velocity spectrum, and reduces the labor intensity of processing staff, thereby providing accurate velocity parameters for the subsequent seismic exploration processing flow.

Inventors

  • WANG XIAOPIN
  • DONG QIANQIAN

Assignees

  • 中国石油化工股份有限公司
  • 中石化石油物探技术研究院有限公司

Dates

Publication Date
20260508
Application Date
20241107

Claims (10)

  1. 1. An automatic correction processing method, comprising: Step 100, inputting CMP gathers of the seismic data, and arranging the CMP gathers according to the sequence of offset distances; step 200, calculating the dynamic matching relation between the seismic channel with the minimum offset and the seismic channels of other channels, and calculating the initial time difference static correction value of each seismic channel; Step 300, constructing an objective function of the time difference correction amount; step 400, searching for an optimal moveout correction amount.
  2. 2. The method of claim 1, wherein the step 100 of inputting CMP gathers of the seismic data and arranging the CMP gathers in the order of offset comprises: the CMP gathers of the seismic data are orderly arranged according to the order of the offset from small to large, and are marked as a two-dimensional signal mathematical expression as follows: S(t,n)=[s 1 ,s 2 ....s N ] (1) Wherein s 1 ....s N is the seismic traces arranged in order of the offset from small to large, and s 1 is the seismic trace with the smallest offset.
  3. 3. The method of claim 1, wherein the step 200 of calculating a dynamic matching relationship between the trace of the minimum offset and the trace of the remaining traces comprises: Extracting a seismic channel s i (i < =n, and i not equal to 1) and a seismic channel s 1 with the minimum offset from a CMP channel set of seismic data, calculating a time sequence matching relation of the seismic channels s 1 and s 1 by adopting a dynamic time sequence planning algorithm, wherein the dynamic matching result of s i and s 1 is as follows: D i ＝DTW(s 1 ,s i ) (2) wherein, DTW is a dynamic time sequence planning algorithm.
  4. 4. A method according to claim 3, wherein the calculating of the initial moveout static correction for each seismic trace in step 200 comprises: And (3) calculating the dynamic matching relation of the seismic channels s i and s 1 according to the formula (2), and further calculating the initial time difference correction quantity of each seismic channel.
  5. 5. The method of claim 1, wherein the objective function of the moveout correction constructed in step 300 is: Wherein, the Wherein σ represents the horizontal spatial distance between two seismic traces, s trj is the seismic trace to be corrected, i and j are the time series of s 1 ,s tr respectively, s 1 is the seismic trace with the minimum offset in the CMP trace set, and s tr is the data of any other seismic trace in the CMP trace set.
  6. 6. The method of claim 5, wherein searching for the optimal moveout correction amount in step 4 comprises: And (3) continuously optimizing the initial time difference intersection correction amount to optimize the objective function, and rewriting the formula (3) to obtain the final product: wherein t2 is the j-moment correction amount of the seismic trace to be corrected, and the optimal time difference correction amount can be optimized by a conjugate gradient method.
  7. 7. An automatic correction processing apparatus, comprising: the input unit is used for inputting the CMP gathers of the seismic data and arranging the CMP gathers according to the sequence of offset distances; the calculating unit is used for calculating the dynamic matching relation between the seismic channel with the minimum offset and the seismic channels of other channels and calculating the initial time difference static correction value of each seismic channel; an objective function construction unit for constructing an objective function of the time difference correction amount; And a time difference correction amount acquisition unit configured to search for an optimal time difference correction amount.
  8. 8. The apparatus according to claim 1, wherein the computing unit performs in particular the following operations: The dynamic time sequence planning algorithm is used for analyzing the CMP gather data, one seismic channel s i (i < =N, and i is not equal to 1) is extracted from the CMP gather of the seismic data, the seismic channel s 1 with the minimum offset is calculated by adopting the dynamic time sequence planning algorithm, the time sequence matching relation between the seismic channels s i and s 1 is calculated, and the dynamic matching result is as follows: D i ＝DTW(s 1 ,s i ) (2) the DTW is a dynamic timing planning algorithm, and is a conventional algorithm, which is not described herein. And (3) calculating the dynamic matching relation of the seismic channels s i and s 1 according to the formula (2), and further calculating the initial time difference correction quantity of each seismic channel.
  9. 9. The apparatus according to claim 7, wherein the objective function of the time difference correction amount constructed by the objective function constructing unit is: Wherein, the Wherein σ represents the horizontal spatial distance between two seismic traces, s trj is the seismic trace to be corrected, i and j are the time series of s 1 ,s tr respectively, s 1 is the seismic trace with the minimum offset in the CMP trace set, and s tr is the data of any other seismic trace in the CMP trace set.
  10. 10. A computer-readable storage medium storing at least one program executable by a computer, which when executed by the computer, causes the computer to perform the steps in the automatic correction processing method according to any one of claims 1 to 6.

Description

Automatic dynamic correction processing method, device and medium Technical Field The invention belongs to the field of oil and gas geophysical exploration, and particularly relates to an automatic dynamic correction processing method, an automatic dynamic correction processing device and a medium. Background With the advance of petroleum exploration towards deep ultra-deep layers, the collected seismic data gradually increase, the picking up of velocity spectrum in the seismic data dynamic correction process is particularly time-consuming and mostly is manual picking up, so that on one hand, too much effort is consumed for interpretation personnel, and on the other hand, the data processing period is prolonged. The traditional motion correction method at present is to divide the energy of the CMP track set to the calculated track set according to time-speed to form a speed spectrum, manually pick up the energy group to obtain the optimal speed according to time sequence, and the speed spectrum track set can be automatically leveled and overlapped to obtain simple overlapped data, so that the speed spectrum pick-up in the motion correction process is most important, the traditional speed spectrum pick-up method can be simply divided into semi-automatic speed spectrum pick-up and full-automatic pick-up, semi-automatic is mainly some nonlinear algorithms, the algorithms need a large amount of constraint information and are not suitable for complex geological scenes, and the error constraint condition can lead to larger result error. With the rise of artificial intelligence technology in recent years, deep learning technology is also used for dynamic correction, but the defect of deep learning is that a large amount of training data is needed, the geological structure of each work area is possibly different greatly for seismic data processing, the existing processing result is difficult to use in other work areas, and a comparison compromise method is to manually analyze the velocity spectrum of the data according to coarse grids such as CMP data extracted from 50×50 grids, prepare the processed velocity spectrum and superposition speed into training data to be brought into deep network training and then predict other CMP gather data, which is effectively degraded into semiautomatic velocity spectrum pickup. Disclosure of Invention The invention aims to solve the problems in the prior art and provide an automatic motion correction processing method, an automatic motion correction processing device and a medium, so that the time for picking up a velocity spectrum is shortened, the precision for picking up the velocity spectrum is improved, and accurate velocity parameters are provided for the subsequent seismic exploration processing flow. The invention is realized by the following technical scheme: In a first aspect of the present invention, there is provided an automatic correction processing method comprising: Step 100, inputting CMP gathers of the seismic data, and arranging the CMP gathers according to the sequence of offset distances; step 200, calculating the dynamic matching relation between the seismic channel with the minimum offset and the seismic channels of other channels, and calculating the initial time difference static correction value of each seismic channel; Step 300, constructing an objective function of the time difference correction amount; step 400, searching for an optimal moveout correction amount. The invention further improves that: In step 100, a CMP gather of seismic data is input, and the CMP gathers are arranged in an offset order, and the specific operations include: the CMP gathers of the seismic data are orderly arranged according to the order of the offset from small to large, and are marked as a two-dimensional signal mathematical expression as follows: S(t,n)=[s1,s2....sN] (1) Wherein s 1....sN is the seismic traces arranged in order of the offset from small to large, and s 1 is the seismic trace with the smallest offset. The invention further improves that: In step 200, the dynamic matching relationship between the seismic trace of the minimum offset and the seismic traces of the rest of the traces is calculated, and the specific operations include: Extracting a seismic channel s i (i < =n, and i not equal to 1) and a seismic channel s 1 with the minimum offset from a CMP channel set of seismic data, calculating a time sequence matching relation of the seismic channels s i and s 1 by adopting a dynamic time sequence planning algorithm, wherein the dynamic matching result of s i and s 1 is as follows: Di＝DTW(s1,si) (2) wherein, DTW is a dynamic time sequence planning algorithm. The invention further improves that: In step 200, an initial moveout static correction value of each seismic trace is calculated, and specific operations include: And (3) calculating the dynamic matching relation of the seismic channels s i and s 1 according to the formula (2), and further calculating