CN103472481A - Reverse-time migration and angle gather extraction method by GPU - Google Patents
Reverse-time migration and angle gather extraction method by GPU Download PDFInfo
- Publication number
- CN103472481A CN103472481A CN2012101845928A CN201210184592A CN103472481A CN 103472481 A CN103472481 A CN 103472481A CN 2012101845928 A CN2012101845928 A CN 2012101845928A CN 201210184592 A CN201210184592 A CN 201210184592A CN 103472481 A CN103472481 A CN 103472481A
- Authority
- CN
- China
- Prior art keywords
- wave field
- reverse
- time
- angular
- gpu
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a reverse-time migration and angle gather extraction method by GPU and belongs to the field of seismic exploration and exploitation in the petroleum industry. The method provided by the invention comprises the following steps: (1) a wave field at a spot point continues forward; (2) concurrent execution of reverse-time migration and angle information statistics is realized; and (3) angle gather extraction is carried out. The invention provides a rapid and stable method for reverse-time migration and angle gather extraction. The method provided by the invention has concurrent execution ability of reverse-time migration and angle information statistics, is more suitable for implementation of a GPU program, and is adopted to solve the problem that reverse-time migration and angle gather extraction must be serially executed formerly and solve disadvantages of low computational efficiency, too much storage requirements, poor fault-tolerant ability and difficulty in trial and error.
Description
Technical field
The invention belongs to seismic prospecting and the development field of petroleum industry, be specifically related to a kind of GPU of utilization and carry out the method that angular-trace gather is extracted in reverse-time migration, the common picture point road collection of the migration imaging software systems that record for artificial earthquake extracts.
Background technology
The English full name of GPU is Graphic Processing Unit, translator of Chinese is " graphic process unit ", it be from CPU (central processing unit), be derived a concept, it is the heart of video card, carry internal memory with the video card plate and form a subsystem, jointly determine the graphics process performance of PC system.Increasing graphical application makes GPU become more and more important in modern computing machine, this chip emerges rapidly in after CPU dominates the PC performance to move towards decades nearly 2 years, its Floating-point Computation ability has reached ten times even tens times of the CPU same period, at high-performance computing sector, almost reached the status on an equal footing with CPU.
Graphic process unit (GPU) has been become to a hot issue at present for the problem that solves general-purpose computations (GPGPU) and related fields thereof.In fact, whole IT industry has all recognized that the GPU general-purpose computations will bring revolutionary variation to PC observantly, and then has influence on the development of CPU.GPU aims at the image Treatment Design, and its storage system is actually the fragmented storage space of a two dimension, comprises a segment number (therefrom reading images) and two-dimensional address (X in image, Y coordinate).In addition, without any indirect write command.Write address output is definite by raster processor, and can not be by routine change.This is great challenge for the algorithm among storer for NATURAL DISTRIBUTION, and the access efficiency of storer is the key issue of boosting algorithm efficiency.Utilized GPU to do high-performance calculation is very difficult work in the past, and the appearance of the senior programming languages such as CUDA and the OpenCL new hand that makes to programme also can grasp the calculated performance advantage of GPU after 2007.
In the application of reverse-time migration, the research work on a basis is exactly the extraction of angular-trace gather.Since entering 21 century, wave equation output is total to the focus that the reflection angle image set just becomes industry research.Its implementation can be divided into two large classes: a class is the extracting method (as shown in Figure 1) before imaging, hereinafter be called method one, the method need to be calculated the wave field angle in real time in the process of reverse-time migration, and real-time output angle road collection, the indivisible realization of three (as shown in Figure 5, Fig. 5 is the refinement flow process of Fig. 1, wherein angular-trace gather extraction and wave field extrapolation are coupling in a module and carry out), and needing the memory headroom of a large amount of (TB level) as guarantee, is the implacable condition of current computer technology; Another kind of is extracting method (as shown in Figure 2) after imaging, hereinafter is called method two, and the method adopts translation or time shift effect condition at first to generate offset distance picture point road collection altogether, then by the methods such as Tau-P convert conversion generation angular-trace gather.Extracting method after imaging is suitable for offset module in the past, only need to do less change and follow-up conversion process just can obtain angular-trace gather, is also the more method of use at present.But, due to common offset distance and the transformational relation of angle can not accomplish corresponding one by one, so equally also there is illusion in the angular-trace gather of generation.
It is even four-dimensional that seismic data acquisition of today has developed into three-dimensional, its data volume is very huge especially, reached the even prestack list big gun data of 1OTB level of TB level, computing time and the IO time of so huge data are also the difficult problems that reverse-time migration angular-trace gather extraction work must face.
Due to the develop rapidly of GPU high-performance calculation, can increase substantially the problem that reverse-time migration and angular-trace gather calculate, but still do not solve storage problem in the implementation process of prior art.
Summary of the invention
The object of the invention is to solve the difficult problem existed in above-mentioned prior art, a kind of method that provides GPU of utilization to carry out the extraction of reverse-time migration angular-trace gather, can extract angular-trace gather in fast and stable ground.
The present invention is achieved by the following technical solutions:
A kind of GPU of utilization carries out the method for reverse-time migration angular-trace gather extraction, said method comprising the steps of:
(1) shot point wave field forward continuation;
(2) executed in parallel reverse-time migration and angle information statistics;
(3) angular-trace gather extracts.
Described step (1) comprises the following steps:
(A1) put into artificial wavelet focus in sp location;
(A2) wave field of the seismic event that utilizes finite difference method to go out the artificial wavelet epicenter excitation in step (A1) carries out the forward continuation along time orientation, obtains not shot point wave field in the same time;
Described step (A2) completes on GPU.
In described step (2), described reverse-time migration comprises the following steps:
(B1) the not shot point wave field in the same time obtained in step (A2) is stored in the wave field Data buffer;
(B2) the shot point wave field that will obtain from step (A2), along the continuation of time orientation forward to constantly maximum, obtains the maximum wave field constantly of shot point; The forward continuation refers to the direction continuation wave field increased progressively along the time;
(B3) read the maximum wave field constantly of geophone station from the geological data of field acquisition, then from maximum constantly, shot point wave field and geophone station wave field, along the time orientation backward extension, are obtained to not wave field in the same time of shot point and geophone station; Backward extension refers to the direction continuation wave field successively decreased along the time;
(B4) by each, shot point wave field constantly and geophone station wave field multiply each other and obtain imaging results;
(B5) the corresponding single big gun of each shot point, the described imaging results using each single big gun in corresponding step (B4) is stored in separately on disk as single big gun migration imaging result;
(B6) each geological data from field acquisition comprises a lot of single big guns, utilizes the described imaging results in the step that each single big gun is corresponding (B4) to do the migration imaging result that many big gun stacks obtain many big gun stacks; Described many big gun stacks refer to that the imaging value for underground same position place is added up a plurality of single big gun imaging results.
In described step (2), described angle information statistics comprises the following steps:
(C1) angle information statistics: the shot point wave field in described wave field Data buffer is stored according to time sequencing, and extracted each amplitude and propagation angle of each locus constantly;
(C2) after completing steps (A2), using the shot point wave field that obtains in step (C1) each locus peak swing propagation angle constantly as angle information data volume store storage on disk;
(C3), on the migration imaging result of the described many big gun stacks that obtain in described step (B6), extract stratigraphic dip and be stored on disk as stratigraphic dip information data body.
Wherein, described step (B2) to (B5) and step (C1) to (C2) executed in parallel, perform step (C3) after step (B6) is finished again.
Described step (3) comprises the following steps:
(D1) from disk read step (C2) and (C3) in the angle information data volume and the stratigraphic dip information data body that obtain;
(D2) calculate incident angle according to described angle information data volume and stratigraphic dip information data body:
Incidence angle θ equals the shot point wave field and reaches half poor of propagation angle β constantly of peak swing and stratigraphic dip α at place, each locus, and equation expression is: θ=β-0.5 α;
(D3) single big gun migration imaging result of read step (B5) output from disk, the road that then single big gun migration imaging result is outputed to the incident angle obtained corresponding to step (D2) is concentrated, and obtains final angular-trace gather data volume.
Data Structures in described Data buffer comprises floating type data, integer data and a floating type array, and wherein floating type data storage current time value, integer data are stored the length of described floating-point array, the wave field of floating type storage of array current time.
Compared with prior art, the invention has the beneficial effects as follows:
(1) utilize the data buffering pool technology in the present invention, realized the parallelization embodiment that reverse-time migration, angle information are added up.Avoid in the past the method that the two can only serial implementation, improve counting yield.
(2), because the angle information of shot point wave field is added up generation before imaging, kept, under the front angle calculation precision conditions of imaging, realizing after reverse-time migration completes carrying out again the extraction work of angular-trace gather.In like this can methods to reduce noises one, reverse-time migration and angular-trace gather extract the memory cost simultaneously carried out, the imaging illusion that angular-trace gather extracts after imaging in simultaneously can removing method two.
(3) this patent is owing to only adding up wave field ceiling capacity angle information constantly, so can leave identical internal memory or disk space in angle calculation result in the same time not, realize the recycling of storage space, can significantly reduce like this demand of internal memory and disk.
The accompanying drawing explanation
Fig. 1 is the step block diagram of method one of the prior art.
Fig. 2 is the step block diagram of method two of the prior art.
Fig. 3 is that the present invention utilizes the GPU reverse-time migration to extract the modular structure figure designed in the method for angular-trace gather.
Fig. 4 is that the present invention utilizes GPU to extract the detailed step block diagram of reverse-time migration angular-trace gather method.
Fig. 5 is the detailed step block diagram that employing method one is extracted angular-trace gather.
Fig. 6 is that the present invention calculates the schematic diagram of incident angle according to shot point wave field propagation angle and stratigraphic dip.
Fig. 7 is the schematic diagram that in the present invention, single big gun descends an incident angle of imaging point accordingly.
Fig. 8 is the sigseeb2a illustraton of model.
Fig. 9-1st, utilize the inventive method to carry out angular-trace gather to the illustraton of model of Fig. 8 and extract the design sketch obtained.
Fig. 9-2nd, utilize the illustraton of model of a pair of Fig. 8 of method to carry out the design sketch that the angular-trace gather extraction obtains.
Fig. 9-3rd, utilize method two to carry out to the illustraton of model of Fig. 8 the design sketch that the angular-trace gather extraction obtains.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail:
The method that the present invention is directed to GPU reverse-time migration extraction angular-trace gather comprises three modules: reverse-time migration module (() in Fig. 4), angle information statistical module ((two) in Fig. 4) and angular-trace gather extraction module ((three) in Fig. 4), the segmentation of module is as Fig. 4, each step can realize with CPU and GPU, but efficiency can improve greatly on GPU, only have the needs that utilize GPU could meet commercial Application.Fig. 3 is the collaborative module map that realizes that the reverse-time migration angular-trace gather extracts of calculating of CPU/GPU.Wherein reverse-time migration module and angle information statistical module can parallel practices, for the angular-trace gather extraction module provides necessary input data; Shot point wave field forward continuation step in the reverse-time migration module need to provide the wave field data for Data buffer, and the angle information statistical module directly reads the wave field data from Data buffer.
Fig. 4 is detail flowchart (() in Fig. 4 of module parallel processing, (2), (3) mean which module is concrete steps belong to), wherein white background means reverse-time migration flow process (belonging to prior art), dark-background means that the implementing procedure of angle information statistical module (comprises the wave field Data buffer, the angle information statistics, peak swing propagation angle calculating constantly and the stratigraphic dip calculating and output of output and imaging section), after completing, these two flow processs utilize single big gun imaging results of reverse-time migration, the stratigraphic dip information extraction angular-trace gather of the peak swing angle information data volume of locus and imaging section.Idiographic flow is as follows:
(1) shot point wave field forward continuation
(2) executed in parallel reverse-time migration and angle information statistics
(3) angular-trace gather extracts.
In described step (1), the continuation of described shot point wave field forward comprises the following steps:
(A1) put into artificial wavelet focus in sp location;
(A2) utilize the continuation of finite difference method forward, can obtain not shot point wave field in the same time.This is prior art, must on GPU, carry out, because the upper execution of GPU just can reach satisfied counting yield, meets the commercial Application requirement; Described finite difference method is one of algorithm of basic wave field extrapolation;
In described step (2), described reverse-time migration comprises the following steps:
(B1) the not shot point wave field in the same time obtained in (A2) step is stored in the wave field Data buffer; In the angle information statistical module, in order to guarantee to implement with the parallelization of reverse-time migration module, the present invention has adopted the intermediate link of wave field Data buffer, this wave field Data buffer is mainly that the not wave field information in the same time that the shot point wave field is just being drilled in process temporarily is stored on internal memory or disk, for the statistical work of angle information.Wherein should to comprise wave field extrapolation on GPU parallel for executed in parallel, and CPU is upper or another GPU on the parallelization of angle information statistical module carry out.Can on the basis that does not affect original reverse-time migration efficiency, realize the extraction work of angular-trace gather like this;
(B2) the shot point wave field that will obtain from step (A2), along time orientation forward continuation (just be called and drill again) to constantly maximum, obtains the maximum wave field constantly of shot point;
(B3) read the maximum wave field constantly of geophone station from the geological data of field acquisition, then from maximum constantly, by shot point wave field and geophone station wave field backward extension, obtain shot point and geophone station not wave field in the same time (because the image-forming condition of reverse-time migration is along the time reversal integration, so need imaging in back-propagation process, if in forward continuation process, the wave field in each moment is stored, need a lot of disk spaces, can reach the demand of tens bases to hundreds of TB during three-dimensional work area, backwards calculation once can be avoided this storage, this is existing technology), backward extension just refers to the direction continuation wave field successively decreased along the time,
(B4) by each, shot point wave field constantly and geophone station wave field multiply each other and obtain imaging results (being the application dependent imaging condition in Fig. 4); Wave field is that its physical meaning of function of room and time is amplitude.
(B5) the corresponding single big gun of each shot point, the described imaging results using each single big gun in corresponding step (B4) is stored in separately on disk as single big gun migration imaging result;
(B6) each geological data from field acquisition comprises a lot of single big guns, utilizes the described imaging results of each single big gun to do the migration imaging result that many big gun stacks obtain many big gun stacks; Described many big guns stack refers to a plurality of single big gun imaging results, for the imaging value of underground same position, added up.
In described step (2), described angle information statistics comprises the following steps:
(C1) angle information statistics: the wave field in described wave field Data buffer is stored according to time sequencing, the angle information extraction module is directly accessed the wave field in the wave field Data buffer, and extract each amplitude and propagation angle (extracting method can be Poynting vector method or local wavelength-division solution, and this is prior art) of each locus constantly, when in the wave field Data buffer, wave field number constantly is greater than 1, at first then wave field amplitude does not more in the same time extract peak swing propagation angle constantly, adding up each locus reaches peak swing propagation angle constantly (in calculating, the time is discretize, with equally spaced moment sequence, embody, in the continuation process, can be temporarily stored in wave field in Buffer Pool, in order to add up angle information, each angle information constantly can be covered by next wave field constantly after having added up, when if the speed of angle information statistics is slower than continuation speed, can there is the wave field in a plurality of moment in Buffer Pool simultaneously, now need more existingly to be in Buffer Pool that not wave field amplitude in the same time is to select the wherein wave field of amplitude maximum to calculate angle information, having added up rear Buffer Pool empties, wait for that back wave field constantly enters Buffer Pool, carry out again identical statistics.)。Wave field is the function of room and time, and the result of statistics is also the function in space.
The implementation method of angle information statistics is also important protection content of the present invention, the Data Structures of wave field Data buffer comprises floating type data, integer data and a floating type array, wherein floating data is stored the current time value, the wave field of floating-point storage of array current time (the direction continuation wave field increased progressively along the time exactly that the continuation of wave field forward refers to, its result is not wave field in the same time, the wave field comprised in the wave field Data buffer so retains in order to calculate angle, each constantly calculates once and can constantly be override by next, so need to indicate calculated angle, it it is that moment, it is described current time.), the length of integer data storages floating-point array.The angle information statistical module directly will be not from the wave field Data buffer wave field value in the same time read, then calculate each amplitude and propagation angle of each locus constantly.It is to be noted, (may there be the wave field in a plurality of moment in the angle information that when wave field is constantly more in the wave field Data buffer, at first relatively then amplitude extracts peak swing simultaneously in Buffer Pool, now need the comparison amplitude), in order to follow-up angular-trace gather extraction work.Can guarantee like this with reverse-time migration in the synchronism of wave field extrapolation, and, without any need for communication, only need time and a wave field in access wave field Data buffer to get final product.
(C2) (after whole continuation constantly) after completing steps (A2), using the shot point wave field that obtains in step (C1) each locus peak swing propagation angle constantly as the storage of angle information data volume (because can not guarantee that in the C1 implementation all moment comprise, whole complete this angle information data volumes that just can obtain of moment continuation.) be stored on disk; The part of the input data that described angle information data volume extracts as follow-up angular-trace gather;
(C3), on the migration imaging result of the many big gun stacks that obtain in described step (B6), extract stratigraphic dip and be stored on disk as stratigraphic dip information data body; The part of the input data that described stratigraphic dip information data body extracts as follow-up angular-trace gather;
Wherein, step (B2) to (B5) can with (C1) to (C2) executed in parallel, (B6) carry out again C3 after being finished;
In described step (3), described angular-trace gather extracts and comprises the following steps:
(D1) from disk read step (C2) and (C3) in the angle information data volume and the stratigraphic dip information data body that obtain;
(D2) calculate incident angle according to described angle information data volume and stratigraphic dip information data body: what angle information statistical module of the present invention finally preserved is that the shot point wave field reaches the wave field propagation angle in the peak swing moment and the stratigraphic dip after many big guns of reverse-time migration stacking image at each place, locus.Utilize this two angle informations, calculate that the incident angle information of seismic event extracts the imaging value generation angular-trace gather in imaging results.Concrete projectional technique is (in Fig. 6, S means that shot point, R mean geophone station) as shown in Figure 6, incidence angle θ equals the shot point wave field and reaches half poor of wave field propagation angle β constantly of peak swing and stratigraphic dip (being called again the reflection line-ups inclination angle) α at place, each locus, and equation expression is: θ=β-0.5 α.
(D3) single big gun migration imaging result of read step (B5) output from disk, the road that then single big gun migration imaging result is outputed to the incident angle obtained corresponding to step (D2) is concentrated, and obtains final angular-trace gather data volume.It is self-assembling formation that Fig. 7 has described the corresponding different incident angle of different single big guns.
In the angular-trace gather abstraction module, each imaging point needs the information of different incidence angles, and the present invention is mapped as different single emplacements to different incidence angles and puts, as shown in Figure 7 (X
smean the sp location coordinate), the corresponding incident angle of the imaging results of single big gun data, preserve respectively the imaging results (being the result that the single big gun migration result in step (B5) is stored separately) of single big gun data in implementation process, as the input data of follow-up angular-trace gather extraction module.
Described step (A2) completes on GPU, and other step completes on GPU or CPU.The present invention mainly emphasizes that the division methods of this flow process is conducive to parallel and raises the efficiency, and (A2) the finite difference calculated amount of step is very large needs use GPU, but angle extraction can be taked to complete on GPU or CPU for distinct methods.
The international SEG master pattern sigseeb2a (as shown in Figure 8) of below usining describes effect of the present invention as instance data of the present invention, and the earthquake work area information of this model is as shown in table 1.
| The big gun number | 500 big guns |
| Sampling interval | 8ms |
| Sampling length | 12s |
| Dan Bao road number | 1067 roads |
| Migration aperture | 35 kilometers of maximums diameter of hole (diameter) |
| Deviation frequency | Maximum frequency 105HZ |
| Sizing grid | Grid 2133 * 1200 |
| Bin size | Bin 37.5 * 25inch |
Table 1
Its hardware environment of using is as shown in table 2.
| |
1 |
| |
1 |
| Intel CPU overall performance | 4 cores * 2.0GHZ=8GHZ |
| The GPU overall performance | 4×240×1.44GHZ＝1382GHZ |
| Host memory | 24GB |
| The GPU video memory | 16GB |
| Hard disk | 2TB |
| Network environment | 1000M b/s |
| Operating system | LIN UX CentOS 4.8 |
Table 2
The angle intervals of output is 1 degree, exports 180 angles, utilizes respectively the inventive method and method one, method two to realize that angular-trace gather extracts, and the time of using is as shown in table 3.
| The inventive method | Method one | Method two | |
| Computing time | 2.5 hour | 5 hours | 4.8 hour |
| The required memory space | 78MB | 13.7GB | 3.8GB |
| Disk space required | 6.5GB | 576GB | 160GB |
Table 3
As can be seen from Table 3, the inventive method calculating and storage resources used is minimum.
Method one is compared with the inventive method, and the rise of computing time is due to reverse-time migration part in method one and angular-trace gather Extraction parts serial computing, often does a step wave field extrapolation and all will wait for that completing just of angular-trace gather extraction work can carry out next step continuation.The rise of method one disk space required and memory headroom is owing to need to be the angular-trace gather opening space in reverse-time migration, and the Output rusults of each big gun is all one and is angular-trace gather, and it increases multiple and equates with the output angle number.Utilize the angular-trace gather of method one calculating and come to the same thing (the seeing accompanying drawing 9-1 and 9-2) that the inventive method is calculated.
Method two is compared with the inventive method, and the rise of computing time is to need the imaging value of multiple offset distance in the reverse-time migration process due to method two, and the increase of its calculated amount is the use of translation image-forming condition.The disk space that method two is required and the rise of memory headroom are owing to needing to export common offset road collection in the reverse-time migration process, it is relevant to exported offset distance number that it increases multiple, the migration result of each single big gun is an offset gather, and its rise multiple is the number of offset distance.The present embodiment is to adopt 50 offset distances of output to calculate.The present invention extracts the illusion of angular-trace gather and wants decrease, Fig. 9-3rd, the angular-trace gather extraction effect that utilizes method two to obtain, Fig. 9-1st, the angular-trace gather extraction effect that utilizes the inventive method to obtain, comparison diagram 9-1 and Fig. 9-3 are visible, the angular-trace gather that the present invention extracts is all evened up substantially at reflecting interface Chu road collection, has hardly illusion.
The innovative point of this patent is as follows:
(1) innovation of the extraction scheme of angular-trace gather: the present invention is the advantage in conjunction with the angular-trace gather extracting method of (method two) after (method one) and imaging before imaging in prior art, under the prerequisite of the computational accuracy of angle Information Statistics, realized the method that after the imaging, angular-trace gather extracts before the succession imaging.
(2) introduce the parallel computation that Data buffer has been realized angle information statistics part and reverse-time migration part, improved counting yield.
(3) the independent storage of stratigraphic dip data volume in the independent storage of angle-data body, angle information statistical module in the independent storage of single big gun migration result, angle information statistical module in reverse-time migration module, guarantee to carry out again the feasibility of angular-trace gather extraction work after reverse-time migration completes, and avoided reverse-time migration in the past and angular-trace gather to extract the memory cost simultaneously carried out.
(4) this patent is owing to only adding up wave field ceiling capacity angle information constantly, so be placed on identical internal memory or disk space in angle calculation result in the same time not, realize the recycling of storage space, can significantly reduce like this demand of internal memory and disk.
Technique scheme is one embodiment of the present invention, for those skilled in the art, on the basis that the invention discloses application process and principle, be easy to make various types of improvement or distortion, and be not limited only to the described method of the above-mentioned embodiment of the present invention, therefore previously described mode is just preferred, and does not have restrictive meaning.
Claims (7)
1. a method of utilizing GPU to carry out the extraction of reverse-time migration angular-trace gather is characterized in that: said method comprising the steps of:
(1) shot point wave field forward continuation;
(2) executed in parallel reverse-time migration and angle information statistics;
(3) angular-trace gather extracts.
2. the GPU of utilization according to claim 1 carries out the method for reverse-time migration angular-trace gather extraction, it is characterized in that: described step (1) comprises the following steps:
(A1) put into artificial wavelet focus in sp location;
(A2) wave field of the seismic event that utilizes finite difference method to go out the artificial wavelet epicenter excitation in step (A1) carries out the forward continuation along time orientation, obtains not shot point wave field in the same time;
Described step (A2) completes on GPU.
3. the GPU of utilization according to claim 2 carries out the method for reverse-time migration angular-trace gather extraction, it is characterized in that: in described step (2), described reverse-time migration comprises the following steps:
(B1) the not shot point wave field in the same time obtained in step (A2) is stored in the wave field Data buffer;
(B2) the shot point wave field that will obtain from step (A2), along the continuation of time orientation forward to constantly maximum, obtains the maximum wave field constantly of shot point; The forward continuation refers to the direction continuation wave field increased progressively along the time;
(B3) read the maximum wave field constantly of geophone station from the geological data of field acquisition, then from maximum constantly, shot point wave field and geophone station wave field, along the time orientation backward extension, are obtained to not wave field in the same time of shot point and geophone station; Backward extension refers to the direction continuation wave field successively decreased along the time;
(B4) by each, shot point wave field constantly and geophone station wave field multiply each other and obtain imaging results;
(B5) the corresponding single big gun of each shot point, the described imaging results using each single big gun in corresponding step (B4) is stored in separately on disk as single big gun migration imaging result;
(B6) each geological data from field acquisition comprises a lot of single big guns, utilizes the described imaging results in the step that each single big gun is corresponding (B4) to do the migration imaging result that many big gun stacks obtain many big gun stacks; Described many big gun stacks refer to that the imaging value for underground same position place is added up a plurality of single big gun imaging results.
4. the GPU of utilization according to claim 3 carries out the method for reverse-time migration angular-trace gather extraction, it is characterized in that: in described step (2), described angle information statistics comprises the following steps:
(C1) angle information statistics: the shot point wave field in described wave field Data buffer is stored according to time sequencing, and extracted each amplitude and propagation angle of each locus constantly;
(C2) after completing steps (A2), using the shot point wave field that obtains in step (C1) each locus peak swing propagation angle constantly as angle information data volume store storage on disk;
(C3), on the migration imaging result of the described many big gun stacks that obtain in described step (B6), extract stratigraphic dip and be stored on disk as stratigraphic dip information data body.
5. the GPU of utilization according to claim 4 carries out the method for reverse-time migration angular-trace gather extraction, it is characterized in that: described step (B2) to (B5) and step (C1) to (C2) executed in parallel perform step (C3) after step (B6) is finished again.
6. the GPU of utilization according to claim 5 carries out the method for reverse-time migration angular-trace gather extraction, it is characterized in that: described step (3) comprising:
(D1) from disk read step (C2) and (C3) in the angle information data volume and the stratigraphic dip information data body that obtain;
(D2) calculate incident angle according to described angle information data volume and stratigraphic dip information data body:
Incidence angle θ equals the shot point wave field and reaches half poor of propagation angle β constantly of peak swing and stratigraphic dip α at place, each locus, and equation expression is: θ=β-0.5 α;
(D3) single big gun migration imaging result of read step (B5) output from disk, the road that then single big gun migration imaging result is outputed to the incident angle obtained corresponding to step (D2) is concentrated, and obtains final angular-trace gather data volume.
7. the GPU of utilization according to claim 3 carries out the method for reverse-time migration angular-trace gather extraction, it is characterized in that: the Data Structures in described Data buffer comprises floating type data, integer data and a floating type array, and wherein floating type data storage current time value, integer data are stored the length of described floating-point array, the wave field of floating type storage of array current time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210184592.8A CN103472481B (en) | 2012-06-06 | 2012-06-06 | A kind of GPU of utilization carries out the method that angular-trace gather is extracted in reverse-time migration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210184592.8A CN103472481B (en) | 2012-06-06 | 2012-06-06 | A kind of GPU of utilization carries out the method that angular-trace gather is extracted in reverse-time migration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103472481A true CN103472481A (en) | 2013-12-25 |
| CN103472481B CN103472481B (en) | 2016-03-23 |
Family
ID=49797398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210184592.8A Active CN103472481B (en) | 2012-06-06 | 2012-06-06 | A kind of GPU of utilization carries out the method that angular-trace gather is extracted in reverse-time migration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103472481B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103969684A (en) * | 2014-05-27 | 2014-08-06 | 中国海洋石油总公司 | Method and device for dipole hypocenter seismic imaging |
| CN104199088A (en) * | 2014-08-15 | 2014-12-10 | 中国石油天然气集团公司 | Incident angle gather extraction method and system |
| CN105891881A (en) * | 2014-05-15 | 2016-08-24 | 蔡晓刚 | CCP (Common Converted Point) trace extraction method of horizontal interface VTI medium |
| CN108072895A (en) * | 2016-11-09 | 2018-05-25 | 中国石油化工股份有限公司 | A kind of anisotropy pre-Stack Reverse optimization method based on GPU |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100565773B1 (en) * | 2005-01-27 | 2006-03-28 | 윤광진 | Seismic imaging method with reverse-time migration using poynting vector |
| US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
| CN102269820A (en) * | 2010-06-01 | 2011-12-07 | 潜能恒信能源技术股份有限公司 | Three-dimensional seismic pre-stack reverse-time migration imaging method based on GPU (graphics processing unit) staggered grid with small memory capacity |
-
2012
- 2012-06-06 CN CN201210184592.8A patent/CN103472481B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100565773B1 (en) * | 2005-01-27 | 2006-03-28 | 윤광진 | Seismic imaging method with reverse-time migration using poynting vector |
| US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
| CN102269820A (en) * | 2010-06-01 | 2011-12-07 | 潜能恒信能源技术股份有限公司 | Three-dimensional seismic pre-stack reverse-time migration imaging method based on GPU (graphics processing unit) staggered grid with small memory capacity |
Non-Patent Citations (3)
| Title |
|---|
| YUZHANG 等: "通过逆时偏移计算角道集研究", 《勘探前沿》, 30 November 2010 (2010-11-30) * |
| 丁亮 等: "逆时偏移成像技术研究进展", 《地球物理学进展》, vol. 26, no. 3, 30 June 2011 (2011-06-30), pages 1085 - 1095 * |
| 刘红伟 等: "起伏地表叠前逆时偏移理论及GPU加速技术", 《地球物理学报》, vol. 54, no. 7, 31 July 2011 (2011-07-31) * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105891881A (en) * | 2014-05-15 | 2016-08-24 | 蔡晓刚 | CCP (Common Converted Point) trace extraction method of horizontal interface VTI medium |
| CN103969684A (en) * | 2014-05-27 | 2014-08-06 | 中国海洋石油总公司 | Method and device for dipole hypocenter seismic imaging |
| CN104199088A (en) * | 2014-08-15 | 2014-12-10 | 中国石油天然气集团公司 | Incident angle gather extraction method and system |
| CN104199088B (en) * | 2014-08-15 | 2017-02-08 | 中国石油天然气集团公司 | Incident angle gather extraction method and system |
| CN108072895A (en) * | 2016-11-09 | 2018-05-25 | 中国石油化工股份有限公司 | A kind of anisotropy pre-Stack Reverse optimization method based on GPU |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103472481B (en) | 2016-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103336959B (en) | A kind of vehicle checking method accelerated based on GPU multi-core parallel concurrent | |
| CN102565854B (en) | Mass data GPU (graphics processing unit) wave equation reverse time migration imaging method | |
| CN102490909B (en) | Multi-body separating analogy method of aircraft | |
| CN103472481B (en) | A kind of GPU of utilization carries out the method that angular-trace gather is extracted in reverse-time migration | |
| CN102609982B (en) | Topology discovery method of space geological data based on unstructured mode | |
| CN102692644A (en) | Depth domain common-image gather generation method | |
| CN103713314B (en) | Pre-stack time migration parallel processing method | |
| CN103049241A (en) | Method for improving computation performance of CPU (Central Processing Unit) +GPU (Graphics Processing Unit) heterogeneous device | |
| Ma et al. | 3DMAX-Net: A multi-scale spatial contextual network for 3D point cloud semantic segmentation | |
| CN103675908A (en) | Wave-equation reverse-time migration imaging method for mass-data graphic processing unit | |
| CN105005072A (en) | PML boundary three-dimensional seismic wave propagation simulation method utilizing CUDA | |
| CN102682476B (en) | The boolean operation method of triangle grid data and system thereof | |
| CN108460195A (en) | Tsunami mathematical calculation model is based on rapid implementation method parallel GPU | |
| CN109272029A (en) | The extensive spectral clustering seismic facies analysis method of well control sparse representation | |
| CN111638551A (en) | Seismic first-motion wave travel time chromatography method and device | |
| Carr et al. | Scalable contour tree computation by data parallel peak pruning | |
| CN104570080A (en) | Multi-GPU-card cooperative and quick calculation method for pre-stack reverse time migration of mass data | |
| CN105894439A (en) | Ocean eddy and Argo buoy intersection data rapid extraction algorithm based on CUDA | |
| CN106484532B (en) | GPGPU parallel calculating method towards SPH fluid simulation | |
| CN102722653B (en) | Ray tracking and accelerating algorithm based on MapReduce | |
| Zellmann et al. | Design and Evaluation of a GPU Streaming Framework for Visualizing Time-Varying AMR Data. | |
| CN109657197A (en) | A kind of pre-stack depth migration calculation method and system | |
| CN104199105A (en) | Pickup method and device of first arrivals in high-density three-dimensional exploration | |
| CN106662665B (en) | The interpolation and convolution of rearrangement for the processing of faster staggered-mesh | |
| CN106646593B (en) | A kind of three-dimensional relief surface Acoustic Forward Modeling method that cross-node is parallel |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |