Method, device and storage medium for suppressing multiple waves

Abstract

The application provides a multiple wave suppression method, device, equipment and storage medium, and belongs to the technical field of seismic data processing. The method comprises the following steps: determining first seismic data of a plurality of first seismic trace sets; determining the curvature and frequency of each seismic trace in the first seismic trace set; determining the maximum curvature and minimum curvature in the curvatures of a plurality of seismic traces, and determining the maximum frequency in the frequencies of a plurality of seismic traces; generating a three-dimensional conical body filter based on the maximum curvature, the minimum curvature and the maximum frequency; performing three-dimensional Radon transform on the first seismic data of the first seismic trace set to obtain three-dimensional Radon transform data; filtering the three-dimensional Radon transform data based on the three-dimensional conical body filter to obtain three-dimensional Radon multiple wave data; and determining target seismic data after suppressing multiple waves based on the difference between the first seismic data of a plurality of first seismic traces and the three-dimensional Radon multiple wave data of a plurality of first seismic trace sets. The method improves the efficiency of suppressing multiple waves.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of seismic data processing, and particularly relates to a multiple wave suppression method and device, equipment and a storage medium. BACKGROUND

[0002] In the process of seismic data processing, seismic waves generally include multiple waves and primary waves. The primary waves are used for imaging, and the multiple waves are usually set as a kind of noise. The existence of the multiple waves will affect the reliability and authenticity of the primary wave imaging. In order to avoid the interference of the multiple waves on the primary wave imaging, the multiple waves in the seismic waves need to be suppressed.

[0003] In the related art, the three-dimensional Radon transform method is generally used to suppress the multiple waves. In the method, the frequency separation boundary of the multiple waves and the primary waves is artificially defined, and the frequencies of the multiple waves and the primary waves are respectively mapped to the two sides of the frequency separation boundary, so as to realize the division of the multiple waves and the primary waves, and then the suppression of the multiple waves is realized based on the divided multiple waves. Since the frequency separation boundary of the multiple waves and the primary waves is artificially defined, a large amount of manpower is consumed when a large amount of seismic data is processed, and the calculation efficiency of suppressing the multiple waves is reduced. SUMMARY

[0004] The embodiments of the present application provide a multiple wave suppression method, device, equipment and storage medium, which can improve the efficiency of suppressing the multiple waves. The technical solution is as follows:

[0005] In one aspect, a multiple wave suppression method is provided, and the method comprises:

[0006] determining first seismic data of a plurality of first seismic trace sets. For each first seismic trace set, the first seismic trace set comprises a plurality of seismic traces, and the first seismic data of the first seismic trace set comprises coordinates of each seismic trace at each sampling time;

[0007] determining the curvature and the frequency of each seismic trace in the first seismic trace set;

[0008] determining the maximum curvature and the minimum curvature in the curvatures of the plurality of seismic traces, and determining the maximum frequency in the frequencies of the plurality of seismic traces;

[0009] generating a three-dimensional conical filter based on the maximum curvature, the minimum curvature and the maximum frequency;

[0010] performing three-dimensional Radon forward transform on the first seismic data of the first seismic trace set to obtain three-dimensional Radon forward transform data;

[0011] performing filtering processing on the three-dimensional Radon forward transform data based on the three-dimensional conical filter to obtain three-dimensional Radon multiple wave data;

[0012] determine target seismic data after suppressing multiples based on a difference between the first seismic data of the plurality of first seismic traces and the three-dimensional Radon multiple wave data of the plurality of first seismic traces.

[0013] In a possible implementation, the generating the three-dimensional pyramid filter based on the maximum curvature, the minimum curvature and the maximum frequency comprises:

[0014] determining a product of the maximum curvature and the maximum frequency to obtain a first frequency maximum value of the three-dimensional pyramid filter;

[0015] determining a product of the minimum curvature and the maximum frequency to obtain a first frequency minimum value of the three-dimensional pyramid filter;

[0016] generating the three-dimensional pyramid filter based on the first frequency maximum value and the first frequency minimum value, and a multiple frequency preservation area of the three-dimensional pyramid filter is an interval composed of the first frequency maximum value and the first frequency minimum value.

[0017] In a possible implementation, the performing three-dimensional Radon forward transformation on the first seismic data of the first seismic trace set to obtain the three-dimensional Radon forward transformation data comprises:

[0018] performing one-dimensional Fourier forward transformation on the first seismic data to obtain to-be-transformed seismic data;

[0019] determining a three-dimensional Radon transformation operator of the first seismic trace set based on the curvature of each seismic trace, the maximum frequency and the to-be-transformed seismic data;

[0020] determining a generalized inverse operator matched with the three-dimensional Radon transformation operator;

[0021] performing three-dimensional Radon forward transformation on the to-be-transformed seismic data based on the generalized inverse operator to obtain the three-dimensional Radon forward transformation data.

[0022] In a possible implementation, the determining the three-dimensional Radon transformation operator of the first seismic trace set based on the curvature of each seismic trace, the maximum frequency and the to-be-transformed seismic data comprises:

[0023] determining a three-dimensional Radon transformation parameter of each seismic trace based on the curvature of each seismic trace and the maximum frequency;

[0024] determining the three-dimensional Radon transformation operator based on the three-dimensional Radon transformation parameter of each seismic trace and the to-be-transformed seismic data through first relationship data, the first relationship data being relationship data among a three-dimensional Radon transformation operator, a three-dimensional Radon transformation parameter and to-be-transformed seismic data.

[0025] In a possible implementation, the target seismic data after suppressing the multiple waves is determined based on a difference between the first seismic data of the plurality of first seismic gathers and the three-dimensional Rayleigh multiple wave data of the plurality of first seismic gathers, including:

[0026] The plurality of three-dimensional Rayleigh multiple wave data are deprocessed to obtain a plurality of dynamic correction multiple wave data;

[0027] The plurality of second seismic data corresponding to the plurality of first seismic data are obtained;

[0028] A difference between the plurality of second seismic data and the plurality of dynamic correction multiple wave data is determined to obtain the target seismic data.

[0029] In a possible implementation, the plurality of three-dimensional Rayleigh multiple wave data are deprocessed to obtain a plurality of dynamic correction multiple wave data, including:

[0030] For each three-dimensional Rayleigh multiple wave data, three-dimensional Rayleigh inverse transformation is performed on the three-dimensional Rayleigh multiple wave data to obtain to-be-transformed multiple wave data;

[0031] One-dimensional Fourier inverse transformation is performed on the to-be-transformed multiple wave data to obtain regularized multiple wave data;

[0032] The regularized multiple wave data are subjected to de-regularization processing and inverse moveout correction processing to obtain the dynamic correction multiple wave data.

[0033] In a possible implementation, the first seismic data of the plurality of first seismic gathers is determined, including:

[0034] Original seismic data are obtained;

[0035] The original seismic data are subjected to dynamic correction preprocessing to obtain second seismic data of a plurality of second seismic gathers;

[0036] For each second seismic gather, grid regularization processing and moveout correction processing are performed on the second seismic data of the second seismic gather to obtain the first seismic data of the first seismic gather.

[0037] In another aspect, a multiple wave suppression device is provided, including:

[0038] A first determination module is configured to determine first seismic data of a plurality of first seismic gathers, for each first seismic gather, the first seismic gather including a plurality of seismic traces, and the first seismic data of the first seismic gather including coordinates of each seismic trace at each sampling time;

[0039] a second determining module, configured to determine a curvature and a frequency of each seismic trace in the first seismic trace set;

[0040] a third determining module, configured to determine a maximum curvature and a minimum curvature in the curvatures of the plurality of seismic traces, and determine a maximum frequency in the frequencies of the plurality of seismic traces;

