CN112130198B - Multi-period three-dimensional seismic data fusion method and device based on time difference correction - Google Patents
Multi-period three-dimensional seismic data fusion method and device based on time difference correction Download PDFInfo
- Publication number
- CN112130198B CN112130198B CN201910547387.5A CN201910547387A CN112130198B CN 112130198 B CN112130198 B CN 112130198B CN 201910547387 A CN201910547387 A CN 201910547387A CN 112130198 B CN112130198 B CN 112130198B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- dimensional seismic
- correction
- time difference
- result
- 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 148
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 16
- 238000011835 investigation Methods 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000004927 fusion Effects 0.000 claims description 15
- 230000003068 static effect Effects 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 19
- 238000007499 fusion processing Methods 0.000 description 11
- 230000008901 benefit Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000003860 storage Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011173 large scale experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002344 surface layer Substances 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
-
- 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/53—Statics correction, e.g. weathering layer or transformation to a datum
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 multi-period three-dimensional seismic data fusion method and device based on time difference correction, wherein the method comprises the following steps: acquiring multi-period three-dimensional seismic data; respectively superposing a plurality of historical three-dimensional seismic data and target acquisition three-dimensional seismic data at the same speed to obtain a plurality of superposition results; the plurality of superposition results are in one-to-one correspondence with the plurality of historical acquisition three-dimensional seismic data; performing cross-correlation gridding time difference investigation on the superposition result to obtain an investigation result; each investigation result comprises a plurality of cross-correlation time difference correction amounts; correcting the historical acquired three-dimensional seismic data according to the plurality of cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data; and fusing the plurality of corrected three-dimensional seismic data to obtain fused three-dimensional seismic data. According to the method, the fused three-dimensional seismic data with higher quality is obtained by improving the precision of time difference correction of the multi-period three-dimensional seismic data.
Description
Technical Field
The invention relates to the technical field of geophysical exploration, in particular to a multi-period three-dimensional seismic data fusion method and device based on time difference correction.
Background
In recent years, three-dimensional seismic data has gained tremendous benefits in oil and gas exploration and development, and three-dimensional seismic data processing technology has been rapidly developed and perfected. With further depth of oil field exploration and development, the defects of single small three-dimensional seismic data are gradually displayed, and the overall evaluation of resources is affected. The foreign related documents and a plurality of domestic large-scale experiments in recent years show that the fusion treatment of the seismic data is an important way for effectively solving the problems on the premise of saving a large amount of acquisition funds.
The fusion processing is to take the omnibearing high-density uniform sampling as a core, utilize the concept of introducing time scale by multi-period acquisition, decompose the high-density acquisition of a single time scale into the high-density acquisition of a plurality of time scales, and finally form a set of omnibearing high-density uniform sampling data body through the high-density acquisition of the plurality of time scales. The multi-period collection refers to collection of data of different years in the same area. Because of the stages of exploration targets and deployment, in the past decades, some old oil fields and old detection areas accumulate monolithic three-dimensional data with different quality of acquisition parameters in different years. The key point of the data fusion processing is to solve the problem of data consistency, because the surface layer structure is continuously changed along with the time, such as the depth change of the diving surface, three-dimensional data volumes acquired in different periods often have time differences, and the investigation of the time differences among blocks in the fusion processing is a very important step.
The current common method for time difference investigation is to survey time differences through overlapping cross lines of different acquisition years, for example, the time difference of a main line direction is 4ms, the time difference of a line connection direction is 6ms, and then the time difference of the cross crossing point is averaged for 5ms. The existing time difference investigation method has lower precision of time difference correction, and based on the time difference investigation result with lower precision, the fusion quality of the three-dimensional seismic data cannot be ensured.
Disclosure of Invention
The invention provides a multi-period three-dimensional seismic data fusion method and device based on time difference correction, which can improve the precision of time difference correction of multi-period three-dimensional seismic data and relieve the problem of inconsistent time difference during three-dimensional seismic data fusion processing, thereby obtaining fused three-dimensional seismic data with higher quality.
In a first aspect, an embodiment of the present invention provides a method for fusion of multi-period three-dimensional seismic data based on time difference correction, where the method includes: acquiring multi-period three-dimensional seismic data; the multi-period three-dimensional seismic data comprises target acquisition three-dimensional seismic data and a plurality of historical acquisition three-dimensional seismic data; superposing the plurality of historical three-dimensional seismic data and the target acquisition three-dimensional seismic data at the same speed to obtain a plurality of superposition results; the plurality of superposition results are in one-to-one correspondence with the plurality of historical acquisition three-dimensional seismic data; performing cross-correlation gridding time difference investigation on the superposition result to obtain an investigation result; each survey result comprises a plurality of cross-correlation time difference correction amounts; correcting the historical acquired three-dimensional seismic data according to the cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data; and fusing the plurality of corrected three-dimensional seismic data to obtain fused three-dimensional seismic data.
In a second aspect, an embodiment of the present invention further provides a multi-period three-dimensional seismic data fusion apparatus based on time difference correction, where the apparatus includes: the acquisition module is used for acquiring multi-period three-dimensional seismic data; the multi-period three-dimensional seismic data comprises target acquisition three-dimensional seismic data and a plurality of historical acquisition three-dimensional seismic data; the superposition module is used for respectively superposing the plurality of historical three-dimensional seismic data and the target acquisition three-dimensional seismic data at the same speed to obtain a plurality of superposition results; the plurality of superposition results are in one-to-one correspondence with the plurality of historical acquisition three-dimensional seismic data; the investigation module is used for carrying out cross-correlation gridding time difference investigation on the superposition results to obtain investigation results; each survey result comprises a plurality of cross-correlation time difference correction amounts; the correction module corrects the historical acquired three-dimensional seismic data according to the plurality of cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data; and the fusion module is used for fusing the plurality of corrected three-dimensional seismic data to obtain fused three-dimensional seismic data.
In a third aspect, an embodiment of the present invention further provides a computer device, including a memory, and a processor, where the memory stores a computer program that can run on the processor, and when the processor executes the computer program, the method for implementing the multi-period three-dimensional seismic data fusion based on time difference correction is implemented.
In a fourth aspect, embodiments of the present invention also provide a computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the above-described time difference correction based multi-term three-dimensional seismic data fusion method.
The embodiment of the invention has the following beneficial effects: the embodiment of the invention provides a multi-period three-dimensional seismic data fusion method and device based on time difference correction. According to the embodiment of the invention, the time difference correction precision of the multi-period three-dimensional seismic data is improved, the problem of inconsistent time difference during the fusion processing of the three-dimensional seismic data is solved, and the fused three-dimensional seismic data with higher quality can be obtained.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the above objects, features and advantages of the present invention more comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the present invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a multi-period three-dimensional seismic data fusion method based on time difference correction according to an embodiment of the invention;
FIG. 2 is a schematic diagram of an implementation flow of generating corrected three-dimensional seismic data of different periods by using a multi-period three-dimensional seismic data fusion method based on time difference correction according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a multi-period acquisition three-dimensional data fusion process according to an embodiment of the present invention;
fig. 4 is a schematic diagram of a meshed cross-correlation time difference investigation result provided by an embodiment of the present invention;
FIG. 5 is a cross-sectional view of an interactive display of a three-dimensional seismic data stack prior to moveout correction in accordance with an embodiment of the present invention;
FIG. 6 is a cross-sectional view of an interactive display of a three-dimensional seismic data stack after moveout correction according to an embodiment of the present invention;
FIG. 7 is a block diagram of a multi-period three-dimensional seismic data fusion apparatus based on time difference correction according to an embodiment of the present invention;
fig. 8 is a block diagram of a computer device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
At present, the conventional time difference investigation method has low precision and accuracy, cannot accurately obtain specific time difference values of each point, and in addition, the condition of inconsistent time difference is easy to occur when three-dimensional data are acquired for multiple times in fusion processing, and the quality of a three-dimensional seismic data fusion result can be influenced due to inaccurate time difference solving among seismic data of different acquisition years.
Based on the above, the multi-period three-dimensional seismic data fusion method and device based on time difference correction provided by the embodiment of the invention can improve the precision of time difference investigation, eliminate the problem of inconsistent time difference during multi-period three-dimensional data acquisition fusion processing, and lay a good data foundation for subsequent consistency processing.
For the convenience of understanding the present embodiment, the method for fusion of multi-period three-dimensional seismic data based on time difference correction disclosed in the present embodiment will be described in detail.
The embodiment of the invention provides a multi-period three-dimensional seismic data fusion method based on time difference correction, which is shown in a flow chart of the multi-period three-dimensional seismic data fusion method based on time difference correction in FIG. 1, and comprises the following steps:
Step S102, multi-period three-dimensional seismic data are acquired.
In an embodiment of the invention, the multi-period three-dimensional seismic data includes a target acquisition three-dimensional seismic data and a plurality of historical acquisition three-dimensional seismic data. The target three-dimensional seismic data is the latest three-dimensional seismic data of the target layer, the historical three-dimensional seismic data is the three-dimensional seismic data acquired by the target layer in a historical period, the historical three-dimensional seismic data can comprise a plurality of data in different periods, the specific number and the using period can be set according to actual requirements, and the embodiment of the invention is not particularly limited.
Step S104, respectively superposing a plurality of historical three-dimensional seismic data and the target acquisition three-dimensional seismic data at the same speed to obtain a plurality of superposition results.
In the embodiment of the invention, for example, if there are N historical three-dimensional seismic data, each historical three-dimensional seismic data is overlapped with the target acquisition three-dimensional seismic data, so that N overlapped results can be obtained, and the N overlapped results are in one-to-one correspondence with the N historical three-dimensional seismic data. The same speed is used for data superposition, so that consistency of superposition results in different historical periods can be ensured.
And S106, performing cross-correlation gridding time difference investigation on the superposition result to obtain an investigation result.
In the embodiment of the invention, the execution process of the cross-correlation meshed time difference investigation can be that a plurality of lines are determined according to the main line direction, a plurality of investigation points are obtained according to the line connection direction, a cross-correlation time difference correction amount is generated for each investigation point, a investigation result can be obtained for each superposition result, and each investigation result comprises a plurality of cross-correlation time difference correction amounts. The cross-correlation time difference correction amount represents the difference between the investigation point in the historical three-dimensional seismic data corresponding to the investigation result and the investigation point in the target three-dimensional seismic data.
Referring to the grid-like cross-correlation moveout survey result diagram shown in fig. 4, the diagram shows a plurality of survey point cross-correlation moveout correction amounts in the form of a grid. The method comprises the steps of determining that 10 lines are different in each grid in the line connecting direction when the table transversely represents the line connecting direction and the longitudinal represents the main line direction according to the superposition result and generating the table in the graph, wherein each line is spaced by 20 meters, determining that 25 points are different in each grid in the main line direction, and each point is spaced by 20 meters. In practical implementation, the number of lines or points of the phase difference between the main line and the line connection direction per cell and the distance between the lines or points may be set according to the requirements, and the embodiment of the present invention is not limited specifically.
And S108, correcting the historical acquired three-dimensional seismic data according to the plurality of cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data.
In the embodiment of the invention, the investigation points in the corresponding historical acquisition three-dimensional seismic data are corrected according to the cross-correlation time difference correction amount, and as the cross-correlation time difference correction amount is multiple, the correction of multiple points in the historical acquisition three-dimensional seismic data can be realized, and the correction accuracy is improved.
Step S110, fusing the plurality of corrected three-dimensional seismic data to obtain fused three-dimensional seismic data.
In embodiments of the present invention, multiple corrected three-dimensional seismic data may be used to improve the resolution of the fused three-dimensional seismic data. Referring to the schematic diagram of multi-period acquisition three-dimensional data fusion processing shown in fig. 3, three-dimensional seismic data acquired in a plurality of historical periods are fused, namely, a plurality of processes for correcting the three-dimensional seismic data and enriching the three-dimensional data body are used, and finally, fused three-dimensional seismic data with higher density can be obtained, so that the quality of the fused three-dimensional seismic data is improved.
The embodiment of the invention provides a multi-period three-dimensional seismic data fusion method based on time difference correction, which comprises the steps of firstly obtaining multi-period three-dimensional seismic data comprising target acquisition three-dimensional seismic data and a plurality of historical acquisition three-dimensional seismic data, then respectively superposing the plurality of historical three-dimensional seismic data and the target acquisition three-dimensional seismic data at the same speed to obtain a plurality of superposition results, carrying out cross-correlation gridding time difference investigation on the plurality of superposition results to obtain a plurality of investigation results, correcting the corresponding historical acquisition three-dimensional seismic data according to the plurality of cross-correlation time difference correction amounts, thereby correcting a plurality of points in the historical acquisition three-dimensional seismic data through the plurality of cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data with higher correction accuracy, and finally fusing the plurality of corrected three-dimensional seismic data to obtain fused three-dimensional seismic data. According to the embodiment of the invention, the time difference correction precision of the multi-period three-dimensional seismic data is improved, the problem of inconsistent time difference during the fusion processing of the three-dimensional seismic data is solved, and the fused three-dimensional seismic data with higher quality can be obtained.
In order to further improve correction accuracy, the correction is performed on the corresponding historical acquired three-dimensional seismic data according to a plurality of cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data, and the method specifically comprises the following steps:
(1) And performing interpolation processing on the cross-correlation time difference correction amount corresponding to the historical acquired three-dimensional seismic data to obtain a three-dimensional static correction amount library.
Each history acquisition three-dimensional seismic data corresponds to a plurality of cross-correlation time difference correction amounts, interpolation processing is respectively carried out on the cross-correlation time difference correction amounts of each history acquisition three-dimensional seismic data, a plurality of three-dimensional static correction amount libraries can be obtained, and the three-dimensional static correction amount libraries corresponding to each history acquisition three-dimensional seismic data together form an integral static correction amount library.
The kriging method is a random interpolation technology based on a general least square algorithm, and uses a variance diagram as a weight function; this technique can be applied to any phenomenon that requires the use of point data to estimate its distribution over the surface. The embodiment of the invention can use the kriging interpolation for processing.
(2) And correcting the historical acquisition three-dimensional seismic data according to the three-dimensional static correction amount library to obtain a correction result.
In the embodiment of the invention, the process of correcting the historical acquisition three-dimensional seismic data, namely the process of increasing or decreasing the numerical value of each check point according to the three-dimensional static correction amount library, can relieve the problem of data consistency caused by the influence of time difference on the surface structure through correction.
(3) And if the correction result meets the correction standard, taking the correction result as the corrected three-dimensional seismic data.
After the correction result is obtained, it is necessary to determine whether the correction result meets a correction standard, where the correction standard may be set according to requirements, which is not specifically limited in the embodiment of the present invention.
And judging whether the correction result meets the correction standard, overlapping the correction result with the target acquired three-dimensional seismic data to obtain a correction overlapping result, judging whether the overlapping result meets the correction standard, and taking the correction result as the correction three-dimensional seismic data if various parameters such as the accuracy of the overlapping result meet the correction standard.
It should be noted that, in the embodiment of the present invention, each correction result is obtained according to a three-dimensional static correction amount library, each three-dimensional static correction amount library corresponds to a survey result, each survey result is obtained by a superposition result, and the superposition result is obtained by superposing a historical acquisition three-dimensional seismic data with the acquisition three-dimensional seismic data, so that a plurality of correction three-dimensional seismic data corresponding to the historical acquisition three-dimensional seismic data one by one can be obtained.
The embodiment of the invention provides a multi-period three-dimensional seismic data fusion method based on time difference correction, which is shown in a cross section of a three-dimensional seismic data superposition section before time difference correction in the figure 5, wherein the cross section represents a line connection direction transversely and represents time longitudinally, the three-dimensional seismic data superposition effect of the three-dimensional seismic data and a target acquisition three-dimensional seismic data is poor before time difference correction is carried out, the cross section of the three-dimensional seismic data superposition section after time difference correction in the figure 6 is shown in the figure, the cross section represents the line connection direction transversely and represents time longitudinally, and the three-dimensional seismic data superposition effect of the three-dimensional seismic data and the target acquisition three-dimensional seismic data is good after time difference correction is carried out. Referring to the implementation flow diagram of the time difference correction-based multi-period three-dimensional seismic data fusion method shown in fig. 2, the implementation flow diagram of the correction three-dimensional seismic data of different periods is generated, and the accuracy of time difference calculation among seismic data of different acquisition years is improved, so that consistency of data time differences in fusion processing is ensured, a good data foundation is laid for subsequent consistency processing, and better three-dimensional seismic data can be obtained through fusion.
The embodiment of the invention also provides a multi-period three-dimensional seismic data fusion device based on time difference correction, referring to a structural block diagram of the multi-period three-dimensional seismic data fusion device based on time difference correction shown in fig. 7, the device comprises:
An acquisition module 71 for acquiring multi-period three-dimensional seismic data; the multi-period three-dimensional seismic data comprises target acquisition three-dimensional seismic data and a plurality of historical acquisition three-dimensional seismic data; the superposition module 72 is configured to respectively superimpose the plurality of historical three-dimensional seismic data and the target acquisition three-dimensional seismic data at the same speed to obtain a plurality of superimposed results; the plurality of superposition results are in one-to-one correspondence with the plurality of historical acquisition three-dimensional seismic data; a survey module 73, configured to perform cross-correlation gridding time difference survey on the plurality of superimposed results, so as to obtain a plurality of survey results; each investigation result comprises a plurality of cross-correlation time difference correction amounts; the correction module 74 corrects the corresponding historical acquired three-dimensional seismic data according to the plurality of cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data; the fusion module 75 is configured to fuse the plurality of corrected three-dimensional seismic data to obtain fused three-dimensional seismic data.
The correction module is specifically used for: interpolation processing is carried out on the cross-correlation time difference correction quantity corresponding to the historical acquired three-dimensional seismic data, so that a three-dimensional static correction quantity library is obtained; correcting the historical acquisition three-dimensional seismic data according to the three-dimensional static correction amount library to obtain a correction result; and if the correction result meets the correction standard, taking the correction result as the corrected three-dimensional seismic data.
The correction module is specifically used for: and performing Kriging interpolation processing on the cross-correlation time difference correction corresponding to the historical acquired three-dimensional seismic data to obtain a three-dimensional static correction library.
The correction module is specifically used for: superposing the correction result and the target acquisition three-dimensional seismic data to obtain a correction superposition result; judging whether the correction superposition result meets the correction standard; if so, the correction result is used as the correction three-dimensional seismic data.
The embodiment of the present invention further provides a computer device, referring to a schematic block diagram of a computer device structure shown in fig. 8, where the computer device includes a memory 81 and a processor 82, and the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, the processor implements steps of any of the methods described above.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described computer device may refer to corresponding procedures in the foregoing method embodiments, which are not repeated here
Embodiments of the present invention also provide a computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the steps of any of the methods described above.
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 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.
Finally, it should be noted that: the above examples are only specific embodiments of the present invention, and are not intended to limit the scope of the present invention, but it should be understood by those skilled in the art that the present invention is not limited thereto, and that the present invention is described in detail with reference to the foregoing examples: any person skilled in the art may modify or easily conceive of the technical solution described in the foregoing embodiments, or perform equivalent substitution of some of the technical features, while remaining within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (6)
1. A multi-period three-dimensional seismic data fusion method based on time difference correction is characterized by comprising the following steps:
Acquiring multi-period three-dimensional seismic data; the multi-period three-dimensional seismic data comprises target acquisition three-dimensional seismic data and a plurality of historical acquisition three-dimensional seismic data;
Superposing the plurality of historical three-dimensional seismic data and the target acquisition three-dimensional seismic data at the same speed to obtain a plurality of superposition results; the plurality of superposition results are in one-to-one correspondence with the plurality of historical acquisition three-dimensional seismic data;
Performing cross-correlation gridding time difference investigation on the superposition result to obtain an investigation result; each survey result comprises a plurality of cross-correlation time difference correction amounts; the execution process of the cross-correlation meshed time difference investigation is to determine a plurality of lines according to the main line direction, determine a plurality of lines according to the line connection direction, obtain a plurality of investigation points, generate a cross-correlation time difference correction amount for each investigation point, obtain a investigation result from each superposition result, and the cross-correlation time difference correction amount represents the difference between the investigation point in the historical acquisition three-dimensional seismic data corresponding to the investigation result and the investigation point in the target acquisition three-dimensional seismic data;
Correcting the historical acquired three-dimensional seismic data according to the cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data;
fusing the corrected three-dimensional seismic data to obtain fused three-dimensional seismic data;
correcting the historical acquired three-dimensional seismic data according to the plurality of cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data, including:
Performing Kriging interpolation processing on the cross-correlation time difference correction amount corresponding to the historical acquired three-dimensional seismic data to obtain a three-dimensional static correction amount library;
Correcting the historical acquired three-dimensional seismic data according to the three-dimensional static correction amount library to obtain a correction result;
and if the correction result meets the correction standard, taking the correction result as corrected three-dimensional seismic data.
2. The method of claim 1, wherein if the correction results meet a correction criterion, using the correction results as corrected three-dimensional seismic data comprises:
superposing the correction result and the target acquisition three-dimensional seismic data to obtain a correction superposition result;
Judging whether the correction superposition result meets a correction standard or not;
And if so, taking the correction result as corrected three-dimensional seismic data.
3. A multi-period three-dimensional seismic data fusion apparatus based on time difference correction, comprising:
The acquisition module is used for acquiring multi-period three-dimensional seismic data; the multi-period three-dimensional seismic data comprises target acquisition three-dimensional seismic data and a plurality of historical acquisition three-dimensional seismic data;
The superposition module is used for respectively superposing the plurality of historical three-dimensional seismic data and the target acquisition three-dimensional seismic data at the same speed to obtain a plurality of superposition results; the plurality of superposition results are in one-to-one correspondence with the plurality of historical acquisition three-dimensional seismic data;
The investigation module is used for carrying out cross-correlation gridding time difference investigation on the superposition results to obtain investigation results; each survey result comprises a plurality of cross-correlation time difference correction amounts; the execution process of the cross-correlation meshed time difference investigation is to determine a plurality of lines according to the main line direction, determine a plurality of lines according to the line connection direction, obtain a plurality of investigation points, generate a cross-correlation time difference correction amount for each investigation point, obtain a investigation result from each superposition result, and the cross-correlation time difference correction amount represents the difference between the investigation point in the historical acquisition three-dimensional seismic data corresponding to the investigation result and the investigation point in the target acquisition three-dimensional seismic data;
The correction module corrects the historical acquired three-dimensional seismic data according to the plurality of cross-correlation time difference correction amounts to obtain corrected three-dimensional seismic data;
the fusion module is used for fusing the plurality of corrected three-dimensional seismic data to obtain fused three-dimensional seismic data;
The correction module is specifically configured to:
Performing Kriging interpolation processing on the cross-correlation time difference correction amount corresponding to the historical acquired three-dimensional seismic data to obtain a three-dimensional static correction amount library;
Correcting the historical acquired three-dimensional seismic data according to the three-dimensional static correction amount library to obtain a correction result;
and if the correction result meets the correction standard, taking the correction result as corrected three-dimensional seismic data.
4. A device according to claim 3, characterized in that the correction module is specifically configured to:
superposing the correction result and the target acquisition three-dimensional seismic data to obtain a correction superposition result;
Judging whether the correction superposition result meets a correction standard or not;
And if so, taking the correction result as corrected three-dimensional seismic data.
5. A computer device comprising a memory, a processor, the memory having stored therein a computer program executable on the processor, characterized in that the processor, when executing the computer program, implements the steps of the method of any of the preceding claims 1-2.
6. A computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the method of any of the preceding claims 1-2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910547387.5A CN112130198B (en) | 2019-06-24 | 2019-06-24 | Multi-period three-dimensional seismic data fusion method and device based on time difference correction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910547387.5A CN112130198B (en) | 2019-06-24 | 2019-06-24 | Multi-period three-dimensional seismic data fusion method and device based on time difference correction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112130198A CN112130198A (en) | 2020-12-25 |
| CN112130198B true CN112130198B (en) | 2024-04-30 |
Family
ID=73849188
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910547387.5A Active CN112130198B (en) | 2019-06-24 | 2019-06-24 | Multi-period three-dimensional seismic data fusion method and device based on time difference correction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112130198B (en) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4759636A (en) * | 1985-12-16 | 1988-07-26 | Amoco Corporation | Method and system for real-time processing of seismic data |
| US5946271A (en) * | 1997-03-21 | 1999-08-31 | Western Atlas International, Inc. | Calibration system for use in time lapse tomography |
| US6289285B1 (en) * | 1999-12-16 | 2001-09-11 | Phillips Petroleum Company | Seismic velocity analysis for class II sands |
| WO2006103929A1 (en) * | 2005-03-29 | 2006-10-05 | National University Corporation The University Of Electro-Communications | Earthquake-prediction system, method and program, and medium recording earthquake prediction program |
| CN101031821A (en) * | 2004-07-23 | 2007-09-05 | 输入输出公司 | Offeset-azimuth binning for migration and velocity analysis |
| CN102707316A (en) * | 2012-06-06 | 2012-10-03 | 中国石油天然气股份有限公司 | 2.5-times three-dimensional seismic exploration method of complex geological target |
| CN106842322A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | A kind of carbon dioxide drive monitors seismic time difference bearing calibration |
| CN107797140A (en) * | 2016-09-07 | 2018-03-13 | 中国石油化工股份有限公司 | Beading seismic response features focus on optimization method and system |
| FR3059107A1 (en) * | 2016-11-18 | 2018-05-25 | Landmark Graphics Corporation | AUTOMATED ANALYSIS OF DIFFERENCES AND CORRECTION BETWEEN SEISMIC BI-DIMENSIONAL (2D) REPORTS |
| CN108196305A (en) * | 2018-01-17 | 2018-06-22 | 东华理工大学 | A kind of mountainous region static correcting method |
| CN108897036A (en) * | 2018-05-15 | 2018-11-27 | 中国石油天然气集团有限公司 | A kind of seismic data processing technique and device |
| CN109884692A (en) * | 2019-04-10 | 2019-06-14 | 中煤科工集团西安研究院有限公司 | Time lapse seismic data processing method and system in a kind of detection of Coal Mine Goaf |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8462583B2 (en) * | 2008-05-29 | 2013-06-11 | Woodside Energy Ltd. | Method of marine seismic data acquisition |
-
2019
- 2019-06-24 CN CN201910547387.5A patent/CN112130198B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4759636A (en) * | 1985-12-16 | 1988-07-26 | Amoco Corporation | Method and system for real-time processing of seismic data |
| US5946271A (en) * | 1997-03-21 | 1999-08-31 | Western Atlas International, Inc. | Calibration system for use in time lapse tomography |
| US6289285B1 (en) * | 1999-12-16 | 2001-09-11 | Phillips Petroleum Company | Seismic velocity analysis for class II sands |
| CN101031821A (en) * | 2004-07-23 | 2007-09-05 | 输入输出公司 | Offeset-azimuth binning for migration and velocity analysis |
| WO2006103929A1 (en) * | 2005-03-29 | 2006-10-05 | National University Corporation The University Of Electro-Communications | Earthquake-prediction system, method and program, and medium recording earthquake prediction program |
| CN102707316A (en) * | 2012-06-06 | 2012-10-03 | 中国石油天然气股份有限公司 | 2.5-times three-dimensional seismic exploration method of complex geological target |
| CN106842322A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | A kind of carbon dioxide drive monitors seismic time difference bearing calibration |
| CN107797140A (en) * | 2016-09-07 | 2018-03-13 | 中国石油化工股份有限公司 | Beading seismic response features focus on optimization method and system |
| FR3059107A1 (en) * | 2016-11-18 | 2018-05-25 | Landmark Graphics Corporation | AUTOMATED ANALYSIS OF DIFFERENCES AND CORRECTION BETWEEN SEISMIC BI-DIMENSIONAL (2D) REPORTS |
| CN108196305A (en) * | 2018-01-17 | 2018-06-22 | 东华理工大学 | A kind of mountainous region static correcting method |
| CN108897036A (en) * | 2018-05-15 | 2018-11-27 | 中国石油天然气集团有限公司 | A kind of seismic data processing technique and device |
| CN109884692A (en) * | 2019-04-10 | 2019-06-14 | 中煤科工集团西安研究院有限公司 | Time lapse seismic data processing method and system in a kind of detection of Coal Mine Goaf |
Non-Patent Citations (8)
| Title |
|---|
| Imaging a near-surface feature using cross-correlation analysis of multi-channel surface wave data;K. Samyn等;Near Surface Geophysics;20120901;第11卷(第01期);第1-10页 * |
| 三维地震勘探技术在多组系复杂构造区的应用;徐敏;刘克难;梅莉;;断块油气田;20110731;第18卷(第04期);第453-456页 * |
| 互相关时差分析技术及其应用――以胜利油田三维地震资料连片处理为例;李继光;耿林;顾庆雷;陈新荣;徐辉;;石油物探;第49卷(第01期);第24-27页 * |
| 克里金算法与多重分形理论在岩土参数随机场分析中的应用;王长虹等;岩土力学;20141031;第35卷(第S2期);第386-392页 * |
| 地震纹理属性在JJD工区断层识别中的应用;崔世凌;张军华;王伟;谭明友;陈海云;物探化探计算技术;20100531;第32卷(第03期);第304-309页 * |
| 地震资料的高保真融合处理技术在冀东油田南堡地区的应用;石文武;黄荣善;王冬娜;张志让;;海相油气地质;第20卷(第02期);第63-71页 * |
| 影响煤层气地区资料成像的因素分析;米信武;程展展;王伟;;中国石油和化工(第S1期);第283，285页 * |
| 非重复性时移地震数据关键处理技术应用研究;郭念民;尚新民;刘欣欣;王淑玲;;石油地球物理勘探;第46卷(第04期);第582-585页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112130198A (en) | 2020-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113740901B (en) | Land seismic data full-waveform inversion method and device based on complex undulating surface | |
| CN103941286A (en) | Rapid three-dimensional horizon interpretation method based on horizontal navigation | |
| CN112464520B (en) | Local gravity anomaly depth inversion method and device | |
| CN106959467B (en) | Seimic wave velocity inversion method and device | |
| CN111596347B (en) | Method and device for rapidly obtaining surface layer longitudinal and transverse wave speed ratio | |
| CN103713322A (en) | Orientation higher-order coherent amplitude change-rate crack detection method | |
| CN104199103A (en) | Static correction method and static correction device for converted wave | |
| CN112130198B (en) | Multi-period three-dimensional seismic data fusion method and device based on time difference correction | |
| CN108828669B (en) | Two-dimensional intersection survey line static correction processing method, device and system | |
| CN102968470B (en) | A kind of accurate calculating has the method for the gully head furrow bank metamorphosis of washing hole | |
| CN113568033A (en) | Design method and device of three-dimensional irregular sampling seismic acquisition observation system | |
| CN107870374B (en) | A kind of strata denudation evolution recovery method | |
| CN112419493B (en) | Shale reservoir three-dimensional attribute model building method and device | |
| CN112505780B (en) | Formation depth data correction method and device | |
| CN113267810B (en) | Seismic exploration full-depth velocity modeling method and device | |
| CN112231878B (en) | Interlayer influence downward permeability coarsening equivalent calculation method and device | |
| CN111045090B (en) | Magnetic anomaly grid rapid stitching method | |
| CN114910966B (en) | Method and device for tomographic inversion during azimuth travel | |
| CN104533405A (en) | Method and device for detecting total organic carbon of shale reservoir | |
| CN112649875A (en) | Method and device for positioning submarine node geophone | |
| CN113156494B (en) | Method and device for calculating equivalent coverage times of surface elements | |
| CN112464521B (en) | Walking and sliding fracture determination method and device | |
| CN112363221B (en) | WALKAWAY VSP line layout method and WALKAWAY VSP line layout device | |
| CN112666607B (en) | Method and device for inverting thickness distribution of yellow soil layer by gravity | |
| CN114076982B (en) | Karst paleo-landform restoration method and device based on waveform characteristic difference |
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 |