CN113031071B - Static correction method and device for long wavelength of earthquake wave - Google Patents
Static correction method and device for long wavelength of earthquake wave Download PDFInfo
- Publication number
- CN113031071B CN113031071B CN201911247971.5A CN201911247971A CN113031071B CN 113031071 B CN113031071 B CN 113031071B CN 201911247971 A CN201911247971 A CN 201911247971A CN 113031071 B CN113031071 B CN 113031071B
- Authority
- CN
- China
- Prior art keywords
- data
- arrival
- arrival data
- long wavelength
- target
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 150
- 230000003068 static effect Effects 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000010219 correlation analysis Methods 0.000 claims abstract description 35
- 238000012545 processing Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 13
- 238000010606 normalization Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 230000008030 elimination Effects 0.000 claims description 5
- 238000003379 elimination reaction Methods 0.000 claims description 5
- 238000012882 sequential analysis Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000003325 tomography 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/612—Previously recorded data, e.g. time-lapse or 4D
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 embodiment of the application provides a static correction method and device for a long wavelength of a seismic wave, wherein the method comprises the following steps: performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and setting the first arrival data with the data similarity higher than a preset threshold value as target first arrival data; determining a long wavelength correction amount of the first arrival data according to the target first arrival data and a relative reference time, wherein the relative reference time is determined according to the horizon data; the application can effectively reduce the error generated by the long wavelength problem, thereby improving the construction interpretation precision.
Description
Technical Field
The application relates to the field of seismic data processing, in particular to a static correction method and device for a long wavelength of seismic waves.
Background
When the seismic wave passes through the near-surface low-speed stratum, if the thickness of the near-surface low-speed stratum is large and the transverse change is quick, the seismic signals received at adjacent positions have large time difference, and the static correction processing can correct the time difference caused by the near-surface low-speed stratum, so that the continuity of a same phase axis is improved, the false image is reduced, and the interpretation precision is improved.
Because the calculation of the static correction is based on the assumption of uniform speed or slow speed in the transverse direction and the assumption that the emergent ray is almost vertical emergent, when the actual situation is not consistent with the assumption, the calculation of the static correction amount can generate errors, and when the errors are expressed as a low-frequency long period, the long-wavelength problem is solved, and the static correction method for solving the long-wavelength problem is long-wavelength static correction. The current static correction method is still based on the two-point assumption, so the long wavelength problem is not solved effectively.
Disclosure of Invention
Aiming at the problems in the prior art, the application provides a static correction method and device for the long wavelength of a seismic wave, which can effectively reduce errors generated by the long wavelength problem, thereby improving the construction interpretation precision.
In order to solve at least one of the problems, the application provides the following technical scheme:
in a first aspect, the present application provides a long wavelength static correction method for seismic waves, comprising:
performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and setting the first arrival data with the data similarity higher than a preset threshold value as target first arrival data;
and determining a long wavelength correction amount of the first arrival data according to the target first arrival data and a relative reference time, wherein the relative reference time is determined according to the horizon data.
Further, the cross-correlation analysis is performed on the first arrival data in the seismic data of the target work area and the horizon data of the preset mark layer to obtain the corresponding data similarity, which comprises the following steps:
dividing the first arrival data into a plurality of groups according to a preset offset interval;
and performing cross-correlation analysis on each group of first arrival data and the horizon data respectively to obtain corresponding data similarity.
Further, the setting the first arrival data with the data similarity higher than the preset threshold as the target first arrival data includes:
normalizing each group of first arrival data according to the maximum value and the minimum value of each group of first arrival data with the data similarity higher than a preset threshold value and the maximum value and the minimum value of the horizon data;
and adding the normalized first arrival data of each group to obtain the target first arrival data.
Further, before the determining of the long wavelength correction amount of the first arrival data based on the target first arrival data and the relative reference time, it includes:
and determining the relative reference time according to any one of the average value of the horizon data, the average value of the maximum value and the minimum value of the horizon data and any value of the horizon data between the maximum value and the minimum value.
Further, the determining the long wavelength correction amount of the first arrival data according to the target first arrival data and the relative reference time includes:
a difference between the target first arrival data and the relative reference time is set as a long wavelength correction amount of the first arrival data.
Further, after the determining the long wavelength correction amount of the first arrival data according to the target first arrival data and the relative reference time, further comprising:
determining shot point correction amount and/or wave detection point correction amount of the corresponding target work area according to the long wavelength correction amount of the first arrival data;
and correcting the seismic data by applying the shot correction amount and/or the wave detection point correction amount to generate corrected seismic data.
In a second aspect, the present application provides a long wavelength static correction device for seismic waves, comprising:
the high coherence determining module is used for performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and setting the first arrival data with the data similarity higher than a preset threshold value as target first arrival data;
and the high coherence elimination module is used for determining a long wavelength correction amount of the first arrival data according to the target first arrival data and a relative reference time, wherein the relative reference time is determined according to the horizon data.
Further, the high coherence determining module includes:
the first arrival data grouping unit is used for dividing the first arrival data into a plurality of groups according to a preset offset interval;
and each group of sequential analysis units is used for performing cross-correlation analysis on each group of first arrival data and the horizon data respectively to obtain corresponding data similarity.
Further, the high coherence determining module includes:
each group of normalization processing units is used for performing normalization processing on each group of first arrival data according to the maximum value and the minimum value of each group of first arrival data with the data similarity higher than a preset threshold value and the maximum value and the minimum value of the horizon data;
and each group of processed summarizing units is used for adding the normalized first arrival data of each group to obtain the target first arrival data.
Further, the method further comprises the following steps:
and the relative reference time determining unit is used for determining the relative reference time according to any one of the average value of the horizon data, the average value of the maximum value and the minimum value of the horizon data and any value of the horizon data between the maximum value and the minimum value.
Further, the high coherence elimination module includes:
and a correction amount determination unit configured to set a difference between the target first-arrival data and the relative reference time as a long-wavelength correction amount of the first-arrival data.
Further, the method further comprises the following steps:
a long wavelength correction amount decomposition unit for determining a shot correction amount and/or a detection point correction amount of the corresponding target work area according to the long wavelength correction amount of the first arrival data;
and the seismic data correction unit is used for correcting the seismic data by applying the shot correction amount and/or the geophone correction amount to generate corrected seismic data.
In a third aspect, the present application provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, said processor implementing the steps of said method for long wavelength static correction of seismic waves when said program is executed by said processor.
In a fourth aspect, the present application provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the long wavelength static correction method of seismic waves.
According to the technical scheme, the application provides the static correction method and the static correction device for the long wavelength of the earthquake waves, the first arrival data with long wavelength problems are determined by carrying out cross-correlation analysis on the first arrival data in the earthquake data of the target work area and the horizon data of the preset mark layer, namely, the first arrival data with the data similarity higher than the preset threshold value obtained by the cross-correlation analysis is set as the target first arrival data, and then the long wavelength correction amount of the first arrival data is determined according to the numerical difference between the target first arrival data and the relative reference time, so that the coherence of the first arrival data and the horizon data is eliminated, the error generated by the long wavelength problems is effectively reduced, and the construction interpretation precision is improved.
Drawings
In order to more clearly illustrate the embodiments of the present application 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, and it is obvious that the drawings in the following description are some embodiments of the present application, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a static correction method for long wavelength of seismic waves according to an embodiment of the application;
FIG. 2 is a second flow chart of a static correction method for long wavelength of seismic waves according to an embodiment of the application;
FIG. 3 is a third flow chart of a static correction method for long wavelength of seismic waves according to an embodiment of the application;
FIG. 4 is a flow chart of a static correction method for long wavelength of seismic waves according to an embodiment of the application;
FIG. 5 is a block diagram of a long wavelength static correction device for seismic waves according to an embodiment of the present application;
FIG. 6 is a second block diagram of a static correction device for long wavelength of seismic waves in an embodiment of the application;
FIG. 7 is a third block diagram of a long wavelength static correction device for seismic waves according to an embodiment of the present application;
FIG. 8 is a fourth block diagram of a static correction device for long wavelength of seismic waves in an embodiment of the application;
FIG. 9 is a fifth block diagram of a static correction device for long wavelength of seismic waves in an embodiment of the application;
fig. 10 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more apparent, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application, and it is apparent that the described embodiments are some embodiments of the present application, but not all embodiments of the present application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
Considering that the calculation of the current static correction is based on the assumption that the velocity is transversely uniform or slowly variable, and the assumption that the emergent ray is almost vertical emergent is simultaneously based, when the actual situation is not consistent with the assumption, the calculation of the static correction amount can generate errors, and becomes a long wavelength problem when the errors are represented as low-frequency long periods, the static correction method for solving the long wavelength problem is long wavelength static correction, and the current static correction method is still based on the two-point assumption, so that the long wavelength problem is not effectively solved.
In order to effectively reduce the error generated by the long wavelength problem and thus improve the construction interpretation precision, the application provides an embodiment of a static correction method for the long wavelength of the seismic wave, referring to fig. 1, the static correction method for the long wavelength of the seismic wave specifically comprises the following contents:
step S101: and performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and setting the first arrival data with the data similarity higher than a preset threshold value as target first arrival data.
It will be appreciated that the first arrival data and the horizon data may be time data, and that in seismic exploration, particularly shallow seismic exploration, tomography and seismic prospecting advanced exploration, the first arrival time pickup accuracy directly affects the inversion results of the geological structure of the exploration area.
Optionally, a phase axis with high relative signal to noise ratio, good continuity and small time variation can be preselected as a preset mark layer in the target work area, and cross-correlation analysis is performed on first arrival data in the seismic data of the target work area and horizon data of the preset mark layer to confirm whether long wavelength problems exist or not and specific first arrival data with long wavelength problems.
It will be appreciated that the cross-correlation analysis, i.e. the cross-correlation calculation of two functions in the prior art, reflects the similarity between the two functions (or signals).
Alternatively, a threshold may be preset for setting the first arrival data corresponding to the similarity higher than the threshold as target first arrival data, that is, first arrival data having a long wavelength problem.
Step S102: and determining a long wavelength correction amount of the first arrival data according to the target first arrival data and a relative reference time, wherein the relative reference time is determined according to the horizon data.
It will be understood that the relative reference time may be a time corresponding to an existing reference plane, after static correction, the seismic data is corrected to an absolute elevation plane, where the elevation plane is called a reference plane, and the starting time of the seismic data is changed from 0 before static correction to a time corresponding to the reference plane, which is called a reference time herein, and is used as a correction reference of the target first arrival data obtained in the step S101, so as to remove the correlation between the first arrival data and horizon data, and further reduce the influence on accuracy caused by long wavelength problem.
Specifically, the relative reference time may be determined by various manners related to the horizon data, for example, an average value of the horizon data as a whole is set as the relative reference time, and for another example, an average value of a maximum value and a minimum value of the horizon data is set as the relative reference time, and any value between the maximum value and the minimum value in the horizon data may be set as the relative reference time, and in other embodiments of the present application, the relative reference time may be determined by other calculation manners related to the horizon data.
As can be seen from the above description, the method for static correction of long wavelength of seismic waves provided by the embodiment of the application can determine the target first-arrival data with long wavelength problem by performing cross-correlation analysis on the first-arrival data in the seismic data of the target work area and the horizon data of the preset mark layer, namely, the first-arrival data with the data similarity higher than the preset threshold value obtained by the cross-correlation analysis is set as the target first-arrival data, and then the long wavelength correction amount of the first-arrival data is determined according to the numerical difference between the target first-arrival data and the relative reference time, so as to eliminate the coherence between the first-arrival data and the horizon data, effectively reduce the error generated by long wavelength problem and improve the construction interpretation precision.
In order to overcome the problem that part of the first arrival data may be missing during the cross-correlation analysis, in an embodiment of the static correction method for long wavelength of seismic waves according to the present application, referring to fig. 2, the following may be further included:
step S201: dividing the first arrival data into a plurality of groups according to a preset offset interval.
Step S202: and performing cross-correlation analysis on each group of first arrival data and the horizon data respectively to obtain corresponding data similarity.
It can be understood that when all the first arrival data of the obtained target work area are used for performing cross-correlation analysis, due to the complex operation environment of the actual production site, the problem of partial first arrival data missing may exist, so as to affect subsequent cross-correlation analysis, therefore, the first arrival data having a corresponding relationship with each offset interval may be set as a group according to a set offset, and the cross-correlation analysis is performed on each group of first arrival data and horizon data respectively through each group of first arrival data, so as to obtain the data similarity corresponding to each group, thereby overcoming the problem that the partial missing of the first arrival data may cause affecting subsequent calculation.
In order to simplify the subsequent calculation and reduce the magnitude, in an embodiment of the static correction method for long wavelength of seismic waves of the present application, referring to fig. 3, the following may be further included:
step S301: and carrying out normalization processing on each group of first arrival data according to the maximum value and the minimum value of each group of first arrival data with the data similarity higher than a preset threshold value and the maximum value and the minimum value of the horizon data.
Step S302: and adding the normalized first arrival data of each group to obtain the target first arrival data.
Optionally, one or more groups of first arrival data with data similarity higher than a preset threshold are first arrival data with long wavelength problem, and the maximum value and the minimum value of the group of first arrival data and the maximum value and the minimum value of horizon data are picked up to perform normalization processing, and specifically, the formula of the normalization processing is as follows:
optionally, overlapping the normalized first arrival data of each group to obtain the target first arrival data.
In order to accurately remove coherence, in an embodiment of the method for static correction of long wavelength of seismic waves of the present application, the method may further specifically include the following:
and determining the relative reference time according to any one of the average value of the horizon data, the average value of the maximum value and the minimum value of the horizon data and any value of the horizon data between the maximum value and the minimum value.
Alternatively, the relative reference time may be determined by various manners related to the horizon data, for example, an average value of the horizon data as a whole is set as the relative reference time, and for another example, an average value of a maximum value and a minimum value of the horizon data is set as the relative reference time, and any value between the maximum value and the minimum value in the horizon data may be set as the relative reference time, and in other embodiments of the present application, the relative reference time may be determined by other calculation manners related to the horizon data.
In order to accurately remove coherence, in an embodiment of the method for static correction of long wavelength of seismic waves of the present application, the method may further specifically include the following:
a difference between the target first arrival data and the relative reference time is set as a long wavelength correction amount of the first arrival data.
Alternatively, the difference between the target first-arrival data and the relative reference time may be calculated and set as a long wavelength correction amount for the first-arrival data, and in other embodiments of the present application, other conversion formulas may be used, such as adding a weight calculation factor.
In order to further determine the correction amounts of other elements related to the seismic data according to the long wavelength correction amounts to correct the whole seismic data, in an embodiment of the seismic wave long wavelength static correction method of the present application, referring to fig. 4, the following may be further included:
step S401: and determining the shot point correction amount and/or the wave detection point correction amount of the corresponding target work area according to the long wavelength correction amount of the first arrival data.
Step S402: and correcting the seismic data by applying the shot correction amount and/or the wave detection point correction amount to generate corrected seismic data.
It will be appreciated that after the long wavelength correction of the first arrival data is obtained, shot decomposition and shot decomposition may be performed on the long wavelength correction to obtain final shot and shot statics to complete the correction for the seismic data, specifically, according to the ground consistency assumption, the delay Ri of the shot at the location i and the delay Sj of the shot at the location j are the same for all corresponding seismic traces, so that the delay Ri of the detector group at the location i and the delay Sj of the source at the location j are the same for all corresponding seismic traces. Let the total time shift amount of one track be Tij, then there is tij=ri, ten Sj, the time shift amount of one track relative to the other track (Tij-Tmn) can be found using cross correlation, and the least squares problem is used to minimize the sum of squares of the errors to obtain Ri and Sj.
In order to effectively reduce errors generated by long wavelength problems and thereby improve construction interpretation accuracy, the present application provides an embodiment of a long wavelength static correction device for seismic waves for implementing all or part of the contents of the long wavelength static correction method for seismic waves, referring to fig. 5, the long wavelength static correction device for seismic waves specifically includes the following contents:
the high coherence determining module 10 is configured to perform cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and set the first arrival data with the data similarity higher than a preset threshold value as target first arrival data.
The high coherence elimination module 20 is configured to determine a long wavelength correction amount of the first arrival data according to the target first arrival data and a relative reference time, wherein the relative reference time is determined according to the horizon data.
As can be seen from the above description, the long-wavelength static correction device for seismic waves provided by the embodiment of the application can determine the target first-arrival data with long-wavelength problems by performing cross-correlation analysis on the first-arrival data in the seismic data of the target work area and the horizon data of the preset mark layer, namely, setting the first-arrival data with the data similarity higher than the preset threshold value obtained by the cross-correlation analysis as the target first-arrival data, and then determining the long-wavelength correction amount of the first-arrival data according to the numerical difference between the target first-arrival data and the relative reference time, so as to eliminate the coherence between the first-arrival data and the horizon data, effectively reduce the errors generated by the long-wavelength problems, and improve the construction interpretation precision.
In order to overcome the problem that part of the first arrival data may be missing during the cross-correlation analysis, in one embodiment of the seismic long wavelength static correction apparatus of the present application, referring to fig. 6, the high coherence determining module 10 includes:
the first arrival data grouping unit 11 is configured to divide the first arrival data into a plurality of groups according to a preset offset interval.
And each group of sequential analysis units 12 is configured to perform cross-correlation analysis on each group of first arrival data and the horizon data respectively, so as to obtain corresponding data similarity.
In order to simplify the subsequent calculation and the reduction of the magnitude, in an embodiment of the seismic long wavelength static correction device of the present application, referring to fig. 7, the high coherence determining module 10 includes:
and each group of normalization processing units 13 is configured to perform normalization processing on each group of first arrival data according to the maximum value and the minimum value of each group of first arrival data, and the maximum value and the minimum value of the horizon data, where the data similarity is higher than a preset threshold.
And each group of post-processing summarizing units 14 is used for adding the normalized first arrival data of each group to obtain the target first arrival data.
In order to accurately remove coherence, in an embodiment of the seismic long wavelength static correction device of the present application, it further comprises:
and the relative reference time determining unit is used for determining the relative reference time according to any one of the average value of the horizon data, the average value of the maximum value and the minimum value of the horizon data and any value of the horizon data between the maximum value and the minimum value.
In order to accurately remove the coherence, in an embodiment of the seismic long wavelength static correction device of the present application, referring to fig. 8, the high coherence removing module 20 includes:
a correction amount determining unit 21 for setting a difference between the target first arrival data and the relative reference time as a long wavelength correction amount of the first arrival data.
In order to be able to further determine the correction amounts of other elements related to the seismic data based on the long wavelength correction amounts to correct the entire seismic data, in one embodiment of the seismic wave long wavelength static correction device of the application, see fig. 9, further comprises:
a long wavelength correction amount decomposition unit 31 for determining a shot correction amount and/or a detector correction amount of the corresponding target work area based on the long wavelength correction amount of the first arrival data.
And a seismic data correction unit 32 for correcting the seismic data by applying the shot correction and/or the geophone correction, and generating corrected seismic data.
In order to further explain the scheme, the application also provides a specific application example for realizing the seismic wave long-wavelength static correction method by using the seismic wave long-wavelength static correction device, which specifically comprises the following contents:
step 1, picking up first arrival data and applying the calculated static correction value;
step 2, picking up the layer position data of the mark layer on the overlapped section, which specifically comprises the following steps:
selecting a same phase axis with high relative signal-to-noise ratio, good continuity and small time variation as a mark layer;
step 3, sorting the first arrival data according to the offset group, including:
dividing the offset distances into a plurality of groups according to the given offset distance intervals and the equal intervals;
adding the first arrival data corresponding to the offset range in the same group to obtain the first arrival data of the group;
step 4, comparing the first arrival data with the horizon data according to the offset distance, comprising:
observing the similarity degree of first arrival data and horizon data of different offset groups by adopting a manual method;
judging the similarity degree by adopting a method for solving the correlation coefficient;
recording offset groups with similar fluctuation variation trends;
step 5, extracting offset first arrival information with high correlation degree for normalization processing, comprising the following steps:
extracting the time maximum and minimum values of the layer position of the picking mark layer;
extracting maximum and minimum values of first arrival data corresponding to different offset groups;
the first arrival data corresponding to each group of offset distances are normalized by adopting the following general expressions:
step 6, adding the first arrival information after normalization processing to obtain first arrival superposition data;
and 7, obtaining relative reference time based on the marker layer horizon data, wherein the method comprises the following steps of:
the relative reference time is the average value of the horizon data or the average value of the maximum and minimum values of the horizon data of the mark layer or a value between the maximum and minimum values;
step 8, calculating the difference between each value of the first arrival superposition data and the relative reference time to obtain a long wavelength correction amount;
and 9, decomposing the shot point and the wave detection point for the long wavelength correction amount, and obtaining a final shot point and wave detection point static correction amount through complete correction work.
From the above description, the following technical effects can be achieved by the present application:
the long wavelength problem can be confirmed whether to exist or not through the comparison of the first arrival information and the horizon information, the error caused by long wavelength can be effectively reduced through eliminating the correlation between the first arrival information and the horizon, and the influence of non-earth surface consistency factors during the first arrival correction can be further reduced through decomposing static correction into shot points and detector point static correction amounts, so that the error of the method is effectively reduced.
The embodiment of the present application further provides a specific implementation manner of an electronic device capable of implementing all the steps in the long wavelength static correction method of seismic waves in the foregoing embodiment, and referring to fig. 10, the electronic device specifically includes the following contents:
a processor (processor) 601, a memory (memory) 602, a communication interface (Communications Interface) 603, and a bus 604;
wherein the processor 601, the memory 602, and the communication interface 603 complete communication with each other through the bus 604; the communication interface 603 is used for implementing information transmission among the seismic wave long wavelength static correction device, the online service system, the client device and other participating mechanisms;
the processor 601 is configured to invoke a computer program in the memory 602, where the processor executes the computer program to implement all the steps in the long-wavelength static correction method for seismic waves in the above embodiment, for example, the processor executes the computer program to implement the following steps:
step S101: and performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and setting the first arrival data with the data similarity higher than a preset threshold value as target first arrival data.
Step S102: and determining a long wavelength correction amount of the first arrival data according to the target first arrival data and a relative reference time, wherein the relative reference time is determined according to the horizon data.
As can be seen from the above description, the electronic device provided by the embodiment of the present application is capable of determining the target first arrival data with long wavelength problem by performing cross-correlation analysis on the first arrival data in the seismic data of the target work area and the horizon data of the preset mark layer, that is, setting the first arrival data with the data similarity higher than the preset threshold value obtained after the cross-correlation analysis as the target first arrival data, and then determining the long wavelength correction amount of the first arrival data according to the numerical difference between the target first arrival data and the relative reference time, so as to eliminate the coherence between the first arrival data and the horizon data, effectively reduce the error generated by the long wavelength problem, and improve the construction interpretation precision.
The embodiment of the present application also provides a computer-readable storage medium capable of implementing all the steps of the long-wavelength static correction method of seismic waves in the above embodiment, the computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements all the steps of the long-wavelength static correction method of seismic waves in the above embodiment, for example, the processor implements the following steps when executing the computer program:
step S101: and performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and setting the first arrival data with the data similarity higher than a preset threshold value as target first arrival data.
Step S102: and determining a long wavelength correction amount of the first arrival data according to the target first arrival data and a relative reference time, wherein the relative reference time is determined according to the horizon data.
As can be seen from the above description, the computer readable storage medium provided by the embodiment of the present application can determine the target first arrival data with long wavelength problem by performing cross-correlation analysis on the first arrival data in the seismic data of the target work area and the horizon data of the preset mark layer, that is, the first arrival data with the data similarity higher than the preset threshold value obtained by the cross-correlation analysis is set as the target first arrival data, and then determine the long wavelength correction amount of the first arrival data according to the numerical difference between the target first arrival data and the relative reference time, so as to eliminate the coherence between the first arrival data and the horizon data, effectively reduce the error generated by the long wavelength problem, and improve the construction interpretation precision.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for a hardware+program class embodiment, the description is relatively simple, as it is substantially similar to the method embodiment, as relevant see the partial description of the method embodiment.
The foregoing describes specific embodiments of the present disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims can be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.
Although the application provides method operational steps as described in the examples or flowcharts, more or fewer operational steps may be included based on conventional or non-inventive labor. The order of steps recited in the embodiments is merely one way of performing the order of steps and does not represent a unique order of execution. When implemented by an actual device or client product, the instructions may be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment) as shown in the embodiments or figures.
The system, apparatus, module or unit set forth in the above embodiments may be implemented in particular by a computer chip or entity, or by a product having a certain function. One typical implementation is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a car-mounted human-computer interaction device, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.
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 one typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of computer-readable media.
Computer readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media for a computer include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device. Computer-readable media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
It will be appreciated by those skilled in the art that embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, the present specification embodiments may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects.
The present embodiments may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The embodiments of the specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for system embodiments, since they are substantially similar to method embodiments, the description is relatively simple, as relevant to see a section of the description of method embodiments. In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," 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 embodiments of the present specification. In this specification, schematic representations of the above terms are not necessarily directed 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. Furthermore, the different embodiments or examples described in this specification and the features of the different embodiments or examples may be combined and combined by those skilled in the art without contradiction.
The foregoing is merely an example of the present specification and is not intended to limit the present specification. Various modifications and variations of the illustrative embodiments will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, or the like, which is within the spirit and principles of the embodiments of the present specification, should be included in the scope of the claims of the embodiments of the present specification.
Claims (6)
1. A method for long wavelength static correction of seismic waves, the method comprising:
performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and setting the first arrival data with the data similarity higher than a preset threshold value as target first arrival data; the first arrival data and the horizon data are time data; a phase shaft with high relative signal-to-noise ratio, good continuity and small time variation is selected in advance in a target work area to serve as a preset mark layer;
determining a long wavelength correction amount of the first arrival data according to the target first arrival data and a relative reference time, wherein the relative reference time is determined according to the horizon data; the relative reference time is the time corresponding to the existing reference plane, after static correction processing, the seismic data is corrected to an absolute elevation plane, the elevation plane is called a reference plane, and the starting time of the seismic data is changed from 0 before static correction to the time corresponding to the reference plane, and is called relative reference time;
performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain corresponding data similarity, wherein the cross-correlation analysis comprises the following steps:
dividing the first arrival data into a plurality of groups according to a preset offset interval; adding the first arrival data corresponding to the offset range in the same group to obtain the first arrival data of the group;
performing cross-correlation analysis on each group of first arrival data and the horizon data respectively to obtain corresponding data similarity;
the setting the first arrival data with the data similarity higher than a preset threshold as target first arrival data comprises the following steps:
normalizing each group of first arrival data according to the maximum value and the minimum value of each group of first arrival data with the data similarity higher than a preset threshold value and the maximum value and the minimum value of the horizon data;
adding the normalized first arrival data of each group to obtain target first arrival data;
before the determining of the long wavelength correction amount of the first arrival data according to the target first arrival data and the relative reference time, the method comprises:
determining the relative reference time according to any one of an average value of the horizon data, an average value of a maximum value and a minimum value of the horizon data and any value of the horizon data between the maximum value and the minimum value;
the determining the long wavelength correction amount of the first arrival data according to the target first arrival data and the relative reference time comprises the following steps:
a difference between the target first arrival data and the relative reference time is set as a long wavelength correction amount of the first arrival data.
2. The method of long wavelength static correction of seismic waves according to claim 1, further comprising, after said determining a long wavelength correction amount of said first arrival data based on said target first arrival data and a relative reference time:
determining shot point correction amount and/or wave detection point correction amount of the corresponding target work area according to the long wavelength correction amount of the first arrival data;
and correcting the seismic data by applying the shot correction amount and/or the wave detection point correction amount to generate corrected seismic data.
3. A long wavelength static correction device for seismic waves, comprising:
the high coherence determining module is used for performing cross-correlation analysis on first arrival data in the seismic data of the target work area and horizon data of a preset mark layer to obtain data similarity, and setting the first arrival data with the data similarity higher than a preset threshold value as target first arrival data; the first arrival data and the horizon data are time data; a phase shaft with high relative signal-to-noise ratio, good continuity and small time variation is selected in advance in a target work area to serve as a preset mark layer;
the high coherence elimination module is used for determining a long wavelength correction amount of the first arrival data according to the target first arrival data and relative reference time, wherein the relative reference time is determined according to the horizon data; the relative reference time is the time corresponding to the existing reference plane, after static correction processing, the seismic data is corrected to an absolute elevation plane, the elevation plane is called a reference plane, and the starting time of the seismic data is changed from 0 before static correction to the time corresponding to the reference plane, and is called relative reference time;
the high coherence determination module includes:
the first arrival data grouping unit is used for dividing the first arrival data into a plurality of groups according to a preset offset interval; adding the first arrival data corresponding to the offset range in the same group to obtain the first arrival data of the group;
each group of sequential analysis units is used for performing cross-correlation analysis on each group of first arrival data and the horizon data respectively to obtain corresponding data similarity;
each group of normalization processing units is used for performing normalization processing on each group of first arrival data according to the maximum value and the minimum value of each group of first arrival data with the data similarity higher than a preset threshold value and the maximum value and the minimum value of the horizon data;
each group of processed summarizing units is used for adding the normalized first arrival data of each group to obtain the target first arrival data;
the static correction device for the long wavelength of the earthquake wave further comprises:
a relative reference time determining unit configured to determine the relative reference time according to any one of an average value of the horizon data, an average value of a maximum value and a minimum value of the horizon data, and an arbitrary value of the horizon data between the maximum value and the minimum value;
the high coherence elimination module includes:
and a correction amount determination unit configured to set a difference between the target first-arrival data and the relative reference time as a long-wavelength correction amount of the first-arrival data.
4. A seismic wave long wavelength static correction apparatus according to claim 3, further comprising:
a long wavelength correction amount decomposition unit for determining a shot correction amount and/or a detection point correction amount of the corresponding target work area according to the long wavelength correction amount of the first arrival data;
and the seismic data correction unit is used for correcting the seismic data by applying the shot correction amount and/or the geophone correction amount to generate corrected seismic data.
5. 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 for long wavelength static correction of seismic waves of any one of claims 1 to 2 when said program is executed by said processor.
6. A computer readable storage medium having stored thereon a computer program, characterized in that the computer program, when being executed by a processor, implements the steps of the seismic long wavelength static correction method according to any of claims 1 to 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911247971.5A CN113031071B (en) | 2019-12-09 | 2019-12-09 | Static correction method and device for long wavelength of earthquake wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911247971.5A CN113031071B (en) | 2019-12-09 | 2019-12-09 | Static correction method and device for long wavelength of earthquake wave |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113031071A CN113031071A (en) | 2021-06-25 |
| CN113031071B true CN113031071B (en) | 2023-09-26 |
Family
ID=76450841
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911247971.5A Active CN113031071B (en) | 2019-12-09 | 2019-12-09 | Static correction method and device for long wavelength of earthquake wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113031071B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103116184A (en) * | 2011-11-16 | 2013-05-22 | 中国石油集团东方地球物理勘探有限责任公司 | Method for calculating transverse wave statics correction value in surface survey |
| CN104181599A (en) * | 2014-08-19 | 2014-12-03 | 中国石油天然气集团公司 | Refracted wave static correction handling method and system based on near-surface layer |
| CN104914469A (en) * | 2014-03-13 | 2015-09-16 | 中国石油化工股份有限公司 | Static correction method and device for converted S-wave |
| CN105319595A (en) * | 2014-07-31 | 2016-02-10 | 中国石油化工股份有限公司 | Long wavelength static correction method in loess plateau seismic data well control |
| CN105911592A (en) * | 2016-06-30 | 2016-08-31 | 中国石油化工股份有限公司 | Practical drilling constraint three-dimensional seismic data long wavelength static correction method |
| CN108279436A (en) * | 2017-12-20 | 2018-07-13 | 中国石油天然气股份有限公司 | Method and device for picking up seismic wave first arrivals |
-
2019
- 2019-12-09 CN CN201911247971.5A patent/CN113031071B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103116184A (en) * | 2011-11-16 | 2013-05-22 | 中国石油集团东方地球物理勘探有限责任公司 | Method for calculating transverse wave statics correction value in surface survey |
| CN104914469A (en) * | 2014-03-13 | 2015-09-16 | 中国石油化工股份有限公司 | Static correction method and device for converted S-wave |
| CN105319595A (en) * | 2014-07-31 | 2016-02-10 | 中国石油化工股份有限公司 | Long wavelength static correction method in loess plateau seismic data well control |
| CN104181599A (en) * | 2014-08-19 | 2014-12-03 | 中国石油天然气集团公司 | Refracted wave static correction handling method and system based on near-surface layer |
| CN105911592A (en) * | 2016-06-30 | 2016-08-31 | 中国石油化工股份有限公司 | Practical drilling constraint three-dimensional seismic data long wavelength static correction method |
| CN108279436A (en) * | 2017-12-20 | 2018-07-13 | 中国石油天然气股份有限公司 | Method and device for picking up seismic wave first arrivals |
Non-Patent Citations (1)
| Title |
|---|
| 长波长静校正问题的识别与解决方法--以鄂尔多斯黄土塬地震资料为例;杨城增 等;地球物理学进展;第31卷(第5期);2212-2218 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113031071A (en) | 2021-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210373187A1 (en) | Focus-stacking imaging method and system based on correlation-based seismic interferometry | |
| CN107329174B (en) | The method and device of well logging information is shown on seismic profile | |
| CN109696705A (en) | Deep-sea OBS uplink and downlink Seismic imaging data correcting method and device | |
| CN111781635B (en) | Seabed four-component elastic wave Gaussian beam depth migration method and device | |
| CN108445532A (en) | A kind of Depth Domain inverse migration method and device | |
| CN106054252B (en) | A kind of method and device of pre-stack time migration | |
| CN111239832B (en) | 3D3C-VSP imaging processing method, device and equipment | |
| CN113031071B (en) | Static correction method and device for long wavelength of earthquake wave | |
| CN112379436A (en) | Method and device for calculating static correction value of seismic data excited by mixed seismic source | |
| CN112764102B (en) | Seismic data acquisition method and device | |
| Midzi et al. | 1-D velocity model for use by the SANSN in earthquake location | |
| CN113009569B (en) | Earthquake migration imaging method and device | |
| CN109471173B (en) | Residual static correction method and device | |
| CN112649875B (en) | Submarine node detector positioning method and device | |
| CN112083490B (en) | Seismic data noise attenuation method and device | |
| CN112255679A (en) | Seismic data diffraction depth migration processing method and device | |
| CN111999767B (en) | Offset imaging method and device for undulating surface | |
| CN111948714B (en) | Noise reduction method and noise reduction device for uplink wave in submarine node seismic data | |
| CN111337973A (en) | Seismic data reconstruction method and system | |
| CN111435172A (en) | Method and device for chromatographic static correction | |
| CN113093279B (en) | Converted wave static correction method and device | |
| CN113568047B (en) | Pre-stack imaging gather generation method and device | |
| CN112558157B (en) | Static correction method and device based on multi-wave combination | |
| CN114063165B (en) | Three-dimensional seismic data splicing method and device | |
| CN112946745B (en) | Converted wave residual static correction method and system based on gravity fitting |
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 |