[0041] a generating module, configured to generate a three-dimensional spheroid filter based on the maximum curvature, the minimum curvature and the maximum frequency;

[0042] a transforming module, configured to perform three-dimensional Radon forward transformation on first seismic data of the first seismic trace set to obtain three-dimensional Radon forward transformation data;

[0043] a filtering module, configured to perform filtering processing on the three-dimensional Radon forward transformation data based on the three-dimensional spheroid filter to obtain three-dimensional Radon multiple wave data;

[0044] a fourth determining module, configured to determine target seismic data after suppressing multiple waves based on a difference between the first seismic data of the plurality of first seismic traces and the three-dimensional Radon multiple wave data of the plurality of first seismic trace sets.

[0045] In a possible implementation, the generating module comprises:

[0046] a first determining unit, configured to determine a product of the maximum curvature and the maximum frequency to obtain a first frequency maximum value of the three-dimensional spheroid filter;

[0047] a second determining unit, configured to determine a product of the minimum curvature and the maximum frequency to obtain a first frequency minimum value of the three-dimensional spheroid filter;

[0048] a generating unit, configured to generate the three-dimensional spheroid filter based on the first frequency maximum value and the first frequency minimum value, and a multiple wave frequency reservation area of the three-dimensional spheroid filter is an interval composed of the first frequency maximum value and the first frequency minimum value.

[0049] In a possible implementation, the transforming module comprises:

[0050] a first transforming unit, configured to perform one-dimensional Fourier forward transformation on the first seismic data to obtain to-be-transformed seismic data;

[0051] a third determining unit, configured to determine a three-dimensional Radon transformation operator of the first seismic trace set based on the curvature of each seismic trace, the maximum frequency and the to-be-transformed seismic data;

[0052] a fourth determining unit, configured to determine a generalized inverse operator matched with the three-dimensional Radon transformation operator;

[0053] a second transformation unit, configured to perform a three-dimensional Radon forward transformation on the to-be-transformed seismic data based on the generalized inverse operator to obtain three-dimensional Radon forward transformation data.

[0054] In a possible implementation, the third determination unit comprises:

[0055] a first determination sub-unit, configured to determine a three-dimensional Radon transformation parameter of each seismic trace based on the curvature of the seismic trace and the maximum frequency;

[0056] a second determination sub-unit, configured to determine the three-dimensional Radon transformation operator based on the three-dimensional Radon transformation parameter of each seismic trace and the to-be-transformed seismic data through first relationship data, the first relationship data being relationship data among the three-dimensional Radon transformation operator, the three-dimensional Radon transformation parameter, and the to-be-transformed seismic data.

[0057] In a possible implementation, the fourth determination module comprises:

[0058] a first processing unit, configured to perform inverse processing on the plurality of three-dimensional Radon multiple wave data to obtain a plurality of moveout multiple wave data;

[0059] a first acquisition unit, configured to acquire a plurality of second seismic data corresponding to a plurality of first seismic data;

[0060] a fifth determination unit, configured to determine a difference between the plurality of second seismic data and the plurality of moveout multiple wave data to obtain the target seismic data.

[0061] In a possible implementation, the processing unit comprises:

[0062] a third transformation sub-unit, configured to perform a three-dimensional Radon inverse transformation on each three-dimensional Radon multiple wave data to obtain to-be-transformed multiple wave data;

[0063] a fourth transformation sub-unit, configured to perform one-dimensional Fourier inverse transformation on the to-be-transformed multiple wave data to obtain regularized multiple wave data;

[0064] a processing sub-unit, configured to perform inverse regularization processing and inverse moveout correction processing on the regularized multiple wave data to obtain the moveout multiple wave data.

[0065] In a possible implementation, the first determination module comprises:

[0066] a second acquisition unit, configured to acquire original seismic data;

[0067] a second processing unit, configured to perform pre-processing on the original seismic data to obtain second seismic data of a plurality of second seismic gathers;

[0068] a third processing unit, configured to, for each second seismic gather, perform grid regularization processing and moveout correction processing on the second seismic data of the second seismic gather to obtain the first seismic data of the first seismic gather.

[0069] In another aspect, a computer device is provided, which includes one or more processors and one or more memories having stored therein at least one instruction, which is loaded and executed by the one or more processors to implement operations performed by the method for suppressing multiple waves according to any of the above implementation manners.

[0070] In another aspect, a computer readable storage medium is provided, which has stored therein at least one instruction, which is loaded and executed by a processor to implement operations performed by the method for suppressing multiple waves according to any of the above implementation manners.

[0071] In another aspect, a computer program product or a computer program is provided, which includes computer program code stored in a computer readable storage medium. A processor of a computer device reads the computer program code from the computer readable storage medium, and the processor executes the computer program code to cause the computer device to perform operations performed by the method for suppressing multiple waves according to the above.

[0072] The technical scheme provided by the embodiments of the present application has at least the following beneficial effects:

[0073] The method for suppressing multiple waves provided by the embodiments of the present application avoids the process of artificially defining a three-dimensional spine filter, improves the efficiency of generating a three-dimensional spine filter, and then obtains three-dimensional Rayleigh multiple wave data through the three-dimensional spine filter, which improves the efficiency of obtaining three-dimensional Rayleigh multiple wave data, and further suppresses multiple waves based on the difference between the first seismic data and the three-dimensional Rayleigh multiple wave data, which improves the efficiency of suppressing multiple waves. BRIEF DESCRIPTION OF DRAWINGS

[0074] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings needed in the embodiment description. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0075] Figure 1 is a flowchart of a method for suppressing multiple waves provided by an embodiment of the present application;

[0076] Figure 2 is a lateral profile of a second seismic gather provided by an embodiment of the present application;

[0077] Figure 3 is a schematic diagram of a three-dimensional vertebral filter provided by an embodiment of the present application;

[0078] Figure 4 is a schematic diagram of three-dimensional Radon transformed seismic data provided by an embodiment of the present application;

[0079] Figure 5 is a schematic diagram of three-dimensional Radon transformed seismic data provided by an embodiment of the present application;

[0080] Figure 6 is a schematic diagram of a three-dimensional Radon transformed multiple wave gather provided by an embodiment of the present application;

[0081] Figure 7 is a schematic diagram of a three-dimensional Radon transformed multiple wave gather provided by an embodiment of the present application;

[0082] Figure 8 is a lateral profile of a seismic gather after suppressing multiple waves provided by an embodiment of the present application;

[0083] Figure 9 is a lateral profile of a seismic gather after suppressing multiple waves provided by an embodiment of the present application;

[0084] Figure 10 is a lateral profile of a second seismic gather provided by an embodiment of the present application;

[0085] Figure 11 is a schematic diagram of energy distribution in a three-dimensional Radon transform domain provided by an embodiment of the present application;

[0086] Figure 12 is a lateral profile of a gather of multiple wave data after dynamic correction provided by an embodiment of the present application;

[0087] Figure 13 is a lateral profile of a gather of target seismic data provided by an embodiment of the present application;

[0088] Figure 14 is a schematic diagram of a multiple wave suppression device provided by an embodiment of the present application;

[0089] Figure 15 is a block diagram of a computer device provided by an embodiment of the present application. DETAILED DESCRIPTION

[0090] For the purpose of making the purpose, technical solutions and advantages of the present application clearer, the embodiments of the present application will be described in further detail below with reference to the drawings.

[0091] The terms "first", "second", "third", and "fourth" and the like in the description and claims of the present application and the accompanying drawings are used to distinguish different objects, and are not used to describe a particular order. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but can optionally include other steps or units not listed, or can optionally include other steps or units inherent to the process, method, product or device.

[0092] The embodiments of the present application provide a multiple wave suppression method, see Figure 1 , the method comprising:

[0093] Step 101: The computer device determines the first seismic data of a plurality of first seismic trace sets.

