CN112255684A - Processing method and device for collecting footprint of seismic data - Google Patents
Processing method and device for collecting footprint of seismic data Download PDFInfo
- Publication number
- CN112255684A CN112255684A CN202011155133.8A CN202011155133A CN112255684A CN 112255684 A CN112255684 A CN 112255684A CN 202011155133 A CN202011155133 A CN 202011155133A CN 112255684 A CN112255684 A CN 112255684A
- Authority
- CN
- China
- Prior art keywords
- data
- post
- seismic data
- stack seismic
- slice
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 76
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000007781 pre-processing Methods 0.000 claims abstract description 24
- 238000003825 pressing Methods 0.000 claims abstract description 21
- 238000005070 sampling Methods 0.000 claims description 76
- 238000004590 computer program Methods 0.000 claims description 18
- 238000002156 mixing Methods 0.000 claims description 17
- 230000010354 integration Effects 0.000 claims description 16
- 238000011084 recovery Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 10
- 230000001629 suppression Effects 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 24
- 230000006870 function Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005056 compaction Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000005520 cutting process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/514—Post-stack
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a processing method and a device for collecting a footprint of seismic data, wherein the method comprises the following steps: pre-stack preprocessing is carried out on the seismic data to obtain post-stack seismic data; performing filling processing on the post-stack seismic data to obtain filled post-stack seismic data; rotating the filled post-stack seismic data to obtain rotated post-stack seismic data; performing acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data; carrying out reverse rotation on the pressed post-stack seismic data to obtain the reverse-rotation post-stack seismic data; and recovering the reversely rotated post-stack seismic data to obtain the recovered post-stack seismic data. The device is used for executing the method. The processing method and device for the acquisition footprint of the seismic data, provided by the embodiment of the invention, improve the signal-to-noise ratio and resolution of the seismic data.
Description
Technical Field
The invention relates to the technical field of seismic exploration, in particular to a processing method and device for a collection footprint of seismic data.
Background
Currently, oil and gas exploration is developing in both breadth and depth. The former aims at discovering new prospect field oil and gas fields; the latter requires to find oil and gas reservoirs with large buried depth and high complexity, and simultaneously requires to solve the problems of fine structure, fine description of reservoir parameters and the like. With the continuous deepening of exploration of lithologic oil and gas reservoirs and hidden oil and gas reservoirs, the geological target also turns to lithologic oil and gas finding from the original structure, and the problem of fine description of oil and gas reservoirs and reservoir parameters is well solved, so that the problem of the signal-to-noise ratio of imaging data is involved.
The acquisition footprint noise is also called acquisition trace noise, is an artificial trace generated by human factors, is an artificial trace left in the process of acquiring and processing seismic data, and is represented by the phenomenon that regular amplitude change false images appear on seismic section depth or time slices, and the energy of part of the acquired footprint noise is very strong, so that the prediction precision and efficiency of an oil and gas reservoir are often seriously influenced. For the compression of the collected footprints, the collected footprints can be partially attenuated by combining and interpolating before stacking; the post-stack acquisition footprint pressing method is more, and comprises an inclination filtering method, an F-Kx-Ky filtering method, a wavelet transformation method, a self-adaptive method and the like. The use of weighting functions in DMO or prestack migration processes can also partially attenuate the acquisition footprint. However, different acquisition footprints have different generation mechanisms, and different acquisition footprints have self-applicability and self-limitation, so that the development of targeted acquisition footprints for pressing has important significance.
Disclosure of Invention
The embodiment of the invention provides a processing method and device for seismic data acquisition footprints, and aims to solve the problem of how to improve the signal-to-noise ratio and the resolution of seismic data for randomly distributed or locally and relatively distributed acquisition footprints.
In one aspect, the present invention provides a processing method for collecting footprints of seismic data, including:
pre-stack preprocessing is carried out on the seismic data to obtain post-stack seismic data;
performing filling processing on the post-stack seismic data to obtain filled post-stack seismic data;
rotating the filled post-stack seismic data to obtain rotated post-stack seismic data;
performing acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data;
carrying out reverse rotation on the pressed post-stack seismic data to obtain the reverse-rotation post-stack seismic data;
and recovering the reversely rotated post-stack seismic data to obtain the recovered post-stack seismic data.
In another aspect, the present invention provides a processing apparatus for acquiring footprints of seismic data, comprising:
the preprocessing module is used for preprocessing the seismic data before the stack to obtain post-stack seismic data;
the filling-up processing module is used for filling up the post-stack seismic data to obtain the post-stack seismic data after filling up;
the rotation module is used for rotating the filled post-stack seismic data to obtain the rotated post-stack seismic data;
the pressing module is used for carrying out acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data;
the anti-rotation module is used for carrying out anti-rotation on the pressed post-stack seismic data to obtain the anti-rotation post-stack seismic data;
and the recovery processing module is used for performing recovery processing on the reversely rotated post-stack seismic data to obtain the recovered post-stack seismic data.
In another aspect, the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor implements the steps of the processing method for acquiring the footprint of the seismic data according to any one of the embodiments described above when executing the computer program.
In yet another aspect, the present invention provides a computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the steps of the method for processing the acquisition footprint of seismic data according to any one of the above embodiments.
The processing method and device for the acquisition footprint of the seismic data, provided by the embodiment of the invention, are used for pre-stack preprocessing the seismic data to obtain post-stack seismic data, performing filling-up processing on the post-stack seismic data to obtain filled-up post-stack seismic data, rotating the filled-up post-stack seismic data to obtain rotated post-stack seismic data, performing acquisition footprint pressing on each slice of the rotated post-stack seismic data to obtain pressed post-stack seismic data, performing reverse rotation on the pressed post-stack seismic data to obtain reverse-rotated post-stack seismic data, and performing recovery processing on the reverse-rotated post-stack seismic data to obtain recovered post-stack seismic data, so that the waveform characteristics of the original seismic data can be maintained, and the signal-to-noise ratio and the resolution of the seismic data are improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
fig. 1 is a schematic flow chart of a processing method for acquiring a footprint of seismic data according to an embodiment of the present invention.
Fig. 2 is a schematic flow chart of a processing method for acquiring footprints of seismic data according to another embodiment of the present invention.
Fig. 3 is a schematic flow chart of a processing method for acquiring footprints of seismic data according to another embodiment of the present invention.
Fig. 4 is a flowchart illustrating a processing method for acquiring footprints of seismic data according to still another embodiment of the present invention.
Fig. 5 is a schematic diagram of a part of data of the 1 st time slice data according to an embodiment of the present invention.
Fig. 6 is a schematic diagram of an integration result in the InLine direction of partial data of the 1 st time slice data according to an embodiment of the present invention.
Fig. 7 is a schematic diagram of integrated time slice data corresponding to partial data of the 1 st time slice data according to an embodiment of the present invention.
Fig. 8 is a schematic diagram of a part of data of new 1 st time slice data provided by an embodiment of the present invention.
FIG. 9 is a schematic diagram of the acquisition of post-stack seismic data prior to footprint suppression provided by an embodiment of the present invention.
FIG. 10 is a schematic diagram of denoised post-stack seismic data according to an embodiment of the invention.
FIG. 11 is a schematic diagram of a frequency amplitude spectrum for acquiring post-stack seismic data prior to footprint compression according to an embodiment of the present invention.
FIG. 12 is a schematic diagram of a frequency-amplitude spectrum of denoised post-stack seismic data according to an embodiment of the invention.
FIG. 13 is a schematic illustration of a slice display for acquiring post-stack seismic data prior to footprint compaction, as provided by an embodiment of the invention;
FIG. 14 is a schematic illustration of a slice display of denoised post-stack seismic data according to an embodiment of the invention.
Fig. 15 is a schematic structural diagram of a processing device for acquiring a footprint of seismic data according to an embodiment of the present invention.
Fig. 16 is a schematic structural diagram of a processing device for acquiring a footprint of seismic data according to another embodiment of the present invention.
Fig. 17 is a schematic structural diagram of a processing apparatus for acquiring a footprint of seismic data according to still another embodiment of the present invention.
Fig. 18 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention. It should be noted that the embodiments and features of the embodiments in the present application may be arbitrarily combined with each other without conflict.
Fig. 1 is a schematic flow chart of a processing method for acquiring footprints of seismic data according to an embodiment of the present invention, and as shown in fig. 1, the processing method for acquiring footprints of seismic data according to an embodiment of the present invention includes:
s101, pre-stack preprocessing is carried out on seismic data to obtain post-stack seismic data;
specifically, seismic data are acquired in the field to obtain seismic data. And pre-stack preprocessing is carried out on the seismic data to obtain post-stack seismic data. The prestack preprocessing includes, but is not limited to, observation system setting, static correction, velocity analysis, energy compensation, noise attenuation, cutting and the like of seismic data, and is set according to actual needs, and the embodiment of the invention is not limited. In the pre-stack pretreatment process, new artificial traces and footprints are avoided from being introduced due to unreasonable pre-stack pretreatment according to a normal treatment flow. The seismic data may be Common midpoint gathers (CMPs). The execution main body of the processing method for the acquisition footprint of the seismic data provided by the embodiment of the invention comprises but is not limited to a computer.
S102, performing filling-up processing on the stacked seismic data to obtain filled-up stacked seismic data;
specifically, after the post-stack seismic data are obtained, filling up processing is performed on the post-stack seismic data to fill up seismic channel data of the post-stack seismic data missing in the directions of an InLine (InLine) and a CrossLine (CrossLine), so as to obtain the filled-up post-stack seismic data, and the filled-up post-stack seismic data form a regular matrix in the directions of the InLine and the CrossLine. Wherein missing seismic trace data may be filled with zero values.
For example, because the observation systems are not uniformly distributed, and the line numbers of the inlines and the crosslines are not consistent with the CMP number, the post-stack seismic data needs to be subjected to filling processing, so that the filled post-stack seismic data form a regular matrix in the whole work area, and the subsequent data processing is facilitated. InLine and CrossLine which lack seismic trace data are recorded while the seismic trace data are filled, so that data recovery can be performed later.
S103, rotating the filled post-stack seismic data to obtain rotated post-stack seismic data;
specifically, after obtaining the filled-in post-stack seismic data, the filled-in post-stack seismic data may be rotated to obtain rotated post-stack seismic data, where the rotated post-stack seismic data includes a plurality of slice data, and the slice data may be time slice data or depth slice data. The rotated post-stack seismic data can be stored in a temporary file, so that subsequent data reading is facilitated.
For example, the filled-in post-stack seismic data is represented as (InLine, CrossLine, Z), Z may be time or depth, the rotated post-stack seismic data obtained by rotating the filled-in post-stack seismic data is represented as (Z, InLine, CrossLine), the rotated post-stack seismic data includes a plurality of slice data, each slice data corresponds to one Z value, data reading is performed in the form of slice data, and the method is suitable for processing mass data and can improve data reading efficiency.
S104, performing acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data;
specifically, after the rotated post-stack seismic data is obtained, acquisition footprint pressing may be performed on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data. In an embodiment of the present invention, the slice data is time slice data or depth slice data.
S105, reversely rotating the pressed post-stack seismic data to obtain reversely rotated post-stack seismic data;
specifically, after obtaining the suppressed post-stack seismic data, the suppressed post-stack seismic data is reversely rotated, the direction of data rotation is opposite to the direction of rotation in step S103, and the reversely rotated post-stack seismic data can be obtained in the reverse process of step S103.
For example, the compressed post-stack seismic data is in the form of slices (Z, InLine, CrossLine), and the anti-rotation post-stack seismic data obtained after the anti-rotation is expressed as (InLine, CrossLine, Z), and becomes data of an array distribution.
S106, performing recovery processing on the reversely rotated post-stack seismic data to obtain recovered post-stack seismic data;
specifically, after the reversely rotated post-stack seismic data is obtained, the reversely rotated post-stack seismic data is subjected to recovery processing, and data corresponding to the seismic channel data complemented in step S102 is removed from the reversely rotated post-stack seismic data, so that recovered post-stack seismic data is obtained.
The processing method for the acquisition footprint of the seismic data, provided by the embodiment of the invention, comprises the steps of preprocessing the seismic data before stacking to obtain post-stack seismic data, filling up the post-stack seismic data to obtain filled-up post-stack seismic data, rotating the filled-up post-stack seismic data to obtain rotated post-stack seismic data, pressing the acquisition footprint of each slice of the rotated post-stack seismic data to obtain pressed post-stack seismic data, reversely rotating the pressed post-stack seismic data to obtain reversely rotated post-stack seismic data, and recovering the reversely rotated post-stack seismic data to obtain recovered post-stack seismic data.
Fig. 2 is a schematic flow chart of a processing method for acquiring footprints of seismic data according to another embodiment of the present invention, and as shown in fig. 2, on the basis of the foregoing embodiments, further, performing acquisition footprinting compaction on each slice data of the rotated post-stack seismic data, and obtaining compacted post-stack seismic data includes:
s1041, integrating the amplitude of the slice data along the direction of an inline line to obtain integrated slice data in the inline direction;
specifically, for each slice data of the rotational post-stack seismic data, the amplitude of the slice data may be integrated along the InLine direction, and the starting point of the integration is the minimum number of inlines, so that the integrated slice data in the InLine direction may be obtained.
S1042, integrating the amplitude of the integral slice data in the main measuring line direction along the direction of the cross measuring line to obtain integral slice data;
specifically, after obtaining the integral slice data in the InLine direction, the amplitude of the integral slice data in the InLine direction is integrated along the crossline direction, and integral slice data is obtained. The integrated slice data may be integrated time offset data or integrated depth offset data.
S1043, filling a slice sample value into the integral slice data in the main measuring line direction according to the set main measuring line filling length, and filling a slice sample value into the integral slice data in the tie measuring line direction according to the set tie measuring line filling length to obtain the filled integral slice data;
specifically, after obtaining the integral slice data, the integral slice data may be filled with slice sample values in a main line direction according to a set main line filling length, that is, a main line filling length of slice sample values along each InLine minimum CMP direction of the integral slice data, the filling length being a main line filling length, a main line filling length of slice sample values also along each InLine maximum CMP direction of the integral slice data, the filling length also being a main line filling length, and the integral slice data may be filled with slice sample values in a tie line direction according to a set tie line filling length, that is, a tie line filling length of slice sample values along each CrossLine minimum line direction of the integral slice data, the filling length being a tie line filling length, a tie line filling length also along each CrossLine maximum line direction of the slice, the filling length is also the cross-line filling length, and the filled integral slice data is obtained. The main measuring line filling length and the cross measuring line filling length are preset and are set according to actual needs, and the embodiment of the invention is not limited.
For example, for the integrated slice data, two spatial variation steps in the InLine direction and the CrossLine direction are set, where the step a in the InLine direction and the step b in the CrossLine direction are respectively set, a and b may be equal or unequal, the direction in which the distribution of the footprint energy is stronger is usually collected, the selected step should be longer, and the collected distribution of the footprint energy strength can be obtained by observing the integrated slice data. The main survey line filling length is set as c, c is half of the InLine direction step length, namely c is a/2, the crosssurvey line filling length is set as d, d is half of the CrossLine direction step length b, namely d is b/2. When the integrated slice data is filled, c slice sample values are filled along each InLine minimum CMP direction of the integrated slice data, the filling length is c, c slice sample values are also filled along each InLine maximum CMP direction of the slice, the filling length is c, d slice sample values are filled along each CrossLine minimum line direction of the slice, the filling length is d, d slice sample values are also filled along each CrossLine maximum line direction of the slice, and the filling length is d, so that the filled integrated slice data is obtained. Wherein the padded sample values may take the value of zero.
And S1044, carrying out diagonal amplitude processing on each sampling point data in the integral slicing data according to the filled integral slicing data to obtain integral sampling point amplitude difference value slicing data.
Specifically, after the filled integral slice data is obtained, diagonal amplitude processing may be performed on each sample point data in the integral slice data according to the filled integral slice data, that is, for each sample point data in the integral slice data, four corresponding rectangular corner point data may be obtained from the filled integral slice data, and the amplitudes of the four rectangular corner point data are subjected to data processing, so that integral sample point amplitude difference slice data may be obtained. And (4) performing the data processing processes of the steps S1041, S1042, S1043 and S1044 on each slice data of the rotational post-stack seismic data, so as to obtain the suppressed post-stack seismic data.
Fig. 3 is a schematic flow chart of a processing method for acquiring footprints of seismic data according to still another embodiment of the present invention, and as shown in fig. 3, on the basis of the foregoing embodiments, further, the performing diagonal amplitude processing on each sample point data in the integrated slice data according to the filled integrated slice data to obtain integrated sample point amplitude difference slice data includes:
s10441, obtaining a sampling point data from the integral slice data, and obtaining a rectangular corner point corresponding to the sampling point data from the filled integral slice data; the sampling point data comprises four rectangular corner points corresponding to the sampling point data to form a rectangle, and the sampling point data is positioned in the center of the rectangle;
specifically, one sampling point data is obtained from the integral slice data, the rectangular corner points corresponding to the sampling point data can be obtained from the filled integral slice data, the number of the rectangular corner points corresponding to the sampling point data is four, the four rectangular corner points are connected to form a rectangle, and the sampling point data is located in the center of the rectangle.
For example, any sample point data is obtained from the integrated slice data, and with this sample point data as the center, a rectangle with a side length a and b respectively can be found from the filled integrated slice data, and this sample point data is exactly the center of the rectangle.
S10442, calculating a difference value between the sum of the amplitude values of the integral sampling points corresponding to the first rectangular corner point and the opposite third rectangular corner point and the sum of the amplitude values of the integral sampling points corresponding to the second rectangular corner point and the opposite fourth rectangular corner point to obtain a middle value; the rectangular corner points corresponding to the sampling point data comprise the first rectangular corner point, the second rectangular corner point, the third rectangular corner point and the fourth rectangular corner point;
specifically, the rectangular corner points corresponding to the sampling point data include a first rectangular corner point, a second rectangular corner point, a third rectangular corner point and a fourth rectangular corner point, the first rectangular corner point is opposite to the third rectangular corner point and is two end points of a rectangular diagonal line, and the second rectangular corner point is opposite to the fourth rectangular corner point and is two end points of another diagonal line of the rectangle. And calculating the sum p of the integral sample point amplitude value of the first rectangular corner point and the integral sample point amplitude value of the third rectangular corner point, calculating the sum q of the integral sample point amplitude value of the second rectangular corner point and the integral sample point amplitude value of the fourth rectangular corner point, calculating the difference value of subtracting q from q, and taking the difference value as an intermediate value. The first rectangular corner point is a relative starting point of the integration InLine direction and the crossLine direction, and is the point with the least integration times in the rectangular corner points corresponding to the sampling point data.
And S10443, calculating and obtaining new sampling point data corresponding to the sampling point data according to the intermediate value, the main measuring line filling length, the contact measuring line filling length and a preset formula.
Specifically, after the intermediate value is obtained, the intermediate value, the main survey line filling length, and the crossline filling length are brought into a preset formula, so that new sampling point data corresponding to the sampling point data can be obtained through calculation. And for each sample point data in the integral slice data, performing the processing procedures of the steps S10441, S10442 and S10443 to obtain new sample point data of each sample point data, wherein the new sample point data of each sample point data form the suppressed post-stack seismic data.
For example, the preset formula is E ═ M/[ (2c-1) (2d-1) ], where E is new sample data corresponding to the sample data, M is a middle value, c is a main measuring line filling length, and d is a tie measuring line filling length.
On the basis of the above embodiments, further, the processing method for the acquisition footprint of the seismic data according to the embodiment of the present invention further includes:
and performing mixed wave combination processing according to the post-stack seismic data and the recovered post-stack seismic data to obtain the post-stack seismic data after denoising.
Specifically, after the recovered post-stack seismic data is obtained, the post-stack seismic data after denoising can be obtained by performing mixed wave combination processing according to the post-stack seismic data and the recovered post-stack seismic data.
Fig. 4 is a schematic flow chart of a processing method for acquiring footprints of seismic data according to still another embodiment of the present invention, and as shown in fig. 4, on the basis of the foregoing embodiments, further, performing a mixing combination process according to the post-stack seismic data and the recovered post-stack seismic data to obtain denoised post-stack seismic data includes:
s401, obtaining the mixing percentage corresponding to each sampling point position;
specifically, as the acquisition footprint of the seismic data is more obvious in the shallow layer and becomes weaker along with the increase of the depth of the stratum, different mixing percentages can be set for the seismic data at different depths, and the mixing percentage is larger as the depth is deeper. The post-stack seismic data is generally expressed in terms of depth or time, the same depth or time corresponds to the same mixing percentage, each sampling point position corresponds to depth or time, and the mixing percentage corresponding to each sampling point position can be obtained.
For example, the sampling point positions with the mixing percentages of 1 and 0 in the post-stack seismic data are set, and then the mixing percentages corresponding to the sampling point positions between the sampling point positions with the mixing percentages of 1 and 0 are obtained in an interpolation mode.
S402, calculating a mixed wave combination result of the stacked seismic data and the recovered stacked seismic data according to the mixed wave percentage corresponding to each sampling point position, and using the mixed wave combination result as the denoised stacked seismic data.
Specifically, after the mixing percentage of each sampling point position is obtained, the mixing combination result of the sampling point data of the post-stack seismic data and the recovered post-stack seismic data is calculated according to the mixing percentage corresponding to each sampling point position, that is, the result of calculating the weighted average value of the sampling point data at the same position in the post-stack seismic data and the recovered post-stack seismic data is used as the sampling point data at the same sampling point position in the post-stack seismic data after denoising.
For example, the mixing percentage corresponding to the sampling point position X is g, the sampling point data of the post-stack seismic data at the sampling point position X is m, the sampling point data of the recovered post-stack seismic data at the sampling point position X is n, and then the sampling point data of the post-stack seismic data at the sampling point position X after denoising is m (1-g) + ng.
The following describes an implementation process of the processing method for acquiring the footprint of the seismic data according to an embodiment of the present invention.
And pre-stack preprocessing is carried out on the acquired seismic data. And carrying out observation system arrangement, static correction, velocity analysis, energy compensation and noise attenuation processing on the seismic data to obtain the post-stack seismic data. The purpose of pre-stack preprocessing is to avoid generating new acquisition footprints, weaken and even suppress the original data acquisition footprints, and finally create conditions for post-stack acquisition footprint suppression.
And performing data filling processing on the stacked seismic data, filling up seismic channel data missing from the stacked seismic data in the InLine and CrossLine directions to obtain filled stacked seismic data, and forming a regular matrix on the filled stacked seismic data in the InLine and CrossLine directions so as to facilitate the pressing of the acquired footprint.
And rotating the filled post-stack seismic data, and converting (InLine, CrossLine, T) arrangement data into data distributed in a (T, InLine, CrossLine) arrangement mode, wherein T represents time. Writing the rotated post-stack seismic data into a temporary file in a time slicing manner,
reading the 1 st time slice data of the rotated post-stack seismic data in the temporary file stored in a slice manner, wherein the partial data of the 1 st time slice data is as shown in fig. 5, firstly integrating the amplitude on the 1 st time slice data along each InLine direction on the 1 st time slice data, and forming the integrated time slice data in the InLine direction, wherein the integration starting point is the minimum InLine number. The integration result in the InLine direction of the partial data of the 1 st time slice data shown in fig. 5 is shown in fig. 6, and zero values on the edges are integrated values after being filled with zero values, and only one row or column is shown in the figure. And then integrating the amplitude of the integrated time slice data in the InLine direction again along the crossLine direction, wherein the integration starting point is the minimum crossLine number, and obtaining the integrated time slice data corresponding to the 1 st time slice data. The integration time slice data in the InLine direction shown in FIG. 6, and the integration results along the CrossLine direction are shown in FIG. 7.
Two space change step lengths in the InLine direction and the CrossLine direction are set and are respectively InLine direction step length 6 and CrossLine direction step length 6. And dividing the two space step lengths by 2 respectively to obtain the filling length of the main measuring line as 3 and the filling length of the contact measuring line as 3. 3 slice sample values are filled along each InLine minimum CMP direction of the integration time slice data corresponding to the 1 st time slice data, the filling length is 3, 3 slice sample values are filled along each InLine maximum CMP direction of the integration time slice data corresponding to the 1 st time slice data, the filling length is 3, and the boundary extension filling is zero filling in consideration of the influence of the collected footprint filtering boundary effect on the collected footprint pressing. Similarly, 3 slice sample values are filled along each CrossLine minimum line direction of the integrated time slice data corresponding to the 1 st time slice data, the filling length is 3, and 3 slice sample values are filled along each CrossLine maximum line direction of the integrated time slice data corresponding to the 1 st time slice data, and the filling length is also 3. Through the data filling processing, the filled integral slice data corresponding to the 1 st time slice data can be obtained.
Reading in the 1 st time slice data again, acquiring any sample data on the 1 st time slice data, and taking the sample data as a center, finding a rectangle with side lengths of 3 and 3 respectively on the filled integral slice data corresponding to the 1 st time slice data, as shown in fig. 5, where the sample data on the 5 th column and the 5 th row is 3.1. The 5 th row and 5 th column of samples in fig. 5 correspond to the centers of the 3 rd row and 3 rd column, the 8 th row and 8 th column, the 8 rd row and 3 rd column and the 3 rd row and 8 th column of the filled integrated slice data corresponding to the 1 st time slice data, as shown in fig. 7, the black dotted square frame at the 5.5 th row and 5.5 th column position, because the position of the sample value of the filled integrated slice data corresponding to the 1 st time slice data is different from the position of the sample value of the 1 st time slice data by 0.5 column and 0.5 row, a rectangle with a length and width of 6 is formed with this black dotted square frame as the rectangle center, the amplitude value of the sample at the 3 rd row and 3 rd column corresponding to the first rectangle is 10.5 plus the amplitude value of the sample at the 8 th row and 8 th column at the opposite corner thereof 126.8, and then the amplitude values of the sample at the 3 rd row and 8 th column at the other two opposite corners are subtracted by 39.2 and 36.9 of the amplitude value of the sample at the 8 th column and 8.8 th row and 8, the sample amplitude difference is integrated with the diagonal of the 6-rectangle in length and width, i.e. the median corresponding to the sample data 3.1, which is 126.8+ 10.5-39.2-36.9-61.2. The sample data 3.1 on the 5 th row and the 5 th column of the 1 st time slice data is replaced with 2.4 by 61.2/(6-1)/(6-1) ═ 2.4, and new sample data 2.4 corresponding to the sample data 3.1 is obtained. All the sample data of the 1 st time slice data are replaced by the above method to form a new 1 st time slice data, and part of the new 1 st time slice data is shown in fig. 8.
This new 1 st time slice data is saved with a new temporary file. And then sequentially processing the 2 nd time slice data, the 3 rd time slice data and the 4 th time slice data by adopting the same processing method for the 1 st time slice data until each time slice data in the rotated post-stack seismic data is processed, thereby obtaining the post-stack seismic data after being suppressed.
And reading the pressed post-stack seismic data from the new temporary file, performing reverse rotation, and rotating the slice mode (T, InLine, CrossLine) data into data distributed in an (InLine, CrossLine, T) arrangement mode to obtain the reverse-rotation post-stack seismic data. And removing data corresponding to the filled seismic data during filling processing from the reversely rotated post-stack seismic data to obtain recovered post-stack seismic data.
And finally, predicting the mixing percentage at 100ms as 10%, 60% at 1000ms and 98% at 5000ms according to the severity of the acquired footprint on the time period, and interpolating the mixing percentages at different time points to ensure that each seismic channel sample point has a corresponding mixing percentage. And reading out the stacked seismic data and the recovered stacked seismic data, and performing mixed wave combination processing to obtain the denoised stacked seismic data, wherein the denoised stacked seismic data is the seismic data from which the acquisition footprint is removed.
The collected footprint effect maps before and after data denoising in fig. 5 and 8 can be compared with fig. 9 to 14. Comparing fig. 9 and fig. 10, it can be found that the in-phase axis is smoother after the acquisition footprint is pressed, and the resolution of the signal-to-noise ratio is significantly improved. Comparing fig. 11 and fig. 12, it can be found that the main energy on the frequency spectrum after the collected footprint is pressed is basically consistent, but the frequency spectrum energy change of 40Hz-80Hz is larger, which indicates that the collected footprint energy is mainly distributed in the range of 40Hz-80Hz, and also indicates that the distribution trend of the collected footprint energy in the InLinehe and CrossLine directions and the distribution range in the frequency domain are accurately judged, and the amplitude-preserving capability after the collected footprint is pressed is stronger. Comparing fig. 13 and fig. 14, it can be found that the acquired footprint slice is smoother after the acquired footprint is pressed by the processing method for acquiring footprints of seismic data provided by the embodiment of the present invention, and it is difficult to recognize the acquired footprint with naked eyes in fig. 14, which illustrates that the processing method for acquiring footprints of seismic data provided by the embodiment of the present invention improves the signal-to-noise ratio and resolution on the seismic data time slice, eliminates the interference of the acquired footprint irrelevant to the prediction of the oil and gas reservoir, and creates favorable conditions for the explanation after the next stack and the successful prediction of the oil and gas reservoir.
Fig. 15 is a schematic structural diagram of a processing apparatus for acquiring a footprint of seismic data according to an embodiment of the present invention, and as shown in fig. 15, the processing apparatus for acquiring a footprint of seismic data according to an embodiment of the present invention includes a preprocessing module 1510, a padding processing module 1520, a rotating module 1530, a suppressing module 1540, a derotating module 1550, and a recovery processing module 1560, where:
the preprocessing module 1510 is configured to perform pre-stack preprocessing on the seismic data to obtain post-stack seismic data; the filling-up processing module 1520 is configured to perform filling-up processing on the post-stack seismic data to obtain filled-up post-stack seismic data; the rotation module 1530 is configured to rotate the filled post-stack seismic data to obtain rotated post-stack seismic data; the pressing module 1540 is configured to perform acquisition footprint pressing on each slice data of the rotated post-stack seismic data, and obtain pressed post-stack seismic data; the reverse rotation module 1550 is configured to perform reverse rotation on the pressed post-stack seismic data to obtain reverse-rotated post-stack seismic data; the recovery processing module 1560 is configured to perform recovery processing on the reversely rotated post-stack seismic data to obtain recovered post-stack seismic data.
Specifically, seismic data are acquired in the field to obtain seismic data. The pre-processing module 1510 performs pre-stack pre-processing on the seismic data to obtain post-stack seismic data. The prestack preprocessing includes, but is not limited to, observation system setting, static correction, velocity analysis, energy compensation, noise attenuation, cutting and the like of seismic data, and is set according to actual needs, and the embodiment of the invention is not limited. In the pre-stack pretreatment process, new artificial traces and footprints are avoided from being introduced due to unreasonable pre-stack pretreatment according to a normal treatment flow.
After the post-stack seismic data are obtained, the filling-up processing module 1520 performs filling-up processing on the post-stack seismic data to fill up seismic channel data of the post-stack seismic data missing in the directions of an InLine (InLine) and a CrossLine (CrossLine), so as to obtain filled-up post-stack seismic data, and the filled-up post-stack seismic data form a regular matrix in the directions of the InLine and CrossLine. Wherein missing seismic trace data may be filled with zero values.
After obtaining the filled-in post-stack seismic data, the rotation module 1530 may rotate the filled-in post-stack seismic data to obtain rotated post-stack seismic data, where the rotated post-stack seismic data includes a plurality of slice data, and the slice data may be time slice data or depth slice data. The rotated post-stack seismic data can be stored in a temporary file, so that subsequent data reading is facilitated.
After obtaining the rotated post-stack seismic data, the suppression module 1540 may perform acquisition footprint suppression on each slice data of the rotated post-stack seismic data to obtain a suppressed post-stack seismic data. In an embodiment of the present invention, the slice data is time slice data or depth slice data.
After obtaining the suppressed post-stack seismic data, a counter rotation module 1550 may counter rotate the suppressed post-stack seismic data to obtain counter-rotated post-stack seismic data.
After obtaining the derotated post-stack seismic data, the recovery processing module 1560 performs recovery processing on the derotated post-stack seismic data, removes data corresponding to the complemented seismic trace data from the derotated post-stack seismic data, and obtains recovered post-stack seismic data.
The processing device for acquiring the footprint of the seismic data, provided by the embodiment of the invention, is used for preprocessing the seismic data before stacking to obtain post-stack seismic data, performing filling-up processing on the post-stack seismic data to obtain filled-up post-stack seismic data, rotating the filled-up post-stack seismic data to obtain rotated post-stack seismic data, performing acquisition footprint pressing on each slice of the rotated post-stack seismic data to obtain pressed post-stack seismic data, performing reverse rotation on the pressed post-stack seismic data to obtain reverse-rotated post-stack seismic data, and performing recovery processing on the reverse-rotated post-stack seismic data to obtain recovered post-stack seismic data.
Fig. 16 is a schematic structural diagram of a processing apparatus for acquiring a footprint of seismic data according to another embodiment of the present invention, and as shown in fig. 16, on the basis of the foregoing embodiments, further, the suppressing module 1540 includes a first integrating unit 1541, a second integrating unit 1542, a filling unit 1543, and a diagonal amplitude processing unit 1544, where:
the first integrating unit 1541 is configured to integrate the amplitude of the slice data along a main line direction to obtain integrated slice data in the main line direction; the second integrating unit 1542 is configured to integrate the amplitude of the integrated slice data in the main line direction along the crossline direction to obtain integrated slice data; the filling unit 1543 is configured to fill a slice sample value into the integral slice data in a main line direction according to a set main line filling length, and fill a slice sample value into the integral slice data in an interconnection line direction according to a set interconnection line filling length, so as to obtain filled integral slice data; the diagonal amplitude processing unit 1544 is configured to perform diagonal amplitude processing on each sampling point data in the integrated slice data according to the filled integrated slice data, so as to obtain post-stack seismic data after compaction.
Specifically, for each slice data of the rotational post-stack seismic data, the first integration unit 1541 may integrate the amplitude of the slice data along the InLine direction, where the starting point of integration is the minimum number of inlines, and may obtain integrated slice data in the InLine direction.
After obtaining the integrated slice data in the InLine direction, the second integrating unit 1542 integrates the amplitude of the integrated slice data in the InLine direction along the crossline direction, obtaining the integrated slice data. The integrated slice data may be integrated time offset data or integrated depth offset data.
After obtaining the integrated slice data, the filling unit 1543 may fill the integrated slice data with slice sample values in a InLine direction according to a set InLine filling length, i.e., fill InLine filling lengths of slice sample values along each InLine minimum CMP direction of the integrated slice data, the filling length being an InLine filling length, also fill InLine filling lengths of slice sample values along each InLine maximum CMP direction of the integrated slice data, the filling length also being a InLine filling length, and fill the integrated slice data with slice sample values in an InLine direction according to a set InLine filling length, i.e., fill InLine filling lengths of slice sample values along each CrossLine minimum line direction of the integrated slice data, the filling length being a CrossLine filling length, also fill CrossLine filling lengths of sample values along each CrossLine maximum line direction of the slice, the filling length is also the cross-line filling length, and the filled integral slice data is obtained. The main measuring line filling length and the cross measuring line filling length are preset and are set according to actual needs, and the embodiment of the invention is not limited.
After obtaining the filled integral slice data, the diagonal amplitude processing unit 1544 may perform diagonal amplitude processing on each sample point data in the integral slice data according to the filled integral slice data, that is, for each sample point data in the integral slice data, four corresponding rectangular corner point data may be obtained from the filled integral slice data, and perform data processing on the amplitudes of the four rectangular corner point data, so as to obtain integral sample point amplitude difference slice data.
Fig. 17 is a schematic structural diagram of a processing apparatus for acquiring a footprint of seismic data according to yet another embodiment of the present invention, and as shown in fig. 17, on the basis of the foregoing embodiments, further, a filling unit 1543 according to an embodiment of the present invention includes an acquiring subunit 15431, a first calculating subunit 15432, and a second calculating subunit 15433, where:
the obtaining subunit 15431 is configured to obtain sample point data from the integrated slice data, and obtain a rectangular corner point corresponding to the sample point data from the filled integrated slice data; the sampling point data comprises four rectangular corner points corresponding to the sampling point data to form a rectangle, and the sampling point data is positioned in the center of the rectangle; the first calculating subunit 15432 is configured to calculate a difference between a sum of the integrated sample amplitude values corresponding to the first rectangular corner point and the third opposing rectangular corner point and a sum of the integrated sample amplitude values corresponding to the second rectangular corner point and the fourth opposing rectangular corner point, so as to obtain an intermediate value; the rectangular corner points corresponding to the sampling point data comprise the first rectangular corner point, the second rectangular corner point, the third rectangular corner point and the fourth rectangular corner point; the second calculating subunit 15433 is configured to calculate and obtain new sampling point data corresponding to the sampling point data according to the intermediate value, the main measurement line filling length, the crossline filling length, and a preset formula.
Specifically, the obtaining subunit 15431 obtains one sampling point data from the integral slice data, where the rectangular corner points corresponding to the sampling point data can be obtained from the filled integral slice data, there are four rectangular corner points corresponding to the sampling point data, and the four rectangular corner points are connected to form a rectangle, where the sampling point data is located in the center of the rectangle.
The rectangular corner points corresponding to the sampling point data comprise a first rectangular corner point, a second rectangular corner point, a third rectangular corner point and a fourth rectangular corner point, wherein the first rectangular corner point is opposite to the third rectangular corner point and is two end points of a rectangular diagonal line, and the second rectangular corner point is opposite to the fourth rectangular corner point and is two end points of another diagonal line of the rectangle. The first calculating subunit 15432 calculates a sum p of the integrated sample amplitude value of the first rectangular corner point and the integrated sample amplitude value of the third rectangular corner point, calculates a sum q of the integrated sample amplitude value of the second rectangular corner point and the integrated sample amplitude value of the fourth rectangular corner point, and calculates a difference obtained by subtracting q from q, and takes the difference as an intermediate value. The first rectangular corner point is a relative starting point of the integration InLine direction and the crossLine direction, and is the point with the least integration times in the rectangular corner points corresponding to the sampling point data.
After obtaining the intermediate value, the second calculating subunit 15433 brings the intermediate value, the main survey line filling length, and the crossline filling length into a preset formula, so as to calculate and obtain new sample point data corresponding to the sample point data.
The embodiment of the apparatus provided in the embodiment of the present invention may be specifically configured to execute the processing flows of the above method embodiments, and the functions of the apparatus are not described herein again, and refer to the detailed description of the above method embodiments.
Fig. 18 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 18, the electronic device may include: a processor (processor)1801, a communication Interface (Communications Interface)1802, a memory (memory)1803, and a communication bus 1804, wherein the processor 1801, the communication Interface 1802, and the memory 1803 communicate with each other via the communication bus 1804. The processor 1801 may call logic instructions in the memory 1803 to perform the following method: pre-stack preprocessing is carried out on the seismic data to obtain post-stack seismic data; performing filling processing on the post-stack seismic data to obtain filled post-stack seismic data; rotating the filled post-stack seismic data to obtain rotated post-stack seismic data; performing acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data; carrying out reverse rotation on the pressed post-stack seismic data to obtain the reverse-rotation post-stack seismic data; and recovering the reversely rotated post-stack seismic data to obtain the recovered post-stack seismic data.
In addition, the logic instructions in the memory 1803 may be implemented in software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as a stand-alone product. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The present embodiment discloses a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions which, when executed by a computer, enable the computer to perform the method provided by the above-mentioned method embodiments, for example, comprising: pre-stack preprocessing is carried out on the seismic data to obtain post-stack seismic data; performing filling processing on the post-stack seismic data to obtain filled post-stack seismic data; rotating the filled post-stack seismic data to obtain rotated post-stack seismic data; performing acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data; carrying out reverse rotation on the pressed post-stack seismic data to obtain the reverse-rotation post-stack seismic data; and recovering the reversely rotated post-stack seismic data to obtain the recovered post-stack seismic data.
The present embodiment provides a computer-readable storage medium, which stores a computer program, where the computer program causes the computer to execute the method provided by the above method embodiments, for example, the method includes: pre-stack preprocessing is carried out on the seismic data to obtain post-stack seismic data; performing filling processing on the post-stack seismic data to obtain filled post-stack seismic data; rotating the filled post-stack seismic data to obtain rotated post-stack seismic data; performing acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data; carrying out reverse rotation on the pressed post-stack seismic data to obtain the reverse-rotation post-stack seismic data; and recovering the reversely rotated post-stack seismic data to obtain the recovered post-stack seismic data.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In the description herein, reference to the description of the terms "one embodiment," "a particular embodiment," "some embodiments," "for example," "an example," "a particular example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A method of processing seismic data acquisition footprints, comprising:
pre-stack preprocessing is carried out on the seismic data to obtain post-stack seismic data;
performing filling processing on the post-stack seismic data to obtain filled post-stack seismic data;
rotating the filled post-stack seismic data to obtain rotated post-stack seismic data;
performing acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data;
carrying out reverse rotation on the pressed post-stack seismic data to obtain the reverse-rotation post-stack seismic data;
and recovering the reversely rotated post-stack seismic data to obtain the recovered post-stack seismic data.
2. The method of claim 1, wherein the performing acquisition footprint suppression on each slice of the rotated post-stack seismic data, obtaining suppressed post-stack seismic data comprises:
integrating the amplitude of the slice data along the direction of the main measuring line to obtain integrated slice data in the direction of the main measuring line;
integrating the amplitude of the integral slice data in the main measuring line direction along the direction of the cross measuring line to obtain integral slice data;
filling a slice sample value into the integral slice data in the main measuring line direction according to the set main measuring line filling length, and filling a slice sample value into the integral slice data in the cross measuring line direction according to the set cross measuring line filling length to obtain the filled integral slice data;
and carrying out diagonal amplitude processing on each sampling point data in the integral slice data according to the filled integral slice data to obtain suppressed post-stack seismic data.
3. The method of claim 2, wherein the performing diagonal amplitude processing on each sample point data in the integrated slice data from the filled integrated slice data to obtain squashed post-stack seismic data comprises:
acquiring sampling point data from the integral slice data, and acquiring rectangular corner points corresponding to the sampling point data from the filled integral slice data; the sampling point data comprises four rectangular corner points corresponding to the sampling point data to form a rectangle, and the sampling point data is positioned in the center of the rectangle;
calculating the difference value between the sum of the amplitude values of the integral sampling points corresponding to the first rectangular corner point and the opposite third rectangular corner point and the sum of the amplitude values of the integral sampling points corresponding to the second rectangular corner point and the opposite fourth rectangular corner point to obtain an intermediate value; the rectangular corner points corresponding to the sampling point data comprise the first rectangular corner point, the second rectangular corner point, the third rectangular corner point and the fourth rectangular corner point;
and calculating to obtain new sampling point data corresponding to the sampling point data according to the intermediate value, the main measuring line filling length, the contact measuring line filling length and a preset formula.
4. The method of any of claims 1 to 3, further comprising:
and performing mixed wave combination processing according to the post-stack seismic data and the recovered post-stack seismic data to obtain the post-stack seismic data after denoising.
5. The method of claim 4, wherein performing a mixture combining process on the post-stack seismic data and the recovered post-stack seismic data to obtain denoised post-stack seismic data comprises:
acquiring the mixing percentage corresponding to each sampling point position;
and calculating the mixed wave combination result of the sample point data of the post-stack seismic data and the restored post-stack seismic data according to the mixed wave percentage corresponding to each sample point position, and using the mixed wave combination result as the sample point data of the post-stack seismic data after denoising.
6. A seismic data acquisition footprint processing apparatus, comprising:
the preprocessing module is used for preprocessing the seismic data before the stack to obtain post-stack seismic data;
the filling-up processing module is used for filling up the post-stack seismic data to obtain the post-stack seismic data after filling up;
the rotation module is used for rotating the filled post-stack seismic data to obtain the rotated post-stack seismic data;
the pressing module is used for carrying out acquisition footprint pressing on each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data;
the anti-rotation module is used for carrying out anti-rotation on the pressed post-stack seismic data to obtain the anti-rotation post-stack seismic data;
and the recovery processing module is used for performing recovery processing on the reversely rotated post-stack seismic data to obtain the recovered post-stack seismic data.
7. The apparatus of claim 6, wherein the compression module comprises:
the first integration unit is used for integrating the amplitude of the slice data along the direction of the main measuring line to obtain integrated slice data in the direction of the main measuring line;
the second integral unit is used for integrating the amplitude of the integral slice data in the main measuring line direction along the direction of the cross measuring line to obtain integral slice data;
the filling unit is used for filling the integral slice data with slice sample values in the main measuring line direction according to the set main measuring line filling length and filling the integral slice data with slice sample values in the tie measuring line direction according to the set tie measuring line filling length to obtain the filled integral slice data;
and the diagonal amplitude processing unit is used for carrying out diagonal amplitude processing on each sampling point data in the integral slice data according to the filled integral slice data to obtain the suppressed post-stack seismic data.
8. The apparatus of claim 7, wherein the filling unit comprises:
the acquisition subunit is configured to acquire sampling point data from the integral slice data, and acquire a rectangular corner point corresponding to the sampling point data from the filled integral slice data; the sampling point data comprises four rectangular corner points corresponding to the sampling point data to form a rectangle, and the sampling point data is positioned in the center of the rectangle;
the first calculating subunit is used for calculating the difference value between the sum of the amplitude values of the integral sampling points corresponding to the first rectangular corner point and the opposite third rectangular corner point and the sum of the amplitude values of the integral sampling points corresponding to the second rectangular corner point and the opposite fourth rectangular corner point to obtain a middle value; the rectangular corner points corresponding to the sampling point data comprise the first rectangular corner point, the second rectangular corner point, the third rectangular corner point and the fourth rectangular corner point;
and the second calculating subunit is used for calculating and obtaining new sampling point data corresponding to the sampling point data according to the intermediate value, the main measuring line filling length, the contact measuring line filling length and a preset formula.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method according to any of claims 1 to 5 are implemented when the computer program is executed by the processor.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011155133.8A CN112255684B (en) | 2020-10-26 | 2020-10-26 | Processing method and device for acquisition footprint of seismic data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011155133.8A CN112255684B (en) | 2020-10-26 | 2020-10-26 | Processing method and device for acquisition footprint of seismic data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112255684A true CN112255684A (en) | 2021-01-22 |
| CN112255684B CN112255684B (en) | 2024-03-01 |
Family
ID=74261572
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011155133.8A Active CN112255684B (en) | 2020-10-26 | 2020-10-26 | Processing method and device for acquisition footprint of seismic data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112255684B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112859169A (en) * | 2021-03-08 | 2021-05-28 | 中国石油天然气集团有限公司 | High-dimensional pre-stack seismic data filtering method and device |
| CN113436138A (en) * | 2021-03-31 | 2021-09-24 | 成都飞机工业（集团）有限责任公司 | Image preprocessing method for aviation structural part identification based on digital image |
| US11835671B2 (en) | 2021-07-29 | 2023-12-05 | Saudi Arabian Oil Company | Method and system for eliminating seismic acquisition footprint through geological guidance |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101852866A (en) * | 2010-05-13 | 2010-10-06 | 中国石油天然气股份有限公司 | Post-stack seismic data filtering method |
| CN103558636A (en) * | 2013-11-08 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for collecting attenuation of footprints from post-stack seismic data |
| US20140288843A1 (en) * | 2013-03-22 | 2014-09-25 | Cgg Services Sa | Vector-dip filtering of seismic data in the time-frequency domain |
| US20150276956A1 (en) * | 2014-03-28 | 2015-10-01 | Cgg Services Sa | Wave-equation based processing and analysis of imaged seismic data |
| WO2016126898A1 (en) * | 2015-02-04 | 2016-08-11 | Schlumberger Technology Corporation | Acquisition footprint attenuation in seismic data |
| CN106353816A (en) * | 2016-08-09 | 2017-01-25 | 中国石油天然气集团公司 | Seismic acquisition footprint noise suppression method and system |
| US20170192118A1 (en) * | 2016-01-05 | 2017-07-06 | Schlumerger Technology Corporation | Amplitude Inversion on Partitioned Depth Image Gathers Using Point Spread Functions |
| CN107247290A (en) * | 2017-07-06 | 2017-10-13 | 西安交通大学 | A kind of seismic data noise drawing method filtered based on space-time fractional order |
| US20180188399A1 (en) * | 2017-01-05 | 2018-07-05 | China Petroleum & Chemical Corporation | Anisotropy matching filtering for attenuation of seismic migration artifacts |
| CN110618449A (en) * | 2018-06-20 | 2019-12-27 | 中国石油化工股份有限公司 | Method and system for processing seismic data |
| CN111751870A (en) * | 2019-03-26 | 2020-10-09 | 中国石油天然气集团有限公司 | Method and device for suppressing multiple waves between stacked layers |
-
2020
- 2020-10-26 CN CN202011155133.8A patent/CN112255684B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101852866A (en) * | 2010-05-13 | 2010-10-06 | 中国石油天然气股份有限公司 | Post-stack seismic data filtering method |
| US20140288843A1 (en) * | 2013-03-22 | 2014-09-25 | Cgg Services Sa | Vector-dip filtering of seismic data in the time-frequency domain |
| CN103558636A (en) * | 2013-11-08 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for collecting attenuation of footprints from post-stack seismic data |
| US20150276956A1 (en) * | 2014-03-28 | 2015-10-01 | Cgg Services Sa | Wave-equation based processing and analysis of imaged seismic data |
| WO2016126898A1 (en) * | 2015-02-04 | 2016-08-11 | Schlumberger Technology Corporation | Acquisition footprint attenuation in seismic data |
| US20170192118A1 (en) * | 2016-01-05 | 2017-07-06 | Schlumerger Technology Corporation | Amplitude Inversion on Partitioned Depth Image Gathers Using Point Spread Functions |
| CN106353816A (en) * | 2016-08-09 | 2017-01-25 | 中国石油天然气集团公司 | Seismic acquisition footprint noise suppression method and system |
| US20180188399A1 (en) * | 2017-01-05 | 2018-07-05 | China Petroleum & Chemical Corporation | Anisotropy matching filtering for attenuation of seismic migration artifacts |
| CN107247290A (en) * | 2017-07-06 | 2017-10-13 | 西安交通大学 | A kind of seismic data noise drawing method filtered based on space-time fractional order |
| CN110618449A (en) * | 2018-06-20 | 2019-12-27 | 中国石油化工股份有限公司 | Method and system for processing seismic data |
| CN111751870A (en) * | 2019-03-26 | 2020-10-09 | 中国石油天然气集团有限公司 | Method and device for suppressing multiple waves between stacked layers |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112859169A (en) * | 2021-03-08 | 2021-05-28 | 中国石油天然气集团有限公司 | High-dimensional pre-stack seismic data filtering method and device |
| CN113436138A (en) * | 2021-03-31 | 2021-09-24 | 成都飞机工业（集团）有限责任公司 | Image preprocessing method for aviation structural part identification based on digital image |
| CN113436138B (en) * | 2021-03-31 | 2023-07-21 | 成都飞机工业（集团）有限责任公司 | Image preprocessing method for aviation structure identification based on digital image |
| US11835671B2 (en) | 2021-07-29 | 2023-12-05 | Saudi Arabian Oil Company | Method and system for eliminating seismic acquisition footprint through geological guidance |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112255684B (en) | 2024-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112255684A (en) | Processing method and device for collecting footprint of seismic data | |
| Zhang et al. | Physical wavelet frame denoising | |
| CN103995289B (en) | Time-varying method of mixed phase wavelet extraction based on time-frequency spectrum analog | |
| Tugume et al. | Precambrian crustal structure in Africa and Arabia: evidence lacking for secular variation | |
| CN104536044B (en) | The interpolation denoising method and system of a kind of geological data | |
| CN101201407B (en) | Relative non-high-frequency leakage equivalent N-drop smooth spectrum analog deconvolution method | |
| DE3149525A1 (en) | F-K FILTERING MULTIPLE REFLECTIONS FROM A SEISMIC SECTION | |
| CN110907995B (en) | Reverse time migration method and device for VSP seismic data in well | |
| US20210311218A1 (en) | Surface wave prediction and removal from seismic data | |
| CN111736224B (en) | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data | |
| CN106646596A (en) | Undulatory sea surface ghost reflection suppression method and device | |
| EP2992361B1 (en) | Device and method for optimization of 4d and 3d seismic data | |
| US20060235620A1 (en) | Method for treating seismic cubes corresponding to obtained for common zone at different times | |
| CN112198547A (en) | Deep or ultra-deep seismic data processing method and device | |
| CN108919345B (en) | Submarine cable land detection noise attenuation method | |
| US7830746B2 (en) | Method for treating seismic data corresponding to acquisitions obtained for a common zone by means of underwater seismic receivers and by means of surface receivers | |
| Wang et al. | Ground roll wave suppression based on wavelet frequency division and radial trace transform | |
| Chen et al. | Adaptive acquisition footprint suppression based on a 3D stationary wavelet transform: A case study from China | |
| CN113589384B (en) | Pre-stack gather amplitude-preserving denoising method based on characteristic that signal changes along with offset distance | |
| CN112255689B (en) | Method for analyzing fidelity velocity of seismic data in multiple wave development area | |
| CN110231651A (en) | A kind of earthquake data before superposition denoising method based on vibration signal resolution of vectors principle | |
| CN112213775A (en) | Fidelity frequency-boosting method for high-coverage-frequency pre-stack seismic data | |
| CN112799132B (en) | Micro-local linear noise suppression method and device | |
| Buzlukov et al. | Prestack data enhancement using local traveltime approximation | |
| CN117233844B (en) | Data processing method and system for high-precision seismic imaging of fault development area |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |