CN112255684B - Processing method and device for acquisition footprint of seismic data - Google Patents
Processing method and device for acquisition footprint of seismic data Download PDFInfo
- Publication number
- CN112255684B CN112255684B CN202011155133.8A CN202011155133A CN112255684B CN 112255684 B CN112255684 B CN 112255684B CN 202011155133 A CN202011155133 A CN 202011155133A CN 112255684 B CN112255684 B CN 112255684B
- 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.)
- Active
Links
- 238000003672 processing method Methods 0.000 title abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 71
- 230000007812 deficiency Effects 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 39
- 230000001502 supplementing effect Effects 0.000 claims abstract description 32
- 230000006835 compression Effects 0.000 claims abstract description 12
- 238000007906 compression Methods 0.000 claims abstract description 12
- 230000010354 integration Effects 0.000 claims description 29
- 238000002156 mixing Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 17
- 238000011084 recovery Methods 0.000 claims description 16
- 238000003825 pressing Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 10
- 238000007781 pre-processing Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims 2
- 230000001629 suppression Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 18
- 238000005056 compaction Methods 0.000 description 11
- 230000007547 defect Effects 0.000 description 10
- 238000009826 distribution Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 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
- 238000009933 burial Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 238000004519 manufacturing process 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
- 238000012805 post-processing Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
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 processing device for acquisition footprints of seismic data, wherein the method comprises the following steps: pre-stacking the seismic data to obtain post-stack seismic data; performing deficiency supplementing treatment on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing; rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data; performing acquisition footprint compression on each slice data of the rotated post-stack seismic data to obtain compressed post-stack seismic data; performing anti-rotation on the pressed post-stack seismic data to obtain anti-rotation post-stack seismic data; and recovering the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data. The device is used for executing the method. The processing method and the processing device for the acquisition footprint of the seismic data improve the signal-to-noise ratio and the 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 a processing device for acquisition footprint of seismic data.
Background
Currently, oil and gas exploration is evolving in both breadth and depth. The former is aimed at finding new field of oil and gas fields for distant exploration; the latter requires the search for reservoirs with great depths of burial and high complexity, as well as the resolution of fine-structure and fine-description of reservoir parameters. Along with the continuous deep exploration of lithologic hydrocarbon reservoirs and hidden hydrocarbon reservoirs, geological targets are changed from original oil and gas searching structures to lithologic oil and gas searching structures, so that the problems of fine description of the parameters of the hydrocarbon reservoirs and reservoirs are well solved, and the problems of signal to noise ratio of imaging data are involved.
The collected footprint noise is also called as collected trace noise, is seismic noise generated by an artificial factor, is an artificial trace left in the process of collecting and processing seismic data, and is expressed by the fact that a regular amplitude variation false image appears on a seismic section depth or a time slice, and the prediction precision and efficiency of an oil and gas reservoir are often seriously influenced by part of the collected footprint noise energy. For such compaction of the acquisition footprint, partial attenuation of the acquisition footprint is typically achieved by combining and interpolation methods, typically prior to lamination; the post-stack acquisition footprint pressing method is more, and comprises an inclination angle filtering method, an F-Kx-Ky filtering method, a wavelet transformation method, a self-adaptive method and the like. The use of weighting functions during DMO or prestack migration can also partially attenuate the acquisition footprint. However, different collection footprint generation mechanisms are different, and different collection footprint pressing methods have the own applicability and the limitations of the method, so that the development of the targeted collection footprint pressing method research has important significance.
Disclosure of Invention
The embodiment of the invention provides a processing method and a processing device for acquisition footprints of seismic data, which aim at the acquisition footprints of random distribution or local related distribution, and solve the problem of how to improve the signal-to-noise ratio and the resolution of the seismic data.
In one aspect, the invention provides a method for processing a collection footprint of seismic data, comprising the following steps:
pre-stacking the seismic data to obtain post-stack seismic data;
performing deficiency supplementing treatment on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing;
rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data;
performing acquisition footprint compression on each slice data of the rotated post-stack seismic data to obtain compressed post-stack seismic data;
performing anti-rotation on the pressed post-stack seismic data to obtain anti-rotation post-stack seismic data;
and recovering the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data.
In another aspect, the present invention provides a processing device for acquiring a footprint of seismic data, comprising:
the preprocessing module is used for pre-stacking the seismic data to obtain post-stack seismic data;
The deficiency supplementing processing module is used for carrying out deficiency supplementing processing on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing;
the rotation module is used for rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data;
the pressing module is used for collecting 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 anti-rotation post-stack seismic data;
and the recovery processing module is used for carrying out recovery processing on the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data.
In yet another aspect, the present invention provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the method for processing a collection footprint of seismic data as described in any of the embodiments above when the program is executed.
In yet another aspect, the present invention provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the method for processing a collection footprint of seismic data as described in any of the embodiments above.
According to the processing method and the processing device for the acquisition footprint of the seismic data, pre-stack pretreatment is carried out on the seismic data to obtain post-stack seismic data, defect filling processing is carried out on the post-stack seismic data to obtain post-stack seismic data after defect filling, rotation is carried out on the post-stack seismic data after defect filling to obtain post-stack seismic data after rotation, acquisition footprint compaction is carried out on each slice data of the post-stack seismic data after rotation to obtain post-stack seismic data after compaction, anti-rotation is carried out on the post-stack seismic data after compaction to obtain anti-rotation post-stack seismic data, recovery processing is carried out on the post-stack seismic data after anti-rotation to obtain post-stack seismic data after recovery, waveform characteristics of the original seismic data can be kept, and signal to noise ratio and resolution ratio of the seismic data are improved.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art. In the drawings:
Fig. 1 is a flow chart of a processing method for acquiring a footprint of seismic data according to an embodiment of the invention.
Fig. 2 is a flow chart of a processing method for acquiring a footprint of seismic data according to another embodiment of the invention.
Fig. 3 is a flow chart of a processing method for acquiring a footprint of seismic data according to still another embodiment of the invention.
Fig. 4 is a flowchart of a processing method for acquiring a footprint of seismic data according to still another embodiment of the invention.
Fig. 5 is a schematic diagram of a portion of the 1 st time slice data according to an embodiment of the present invention.
Fig. 6 is a schematic diagram of integration results in an 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 integration time slice data corresponding to a portion of the 1 st time slice data according to an embodiment of the present invention.
Fig. 8 is a schematic diagram of a portion of the data of the new 1 st time slice data according to an embodiment of the present invention.
FIG. 9 is a schematic representation of post-stack seismic data acquisition prior to footprint compaction in accordance with one embodiment of the present invention.
FIG. 10 is a schematic representation of post-stack seismic data after denoising according to one embodiment of the present invention.
FIG. 11 is a schematic diagram of frequency amplitude spectra of post-stack seismic data before acquisition footprint compaction in accordance with an embodiment of the invention.
FIG. 12 is a schematic diagram of the frequency amplitude spectrum of denoised post-stack seismic data provided by an embodiment of the invention.
FIG. 13 is a schematic illustration of a slice display of post-stack seismic data before acquisition footprint compaction in accordance with an embodiment of the invention;
FIG. 14 is a schematic view of a slice display of denoised post-stack seismic data according to one 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 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 invention.
Fig. 17 is a schematic structural diagram of a processing device for acquiring a footprint of seismic data according to still another embodiment of the invention.
Fig. 18 is a schematic physical structure of an electronic device according to an embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention will be described in further detail with reference to the accompanying drawings. The exemplary embodiments of the present invention and their descriptions herein are for the purpose of explaining the present invention, but are not to be construed as limiting the invention. It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be arbitrarily combined with each other.
Fig. 1 is a flow chart of a processing method of a collection footprint of seismic data according to an embodiment of the invention, as shown in fig. 1, the processing method of a collection footprint of seismic data according to an embodiment of the invention includes:
s101, pre-stack pretreatment is carried out on seismic data to obtain post-stack seismic data;
specifically, the seismic data is acquired in the field to obtain the seismic data. And pre-stacking the seismic data to obtain post-stack seismic data. The pre-stack preprocessing includes, but is not limited to, performing a system for observing seismic data, static correction, velocity analysis, energy compensation, noise attenuation, excision and other modes, and is set according to actual needs, and the embodiment of the invention is not limited. In the pre-stack pretreatment process, new artifacts and footprints are avoided from being unreasonably introduced due to the pre-stack pretreatment according to the normal treatment flow. The seismic data may be a common-center point gather (Common Middle Point, CMP for short). 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, carrying out defect supplementing treatment on the post-stack seismic data to obtain post-stack seismic data after defect supplementing;
Specifically, after the post-stack seismic data is obtained, the post-stack seismic data is subjected to a gap filling process to fill up the seismic channel data of the post-stack seismic data, which is missing in the directions of a main line (Inline) and a cross line (cross line), so as to obtain the post-stack seismic data after the gap filling, and the post-stack seismic data after the gap filling forms a regular matrix in the directions of the Inline and the cross line. Wherein missing seismic trace data may be padded with zero values.
For example, because the observation system is unevenly distributed, the Inline direction line number is inconsistent with the cross line direction line number and the CMP number, the post-stack seismic data needs to be subjected to the deficiency treatment, so that the post-stack seismic data after deficiency forms a regular matrix in the whole work area, and the subsequent data treatment is convenient. Inline and cross line of the missing seismic trace data are recorded while the seismic trace data are supplemented, so that the subsequent recovery of the data is facilitated.
S103, rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data;
specifically, after the post-stack seismic data after the deficiency is obtained, the post-stack seismic data after the deficiency can be rotated to obtain rotated post-stack seismic data, the rotated post-stack seismic data includes a plurality of slice data, and the slice data can be time slice data or depth slice data. The rotated post-stack seismic data can be stored in a temporary file, so that the subsequent data reading is facilitated.
For example, the post-stack seismic data after the deficiency is represented as (InLine, crossLine, Z), Z may be time or depth, the post-stack seismic data after the deficiency is rotated is represented as (Z, inLine, crossLine), the post-stack seismic data after the rotation includes a plurality of slice data, each slice data corresponds to a Z value, and the data reading is performed in the form of slice data, which is suitable for processing mass data and can improve the data reading efficiency.
S104, collecting footprint pressing is carried out on each slice data of the rotated post-stack seismic data, and pressed post-stack seismic data are obtained;
specifically, after the rotated post-stack seismic data is obtained, each slice data of the rotated post-stack seismic data may be subjected to acquisition footprint compaction to obtain compacted post-stack seismic data. In the embodiment of the present invention, the slice data is time slice data or depth slice data.
S105, performing anti-rotation on the pressed post-stack seismic data to obtain anti-rotation post-stack seismic data;
specifically, after the post-stack seismic data after the compacting is obtained, the post-stack seismic data after compacting is rotated reversely, the direction of data rotation is opposite to the direction of rotation in the step S103, and the reverse process in the step S103 is performed, so that the post-stack seismic data after the reverse rotation can be obtained.
For example, the compressed post-stack seismic data is in the form of slices (Z, inLine, crossLine), and the anti-rotated post-stack seismic data obtained after anti-rotation is represented as (InLine, crossLine, Z) as data in a spread distribution.
S106, recovering the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data;
specifically, after the anti-rotation post-stack seismic data is obtained, recovering the anti-rotation post-stack seismic data, and removing data corresponding to the seismic channel data supplemented in the step S102 from the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data.
According to the processing method for the acquisition footprint of the seismic data, pre-stack pretreatment is carried out on the seismic data, post-stack seismic data are obtained, defect-filling treatment is carried out on the post-stack seismic data, post-stack seismic data after defect filling are obtained, rotation is carried out on the post-stack seismic data after defect filling, post-stack seismic data after rotation are obtained, acquisition footprint compaction is carried out on each slice data of the post-stack seismic data after rotation, pressed post-stack seismic data are obtained, anti-rotation post-stack seismic data after compaction are obtained, recovery treatment is carried out on the post-stack seismic data after anti-rotation, the waveform characteristics of the original seismic data can be kept, the acquisition footprint in the original seismic data is compacted, and the signal to noise ratio and the resolution ratio of the seismic data are improved.
Fig. 2 is a flow chart of a processing method for acquiring a footprint of seismic data according to another embodiment of the present invention, as shown in fig. 2, further, based on the above embodiments, the step of performing acquisition footprint compaction on each slice data of the rotated post-stack seismic data, to obtain the compacted post-stack seismic data includes:
s1041, integrating the amplitude of the slice data along the direction of the main line to obtain integrated slice data along the direction of the main line;
specifically, for each slice data of the rotated post-stack seismic data, the amplitude of the slice data may be integrated along the InLine direction, the starting point of the integration is a minimum InLine number, and integrated slice data in the InLine direction may be obtained.
S1042, integrating the amplitude of the integral slice data in the direction of the main line along the direction of the cross-line to obtain integral slice data;
specifically, after the integration slice data in the InLine direction is obtained, the amplitude of the integration slice data in the InLine direction is integrated in the crossline direction to obtain the integration slice data. The integrated slice data may be integrated time offset data or integrated depth offset data.
S1043, filling slice sample point values for the integral slice data in the direction of the main line according to the set main line filling length, and filling slice sample point values for the integral slice data in the direction of the cross line according to the set cross line filling length, so as to obtain filled integral slice data;
specifically, after the integrated slice data is obtained, slice sample values of the integrated slice data may be filled in a main line direction according to a set main line filling length, that is, a plurality of slice sample values of a main line filling length are filled in a minimum InLine CMP direction of the integrated slice data, the filling length is a main line filling length, a plurality of slice sample values of a main line filling length are also filled in a maximum InLine CMP direction of the integrated slice data, the filling length is also a main line filling length, and the integrated slice data is filled with slice sample values of a cross line filling length in a cross line direction according to the set cross line filling length, that is, a plurality of slice sample values of a cross line filling length are filled in a minimum InLine direction of the integrated slice data, the filling length is a cross line filling length, the cross line filling length sample values are also filled in a maximum InLine direction of a slice, the filling length is also a cross line filling length, and the filling length is also a cross line filling length, so as to obtain the integrated slice data after filling. The filling length of the main line and the filling length of the cross-connection line are preset and set according to actual needs, and the embodiment of the invention is not limited.
For example, for the integral slice data, two spatial variation step sizes of the Inline direction and the cross line direction are set, the step sizes of the Inline direction and the cross line direction are respectively equal or unequal, the direction of the stronger the collected footprint energy distribution is usually, the selected step size is longer, and the collected footprint energy intensity distribution can be obtained through integral slice data observation. The CrossLine fill length is set to c, which is half the InLine direction step, i.e., c=a/2, and d, which is half the CrossLine direction step, b, i.e., d=b/2. And when the integral slice data is filled, filling c slice sample point values along the minimum CMP direction of each Inline of the integral slice data, wherein the filling length is c, filling c slice sample point values along the maximum CMP direction of each Inline of the slice, the filling length is c, filling d slice sample point values along the minimum line direction of each cross line of the slice, the filling length is d, filling d slice sample point values along the maximum line direction of each cross line of the slice, and the filling length is d, thereby obtaining the filled integral slice data. Wherein the filled sample values may take on zero values.
S1044, performing diagonal amplitude processing on each sample point data in the integral slice data according to the filled integral slice data to obtain integral sample point amplitude difference value slice 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 data processing may be performed on the amplitudes of the four rectangular corner point data, so as to obtain integral sample point amplitude difference slice data. And (3) carrying out the data processing process of the steps S1041, S1042, S1043 and S1044 on each slice data of the rotated post-stack seismic data, and obtaining the pressed post-stack seismic data.
Fig. 3 is a flowchart of a processing method for acquiring a footprint of seismic data according to another embodiment of the present invention, as shown in fig. 3, further, based on the above embodiments, the performing diagonal amplitude processing on each sample data in the integral slice data according to the filled integral slice data, to obtain integral sample amplitude difference slice data includes:
S10441, acquiring one sample point data from the integral slice data, and acquiring a rectangular corner point corresponding to the sample point data from the filled integral slice data; four rectangular corner points corresponding to the sample data form a rectangle, and the sample data are positioned in the center of the rectangle;
specifically, one sample point data is obtained from the integral slice data, rectangular corner points corresponding to the sample point data can be obtained from the filled integral slice data, four rectangular corner points corresponding to the sample point data are connected to form a rectangle, and the sample point data is located at the center of the rectangle.
For example, any sample data is acquired in the integral slice data, and a rectangle with a side length of a and b can be found in the filled integral slice data by taking the sample data as a center, and the sample data is exactly the center of the rectangle.
S10442, calculating a difference value of the sum of integrated sample point amplitude values corresponding to the first rectangular corner point and the opposite third rectangular corner point and the sum of integrated sample point amplitude values corresponding to the second rectangular corner point and the opposite fourth rectangular corner point, and obtaining an intermediate value; the rectangular corner points corresponding to the sample points 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 sample points comprise 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 diagonal line of a rectangle, 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. 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, and then calculating the difference value of subtracting q from q, wherein the difference value is taken as an intermediate value. The first rectangular corner point is a relative starting point of the integral Inline direction and the cross line direction, and is the point with the least integral times in the rectangular corner points corresponding to the sample points.
S10443, calculating to obtain new sample point data corresponding to the sample point data according to the intermediate value, the main line filling length, the cross line filling length and a preset formula.
Specifically, after the intermediate value is obtained, the intermediate value, the main line filling length and the cross-line filling length are brought into a preset formula, so that new sample point data corresponding to the sample point data can be obtained through calculation. For each sample data in the integral slice data, the processing procedures of steps S10441, S10442 and S10443 are performed to obtain new sample data of each sample data, where the new sample data of each sample data constitutes post-stack seismic data after compacting.
For example, the preset formula is e=m/[ (2 c-1) (2 d-1) ], where E is new sample data corresponding to the sample data, M is an intermediate value, c is a filling length of the main line, and d is a filling length of the cross line.
On the basis of the above embodiments, the processing method for acquiring the footprint of the seismic data provided by the embodiment of the invention further includes:
and carrying out mixed wave combination processing according to the post-stack seismic data and the recovered post-stack seismic data to obtain denoised post-stack seismic data.
Specifically, after the recovered post-stack seismic data is obtained, performing mixed wave combination processing according to the post-stack seismic data and the recovered post-stack seismic data, so that denoised post-stack seismic data can be obtained.
Fig. 4 is a flow chart of a processing method for acquiring a footprint of seismic data according to still another embodiment of the present invention, as shown in fig. 4, further, based on the above embodiments, the performing a wave mixing and combining 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 a mixing percentage corresponding to each sample point position;
specifically, because the acquisition footprint of the seismic data is more obvious in the shallow layer, the acquisition footprint becomes weaker as the depth of the stratum increases, so that different mixing percentages can be set for the seismic data with different depths, and the deeper the depth, the greater the mixing percentage. The post-stack seismic data is typically represented in depth or time, where the same depth or time corresponds to the same mixing percentage, and each sample point corresponds to the depth or time, and the mixing percentage corresponding to each sample point may be obtained.
For example, sample positions with the mixing percentages of 1 and 0 in the post-stack seismic data are set, and then the mixing percentages corresponding to each sample position between the sample positions with the mixing percentages of 1 and 0 are obtained by interpolation.
S402, calculating a mixed wave combination result of the post-stack seismic data and the sample point data of the restored post-stack seismic data according to the mixed wave percentage corresponding to each sample point position, and taking the mixed wave combination result as the sample point data of the post-stack seismic data after denoising.
Specifically, after obtaining the mixing percentage of each sample point position, calculating a mixing combination result of the post-stack seismic data and the sample point data of the post-stack seismic data after recovery according to the mixing percentage corresponding to each sample point position, namely calculating a weighted average value of the sample point data of the same sample point position in the post-stack seismic data after recovery and the post-stack seismic data after recovery, and taking the weighted average value as the sample point data of the same sample point position of the post-stack seismic data after denoising.
For example, the mixing percentage corresponding to the sample point position X is g, the sample point data of the post-stack seismic data at the sample point position X is m, the sample point data of the post-stack seismic data at the sample point position X after recovery is n, and then the sample point data of the post-stack seismic data at the sample point position X after denoising is m (1-g) +ng.
The implementation process of the processing method for acquiring the footprint of the seismic data provided by the embodiment of the invention is described in the following with a specific embodiment.
Pre-stack pretreatment is carried out on the collected seismic data. And carrying out observation system setting, static correction, speed analysis, energy compensation and noise attenuation treatment on the seismic data to obtain post-stack seismic data. The aim of pre-stack pretreatment is to avoid generating new acquisition footprints, and also to weaken and even press original data acquisition footprints, and finally to create conditions for post-stack acquisition footprint pressing.
And carrying out data deficiency processing on the post-stack seismic data, and supplementing the seismic channel data of which the post-stack seismic data is missing in the Inline and cross line directions to obtain the post-stack seismic data after deficiency, wherein the post-stack seismic data after deficiency forms a rule matrix in the Inline and cross line directions so as to facilitate the subsequent compression of the acquisition footprint.
The post-stack seismic data after the deficiency is rotated, the (InLine, crossLine, T) arrangement data are converted into data distributed in the (T, inLine, crossLine) arrangement mode, and T represents time, so that the main purpose of the post-processing is to facilitate the slicing processing of mass data and prevent the problems of overlarge data quantity and insufficient program memory. Writing the rotated post-stack seismic data into a temporary file in a time slicing manner,
The 1 st time slice data of the post-stack seismic data after rotation in the temporary file stored in a slice mode is read in, the partial data of the 1 st time slice data is shown in fig. 5, the amplitude of the 1 st time slice data is integrated along each Inline direction on the 1 st time slice data, the integration starting point is the minimum Inline number, and the integration time slice data in the Inline direction is formed. 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, where the zero value on the edge is the integrated value after filling the zero value, and only one line or one column is shown in the figure. And then integrating the amplitude on the integration time slice data in the Inline direction along the cross line direction again, wherein the integration starting point is the minimum cross line number, and the integration time slice data corresponding to the 1 st time slice data is obtained. Integration time slice data in the InLine direction shown in fig. 6, and integration results in the CrossLine direction are shown in fig. 7.
Setting two space change step sizes of the Inline direction and the cross line direction, wherein the two space change step sizes are respectively the Inline direction step size 6 and the cross line direction step size 6. The two spatial steps are divided by 2 to obtain a crossline filling length of 3 and a crossline filling length of 3. 3 slice sample point values are filled along the minimum CMP direction of each Inline of the integration time slice data corresponding to the 1 st time slice data, the filling length is 3, 3 slice sample point values are filled along the maximum CMP direction of each Inline of the integration time slice data corresponding to the 1 st time slice data, the filling length is 3, and the boundary expansion filling is zero filling considering the influence of the collected footprint filtering boundary effect on the collected footprint pressing in the embodiment of the invention. And filling 3 slice sample point values along the minimum line direction of each cross line of the integration time slice data corresponding to the 1 st time slice data, wherein the filling length is 3, and filling 3 slice sample point values along the maximum line direction of each cross line of the integration time slice data corresponding to the 1 st time slice data, wherein the filling length is 3. After the data filling processing, the filled integral slice data corresponding to the 1 st time slice data can be obtained.
And reading in the 1 st time slice data again, acquiring any sample point data on the 1 st time slice data, and taking the sample point data as a center, and finding a rectangle with side lengths of 3 and 3 on the filled integral slice data corresponding to the 1 st time slice data, wherein the sample point data of the 5 th column and the 5 th row in the figure 5 is 3.1. In fig. 5, the sample points in column 5 and column 5 correspond to the centers of column 3, column 8, column 3 and column 8 in column 3 on the filled integral slice data corresponding to the 1 st time slice data, as shown in fig. 7, the black dotted line box in column 5.5 of column 5.5 is obtained by subtracting the sample point amplitude value 39.2 of column 3 and the sample point amplitude value 36.9 of column 3 on the other two diagonal angles from the center of the rectangle with the black dotted line box as the center, forming a rectangle with the length and width of 6, using the sample point amplitude value 10.5 of column 3 and column 3 corresponding to the corner point of the first rectangle of the rectangle to add the sample point amplitude value 126.8 of column 8 on the diagonal angle, and subtracting the sample point amplitude value 39.9 of column 3 and the sample point amplitude value 36.9 of column 3 on the other two diagonal angles to obtain a sample point amplitude value of between the rectangle and the sample point amplitude value of column 3 and the rectangle with the length of 6 = 2.5 and the sample point amplitude value of 6, namely, and the sample amplitude value between the sample points between the two diagonal values of column 3 and column 3 = 6.5.5 and 6.5 = 2.6 are obtained. The new sample data 2.4 corresponding to the sample data 3.1 is obtained by replacing the sample data 3.1 on the 5 th row and 5 th column of the 1 st time slice data with 61.2/(6-1)/(6-1) =2.4 and 2.4. All sample data of the 1 st time-sliced data are replaced by the above method to form a new 1 st time-sliced data, and part of the new 1 st time-sliced data is shown in fig. 8.
The new 1 st time slice data is saved with a new temporary file. And then adopting the same processing method for the 1 st time slice data to sequentially process the 2 nd time slice data, the 3 rd time slice data and the 4 th time slice data until each time slice data in the rotated post-stack seismic data is processed, thereby obtaining the pressed post-stack seismic data.
And reading out the pressed post-stack seismic data from the new temporary file, performing anti-rotation, and rotating the data in a slicing mode (T, inLine, crossLine) into data distributed in a (InLine, crossLine, T) arrangement mode to obtain the anti-rotated post-stack seismic data. And removing the data corresponding to the seismic data which are supplemented in the supplementing process from the anti-rotation post-stack seismic data, and obtaining the recovered post-stack seismic data.
Finally, predicting that the mixing percentage at 100ms is 10%, the mixing percentage at 1000ms is 60% and the mixing percentage at 5000ms is 98% according to the severity degree of the acquired footprint on the time period, and interpolating the mixing percentages at different time points to ensure that each seismic trace sample point has a corresponding mixing percentage. And reading out the post-stack seismic data and the recovered post-stack seismic data, and performing mixed wave combination processing to obtain denoised post-stack seismic data, wherein the denoised post-stack seismic data is the seismic data after the acquisition footprint is removed.
The capture footprint effect graphs before and after denoising the data in fig. 5 and 8 can be compared with fig. 9 to 14. Comparing fig. 9 and 10, it can be found that the same phase axis is smoother after the acquisition footprint is pressed, and the signal to noise ratio resolution is significantly improved. Comparing fig. 11 and fig. 12, it can be found that the main energy on the spectrum after the collected footprint is pressed is basically consistent, but the spectrum energy is changed greatly from 40Hz to 80Hz, which means that the collected footprint energy is mainly distributed in the range from 40Hz to 80Hz, and meanwhile, the invention also means that the distribution trend of the collected footprint energy in the directions of InLinehe and CrossLine and the distribution range of the frequency domain are accurately judged, and the amplitude-preserving capability after the collected footprint is pressed is strong. Comparing fig. 13 and fig. 14 can find that after the processing method of the acquisition footprint of the seismic data provided by the embodiment of the invention suppresses the acquisition footprint, the acquisition footprint slice is smoother, and the acquisition footprint is difficult to identify by naked eyes in fig. 14, which illustrates that the processing method of the acquisition footprint of the seismic data provided by the embodiment of the invention improves the signal-to-noise ratio and the resolution on the time slice of the seismic data, eliminates the acquisition footprint interference irrelevant to the reservoir prediction, and creates an advantage for the next post-stack interpretation of the reservoir success prediction.
Fig. 15 is a schematic structural diagram of a processing device for acquiring a footprint of seismic data according to an embodiment of the invention, as shown in fig. 15, where the processing device for acquiring a footprint of seismic data according to an embodiment of the invention includes a preprocessing module 1510, a deficiency processing module 1520, a rotation module 1530, a pressing module 1540, a counter-rotation module 1550 and a recovery processing module 1560, wherein:
the preprocessing module 1510 is used for pre-stacking the seismic data to obtain post-stack seismic data; the deficiency supplementing processing module 1520 is configured to perform deficiency supplementing processing on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing; the rotation module 1530 is configured to rotate the post-stack seismic data after the deficiency, to obtain rotated post-stack seismic data; the compression module 1540 is configured to perform acquisition footprint compression on each slice data of the rotated post-stack seismic data, to obtain compressed post-stack seismic data; the anti-rotation module 1550 is configured to perform anti-rotation on the pressed post-stack seismic data to obtain anti-rotated post-stack seismic data; the recovery processing module 1560 is configured to perform recovery processing on the anti-rotated post-stack seismic data, and obtain recovered post-stack seismic data.
Specifically, the seismic data is acquired in the field to obtain the seismic data. The pre-processing module 1510 performs pre-stack pre-processing on the seismic data to obtain post-stack seismic data. The pre-stack preprocessing includes, but is not limited to, performing a system for observing seismic data, static correction, velocity analysis, energy compensation, noise attenuation, excision and other modes, and is set according to actual needs, and the embodiment of the invention is not limited. In the pre-stack pretreatment process, new artifacts and footprints are avoided from being unreasonably introduced due to the pre-stack pretreatment according to the normal treatment flow.
After the post-stack seismic data is obtained, a deficiency supplementing processing module 1520 performs deficiency supplementing processing on the post-stack seismic data to supplement the seismic trace data of the post-stack seismic data missing in the directions of a main line (InLine) and a cross line (CrossLine), so as to obtain the post-stack seismic data after deficiency supplementing, and the post-stack seismic data after deficiency supplementing forms a regular matrix in the directions of the InLine and the cross line. Wherein missing seismic trace data may be padded with zero values.
After obtaining the post-stack seismic data after the deficiency, the rotation module 1530 may rotate the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data, where the rotated post-stack seismic data includes a plurality of slice data, which may be time slice data or depth slice data. The rotated post-stack seismic data can be stored in a temporary file, so that the subsequent data reading is facilitated.
After obtaining the rotated post-stack seismic data, the compression module 1540 may perform acquisition footprint compression on each slice data of the rotated post-stack seismic data to obtain compressed post-stack seismic data. In the embodiment of the present invention, the slice data is time slice data or depth slice data.
After the compressed post-stack seismic data is obtained, the anti-rotation module 1550 performs anti-rotation on the compressed post-stack seismic data, and may obtain anti-rotated post-stack seismic data.
After the anti-rotation post-stack seismic data is obtained, the restoration processing module 1560 performs restoration processing on the anti-rotation post-stack seismic data, and removes data corresponding to the supplemented seismic trace data from the anti-rotation post-stack seismic data to obtain restored post-stack seismic data.
The processing device for acquiring the footprint of the seismic data provided by the embodiment of the invention performs pre-stack pretreatment on the seismic data to obtain post-stack seismic data, performs defect supplementing processing on the post-stack seismic data to obtain post-stack seismic data after defect supplementing, rotates the post-stack seismic data after defect supplementing to obtain rotated post-stack seismic data, acquires the footprint of each slice data of the rotated post-stack seismic data to obtain pressed post-stack seismic data, performs anti-rotation on the pressed post-stack seismic data to obtain anti-rotated post-stack seismic data, performs recovery processing on the anti-rotated post-stack seismic data to obtain recovered post-stack seismic data, can maintain waveform characteristics of the original seismic data, presses the acquired footprint in the original seismic data, and improves the signal-to-noise ratio and resolution of the seismic data.
Fig. 16 is a schematic structural diagram of a processing device for acquiring a footprint of seismic data according to another embodiment of the invention, as shown in fig. 16, further, based on the above embodiments, a compacting 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 integration unit 1541 is configured to integrate the amplitude of the slice data along the inline direction, to obtain integrated slice data along the inline direction; the second integrating unit 1542 is configured to integrate the amplitude of the integrated slice data in the crossline direction along the crossline direction, so as to obtain integrated slice data; the filling unit 1543 is configured to fill slice sample values into the integrated slice data in the direction of the main line according to the set filling length of the main line, and fill slice sample values into the integrated slice data in the direction of the cross line according to the set filling length of the cross line, so as to obtain filled integrated slice data; the diagonal amplitude processing unit 1544 is configured to perform diagonal amplitude processing on each sample data in the integrated slice data according to the filled integrated slice data, so as to obtain post-stack seismic data after compacting.
Specifically, for each slice data of the rotated post-stack seismic data, the first integration unit 1541 may integrate the amplitude of the slice data along the InLine direction, where the start point of integration is the minimum InLine number, and may obtain integrated slice data in the InLine direction.
After obtaining the integration slice data in the InLine direction, the second integration unit 1542 integrates the amplitude of the integration slice data in the InLine direction along the crossline direction to obtain the integration slice data. The integrated slice data may be integrated time offset data or integrated depth offset data.
After the integrated slice data is obtained, the filling unit 1543 may fill the integrated slice data with slice sample values in the InLine direction according to a set InLine fill length, that is, fill the InLine fill length with slice sample values in the InLine direction along each InLine minimum CMP direction of the integrated slice data, fill the InLine fill length with slice sample values in the InLine fill length along each InLine maximum CMP direction of the integrated slice data, fill the InLine fill length with slice sample values in the InLine fill length according to the set InLine fill length, and fill the integrated slice data with slice sample values in the CrossLine direction according to the set CrossLine fill length, that is, fill the CrossLine fill length with slice sample values in the CrossLine minimum line direction of the integrated slice data along each CrossLine minimum line direction, fill the CrossLine fill length with slice sample values in the CrossLine fill length along each CrossLine maximum line direction of the slice, and fill the CrossLine fill length with slice sample values in the CrossLine fill length according to the CrossLine fill length. The filling length of the main line and the filling length of the cross-connection line are preset and 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, may obtain, for each sample point data in the integral slice data, four corresponding rectangular corner point data from the filled integral slice data, and perform data processing on 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 device for acquiring a footprint of seismic data according to another embodiment of the invention, as shown in fig. 17, further, based on the above embodiments, a filling unit 1543 provided by the embodiment of the invention includes an acquiring subunit 15431, a first calculating subunit 15432, and a second calculating subunit 15433, where:
the acquiring subunit 15431 is configured to acquire one sample point data from the integral slice data, and acquire a rectangular corner point corresponding to the sample point data from the filled integral slice data; four rectangular corner points corresponding to the sample data form a rectangle, and the sample data are positioned in the center of the rectangle; the first calculating subunit 15432 is configured to calculate a difference value between a sum of integrated sample point amplitude values corresponding to the first rectangular corner and the opposite third rectangular corner and a sum of integrated sample point amplitude values corresponding to the second rectangular corner and the opposite fourth rectangular corner, to obtain an intermediate value; the rectangular corner points corresponding to the sample points 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 sample data corresponding to the sample data according to the intermediate value, the main line filling length, the cross-line filling length, and a preset formula.
Specifically, the acquiring subunit 15431 acquires one sample point data from the integral slice data, and may obtain rectangular corner points corresponding to the sample point data in the filled integral slice data, where four rectangular corner points corresponding to the sample point data are connected to form a rectangle, and the sample point data is located in the center of the rectangle.
The rectangular corner points corresponding to the sample points 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 diagonal line of a rectangle, 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 point amplitude value of the first rectangular corner point and the integrated sample point amplitude value of the third rectangular corner point, calculates a sum q of the integrated sample point amplitude value of the second rectangular corner point and the integrated sample point amplitude value of the fourth rectangular corner point, calculates a difference value of q minus q, and uses the difference value as an intermediate value. The first rectangular corner point is a relative starting point of the integral Inline direction and the cross line direction, and is the point with the least integral times in the rectangular corner points corresponding to the sample points.
After obtaining the intermediate value, the second calculating subunit 15433 brings the intermediate value, the inline filling length, and the crossline filling length into a preset formula, and may calculate and obtain new sample data corresponding to the sample data.
The embodiment of the apparatus provided in the embodiment of the present invention may be specifically used to execute the processing flow of each method embodiment, and the functions thereof are not described herein again, and may refer to the detailed description of the method embodiments.
Fig. 18 is a schematic physical structure of an electronic device according to an embodiment of the present invention, 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 perform communication with each other through the communication bus 1804. The processor 1801 may call logic instructions in the memory 1803 to perform the following method: pre-stacking the seismic data to obtain post-stack seismic data; performing deficiency supplementing treatment on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing; rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data; performing acquisition footprint compression on each slice data of the rotated post-stack seismic data to obtain compressed post-stack seismic data; performing anti-rotation on the pressed post-stack seismic data to obtain anti-rotation post-stack seismic data; and recovering the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data.
Further, the logic instructions in the memory 1803 may be implemented in the form of software functional units and may be stored in a computer readable storage medium when sold or used as a stand alone product. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform 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, random Access Memory), a magnetic disk, or an optical disk, or 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, are capable of performing the methods provided by the above-described method embodiments, for example comprising: pre-stacking the seismic data to obtain post-stack seismic data; performing deficiency supplementing treatment on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing; rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data; performing acquisition footprint compression on each slice data of the rotated post-stack seismic data to obtain compressed post-stack seismic data; performing anti-rotation on the pressed post-stack seismic data to obtain anti-rotation post-stack seismic data; and recovering the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data.
The present embodiment provides a computer-readable storage medium storing a computer program that causes the computer to execute the methods provided by the above-described method embodiments, for example, including: pre-stacking the seismic data to obtain post-stack seismic data; performing deficiency supplementing treatment on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing; rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data; performing acquisition footprint compression on each slice data of the rotated post-stack seismic data to obtain compressed post-stack seismic data; performing anti-rotation on the pressed post-stack seismic data to obtain anti-rotation post-stack seismic data; and recovering the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data.
It will be appreciated by those skilled in the art that 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 flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations 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 of the present specification, reference to the terms "one embodiment," "one 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, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. 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 foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the invention, and is not meant to limit the scope of the invention, but to limit the invention to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (8)
1. A method for processing an acquisition footprint of seismic data, comprising:
pre-stacking the seismic data to obtain post-stack seismic data;
performing deficiency supplementing treatment on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing;
rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data;
performing acquisition footprint compression on each slice data of the rotated post-stack seismic data to obtain compressed post-stack seismic data;
performing anti-rotation on the pressed post-stack seismic data to obtain anti-rotation post-stack seismic data;
recovering the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data;
the step of acquiring footprint compression of each slice data of the rotated post-stack seismic data, the step of obtaining compressed post-stack seismic data comprises the following steps:
integrating the amplitude of the slice data along the direction of the main line to obtain integrated slice data along the direction of the main line;
integrating the amplitude of the integral slice data in the direction of the main line along the direction of the cross-line to obtain secondary integral slice data;
Filling slice sample values of the secondary integral slice data in the direction of the main line according to the set main line filling length, and filling slice sample values of the secondary integral slice data in the direction of the cross line according to the set cross line filling length, so as to obtain filled integral slice data;
and carrying out diagonal amplitude processing on each sample data in the quadratic integral slice data according to the filled integral slice data to obtain pressed post-stack seismic data.
2. The method of claim 1, wherein said performing diagonal amplitude processing on each sample data in said quadratic integrator slice data from said filled integrator slice data to obtain consolidated post-stack seismic data comprises:
acquiring one sample point data from the secondary integral slice data, and acquiring a rectangular corner point corresponding to the sample point data from the filled integral slice data; four rectangular corner points corresponding to the sample data form a rectangle, and the sample data are positioned in the center of the rectangle;
calculating the difference value of the sum of the integrated sample point amplitude values corresponding to the first rectangular corner point and the opposite third rectangular corner point and the sum of the integrated sample point amplitude values 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 sample points 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 sample point data corresponding to the sample point data according to the intermediate value, the main line filling length, the cross-line filling length and a preset formula.
3. The method according to claim 1 or 2, further comprising:
and carrying out mixed wave combination processing according to the post-stack seismic data and the recovered post-stack seismic data to obtain denoised post-stack seismic data.
4. The method of claim 3, wherein performing a wave mixing combining process based 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 a mixed wave combination result of the post-stack seismic data and the restored post-stack seismic data according to the mixed wave percentage corresponding to each sampling point, and taking the mixed wave combination result as the denoised post-stack seismic data.
5. A processing device for acquiring a footprint of seismic data, comprising:
the preprocessing module is used for pre-stacking the seismic data to obtain post-stack seismic data;
the deficiency supplementing processing module is used for carrying out deficiency supplementing processing on the post-stack seismic data to obtain post-stack seismic data after deficiency supplementing;
The rotation module is used for rotating the post-stack seismic data after the deficiency to obtain rotated post-stack seismic data;
the pressing module is used for collecting 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 anti-rotation post-stack seismic data;
the recovery processing module is used for carrying out recovery processing on the anti-rotation post-stack seismic data to obtain recovered post-stack seismic data;
wherein, the suppression module includes:
the first integration unit is used for integrating the amplitude of the slice data along the direction of the main line to obtain integrated slice data along the direction of the main line;
the second integration unit is used for integrating the amplitude of the integration slice data in the direction of the main line along the direction of the cross-line to obtain secondary integration slice data;
the filling unit is used for filling slice sample point values into the secondary integral slice data in the direction of the main line according to the set filling length of the main line, and filling slice sample point values into the secondary integral slice data in the direction of the cross line according to the set filling length of the cross line, so as to obtain filled integral slice data;
And the diagonal amplitude processing unit is used for performing diagonal amplitude processing on each sample point data in the secondary integral slice data according to the filled integral slice data to obtain pressed post-stack seismic data.
6. The apparatus of claim 5, wherein the filling unit comprises:
the acquisition subunit is used for acquiring one sample point data from the secondary integral slice data and acquiring a rectangular corner point corresponding to the sample point data from the filled integral slice data; four rectangular corner points corresponding to the sample data form a rectangle, and the sample data are positioned in the center of the rectangle;
a first calculating subunit, configured to calculate a difference value between a sum of integrated sample point amplitude values corresponding to the first rectangular corner and the opposite third rectangular corner and a sum of integrated sample point amplitude values corresponding to the second rectangular corner and the opposite fourth rectangular corner, to obtain an intermediate value; the rectangular corner points corresponding to the sample points 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 sample point data corresponding to the sample point data according to the intermediate value, the main line filling length, the cross line filling length and a preset formula.
7. 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 processor implements the steps of the method of any of claims 1 to 4 when the computer program is executed.
8. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method according to any one of claims 1 to 4.
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 CN112255684A (en) | 2021-01-22 |
| CN112255684B true 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) |
Families Citing this family (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 |
| 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 |
Citations (7)
| 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 |
| 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 |
| 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 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2781937B1 (en) * | 2013-03-22 | 2022-09-28 | CGG Services SAS | 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 |
| US10295683B2 (en) * | 2016-01-05 | 2019-05-21 | Schlumberger Technology Corporation | Amplitude inversion on partitioned depth image gathers using point spread functions |
| US10677949B2 (en) * | 2017-01-05 | 2020-06-09 | China Petroleum & Chemical Corporation | Anisotropy matching filtering for attenuation of seismic migration artifacts |
-
2020
- 2020-10-26 CN CN202011155133.8A patent/CN112255684B/en active Active
Patent Citations (7)
| 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 |
| 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 |
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112255684A (en) | 2021-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112255684B (en) | Processing method and device for acquisition footprint of seismic data | |
| DE69728411T2 (en) | SIGNAL-RESISTANT VIBRATION SOURCE AND SEISMIC PROCESS WITH SOURCE SEPARATION | |
| CN102998706B (en) | Method and system for attenuating seismic data random noise | |
| CN112083495B (en) | Method for improving resolution ratio based on variation mode decomposition synchronous compression wavelet transformation | |
| CN110907995B (en) | Reverse time migration method and device for VSP seismic data in well | |
| CN107388963B (en) | The digital Moiré patterns phase extraction method that wavelet analysis and low-pass filtering are combined | |
| CN111736224B (en) | Method, storage medium and equipment for suppressing linear interference of pre-stack seismic data | |
| Chowdhury et al. | An atomic technique for removal of Gaussian noise from a noisy gray scale image using lowpass-convoluted Gaussian filter | |
| Lin et al. | Random noise suppression of magnetic resonance sounding oscillating signal by combining empirical mode decomposition and time-frequency peak filtering | |
| CN105911590B (en) | A kind of compression process method and device of harmonic noise | |
| CN110007346A (en) | Adaptive intermediate value very noisy based on wavelet decomposition interferes seismic data noise attenuation method | |
| CN112198547A (en) | Deep or ultra-deep seismic data processing method and device | |
| CN109255764A (en) | A kind of method and device removing FIB-SEM image curtain noise | |
| Da Wang et al. | Super resolution convolutional neural network models for enhancing resolution of rock micro-ct images | |
| CN116413805A (en) | Time-varying correction method, device, equipment and related system for wavelet consistency | |
| Wang et al. | Ground roll wave suppression based on wavelet frequency division and radial trace transform | |
| CN113743193A (en) | Pre-stack seismic data linear interference suppression method and system | |
| 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 | |
| CN112799132B (en) | Micro-local linear noise suppression method and device | |
| Trojanowski et al. | A multichannel convolution filter for correlated noise: Microseismic data application | |
| CN111929726A (en) | Seismic coherent data volume processing method and device | |
| CN117233844B (en) | Data processing method and system for high-precision seismic imaging of fault development area | |
| CN118011473A (en) | High-resolution seismic imaging processing method, device and equipment | |
| Feng et al. | Adaptive hybrid diffusion model using variational mode decomposition for edge preserving noise attenuation |
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 |