[0094] Wherein, for each first seismic trace set, the first seismic trace set includes a plurality of seismic traces, and the first seismic data of the first seismic trace set includes the coordinates of each seismic trace at each sampling time. The sampling time is the travel time of the seismic wave in the stratum, that is, the time from the shot point emitting the seismic wave to the seismic wave reaching the geophone.

[0095] This step can be realized by the following steps (1)-(3):

[0096] (1) The computer device acquires the original seismic data.

[0097] Wherein, the original seismic data is the seismic data to be suppressed multiple waves. The original seismic data includes the original seismic data of a plurality of original seismic traces, and the original seismic data of the plurality of original seismic traces includes the coordinates at each sampling time, which includes the horizontal coordinate and the vertical coordinate.

[0098] (2) The computer device performs moveout correction preprocessing on the original seismic data to obtain the second seismic data of a plurality of second seismic trace sets.

[0099] Wherein, for each second seismic trace set, the second seismic trace set includes a plurality of seismic traces, and the second seismic data of the second seismic trace set includes the coordinates of each seismic trace at each sampling time.

[0100] It should be noted that the computer device also performs the cutting processing on the stretched and distorted seismic trace at the far offset position in each second seismic trace set, and the second seismic data in the second seismic trace set is represented as d0(x, y, t). x and y are the horizontal coordinate and the vertical coordinate of any seismic trace in the second seismic trace set respectively, x is the coordinate in the in-line direction of the original seismic data, y is the coordinate perpendicular to the in-line direction of the original seismic data; and t is the sampling time corresponding to the seismic trace.

[0101] Referring to Figure 2 , Figure 2 is a lateral profile of a second seismic trace set, the second seismic trace set cuts the stretched and distorted seismic trace at the near offset position; the lateral profile extracts four slices for arrangement and display, and it can be seen from the figure that the primary wave event in the second seismic trace set is corrected, and the first-order and second-order multiples are in a parabolic shape.

[0102] (3) The computer device performs the grid regularization processing and the time difference correction processing on the second seismic data in each second seismic trace set to obtain the first seismic data in the first seismic trace set.

[0103] In this embodiment of the present application, the grid regularization processing and the time difference correction processing are performed on the first seismic data in the second seismic trace set, and a regular and complete first seismic trace set and first seismic data can be obtained.

[0104] It should be noted that the time difference correction is needed for each seismic trace in the grid regularization processing, and the correction formula is as follows:

[0105] In this embodiment of the present application, the grid regularization processing and the time difference correction processing are performed on the first seismic data in the second seismic trace set, and a regular and complete first seismic trace set and first seismic data can be obtained.

[0106] In this embodiment of the present application, the grid regularization processing and the time difference correction processing are performed on the first seismic data in the second seismic trace set, and a regular and complete first seismic trace set and first seismic data can be obtained.

[0107] Step 102: The computer device determines the curvature and the frequency of each seismic trace in the first seismic trace set.

[0108] It should be noted that the curvature of each seismic trace includes a lateral curvature and a longitudinal curvature. Continuing to refer to Figure 2 , Figure 2 is a slice of a lateral profile corresponding to a lateral curvature, which includes a flat axis of primary waves and a parabolic axis of multiple waves, the flat axis and the parabolic axis are respectively composed of a plurality of seismic traces, the curvature of the flat axis is substantially close to zero, and the curvature of any point of the parabolic axis is different, i.e., each seismic trace corresponds to a curvature. The frequency of each seismic trace is the inverse of the sampling time.

[0109] wherein the computer device can directly obtain the apparent velocity of the parabolic axis of each seismic trace from the first set of seismic traces, and take the inverse of the square of the apparent velocity as the curvature of the seismic trace, see Formula I.

[0110] Formula I:

[0111] wherein q is the curvature, and v is the apparent velocity.

[0112] Step 103: The computer device determines the maximum curvature and the minimum curvature in the curvatures of the plurality of seismic traces, and determines the maximum frequency in the frequencies of the plurality of seismic traces.

[0113] It should be noted that each seismic trace corresponds to a curvature and a frequency; the computer device determines the maximum value and the minimum value in the plurality of curvatures as the maximum curvature and the minimum curvature, respectively. The computer device determines the maximum value in the plurality of frequencies as the maximum frequency.

[0114] wherein the maximum value of the curvature is the maximum value corresponding to the apparent velocity, and the minimum value of the curvature is the minimum value corresponding to the apparent velocity, see Formula II.

[0115] Formula II:

[0116] wherein q max is the maximum curvature, v max is the maximum apparent velocity, q min is the minimum curvature, and v min is the minimum apparent velocity.

[0117] Step 104: The computer device generates a three-dimensional pyramid filter based on the maximum curvature, the minimum curvature, and the maximum frequency.

[0118] This step can be implemented through the following steps (1)-(3):

[0119] (1) The computer device determines the product of the maximum curvature and the maximum frequency to obtain a first frequency maximum value of the three-dimensional pyramid filter.

[0120] (2) The computer device determines the product of the minimum curvature and the maximum frequency to obtain a first frequency minimum value of the three-dimensional pyramid filter.

[0121] (3) The computer device generates a three-dimensional cone filter based on the first maximum value and the first minimum value of the first frequency. The multiple frequency retention area of ​​the three-dimensional cone filter is the interval composed of the first maximum value and the first minimum value of the first frequency.

[0122] The expression for the three-dimensional cone filter is:

[0123] Where, q max For maximum curvature, q min Let f be the minimum curvature, f be the maximum frequency, and q be the minimum curvature. max f is the maximum value of the first frequency, q min f is the minimum value of the first frequency, and r is the radius of the base of the vertebral body.

[0124] See Figure 3 , Figure 3 This is a schematic diagram of a three-dimensional cone filter. The generatrix length of the outer cone of the three-dimensional cone filter is the maximum curvature q. max The generatrix length of the inner cone is the minimum curvature q. min The length of the cone's axis is the maximum frequency.

[0125] In this design, the multiple frequency retention region of the 3D cone filter is the area inside the cone, i.e., the retention region and transition region shown in the figure. The primary frequency cutoff region of the 3D cone filter is the area outside the cone, i.e., the cutoff region shown in the figure. Based on the spatial distribution characteristics of the primary and multiple frequencies mapped to the fast 3D Radon transform domain, the frequencies of the seismic traces mapped to the retention and transition regions of the 3D cone filter are the frequencies of the multiple waves, while the frequencies of the seismic traces mapped to the cutoff region of the 3D cone filter are the frequencies of the primary waves.

[0126] In this embodiment, a three-dimensional cone filter is generated for multiple wave separation. By utilizing the range of three-dimensional Radon transform parameters within each frequency and the apparent velocity of the multiple waves, the multiple wave retention region and the primary wave cut-off region of the three-dimensional cone filter are defined, avoiding the step of manually defining each region of the three-dimensional cone filter, which is more convenient for industrial production applications.

[0127] Step 105: The computer equipment performs a three-dimensional Radon transform on the first seismic data of the first seismic gather to obtain three-dimensional Radon transform data.

[0128] This step can be achieved through the following steps (1)-(4):

[0129] (1) The computer equipment performs a one-dimensional Fourier forward transform on the first seismic data to obtain the seismic data to be transformed.

[0130] It should be noted that the one-dimensional Fourier forward transform is a transform along the time direction, that is, the first seismic data in the time domain is transformed into the to-be-transformed seismic data in the frequency domain. The first seismic data is represented as d1(x, y, t), and the to-be-transformed seismic data is represented as D(x, y, f), and f is the frequency.

[0131] (2) The computer device determines the three-dimensional Radon transform operator of the first seismic trace set based on the curvature, the maximum frequency and the to-be-transformed seismic data of each seismic trace.

[0132] This step can be implemented through the following steps A1-A2:

[0133] A1: The computer device determines the three-dimensional Radon transform parameter of each seismic trace based on the curvature and the maximum frequency of each seismic trace.

[0134] The three-dimensional Radon transform parameter includes a lateral transform parameter and a longitudinal transform parameter. The computer device determines the product of the lateral curvature and the maximum frequency to obtain the lateral transform parameter; and determines the product of the longitudinal curvature and the maximum frequency to obtain the longitudinal transform parameter. The three-dimensional Radon transform parameter is obtained through the following formula three.

[0135] Formula three: λ x = q x f, λ y = q y f

[0136] Wherein, λ x is the lateral transform parameter, λ y is the longitudinal transform parameter, q x is the lateral curvature, q y is the longitudinal curvature, and f is the maximum frequency.

[0137] In the embodiments of the present application, since the curvature is inversely proportional to the frequency, that is, the lower the frequency, the larger the curvature scanning range, so as to increase the curvature scanning range by reducing the frequency, so as to alleviate the spatial truncation effect caused by the too small curvature scanning range, thereby avoiding the energy group tailing phenomenon caused by the spatial truncation effect.

[0138] A2: The computer device determines the three-dimensional Radon transform operator based on the three-dimensional Radon transform parameter of each seismic trace and the to-be-transformed seismic data through the first relationship data.

[0139] The first relationship data is the relationship data between the three-dimensional Radon transform operator, the three-dimensional Radon transform parameter and the to-be-transformed seismic data.

[0140] It should be noted that the three-dimensional Radon transform operator includes a lateral transform operator and a longitudinal transform operator, and the first relationship data includes lateral relationship data and longitudinal relationship data. The to-be-transformed seismic data includes the coordinates corresponding to each frequency of each seismic trace.

[0141] The lateral relationship data is:

[0142] The longitudinal relationship data is:

[0143] wherein, is a lateral transform operator, λ1, λ2…λ n is a lateral transform parameter of each seismic trace, x1, x2…x n is the abscissa of each seismic trace; is a longitudinal transform operator, λ1, λ2…λ k is a longitudinal transform parameter of each seismic trace, y1, y2…y n is the ordinate of each seismic trace; m is the number of traces in the lateral direction of the first seismic trace set, p is the number of traces in the longitudinal direction of the first seismic trace set, n is the lateral dimension of the three-dimensional Radon transform domain model space, k is the longitudinal dimension of the three-dimensional Radon transform domain model space, and the sizes of the lateral dimension n and the longitudinal dimension k can be set and changed as needed.

[0144] wherein the computer device substitutes the lateral transform parameter and the abscissa of each seismic trace of the seismic data to be transformed into the lateral relationship data to obtain the lateral transform operator, and substitutes the longitudinal transform parameter and the ordinate of each seismic trace of the seismic data to be transformed into the longitudinal relationship data to obtain the longitudinal transform operator.

[0145] In the embodiment of the present application, since the expression of the conventional three-dimensional Radon transform operator is and It can be seen that the conventional three-dimensional Radon transform operator is closely related to frequency, and for any different frequency, the three-dimensional Radon transform operator needs to be recalculated, which is time-consuming and laborious. In the embodiment of the present application, the three-dimensional Radon transform operator is redefined, and the frequency is not involved in the three-dimensional Radon transform operator, so that the three-dimensional Radon transform operator only needs to be calculated once within all frequencies, which saves time and effort, thereby improving the calculation efficiency.

[0146] In the embodiment of the present application, since the first seismic data is data after regular processing, the offset distances in the lateral and longitudinal directions are constant, and then the three-dimensional Radon operator determined based on the first seismic data is relatively stable, and the false frequency phenomenon is not introduced due to the non-uniformity of the spatial position, so that the energy group of the three-dimensional Radon transform domain is more convergent.

[0147] (3) The computer device determines a generalized inverse operator matched with the three-dimensional Radon transform operator.

[0148] The computer device performs conjugate transformation on the matrix of the three-dimensional Radon transformation operator to obtain a conjugate matrix, the conjugate matrix is a conjugate operator, and inverse transformation is performed on the conjugate operator and the matrix of the three-dimensional Radon transformation operator to obtain an inverse matrix, the inverse matrix is a generalized inverse operator.

[0149] It should be noted that the three-dimensional Radon transformation operator includes a transverse transformation operator and a longitudinal transformation operator The conjugate operator includes a transverse conjugate operator and a longitudinal conjugate operator The generalized inverse operator includes a transverse inverse operator and a longitudinal inverse operator

[0150] The computer device performs conjugate transformation on the matrix of the transverse transformation operator to obtain a transverse conjugate operator, and performs inverse transformation on the matrix of the transverse transformation operator and the transverse conjugate operator to obtain a transverse inverse operator. The computer device performs conjugate transformation on the matrix of the longitudinal transformation operator to obtain a longitudinal conjugate operator, and performs inverse transformation on the matrix of the longitudinal transformation operator and the longitudinal conjugate operator to obtain a longitudinal inverse operator.

[0151] In the embodiment of the present application, the three-dimensional Radon operator is reconstructed by using the three-dimensional Radon transformation parameters, the dependence of the three-dimensional Radon operator on frequency is eliminated, the same three-dimensional Radon operator can be used when calculating at all frequencies, and the conjugate operator and the generalized inverse operator of the three-dimensional Radon operator only need to be calculated once, a large number of matrix inverse operation processes are reduced, and the solution of the inverse operation is the most time-consuming part of the entire three-dimensional Radon transformation, so that the calculation efficiency can be greatly improved to achieve the purpose of high calculation efficiency.

[0152] (4) The computer device performs three-dimensional Radon forward transformation on the seismic data to be transformed based on the generalized inverse operator to obtain three-dimensional Radon forward transformation data.

[0153] In this step, the computer device generates a transverse correction inverse operator and a longitudinal correction inverse operator based on the transverse inverse operator and the longitudinal inverse operator; and generates second relationship data based on the transverse correction inverse operator and the longitudinal correction inverse operator, two parameters of the second relationship data are the transverse correction inverse operator and the longitudinal correction inverse operator respectively, the independent variable is the seismic data to be transformed, and the dependent variable is the three-dimensional Radon forward transformation data.

[0154] The second relationship data is:

[0155]

[0156] wherein, M(λ x ,λ yf) represents the three-dimensional Radon transform data, which represents the transverse transform parameters, longitudinal transform parameters, and frequency of each seismic trace. D(x,y,f) represents the seismic data to be transformed, which represents the abscissa, ordinate, and frequency of each seismic trace. For the transverse correction inverse operator, The longitudinal correction inverse operator, μ, is the damping factor ranging from 0.01 to 1, and I is the identity matrix. μI is mainly used to enhance the stability of matrix inversion. Figure 2 Taking the second seismic gather as an example, μ is set to 0.01.

[0157] See Figure 4 , Figure 4 This is a profile of seismic data after conventional 3D Lardon positive transformation, from which four longitudinal curvatures q are extracted. y The same slice, from Figure 4 As can be seen, the seismic data D(x,y,f) to be transformed on the same phase axis is mapped to M(q) in the transform domain. x ,q y After f), the curves are separated into "points" in the transform domain, but due to the limitation of the three-dimensional acquisition aperture, i.e. the curvature scanning range is small, a spatial truncation effect is generated, which in turn produces the energy tailing phenomenon.

