CN110426736A - A kind of acquisition methods and device of offset gather - Google Patents
A kind of acquisition methods and device of offset gather Download PDFInfo
- Publication number
- CN110426736A CN110426736A CN201910711930.0A CN201910711930A CN110426736A CN 110426736 A CN110426736 A CN 110426736A CN 201910711930 A CN201910711930 A CN 201910711930A CN 110426736 A CN110426736 A CN 110426736A
- Authority
- CN
- China
- Prior art keywords
- shot
- imaging results
- offset
- data
- moment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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. analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
Abstract
The embodiment of the present invention provides the acquisition methods and device of a kind of offset gather.Method includes: that each single-shot in single-shots multiple in area to be measured is based on wave equation migration algorithm, obtains corresponding first imaging results of single-shot data of each single-shot and corresponding second imaging results of modulation single-shot data of each single-shot；Wherein, the modulation single-shot data are obtained after being modulated according to offset distance attribute to the single-shot data；According to corresponding first imaging results of each single-shot and the second imaging results, the corresponding earth's surface offset gather of each single-shot is obtained；Summarize the corresponding earth's surface earth's surface offset gather of each single-shot, obtains the offset gather in the area to be measured.Method and device provided in an embodiment of the present invention, calculating 30% calculation amount of increase using GPU can be obtained the offset gather in area to be measured, computing resource be saved, so that computing resource utilization rate is higher.
Description
Technical field
The present invention relates to seismic exploration technique field more particularly to the acquisition methods and device of a kind of offset gather.
Background technique
Seismic exploration data depth imaging method includes Kirchhoff offset, one-way wave offset and reverse-time migration etc., one way
Wave offset and reverse-time migration belong to wave equation migration.In one timing of image taking speed, Kirchhoff offset method computational accuracy
Minimum, one-way wave and reverse-time migration method imaging effect are more preferable.Kirchhoff calculations of offset can be with output offset away from trace gather, partially
Moving can be used to carry out AVO analysis, speed update etc. away from trace gather, be depth migration method indispensability output data.But when construction is multiple
When miscellaneous area is imaged, the higher wave equation migration of precision need to be applied, but because its input data is single shot record, cannot be exported partially
It moves away from trace gather, updated which limits the information excavating of such data depth and speed etc..It is inclined that wave equation is calculated at present
The method away from trace gather of shifting is mainly that single-shot data are divided into different offset distance sections to carry out substep offset, this is equivalent to one group
Data are divided into several groups and calculate again, and high-intensitive calculating limits its application.
Summary of the invention
In view of the problems of the existing technology, the embodiment of the present invention provides the acquisition methods and dress of a kind of offset gather
It sets.
In a first aspect, the embodiment of the present invention provides a kind of acquisition methods of offset gather, comprising:
For each single-shot in multiple single-shots in area to be measured, it is based on wave equation migration algorithm, obtains each single-shot
Corresponding second imaging results of modulation single-shot data of corresponding first imaging results of single-shot data and each single-shot；Wherein,
The modulation single-shot data are obtained after being modulated according to offset distance attribute to the single-shot data；
According to corresponding first imaging results of each single-shot and the second imaging results, it is inclined to obtain the corresponding earth's surface of each single-shot
It moves away from trace gather；
Summarize the corresponding earth's surface offset gather of each single-shot, obtains the offset gather in the area to be measured.
Further, the wave equation migration algorithm is reverse-time migration algorithm；Correspondingly,
Based on wave equation migration algorithm, corresponding first imaging results of single-shot data of each single-shot are obtained, comprising:
Geophone station wave field of each single-shot at each moment and the shot point wave field at each moment are obtained, and as each single-shot
Single-shot data；
Corresponding first imaging results are obtained according to the single-shot data of each single-shot based on the reverse-time migration algorithm.
Further, it is based on the reverse-time migration algorithm, according to the single-shot data of each single-shot, obtains corresponding the first one-tenth
As the calculation formula of result I (x, s) are as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and n is positive integer, and u (x, s, t) is single-shot s in t
The geophone station wave field at quarter, s (x, s, t) are shot point wave field of the single-shot s in t moment.
Further, it is based on wave equation migration algorithm, the modulation single-shot data for obtaining each single-shot are the second one-tenth corresponding
As result, comprising:
The calculation formula of the first imaging results I (x, s) is deformed, deformed first imaging results I is obtained
The calculation formula of (x, s) is as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and r represents geophone station, and n is positive integer, s (x, s, t)
It is single-shot s in the shot point wave field of t moment, u (x, s, r, t) is single track geophone station wave field of the single-shot s in t moment；
According to offset distance attribute, the single track geophone station wave field to each single-shot at each moment is modulated, after obtaining modulation
Single track geophone station wave field；
By each single-shot in the modulated single track geophone station wave field at each moment and in the shot point wave field at each moment, as every
The modulation single-shot data of a single-shot；
Corresponding second imaging knot is obtained according to the modulation single-shot data of each single-shot based on the reverse-time migration algorithm
Fruit.
Further, it is based on the reverse-time migration algorithm, according to the modulation single-shot data of each single-shot, obtains corresponding the
Two imaging results IhThe calculation formula of (x, s) are as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and r represents geophone station, and n is positive integer, s (x, s, t)
It is single-shot s in the shot point wave field of t moment, u (x, s, r, t) is single track geophone station wave field of the single-shot s in t moment, and offset is single
Distance of the big gun s to geophone station r.
Further, according to corresponding first imaging results of each single-shot and the second imaging results, each single-shot pair is obtained
Answer earth's surface offset gather I (x, h | calculation formula s) are as follows:
I (x, h | s)=I (x, s) δ (h-hA(x,s))
Wherein,
Wherein, x represents imaging position, behalf single-shot, and I (x, s) is the first imaging results, Ih(x, s) is the second imaging knot
Fruit, h represent offset distance, hA(x, s) is that the second imaging results of single-shot s are corresponding relative to each imaging point of the first imaging results
Offset distance.
Further, summarize the corresponding earth's surface offset gather of each single-shot, obtain the offset gather in the area to be measured
Calculation formula are as follows:
Wherein, x represents imaging position, behalf single-shot, and I (x, h | it s) is the corresponding earth's surface offset gather of single-shot s.
Second aspect, the embodiment of the present invention provide a kind of acquisition device of offset gather, comprising:
Imaging results obtain module, for it is inclined to be based on wave equation for each single-shot in multiple single-shots in area to be measured
Algorithm is moved, the modulation single-shot data of corresponding first imaging results of single-shot data and each single-shot that obtain each single-shot are corresponding
The second imaging results；Wherein, the modulation single-shot data are after being modulated according to offset distance attribute to the single-shot data
It obtains；
Earth's surface offset gather obtains module, for being tied according to corresponding first imaging results of each single-shot and the second imaging
Fruit obtains the corresponding earth's surface offset gather of each single-shot；
Offset gather obtains module, for summarizing the corresponding earth's surface offset gather of each single-shot, obtains described to be measured
The offset gather in area
The third aspect, the embodiment of the present invention provides a kind of electronic equipment, including memory, processor and is stored in memory
Computer program that is upper and can running on a processor, is realized when the processor executes described program as first aspect provides
Method the step of.
Fourth aspect, the embodiment of the present invention provide a kind of non-transient computer readable storage medium, are stored thereon with calculating
Machine program is realized as provided by first aspect when the computer program is executed by processor the step of method.
The acquisition methods and device of a kind of offset gather provided in an embodiment of the present invention, for each list in area to be measured
Big gun is modulated single-shot data according to offset distance attribute, and is based on wave equation migration algorithm, according to the single-shot before modulation
Data and modulated single-shot data, corresponding first imaging results of single-shot data and modulated single-shot number before being modulated
According to corresponding second imaging results, it is then based on two imaging results and obtains the corresponding earth's surface offset gather of single-shot, last root
According to the corresponding earth's surface offset gather of each single-shot, the offset gather in area to be measured is obtained.Method provided in an embodiment of the present invention,
Different GPU is utilized respectively it is not necessary that single-shot data are divided into several groups to calculate, can be obtained using only a GPU calculating
The offset gather in area to be measured, saves computing resource, so that computing resource utilization rate is higher.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is this hair
Bright some embodiments for those of ordinary skill in the art without creative efforts, can be with root
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is a kind of acquisition methods flow chart of offset gather provided in an embodiment of the present invention；
Fig. 2 is a kind of calculation flow chart of the earth's surface offset gather based on reverse-time migration provided in an embodiment of the present invention；
Fig. 3 is a kind of offset Marmousi trace gather signal obtained based on different offset methods provided in an embodiment of the present invention
Figure；
Fig. 4 is a kind of schematic diagram of offset gather obtained based on different offset methods provided in an embodiment of the present invention；
Fig. 5 is a kind of structural schematic diagram of the acquisition device of offset gather provided in an embodiment of the present invention；
Fig. 6 is the entity structure schematic diagram of a kind of electronic equipment provided in an embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
Fig. 1 is a kind of acquisition methods flow chart of offset gather provided in an embodiment of the present invention, the execution master of this method
Body is the acquisition device of offset gather, such as server, and the present invention is not especially limit this, as shown in Figure 1, should
Method includes:
Step 101, for each single-shot in multiple single-shots in area to be measured, it is based on wave equation migration algorithm, is obtained every
The corresponding second imaging knot of the modulation single-shot data of corresponding first imaging results of the single-shot data of a single-shot and each single-shot
Fruit；Wherein, the modulation single-shot data are obtained after being modulated according to offset distance attribute to the single-shot data.
Step 102, according to corresponding first imaging results of each single-shot and the second imaging results, it is corresponding to obtain each single-shot
Earth's surface offset gather.
Step 103, summarize the corresponding earth's surface offset gather of each single-shot, obtain the offset gather in the area to be measured.
Specifically, it is necessary first to which explanation, seismic prospecting refer to that elastic wave caused by artificial excitation utilizes underground
The difference of dielectric resilient and density is inferred by observing and analyzing propagation law of the seismic wave in underground of artificial earthquake generation
The property of subterranean strata and the geophysical exploration method of form.
The earth's surface range for needing to carry out seismic prospecting is known as area to be measured, in order to obtain the seismic exploration data in area to be measured,
It needs to arrange multiple single-shots (shot point) in area to be measured, and multiple wave detectors (geophone station) is set.It should be noted that of the invention
Embodiment is not especially limited the setting of single-shot and the setting of wave detector.
For each single-shot in area to be measured, the wave field of the elastic wave of single-shot excitation is known as shot point wave field, by each inspection
The wave field that wave device detects is known as geophone station wave field, by geophone station wave field of the single-shot at each moment and the shot point wave field at each moment
The referred to as single-shot data of single-shot.
Based on wave equation algorithm, corresponding first imaging results of single-shot data of each single-shot are obtained.Wherein, fluctuation side
Journey algorithm can be one-way wave migration algorithm or reverse-time migration algorithm etc., and the embodiment of the present invention is not construed as limiting this.
For each single-shot in area to be measured, single-shot data are modulated by offset distance attribute to obtain modulation single-shot number
According to, wherein offset distance attribute refers to single-shot to the distance value of each wave detector, by it to the geophone station wave field in single-shot data
It is modulated, single-shot is known as to the modulation of single-shot in the modulated geophone station wave field at each moment and the shot point wave field at each moment
Single-shot data.
Based on wave equation algorithm, corresponding second imaging results of modulation single-shot data of each single-shot are obtained.Wherein, wave
Dynamic equation algorithm can be one-way wave migration algorithm or reverse-time migration algorithm etc., and the embodiment of the present invention is not construed as limiting this.It needs
Illustrate, it is consistent with the wave equation algorithm of the second imaging results is obtained to obtain the first imaging results.
According to corresponding first imaging results of each single-shot and the second imaging results, it is inclined to obtain the corresponding earth's surface of each single-shot
It moves away from trace gather.Summarize the corresponding earth's surface offset gather of each single-shot, obtains the offset gather in the area to be measured.
Method provided in an embodiment of the present invention, for each single-shot in area to be measured, according to offset distance attribute, to single-shot number
It is obtained according to being modulated, and based on wave equation migration algorithm according to the single-shot data and modulated single-shot data before modulation
Corresponding first imaging results of single-shot data and corresponding second imaging results of modulated single-shot data before modulation, then base
The corresponding earth's surface offset gather of single-shot is obtained in two imaging results, finally according to each single-shot corresponding earth's surface offset distance road
Collection, obtains the offset gather in area to be measured.Method provided in an embodiment of the present invention calculates earth's surface offset gather, meter using GPU
Calculation amount only increases by 30%, saves computing resource, so that computing resource utilization rate is higher.
Based on any of the above-described embodiment, the wave equation migration algorithm in the embodiment of the present invention is reverse-time migration calculation
Method；Correspondingly,
Based on wave equation migration algorithm, corresponding first imaging results of single-shot data of each single-shot are obtained, comprising:
Geophone station wave field of each single-shot at each moment and the shot point wave field at each moment are obtained, and as each single-shot
Single-shot data；
Corresponding first imaging results are obtained according to the single-shot data of each single-shot based on the reverse-time migration algorithm.
Wherein, it is based on the reverse-time migration algorithm, according to the single-shot data of each single-shot, obtains corresponding first imaging knot
The calculation formula of fruit I (x, s) are as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and n is positive integer, and u (x, s, t) is single-shot s in t
The geophone station wave field at quarter, s (x, s, t) are shot point wave field of the single-shot s in t moment.
Specifically, the embodiment of the present invention mainly illustrates how to calculate single-shot according to the single-shot data (single-shot data) before modulation
Corresponding first imaging results.
In embodiments of the present invention, it is (single that the geophone station wave field of single-shot and shot point wave field are known as to the single-shot data before modulating
Big gun data), the geophone station wave field of each single-shot is multiplied with shot point wave field, it is after being multiplied as a result, and by all single-shots
Results added after multiplication obtains the first imaging results.Wherein, geophone station wave field and shot point wave field, which are propagated by wave field, calculates
It obtains.
Based on any of the above-described embodiment, it is based on wave equation migration algorithm, obtains the modulation single-shot data pair of each single-shot
The second imaging results answered, comprising:
The calculation formula of the first imaging results I (x, s) is deformed, deformed first imaging results I is obtained
The calculation formula of (x, s) is as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and r represents geophone station, and n is positive integer, s (x, s, t)
It is single-shot s in the shot point wave field of t moment, u (x, s, r, t) is single track geophone station wave field of the single-shot s in t moment；
According to offset distance attribute, the single track geophone station wave field to each single-shot at each moment is modulated, after obtaining modulation
Single track geophone station wave field；
By each single-shot in the modulated single track geophone station wave field at each moment and in the shot point wave field at each moment, as every
The modulation single-shot data of a single-shot；
Corresponding second imaging knot is obtained according to the modulation single-shot data of each single-shot based on the reverse-time migration algorithm
Fruit.
Wherein, it is based on the reverse-time migration algorithm, according to the modulation single-shot data of each single-shot, obtains corresponding the second one-tenth
As result IhThe calculation formula of (x, s) are as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and r represents geophone station, and n is positive integer, s (x, s, t)
It is single-shot s in the shot point wave field of t moment, u (x, s, r, t) is single track geophone station wave field of the single-shot s in t moment, and offset is single
Distance of the big gun s to geophone station r.
Specifically, the embodiment of the present invention mainly illustrates how to be calculated according to modulated single-shot data (modulation single-shot data)
Corresponding second imaging results of single-shot.
In embodiments of the present invention, the calculation formula of the first imaging results in above-described embodiment is deformed first, that is,
Above-mentioned formula is rewritten into the formula of single track form of calculation, as follows:
Then the modulation for carrying out single track offset distance to the wave field of geophone station according to the distance of single-shot to each geophone station obtains the
Two imaging results IhThe calculation formula of (x, s) are as follows:
It is obtained based on any of the above-described embodiment according to corresponding first imaging results of each single-shot and the second imaging results
The corresponding earth's surface offset gather I of each single-shot (x, h | calculation formula s) are as follows:
I (x, h | s)=I (x, s) δ (h-hA(x,s))
Wherein,
Wherein, x represents imaging position, behalf single-shot, and I (x, s) is the first imaging results, Ih(x, s) is the second imaging knot
Fruit, h represent offset distance, hA(x, s) is that the second imaging results of single-shot s are corresponding relative to each imaging point of the first imaging results
Offset distance.
Wherein, summarize the corresponding earth's surface offset gather of each single-shot, obtain the meter of the offset gather in the area to be measured
Calculate formula are as follows:
Wherein, x represents imaging position, behalf single-shot, and I (x, h | it s) is the corresponding earth's surface offset gather of single-shot s.
Specifically, the embodiment of the present invention is mainly illustrated how according to corresponding first imaging results of each single-shot and the second one-tenth
As a result, obtain the corresponding earth's surface offset gather of each single-shot, and how according to the corresponding earth's surface offset distance of each single-shot
Trace gather obtains the offset gather in area to be measured.
As a preferred embodiment, combination Fig. 2 of the embodiment of the present invention carries out into one the method in above-described embodiment
Step illustrates, and Fig. 2 is a kind of calculation process of the earth's surface offset gather based on reverse-time migration provided in an embodiment of the present invention
Figure, as shown in Fig. 2, calculation process is as follows:
1, in CPU, offset parameter prepares；
2, in CPU, migration velocity is read；
3, in CPU, random boundary condition setting；
4, by the above-mentioned data copy in CPU into GPU；
5, in GPU, it=0 is enabled, according to the data received, carries out shot point main story；
6, in GPU, it is inserted into shot point wavefield data；
7, in GPU, it=it+1 is enabled；
8, in GPU, judgement: whether it is less than itmax；
9, in GPU, if so, carrying out the anti-pass of shot point wave field；
10, in GPU, former detection wave field is propagated, and is inserted into geophone station wave field；It modulates detection field to propagate, is inserted into geophone station wave field；
11, in GPU, offset distance is extracted；
12, it=it-1 is enabled；
13, judge whether it is greater than or equal to 0；
14, if it is not, then re-executeing the steps 8；If it is not, then copying data to CPU；
15, in CPU, according to the data of copy, offset gather is extracted.
It should be noted that, if removing modulation detection wave field part of propagation, be exactly reverse-time migration calculation process in Fig. 2,
Therefore from figure 2 it can be seen that mainly including the positive and negative biography of shot point wave field, geophone station anti-pass modulates geophone station anti-pass, because only
Modulation geophone station anti-pass is increased, therefore under shot point and geophone station aperture unanimous circumstances, calculation amount merely adds 1/3, because
This this method is a kind of efficient calculation method, saves computing resource, so that computing resource utilization rate is higher.
Fig. 3 is a kind of offset Marmousi trace gather signal obtained based on different offset methods provided in an embodiment of the present invention
Figure, wherein the offset Marmousi trace gather of (a) partially to be obtained using the method provided in the embodiment of the present invention is illustrated, the portion (b)
It is divided into and is illustrated using the offset Marmousi trace gather that Kirchhoff offset method obtains, it is clear that is provided in the embodiment of the present invention
Method is more preferable relative to Kirchhoff migration imaging effect.Fig. 4 is provided in an embodiment of the present invention a kind of based on different offset sides
The schematic diagram for the offset gather that method obtains, wherein (a) is partially inclined to be obtained using the method provided in the embodiment of the present invention
The signal away from trace gather is moved, is partially (b) signal of the offset gather obtained using Kirchhoff offset method, it is clear that this hair
The method provided in bright embodiment is smoother relative to the trace gather that Kirchhoff offset generates, and resolution ratio is higher.
Based on any of the above-described embodiment, Fig. 5 is a kind of acquisition device of offset gather provided in an embodiment of the present invention
Structural schematic diagram, as shown in figure 5, the device includes:
Imaging results obtain module 501, for being based on wave equation for each single-shot in multiple single-shots in area to be measured
Migration algorithm obtains corresponding first imaging results of single-shot data of each single-shot and the modulation single-shot data pair of each single-shot
The second imaging results answered；Wherein, the modulation single-shot data are to be modulated according to offset distance attribute to the single-shot data
It obtains afterwards；Earth's surface offset gather obtains module 502, for according to corresponding first imaging results of each single-shot and the second one-tenth
As a result, obtaining the corresponding earth's surface offset gather of each single-shot；Offset gather obtains module 503, for summarizing each list
The corresponding earth's surface offset gather of big gun obtains the offset gather in the area to be measured.
Specifically, device provided in an embodiment of the present invention is specifically used for executing above method embodiment, the embodiment of the present invention
This is no longer repeated.Device provided in an embodiment of the present invention, for each single-shot in area to be measured, according to offset distance category
Property, single-shot data are modulated, and be based on wave equation migration algorithm, according to before modulation single-shot data and modulated list
Big gun data, corresponding second imaging of corresponding first imaging results of single-shot data and modulated single-shot data before being modulated
As a result, being then based on two imaging results obtains the corresponding earth's surface offset gather of single-shot, it is finally corresponding according to each single-shot
Earth's surface offset gather obtains the offset gather in area to be measured.Method provided in an embodiment of the present invention, it is not necessary that single-shot data are equal
It is divided into several groups and is utilized respectively different GPU and calculated, the offset distance road that area to be measured can be obtained is calculated using only a GPU
Collection, saves computing resource, so that computing resource utilization rate is higher.
Fig. 6 is the entity structure schematic diagram of a kind of electronic equipment provided in an embodiment of the present invention, as shown in fig. 6, the electronics
Equipment may include: processor (processor) 601, communication interface (Communications Interface) 602, storage
Device (memory) 603 and communication bus 604, wherein processor 601, communication interface 602, memory 603 pass through communication bus
604 complete mutual communication.Processor 601, which can call, to be stored on memory 603 and can run on processor 601
Computer program, the method to execute the various embodiments described above offer, for example, for each of multiple single-shots in area to be measured
Single-shot, be based on wave equation migration algorithm, obtain each single-shot corresponding first imaging results of single-shot data and each list
Corresponding second imaging results of modulation single-shot data of big gun；Wherein, the modulation single-shot data are according to offset distance attribute to institute
It states and obtains after single-shot data are modulated；According to corresponding first imaging results of each single-shot and the second imaging results, obtain
The corresponding earth's surface offset gather of each single-shot；Summarize the corresponding earth's surface offset gather of each single-shot, obtains the area to be measured
Offset gather.
In addition, the logical order in above-mentioned memory 603 can be realized by way of SFU software functional unit and conduct
Independent product when selling or using, can store in a computer readable storage medium.Based on this understanding, originally
The technical solution of the inventive embodiments substantially part of the part that contributes to existing technology or the technical solution in other words
It can be embodied in the form of software products, which is stored in a storage medium, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes the present invention respectively
The all or part of the steps of a embodiment the method.And storage medium above-mentioned includes: USB flash disk, mobile hard disk, read-only memory
(ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic or disk
Etc. the various media that can store program code.
The embodiment of the present invention also provides a kind of non-transient computer readable storage medium, is stored thereon with computer program,
The computer program is implemented to carry out the transmission method of the various embodiments described above offer when being executed by processor, for example, for
Each single-shot in area to be measured in multiple single-shots is based on wave equation migration algorithm, and the single-shot data for obtaining each single-shot are corresponding
The first imaging results and each single-shot corresponding second imaging results of modulation single-shot data；Wherein, the modulation single-shot
Data are obtained after being modulated according to offset distance attribute to the single-shot data；According to corresponding first imaging of each single-shot
As a result with the second imaging results, the corresponding earth's surface offset gather of each single-shot is obtained；It is inclined to summarize the corresponding earth's surface of each single-shot
It moves away from trace gather, obtains the offset gather in the area to be measured.
The apparatus embodiments described above are merely exemplary, wherein described, unit can as illustrated by the separation member
It is physically separated with being or may not be, component shown as a unit may or may not be physics list
Member, it can it is in one place, or may be distributed over multiple network units.It can be selected according to the actual needs
In some or all of the modules achieve the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying creativeness
Labour in the case where, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Based on this understanding, on
Stating technical solution, substantially the part that contributes to existing technology can be embodied in the form of software products in other words, should
Computer software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes each implementation
Method described in certain parts of example or embodiment.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features；
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (10)
1. a kind of acquisition methods of offset gather characterized by comprising
For each single-shot in multiple single-shots in area to be measured, it is based on wave equation migration algorithm, obtains the single-shot of each single-shot
Corresponding second imaging results of modulation single-shot data of corresponding first imaging results of data and each single-shot；Wherein, described
Modulation single-shot data are obtained after being modulated according to offset distance attribute to the single-shot data；
According to corresponding first imaging results of each single-shot and the second imaging results, the corresponding earth's surface offset distance of each single-shot is obtained
Trace gather；
Summarize the corresponding earth's surface offset gather of each single-shot, obtains the offset gather in the area to be measured.
2. the acquisition methods of offset gather according to claim 1, which is characterized in that the wave equation migration algorithm
For reverse-time migration algorithm；Correspondingly,
Based on wave equation migration algorithm, corresponding first imaging results of single-shot data of each single-shot are obtained, comprising:
Obtain geophone station wave field of each single-shot at each moment and the shot point wave field at each moment, and the single-shot as each single-shot
Data；
Corresponding first imaging results are obtained according to the single-shot data of each single-shot based on the reverse-time migration algorithm.
3. the acquisition methods of offset gather according to claim 2, which is characterized in that calculated based on the reverse-time migration
Method obtains the calculation formula of corresponding first imaging results I (x, s) according to the single-shot data of each single-shot are as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and n is positive integer, and u (x, s, t) is single-shot s in t moment
Geophone station wave field, s (x, s, t) are shot point wave field of the single-shot s in t moment.
4. the acquisition methods of offset gather according to claim 3, which is characterized in that calculated based on wave equation migration
Method obtains corresponding second imaging results of modulation single-shot data of each single-shot, comprising:
The calculation formula of the first imaging results I (x, s) is deformed, deformed first imaging results I (x, s) is obtained
Calculation formula it is as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and r represents geophone station, and n is positive integer, and s (x, s, t) is single
For big gun s in the shot point wave field of t moment, u (x, s, r, t) is single track geophone station wave field of the single-shot s in t moment；
According to offset distance attribute, the single track geophone station wave field to each single-shot at each moment is modulated, and obtains modulated list
Road geophone station wave field；
By each single-shot in the modulated single track geophone station wave field at each moment and in the shot point wave field at each moment, as each list
The modulation single-shot data of big gun；
Corresponding second imaging results are obtained according to the modulation single-shot data of each single-shot based on the reverse-time migration algorithm.
5. the acquisition methods of offset gather according to claim 4, which is characterized in that calculated based on the reverse-time migration
Method obtains corresponding second imaging results I according to the modulation single-shot data of each single-shothThe calculation formula of (x, s) are as follows:
Wherein, x represents imaging position, behalf single-shot, and t represents the moment, and r represents geophone station, and n is positive integer, and s (x, s, t) is single
Big gun s is in the shot point wave field of t moment, and u (x, s, r, t) is single track geophone station wave field of the single-shot s in t moment, and offset is that single-shot s is arrived
The distance of geophone station r.
6. the acquisition methods of offset gather according to claim 5, which is characterized in that according to each single-shot corresponding
One imaging results and the second imaging results, obtain the corresponding earth's surface offset gather I of each single-shot (x, h | calculation formula s)
Are as follows:
I (x, h | s)=I (x, s) δ (h-hA(x,s))
Wherein,
Wherein, x represents imaging position, behalf single-shot, and I (x, s) is the first imaging results, Ih(x, s) is the second imaging results, h
Represent offset distance, hA(x, s) is each imaging point corresponding offset of the second imaging results of single-shot s relative to the first imaging results
Away from.
7. the acquisition methods of offset gather according to claim 6, which is characterized in that summarize each single-shot correspondingly
Table offset gather obtains the calculation formula of the offset gather in the area to be measured are as follows:
Wherein, x represents imaging position, behalf single-shot, and I (x, h | it s) is the corresponding earth's surface offset gather of single-shot s.
8. a kind of acquisition device of offset gather characterized by comprising
Imaging results obtain module, for being calculated based on wave equation migration for each single-shot in multiple single-shots in area to be measured
Method, the modulation single-shot data of corresponding first imaging results of single-shot data and each single-shot that obtain each single-shot corresponding
Two imaging results；Wherein, the modulation single-shot data are obtained after being modulated according to offset distance attribute to the single-shot data
's；
Earth's surface offset gather obtains module, is used for according to corresponding first imaging results of each single-shot and the second imaging results,
Obtain the corresponding earth's surface offset gather of each single-shot；
Offset gather obtains module, for summarizing the corresponding earth's surface offset gather of each single-shot, obtains the area to be measured
Offset gather.
9. a kind of electronic equipment including memory, processor and stores the calculating that can be run on a memory and on a processor
Machine program, which is characterized in that the processor is realized when executing described program such as any one of claim 1 to 7 the method
Step.
10. a kind of non-transient computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer
It is realized when program is executed by processor such as the step of any one of claim 1 to 7 the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910711930.0A CN110426736B (en) | 2019-08-02 | 2019-08-02 | Offset gather acquisition method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910711930.0A CN110426736B (en) | 2019-08-02 | 2019-08-02 | Offset gather acquisition method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110426736A true CN110426736A (en) | 2019-11-08 |
| CN110426736B CN110426736B (en) | 2020-04-07 |
Family
ID=68412328
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910711930.0A Active CN110426736B (en) | 2019-08-02 | 2019-08-02 | Offset gather acquisition method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110426736B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140149046A1 (en) * | 2011-05-23 | 2014-05-29 | Total Sa | Method of processing seismic data by providing surface offset common image gathers |
| CN104635264A (en) * | 2014-08-28 | 2015-05-20 | 中国石油天然气股份有限公司 | Processing method and apparatus of pre-stack seismic data |
| CN106250101A (en) * | 2015-06-12 | 2016-12-21 | 中国石油化工股份有限公司 | Migration before stack method for parallel processing based on MapReduce and device |
| CN107390270A (en) * | 2017-07-13 | 2017-11-24 | 西安交通大学 | A kind of AVA analysis methods based on elastic wave reverse-time migration ADCIGs |
| CN107765298A (en) * | 2017-07-14 | 2018-03-06 | 中国石油天然气集团公司 | A kind of method and device for determining near-surface velocity model |
| CN108802813A (en) * | 2018-06-13 | 2018-11-13 | 中国石油大学(华东) | A kind of multi-component seismic data offset imaging method and system |
-
2019
- 2019-08-02 CN CN201910711930.0A patent/CN110426736B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140149046A1 (en) * | 2011-05-23 | 2014-05-29 | Total Sa | Method of processing seismic data by providing surface offset common image gathers |
| CN104635264A (en) * | 2014-08-28 | 2015-05-20 | 中国石油天然气股份有限公司 | Processing method and apparatus of pre-stack seismic data |
| CN106250101A (en) * | 2015-06-12 | 2016-12-21 | 中国石油化工股份有限公司 | Migration before stack method for parallel processing based on MapReduce and device |
| CN107390270A (en) * | 2017-07-13 | 2017-11-24 | 西安交通大学 | A kind of AVA analysis methods based on elastic wave reverse-time migration ADCIGs |
| CN107765298A (en) * | 2017-07-14 | 2018-03-06 | 中国石油天然气集团公司 | A kind of method and device for determining near-surface velocity model |
| CN108802813A (en) * | 2018-06-13 | 2018-11-13 | 中国石油大学(华东) | A kind of multi-component seismic data offset imaging method and system |
Non-Patent Citations (1)
| Title |
|---|
| 王江斌: "基于逆时偏移的偏移距域共成像点道集计算", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110426736B (en) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772583B (en) | A kind of earthquake diffracted wave separation method and device | |
| US9632192B2 (en) | Method of processing seismic data by providing surface offset common image gathers | |
| Shi et al. | Seismic elastic RTM with vector-wavefield decomposition | |
| CN113740901B (en) | Land seismic data full-waveform inversion method and device based on complex undulating surface | |
| WO2017162731A1 (en) | Method of operating a data-processing system for the simulation of the acoustic wave propagation in the transversely isotropic media comprising an hydrocarbon reservoir | |
| CN109633742B (en) | Full waveform inversion method and device | |
| RU2570827C2 (en) | Hybrid method for full-waveform inversion using simultaneous and sequential source method | |
| CN108363096A (en) | A kind of land and water detector seismic data uplink and downlink wave field separation method and device | |
| AU2017366757A1 (en) | Seismic acquisition geometry evaluation using full-waveform inversion | |
| Faucher et al. | Eigenvector models for solving the seismic inverse problem for the Helmholtz equation | |
| US10215866B2 (en) | Seismic modeling system providing seismic survey data frequency domain inpainting and related methods | |
| US10429527B2 (en) | Seismic modeling system providing seismic survey data inpainting based upon suspect region boundary comparisons and related methods | |
| Birnie et al. | An introduction to distributed training of deep neural networks for segmentation tasks with large seismic data sets | |
| US20240134075A1 (en) | Method, device and computer device for decoupling anisotropic elastic wave | |
| CN105204063A (en) | Seismic data velocity model establishing method and device | |
| US20200132871A1 (en) | Post-Stack Kirchhoff Depth De-Migration Method for Tilted Transverse Isotropic (TTI) and Heterogeneous Media Based on Ray Tracing on Migrated Data | |
| CN112231974B (en) | Deep learning-based method and system for recovering seismic wave field characteristics of rock breaking seismic source of TBM (Tunnel boring machine) | |
| Wang et al. | Waveform inversion including well constraints, anisotropy, and attenuation | |
| US11143769B2 (en) | Seismic modeling system providing seismic survey data spatial domain exemplar inpainting and related methods | |
| Jing et al. | Fault detection using a convolutional neural network trained with point-spread function-convolution-based samples | |
| CN110426736A (en) | A kind of acquisition methods and device of offset gather | |
| Roncoroni et al. | Efficient extraction of seismic reflection with Deep Learning | |
| CN105403919B (en) | A kind of reverse-time migration imaging method and device | |
| EP3859405A1 (en) | Computer-implemented method for determining a velocity image of a domain of the subsurface structural geology in an oil and gas reservoir | |
| US9720117B1 (en) | Imaging subsurface properties using a parallel processing database system |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210702 Address after: Room 421, unit 1, 99 Yuexiu Road, Haidian District, Beijing 100096 Patentee after: BEIJING GEOSTAR SCIENCE & TECHNOLOGY Co.,Ltd. Address before: 100083 No. 29, Haidian District, Beijing, Xueyuan Road Patentee before: China University of Geosciences (Beijing) |