[0158] See Figure 5 , Figure 5 This is a profile of seismic data after three-dimensional Radon transform using the method provided in this application embodiment. Four longitudinal transform parameters λ are extracted from this profile. y The same slice, from Figure 5 As can be seen, the seismic data D(x,y,f) to be transformed along the same phase axis is mapped to the transform domain M(λ). x ,λ y After f), the energy is distributed in a straight line with different curvatures in the transform domain. The vertical distribution is the primary wave energy, and the inclined distribution is the secondary wave energy. The boundary between the primary wave energy and the secondary wave energy is obvious, the energy is relatively convergent, and the separability between the primary wave energy and the secondary wave energy is stronger. This facilitates the separation of the primary wave energy and the secondary wave energy, and allows the secondary wave energy to be attenuated to a greater extent.

[0159] It should be noted that, Figure 4 The dashed lines in the diagram are artificially defined cut lines in the conventional three-dimensional Lardon transform domain. In processing massive amounts of seismic data, this consumes a significant amount of manpower and has low computational efficiency. Meanwhile... Figure 5 The dashed lines in the diagram are cut-off lines defined based on a three-dimensional cone filter. This allows for the automatic division of primary and secondary wave energies, avoiding the need for manual definition of the cut-off area and making it more convenient for industrial production applications.

[0160] Step 106: The computer device filters the three-dimensional Radon transformed data based on the three-dimensional spheroid filter to obtain three-dimensional Radon multiple wave data.

[0161] wherein the three-dimensional Radon transformed data is M(λ x ,λ y ,f); the computer device inputs the three-dimensional Radon transformed data into the three-dimensional spheroid filter, and for each seismic trace in the three-dimensional Radon transformed data, determines the mapping radius of the seismic trace in the three-dimensional spheroid filter by the formula If the mapping radius of the seismic trace is located between the first frequency maximum value and the first frequency minimum value, the seismic trace is mapped to the inside of the three-dimensional spheroid filter, i.e. the multiple wave reservation zone. If the mapping radius is greater than the first frequency maximum value or less than the first frequency minimum value, the seismic trace is mapped to the outside of the unit spheroid filter, i.e. the primary wave frequency cut-off zone. In this way, the primary wave data in the three-dimensional Radon transformed data is mapped to the primary wave frequency cut-off zone, and the multiple wave data is mapped to the multiple wave reservation zone, so as to separate the primary wave data and the multiple wave data, and then the multiple wave data in the multiple wave reservation zone is taken as the three-dimensional Radon multiple wave data, and the expression of the three-dimensional Radon multiple wave data is D m (λ x ,λ y ,f).

[0162] Step 107: The computer device determines target seismic data after suppressing multiple waves based on the difference between the first seismic data of the plurality of first seismic trace sets and the three-dimensional Radon multiple wave data of the plurality of first seismic trace sets.

[0163] This step can be implemented through the following steps (1)-(3):

[0164] (1) The computer device inversely processes the plurality of three-dimensional Radon multiple wave data to obtain a plurality of moveout multiple wave data.

[0165] This step can be implemented through the following steps A1-A3:

[0166] A1: The computer device inversely processes the three-dimensional Radon multiple wave data to obtain to-be-transformed multiple wave data by performing three-dimensional Radon inverse transformation on the three-dimensional Radon multiple wave data for each three-dimensional Radon multiple wave data.

[0167] It should be noted that the computer device inversely processes the three-dimensional Radon multiple wave data to obtain to-be-transformed multiple wave data in the frequency space domain. The expression of the to-be-transformed multiple wave data is D m (x,y,f).

[0168] A2: The computer device inversely processes the to-be-transformed multiple wave data to obtain regularized multiple wave data by performing one-dimensional Fourier inverse transformation on the to-be-transformed multiple wave data.

[0169] It should be noted that the one-dimensional inverse Fourier transform is a transformation along the time direction, that is, the frequency domain to-be-transformed multiple wave data is transformed into time domain regularized multiple wave data, and the regularized multiple wave data is represented as d m0 (x,y,t).

[0170] A3: The computer device performs inverse regularization processing and inverse time difference correction processing on the regularized multiple wave data to obtain the moveout-corrected multiple wave data.

[0171] Wherein, the computer device searches, based on the second seismic data d0(x,y,t) of the second seismic trace set, for each seismic trace in the second seismic data, the closest node in the regular grid of the regularized multiple wave data to the seismic trace, and restores the coordinates of the node to the seismic trace.

[0172] It should be noted that the time of the node needs to be subjected to inverse time difference correction processing to obtain the sampling time corresponding to the second seismic data every time a node is restored.

[0173] Referring to Figure 6 , Figure 6 The moveout-corrected multiple wave data obtained by the method provided in the embodiments of the present application. Figure 7 It can be seen that the multiple wave data obtained by the conventional three-dimensional Radon transform contains the shadow of multiple wave residues, indicating that the multiple wave data is not completely separated; and Figure 6 and Figure 7 The moveout-corrected multiple wave data in the present application does not contain the shadow of residues, indicating that the moveout-corrected multiple wave data is completely separated. Figure 7

[0174] (2) The computer device obtains a plurality of second seismic data corresponding to a plurality of first seismic data.

[0175] It should be noted that the plurality of second seismic data in this step is the second seismic data of the plurality of second seismic trace sets in step 101; each first seismic data corresponds to one second seismic data.

[0176] (3) The computer device determines the difference between the second seismic data and the moveout-corrected multiple wave data to obtain target seismic data.

[0177] Wherein, the computer device uses an adaptive subtraction method to subtract a plurality of moveout-corrected multiple wave data in a plurality of second seismic data, and the difference obtained is the target seismic data for suppressing multiple waves.

[0178] Referring to Figure 8 , Figure 8 ​For suppressing the multiple wave of the gather after the conventional 3D Radon transform, it can be seen from the figure that due to the influence of the spatial truncation effect, there is a residual of part of the multiple wave energy, and the multiple wave separation is not complete. See Figure 9 , Figure 9 For suppressing the multiple wave of the gather after the method provided by the embodiment of the present application, because the energy group is more convergent, the separability of the primary wave energy and the multiple wave energy is stronger, and the residual multiple wave energy is less, thereby making the multiple wave be attenuated to a greater extent.

[0179] Take the second seismic gather of Figure 2 as an example, test the calculation running time of the conventional 3D Radon transform and the 3D Radon transform provided by the embodiment of the present application, and the size of the original seismic data is 512 MB; under the same test platform environment, the conventional 3D Radon transform takes 782 seconds, and the 3D Radon transform provided by the embodiment of the present application takes 94 seconds, which is about 8 times faster than the conventional 3D Radon transform, thereby improving the calculation efficiency of the 3D Radon transform.

[0180] See Figure 10 , Figure 10 The second seismic gather is for a certain marine area, and because the sea surface is a strong reflection interface, a strong energy first-order and second-order free surface multiple wave is formed. See Figure 11 , Figure 11 The energy distribution of the first seismic data transformed to the 3D Radon transform domain by the method provided by the embodiment of the present application can be seen from the figure, and the primary wave energy is vertically distributed near the region of λ x = 0, λ y = 0, and the multiple wave energy is distributed in the non-zero region with different slopes. See Figure 12 , Figure 12 After the 3D Radon forward transform data is filtered by the 3D cone filter and after the inverse processing, the gather of the moveout multiple wave data obtained can be seen, and the strong energy first-order and second-order surface multiple wave is effectively separated. See Figure 13 , Figure 13 The second seismic data and the moveout multiple wave data corresponding to Figure 10 and Figure 12 respectively are processed by the adaptive subtraction method to obtain the gather of the target seismic data after suppressing the multiple wave; it should be noted that Figure 10 and Figure 12 The gathers are preprocessed in amplitude, phase and moveout. From Figure 13 , it can be seen that the multiple wave in the gather of the target seismic data obtained by the method provided by the embodiment of the present application is effectively suppressed.

[0181] The embodiment of the present application provides a multiple wave suppression method. Since the three-dimensional spheroid filter is determined based on the maximum curvature, the minimum curvature and the maximum frequency, the process of artificially defining the three-dimensional spheroid filter is avoided, and the efficiency of generating the three-dimensional spheroid filter is improved. Then, the three-dimensional Rayleigh multiple wave data is obtained through the three-dimensional spheroid filter, and the efficiency of obtaining the three-dimensional Rayleigh multiple wave data is improved. Then, the multiple wave is suppressed based on the difference between the first seismic data and the three-dimensional Rayleigh multiple wave data, and the efficiency of suppressing the multiple wave is improved.

[0182] The embodiment of the present application provides a multiple wave suppression device. Referring to Figure 14 , the device comprises:

[0183] The first determination module 1401 is configured to determine first seismic data of a plurality of first seismic gathers. For each first seismic gather, the first seismic gather comprises a plurality of seismic traces, and the first seismic data of the first seismic gather comprises coordinates of each seismic trace at each sampling time.

[0184] The second determination module 1402 is configured to determine the curvature and the frequency of each seismic trace in the first seismic gather.

[0185] The third determination module 1403 is configured to determine the maximum curvature and the minimum curvature in the curvatures of the plurality of seismic traces, and determine the maximum frequency in the frequencies of the plurality of seismic traces.

[0186] The generation module 1404 is configured to generate a three-dimensional spheroid filter based on the maximum curvature, the minimum curvature and the maximum frequency.

[0187] The transformation module 1405 is configured to perform three-dimensional Rayleigh forward transformation on the first seismic data of the first seismic gather to obtain three-dimensional Rayleigh forward transformation data.

[0188] The filtering module 1406 is configured to perform filtering processing on the three-dimensional Rayleigh forward transformation data based on the three-dimensional spheroid filter to obtain three-dimensional Rayleigh multiple wave data.

[0189] The fourth determination module 1407 is configured to determine target seismic data after suppressing the multiple wave based on the difference between the first seismic data of the plurality of first seismic traces and the three-dimensional Rayleigh multiple wave data of the plurality of first seismic gathers.

[0190] In a possible implementation manner, the generation module 1404 comprises:

[0191] The first determination unit is configured to determine the product of the maximum curvature and the maximum frequency to obtain a first frequency maximum value of the three-dimensional spheroid filter.

[0192] The second determination unit is configured to determine the product of the minimum curvature and the maximum frequency to obtain a first frequency minimum value of the three-dimensional spheroid filter.

[0193] The generating unit is configured to generate a three-dimensional pyramid filter based on the first frequency maximum value and the first frequency minimum value, and a multiple wave frequency reservation area of the three-dimensional pyramid filter is an interval composed of the first frequency maximum value and the first frequency minimum value.

[0194] In a possible implementation, the transforming module 1405 includes:

[0195] The first transforming unit is configured to perform one-dimensional Fourier forward transform on the first seismic data to obtain the to-be-transformed seismic data.

[0196] The third determining unit is configured to determine a three-dimensional Radon transform operator of the first seismic trace set based on the curvature, the maximum frequency of each seismic trace, and the to-be-transformed seismic data.

[0197] The fourth determining unit is configured to determine a generalized inverse operator matched with the three-dimensional Radon transform operator.

[0198] The second transforming unit is configured to perform three-dimensional Radon forward transform on the to-be-transformed seismic data based on the generalized inverse operator to obtain three-dimensional Radon forward transform data.

[0199] In a possible implementation, the third determining unit includes:

[0200] The first determining subunit is configured to determine a three-dimensional Radon transform parameter of each seismic trace based on the curvature and the maximum frequency of each seismic trace.

[0201] The second determining subunit is configured to determine the three-dimensional Radon transform operator based on the three-dimensional Radon transform parameter of each seismic trace and the to-be-transformed seismic data through first relationship data, the first relationship data being relationship data among the three-dimensional Radon transform operator, the three-dimensional Radon transform parameter, and the to-be-transformed seismic data.

[0202] In a possible implementation, the fourth determining module 1407 includes:

[0203] The first processing unit is configured to perform inverse processing on the plurality of three-dimensional Radon multiple wave data to obtain a plurality of moveout multiple wave data.

[0204] The first obtaining unit is configured to obtain a plurality of second seismic data corresponding to the plurality of first seismic data.

[0205] The fifth determining unit is configured to determine a difference between the plurality of second seismic data and the plurality of moveout multiple wave data to obtain target seismic data.

[0206] In a possible implementation, the processing unit includes:

[0207] The third transform subunit is configured to perform three-dimensional Radon inverse transform on the three-dimensional Radon multiple wave data to obtain the to-be-transformed multiple wave data.

[0208] The fourth transform subunit is configured to perform one-dimensional Fourier inverse transform on the to-be-transformed multiple wave data to obtain the regularized multiple wave data.

[0209] The processing subunit is configured to perform de-regularization processing and inverse moveout correction processing on the regularized multiple wave data to obtain the move-corrected multiple wave data.

[0210] In a possible implementation, the first determination module 1401 comprises:

[0211] The second acquisition unit is configured to acquire the original seismic data.

[0212] The second processing unit is configured to perform move correction preprocessing on the original seismic data to obtain second seismic data of a plurality of second seismic gathers.

[0213] The third processing unit is configured to perform grid regularization processing and moveout correction processing on the second seismic data of each second seismic gather to obtain first seismic data of a first seismic gather.

[0214] Figure 15 A structural block diagram of a computer device 1500 provided by an example embodiment of the present application is shown. The computer device 1500 can be a portable mobile computer device, such as a smart phone, a tablet computer, an MP3 (Moving Picture Experts Group Audio Layer III) player, an MP4 (Moving Picture Experts Group Audio Layer IV) player, a notebook computer, or a desktop computer. The computer device 1500 can also be referred to as a user device, a portable computer device, a laptop computer device, a desktop computer device, and other names.

[0215] Generally, the computer device 1500 comprises a processor 1501 and a memory 1502.

[0216] The processor 1501 can include one or more processing cores, such as a 4-core processor, an 8-core processor, and the like. The processor 1501 can be implemented in at least one of a hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), a PLA (Programmable Logic Array). The processor 1501 can also include a main processor and a coprocessor, the main processor being a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit), and the coprocessor being a low-power processor for processing data in a standby state. In some embodiments, the processor 1501 can be integrated with a GPU (Graphics Processing Unit) for rendering and drawing content required to be displayed by the display screen. In some embodiments, the processor 1501 can further include an AI (Artificial Intelligence) processor for processing computing operations related to machine learning.

[0217] The memory 1502 can include one or more computer-readable storage media that can be non-transitory. The memory 1502 can also include high-speed random access memory and nonvolatile memory such as one or more disk storage devices, flash storage devices. In some embodiments, the non-transitory computer-readable storage medium in the memory 1502 is used to store at least one instruction for being executed by the processor 1501 to implement the multiple wave suppression method provided by the method embodiment in the present application.

[0218] In some embodiments, the computer device 1500 can also optionally include a peripheral device interface 1503 and at least one peripheral device. The processor 1501, the memory 1502, and the peripheral device interface 1503 can be connected through a bus or a signal line. Each peripheral device can be connected to the peripheral device interface 1503 through a bus, a signal line, or a circuit board. Specifically, the peripheral device includes at least one of a radio frequency circuit 1504, a display screen 1505, a camera assembly 1506, an audio circuit 1507, a positioning assembly 1508, and a power supply 1509.

[0219] The peripheral interface 1503 can be used to connect at least one I / O (Input / Output) related peripheral device to the processor 1501 and the memory 1502. In some embodiments, the processor 1501, the memory 1502 and the peripheral interface 1503 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 1501, the memory 1502 and the peripheral interface 1503 can be implemented on a separate chip or circuit board, and the present embodiments are not limited in this regard.

[0220] The radio frequency circuit 1504 is used to receive and send RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuit 1504 communicates with a communication network and other communication devices through electromagnetic signals. The radio frequency circuit 1504 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals into electrical signals. Optionally, the radio frequency circuit 1504 includes an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and the like. The radio frequency circuit 1504 can communicate with other computer devices through at least one wireless communication protocol. The wireless communication protocol includes but is not limited to the World Wide Web, a metropolitan area network, an intranet, various generations of mobile communication networks (2G, 3G, 4G and 5G), a wireless local area network and / or a WiFi (Wireless Fidelity) network. In some embodiments, the radio frequency circuit 1504 can also include NFC (Near Field Communication) related circuitry, and the present application is not limited in this regard.

[0221] The display screen 1505 is configured to display a UI (User Interface). The UI can include graphics, text, icons, video, and any combination thereof. When the display screen 1505 is a touch display screen, the display screen 1505 is further configured to capture touch signals on or above the surface of the display screen 1505. The touch signals can be input to the processor 1501 as control signals for processing. In this case, the display screen 1505 can also be configured to provide virtual buttons and / or virtual keyboard, also known as soft buttons and / or soft keyboard. In some embodiments, the display screen 1505 can be one, disposed on the front panel of the computer device 1500; in other embodiments, the display screen 1505 can be at least two, respectively disposed on different surfaces of the computer device 1500 or in a folding design; in other embodiments, the display screen 1505 can be a flexible display screen, disposed on a curved surface or a folding surface of the computer device 1500. Even, the display screen 1505 can also be disposed in an irregular shape other than a rectangle, i.e., a special-shaped screen. The display screen 1505 can be made of materials such as LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), etc.

[0222] The camera assembly 1506 is configured to capture images or videos. Optionally, the camera assembly 1506 includes a front-facing camera and a rear-facing camera. Typically, the front-facing camera is disposed on the front panel of the computer device, and the rear-facing camera is disposed on the back of the computer device. In some embodiments, the rear-facing camera is at least two, which are any one of a main camera, a depth-of-field camera, a wide-angle camera, and a telephoto camera, to realize the background blur function by fusing the main camera and the depth-of-field camera, the panoramic shooting and VR (Virtual Reality) shooting function by fusing the main camera and the wide-angle camera, or other fusion shooting functions. In some embodiments, the camera assembly 1506 can further include a flash. The flash can be a single-color-temperature flash or a dual-color-temperature flash. The dual-color-temperature flash refers to a combination of a warm light flash and a cold light flash, which can be used for light compensation under different color temperatures.

[0223] The audio circuit 1507 can include a microphone and a speaker. The microphone is used to collect sound waves of a user and an environment, and convert the sound waves into an electrical signal input to the processor 1501 for processing, or input to the radio frequency circuit 1504 to realize voice communication. For the purpose of stereo sound collection or noise reduction, the microphone can be multiple, respectively arranged at different parts of the computer device 1500. The microphone can also be an array microphone or an omnidirectional collection type microphone. The speaker is used to convert the electrical signal from the processor 1501 or the radio frequency circuit 1504 into sound waves. The speaker can be a traditional diaphragm speaker, or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only can the electrical signal be converted into a sound wave that a human can hear, but also can be converted into a sound wave that a human cannot hear for ranging purposes. In some embodiments, the audio circuit 1507 can also include a headphone jack.

[0224] The positioning component 1508 is used to position the current geographic location of the computer device 1500 to realize navigation or LBS (Location Based Service). The positioning component 1508 can be a positioning component based on the GPS (Global Positioning System) of the United States, the Beidou system of China or the Galileo system of Russia.

[0225] The power supply 1509 is used to supply power to various components in the computer device 1500. The power supply 1509 can be alternating current, direct current, disposable battery or rechargeable battery. When the power supply 1509 includes a rechargeable battery, the rechargeable battery can be a wired charging battery or a wireless charging battery. The wired charging battery is a battery charged through a wired line, and the wireless charging battery is a battery charged through a wireless coil. The rechargeable battery can also be used to support fast charging technology.

[0226] In some embodiments, the computer device 1500 further includes one or more sensors 1510. The one or more sensors 1510 include but are not limited to: an acceleration sensor 1511, a gyroscope sensor 1512, a pressure sensor 1513, a fingerprint sensor 1514, an optical sensor 1515 and a proximity sensor 1516.

[0227] The acceleration sensor 1511 can detect the acceleration magnitude in three coordinate axes of the coordinate system established by the computer device 1500. For example, the acceleration sensor 1511 can be used to detect the components of the gravitational acceleration in three coordinate axes. The processor 1501 can control the display screen 1505 to display the user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 1511. The acceleration sensor 1511 can also be used for gaming or user motion data collection.

[0228] The gyroscope sensor 1512 can detect the body direction and rotation angle of the computer device 1500, and the gyroscope sensor 1512 can cooperate with the acceleration sensor 1511 to collect the 3D action of the user on the computer device 1500. According to the data collected by the gyroscope sensor 1512, the processor 1501 can realize the following functions: motion sensing (such as changing the UI according to the user's tilt operation), image stabilization when shooting, game control, and inertial navigation.

[0229] The pressure sensor 1513 can be arranged on the side frame of the computer device 1500 and / or the lower layer of the display screen 1505. When the pressure sensor 1513 is arranged on the side frame of the computer device 1500, the user's holding signal on the computer device 1500 can be detected, and the left and right hand recognition or shortcut operation can be performed by the processor 1501 according to the holding signal collected by the pressure sensor 1513. When the pressure sensor 1513 is arranged on the lower layer of the display screen 1505, the controllable control on the UI interface can be controlled by the processor 1501 according to the pressure operation of the user on the display screen 1505. The controllable control includes at least one of a button control, a scroll bar control, an icon control, and a menu control.

[0230] The fingerprint sensor 1514 is used to collect the fingerprint of the user, and the identity of the user can be recognized by the processor 1501 according to the fingerprint collected by the fingerprint sensor 1514, or by the fingerprint sensor 1514 according to the collected fingerprint. When the identity of the user is recognized as a trusted identity, the user is authorized by the processor 1501 to perform related sensitive operations, including unlocking the screen, viewing encrypted information, downloading software, payment, and changing settings. The fingerprint sensor 1514 can be arranged on the front, back or side of the computer device 1500. When the computer device 1500 is provided with a physical button or a manufacturer's logo, the fingerprint sensor 1514 can be integrated with the physical button or the manufacturer's logo.

[0231] The optical sensor 1515 is used to collect the ambient light intensity. In one embodiment, the processor 1501 can control the display brightness of the display screen 1505 according to the ambient light intensity collected by the optical sensor 1515. Specifically, when the ambient light intensity is high, the display brightness of the display screen 1505 is increased; when the ambient light intensity is low, the display brightness of the display screen 1505 is decreased. In another embodiment, the processor 1501 can also dynamically adjust the shooting parameters of the camera assembly 1506 according to the ambient light intensity collected by the optical sensor 1515.

[0232] The proximity sensor 1516, also referred to as a distance sensor, is usually arranged on the front panel of the computer device 1500. The proximity sensor 1516 is used to collect the distance between the user and the front face of the computer device 1500. In one embodiment, when the proximity sensor 1516 detects that the distance between the user and the front face of the computer device 1500 gradually decreases, the display screen 1505 is switched from the bright screen state to the screen-off state under the control of the processor 1501; when the proximity sensor 1516 detects that the distance between the user and the front face of the computer device 1500 gradually increases, the display screen 1505 is switched from the screen-off state to the bright screen state under the control of the processor 1501.

[0233] Those skilled in the art can understand that, Figure 15 The structure shown in the figure does not constitute a limitation on the computer device 1500, and can include more or fewer components than shown, or combine certain components, or use different component arrangements.

[0234] In another aspect, a computer program product or computer program is provided, which includes computer program code stored in a computer readable storage medium. The computer program code is read by the processor of the computer device from the computer readable storage medium, and the processor executes the computer program code, so that the computer device performs the operations performed by the above-mentioned determination method for identifying water-invasion layers.

[0235] In some embodiments, the computer program related to the embodiments of the present application can be deployed to execute on one computer device, or on multiple computer devices located in one place, or on multiple computer devices distributed in multiple places and interconnected through a communication network, which can constitute a blockchain system.

[0236] The embodiments of the present application provide a multiple wave suppression method. Since the three-dimensional spines filter is determined based on the maximum curvature, the minimum curvature and the maximum frequency, the process of artificially defining the three-dimensional spines filter is avoided, and the efficiency of generating the three-dimensional spines filter is improved. Then, the three-dimensional Rayleigh multiple wave data is obtained through the three-dimensional spines filter, and the efficiency of obtaining the three-dimensional Rayleigh multiple wave data is improved. Further, the multiple wave is suppressed based on the difference between the first seismic data and the three-dimensional Rayleigh multiple wave data, and the efficiency of suppressing the multiple wave is improved.

[0237] The above only describes optional embodiments of the present application and does not limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method of suppressing multiple waves, characterized by, The method comprises: determining first seismic data of a plurality of first seismic gathers, for each first seismic gather, the first seismic gather comprising a plurality of seismic traces, the first seismic data of the first seismic gather comprising coordinates of each seismic trace at each sampling time; determining a curvature and a frequency of each seismic trace in the first seismic gather; determining a maximum curvature and a minimum curvature in the curvatures of the plurality of seismic traces, and determining a maximum frequency in the frequencies of the plurality of seismic traces; generating a three-dimensional spheroid filter based on the maximum curvature, the minimum curvature and the maximum frequency; performing three-dimensional Radon forward transformation on the first seismic data of the first seismic gather to obtain three-dimensional Radon forward transformation data; filtering the three-dimensional Radon forward transformation data based on the three-dimensional spheroid filter to obtain three-dimensional Radon multiple wave data; determining target seismic data after suppressing multiple waves based on a difference between the first seismic data of the plurality of first seismic gathers and the three-dimensional Radon multiple wave data of the plurality of first seismic gathers; wherein the generating of the three-dimensional spheroid filter based on the maximum curvature, the minimum curvature and the maximum frequency comprises: determining a product of the maximum curvature and the maximum frequency to obtain a first frequency maximum value of the three-dimensional spheroid filter; determining a product of the minimum curvature and the maximum frequency to obtain a first frequency minimum value of the three-dimensional spheroid filter; and generating the three-dimensional spheroid filter based on the first frequency maximum value and the first frequency minimum value, a multiple wave frequency retention area of the three-dimensional spheroid filter being an interval composed of the first frequency maximum value and the first frequency minimum value.

2. The method of suppressing multiple waves according to claim 1, wherein, The performing of the three-dimensional Radon forward transformation on the first seismic data of the first seismic gather to obtain three-dimensional Radon forward transformation data comprises: performing one-dimensional Fourier forward transformation on the first seismic data to obtain to-be-transformed seismic data; determining a three-dimensional Radon transformation operator of the first seismic gather based on the curvature of each seismic trace, the maximum frequency and the to-be-transformed seismic data; determining a generalized inverse operator matched with the three-dimensional Radon transformation operator; performing three-dimensional Radon forward transformation on the to-be-transformed seismic data based on the generalized inverse operator to obtain the three-dimensional Radon forward transformation data.

3. The method of suppressing multiple waves according to claim 2, wherein, The determining of the three-dimensional Radon transformation operator of the first seismic gather based on the curvature of each seismic trace, the maximum frequency and the to-be-transformed seismic data comprises: determining a three-dimensional Radon transformation parameter of each seismic trace based on the curvature of each seismic trace and the maximum frequency; determining the three-dimensional Radon transformation operator based on the three-dimensional Radon transformation parameter of each seismic trace and the to-be-transformed seismic data through first relationship data, the first relationship data being relationship data among the three-dimensional Radon transformation operator, the three-dimensional Radon transformation parameter and the to-be-transformed seismic data.

4. The method of suppressing multiple waves according to claim 1, wherein, The determining of target seismic data after suppressing multiple waves based on a difference between the first seismic data of the plurality of first seismic gathers and the three-dimensional Radon multiple wave data of the plurality of first seismic gathers comprises: performing inverse processing on a plurality of three-dimensional Radon multiple wave data to obtain a plurality of moveout multiple wave data; Obtaining a plurality of second seismic data corresponding to a plurality of first seismic data; Determining a difference between the plurality of second seismic data and the plurality of moveout multiple wave data to obtain the target seismic data.

5. The method of suppressing multiple waves according to claim 4, wherein, The anti-processing of the plurality of three-dimensional Radon multiple wave data to obtain the plurality of moveout multiple wave data comprises: For each three-dimensional Radon multiple wave data, performing three-dimensional Radon inverse transformation on the three-dimensional Radon multiple wave data to obtain to-be-transformed multiple wave data; Performing one-dimensional Fourier inverse transformation on the to-be-transformed multiple wave data to obtain regularized multiple wave data; Performing anti-regularization processing and anti-time difference correction processing on the regularized multiple wave data to obtain the moveout multiple wave data.

6. The method of suppressing multiple waves according to claim 1, wherein, The first seismic data of the plurality of first seismic traces comprises: Obtaining original seismic data; Performing moveout correction preprocessing on the original seismic data to obtain second seismic data of a plurality of second seismic traces; For each second seismic trace, performing grid regularization processing and time difference correction processing on the second seismic data of the second seismic trace to obtain the first seismic data of the first seismic trace.

7. A device for suppressing multiple waves, characterized by comprising: The device comprises: A first determination module configured to determine first seismic data of a plurality of first seismic traces, wherein each first seismic trace comprises a plurality of seismic traces, and the first seismic data of each first seismic trace comprises coordinates of each seismic trace at each sampling time; A second determination module configured to determine curvature and frequency of each seismic trace in the first seismic trace; A third determination module configured to determine maximum curvature and minimum curvature in the curvature of the plurality of seismic traces, and determine maximum frequency in the frequency of the plurality of seismic traces; A generation module configured to generate a three-dimensional cone filter based on the maximum curvature, the minimum curvature and the maximum frequency; A transformation module configured to perform three-dimensional Radon forward transformation on the first seismic data of the first seismic trace to obtain three-dimensional Radon forward transformation data; A filtering module configured to perform filtering processing on the three-dimensional Radon forward transformation data based on the three-dimensional cone filter to obtain three-dimensional Radon multiple wave data; A fourth determination module configured to determine target seismic data after multiple wave suppression based on a difference between the first seismic data of the plurality of first seismic traces and the three-dimensional Radon multiple wave data of the plurality of first seismic traces; The generation module is configured to: determine a product of the maximum curvature and the maximum frequency to obtain a first frequency maximum value of the three-dimensional cone filter; determine a product of the minimum curvature and the maximum frequency to obtain a first frequency minimum value of the three-dimensional cone filter; and generate the three-dimensional cone filter based on the first frequency maximum value and the first frequency minimum value, wherein a multiple wave frequency retention area of the three-dimensional cone filter is an interval composed of the first frequency maximum value and the first frequency minimum value.

8. A computer device, comprising: The computer device comprises one or more processors and one or more memories, and the one or more memories store at least one instruction, which is loaded and executed by the one or more processors to implement the operation performed by the method for suppressing multiple waves according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The storage medium stores at least one instruction, which is loaded and executed by the processor to implement the operation performed by the method for suppressing multiple waves according to any one of claims 1 to 6.

Images