CN104765067B - A kind of efficient pre-stack time migration velocity analysis method - Google Patents
A kind of efficient pre-stack time migration velocity analysis method Download PDFInfo
- Publication number
- CN104765067B CN104765067B CN201410003240.7A CN201410003240A CN104765067B CN 104765067 B CN104765067 B CN 104765067B CN 201410003240 A CN201410003240 A CN 201410003240A CN 104765067 B CN104765067 B CN 104765067B
- Authority
- CN
- China
- Prior art keywords
- subprocess
- time migration
- stack time
- migration velocity
- velocity analysis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The present invention provides a kind of efficient pre-stack time migration velocity analysis method, is applied to oil exploration technology field, and the method includes：Apply for multiple processes；Determine distribution information；For each subprocess storage allocation space；Host process is each subprocess distribution pre-stack time migration velocity analysis position and earthquake data before superposition road according to the distribution information；Each subprocess carries out pre-stack time migration using the pre-stack time migration velocity analysis position and the earthquake data before superposition road, generates pre-stack time migration result and is sent to the host process；Host process obtains pre-stack time migration velocity analysis result according to the pre-stack time migration result that all subprocess send.The present invention takes full advantage of computer cluster technology, pre-stack time migration velocity analysis is calculated and is assigned to multiple subprocess up, calculating speed can be significantly improved, shortens the activity duration, the 3-d seismic exploration work area being particularly suited for more velocity analysis position.
Description
Technical field
A kind of the present invention relates to oil exploration technology field, in particular it relates to efficient pre-stack time migration velocity analysis side
Method.
Background technology
During oil-gas resource Seismic Exploration Data Processing, pre-stack time migration is conventional imaging method,
Conventional treatment means are had become as, it is required parameter to carry out pre-stack time migration hour offset velocity field.So, prestack time is inclined
Move velocity analysis and build an instrument whether efficiently, convenient, practical, be pre-stack time migration can popularize decisive because
Element.
Existing pre-stack time migration velocity analysis method mainly has constant-velocity scan method, progressively revised law and divide along interval velocity
Three kinds of analysis method, but these three methods have the shortcomings that computationally intensive, process step is many, it is long to expend the cycle.In order to overcome this lacking
Point, occurs in that other pre-stack time migration velocity analysis methods, such as Application No. in prior art again
200910080622.9 Chinese patent proposes a kind of method for determining pre-stack time migration speed, and the method is only to needing
The information exported on velocity analysis point is calculated, and is implemented simple and is significantly reduced amount of calculation, but due to its EMS memory occupation
Amount and calculating job run time are all directly proportional to the number of velocity analysis location point, for more velocity analysis position
Large-area three-dimensional seismic prospecting work area, its application will be restricted.For example, if there is hundreds of normal-moveout spectrum in operation, may
The internal memory of nearly 10G byte is taken, that is out the configuration of most of active computers；And such operation may run hundreds of
Hour even tens days.If occurring the fortuitous events such as computer glitch during job run, operation will all that has been achieved is spoiled, makes
The risk for becoming process contract be timely completed can not be ignored.
Analyze aforementioned various pre-stack time migration velocity analysis methods and understand that aforementioned various methods are all to be applied to single meter
Calculate machine equipment to process on basis, limited by the disposal ability of single computer equipment, therefore aforementioned various methods are also all
There is the defect of calculating cycle length to varying degrees, become the unfavorable factor for smoothly completing pre-stack time migration operation.
Content of the invention
The main purpose of the embodiment of the present invention is to provide a kind of efficient pre-stack time migration velocity analysis method, to solve
The problem of the calculating cycle length that existing pre-stack time migration velocity analysis method is present.
To achieve these goals, the embodiment of the present invention provides a kind of efficient pre-stack time migration velocity analysis method, bag
Include：
Apply for multiple processes, in the plurality of process is defined as host process, remaining is defined as subprocess；
Determine distribution information, the distribution information includes each subprocess and will distribute to the prestack time of the subprocess
Corresponding relation between migration velocity analysis position；
For each subprocess storage allocation space；
Host process is each subprocess distribution pre-stack time migration velocity analysis position and prestack according to the distribution information
Seismic-data traces；
Each subprocess is folded using the pre-stack time migration velocity analysis position and the earthquake data before superposition road
Front time migration, generates pre-stack time migration result and is sent to the host process；
The host process obtains pre-stack time migration speed according to the pre-stack time migration result that all subprocess send and divides
Analysis result, the pre-stack time migration velocity analysis result include whole pre-stack time migration velocity analysis positions and whole prestacks
The corresponding common reflection point CRP road collection of seismic-data traces, migrated section section and Migration velocity spectra.
By means of technique scheme, pre-stack time migration velocity analysis is calculated and is assigned on multiple subprocess by the present invention
Go, computer cluster technology is taken full advantage of, the speed of pre-stack time migration velocity analysis calculating is significantly improved, is significantly shortened
The activity duration that pre-stack time migration velocity analysis is calculated, the large area three that is particularly suited for more velocity analysis position
Dimension seismic prospecting work area, with preferable practicality.
Description of the drawings
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, embodiment will be described below
Needed for accompanying drawing to be used be briefly described, it should be apparent that, drawings in the following description be only the present invention some
Embodiment, for those of ordinary skill in the art, on the premise of not paying creative work, can also be attached according to these
Figure obtains other accompanying drawings.
Fig. 1 is the flow chart of the efficient pre-stack time migration velocity analysis method that the present invention is provided；
Fig. 2 is the work matching relationship schematic diagram between the host process that the present invention is provided and subprocess；
Fig. 3 be the embodiment of the present invention one provide when m- velocity plane on Migration velocity spectra Local map；
Fig. 4 is the friction speed of the offer of the embodiment of the present invention one to difference geophone offset Chu CRP road on reflection line-ups
Collection Local map；
Fig. 5 is skew offset segment Local map of the friction speed of the offer of the embodiment of the present invention one to reflection line-ups；
Fig. 6 is the Migration velocity spectra Local map calculated using partial data that the embodiment of the present invention one is provided；
Fig. 7 is the CRP road collection Local map calculated using partial data that the embodiment of the present invention one is provided；
Fig. 8 is the migrated section section Local map calculated using partial data that the embodiment of the present invention one is provided；
Fig. 9 is the pre-stack time migration velocity analysis and the process chart for building field that the embodiment of the present invention one is provided.
Specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the present invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
With the development of computer technology, computer cluster has become as most of petroleum exploration seismic data processing centers
The workbench of common configuration, it have numerous nodes, can undertake large-scale calculating task in a parallel fashion.The present invention
Using computer cluster technology, the calculating process to pre-stack time migration velocity analysis adopts parallelization operating type, improves
The speed of pre-stack time migration velocity analysis, shortens duty cycle.
As shown in figure 1, the present invention provides a kind of efficient pre-stack time migration velocity analysis method, including：
Step S11, applies for multiple processes, in the plurality of process is defined as host process, remaining is defined as son
Process.
Specifically, the step is applied multiple in computer cluster（At least three）The resource of process, one of conduct
Host process, remaining is used as subprocess；Wherein, host process is used for being responsible for distributing data to subprocess, reclaims the calculating knot of subprocess
Really, the data that subprocess is used for issuing using host process carry out pre-stack time migration calculating.It is illustrated in figure 2 master in the present invention to enter
Work matching relationship between journey and subprocess.
The step can estimate Activity Calculation amount in advance, and consider the process use of current computer cluster device
Situation, then determine altogether and apply for how many processes.
Step S12, determines distribution information, and the distribution information includes each subprocess and will distribute to the subprocess
Corresponding relation between pre-stack time migration velocity analysis position.
In order to the processing speed for ensureing each subprocess is suitable, to be rapidly completed calculating task and not waste calculating money
Source, in the step, can be to the suitable pre-stack time migration velocity analysis position of each subprocess distribution number.If prestack time
The quantity of migration velocity analysis position is the integral multiple of subprocess quantity, and the step is can determine that whole pre-stack time migrations speed
Degree analysis position is evenly distributed to each subprocess, and is determined each subprocess and will be distributed to the prestack time of the subprocess
Corresponding relation between migration velocity analysis position.If the quantity of pre-stack time migration velocity analysis position is not subprocess number
The integral multiple of amount, then the velocity analysis number of positions difference on different subprocess is less than one.
Step S13, is each subprocess storage allocation space.
Specifically, the step is directed to each subprocess, according to the corresponding pre-stack time migration velocity analysis position of the subprocess
Put quantity, the sampling point quantity in earthquake data before superposition road and default pre-stack time migration calculating parameter and memory headroom is calculated, and
Distribute to the subprocess.
When the pre-stack time migration velocity analysis position data output for distributing to each subprocess is roughly the same, enter per height
The pre-stack time migration that carries out in journey calculate needed for memory headroom mainly with the sampling point quantity in earthquake data before superposition road and preset
Pre-stack time migration calculating parameter related.
Step S14, host process distribute pre-stack time migration velocity analysis position according to the distribution information for each subprocess
Put and earthquake data before superposition road.
In the step, host process can be divided into the continuous and two kinds of shapes in interval for each subprocess distribution earthquake data before superposition road
Formula, wherein, conitnuous forms refer to that whole earthquake data before superposition roads are distributed on the subprocess by host process, and alternate form is referred to
The subprocess is distributed in part earthquake data before superposition road by host process.
During using conitnuous forms, the data input amount that host process is sent to each subprocess is bigger, and the consuming time is longer, but
Pre-stack time migration velocity analysis result obtained by subprocess can be truer；During using alternate form, host process is sent to every
The data input amount of individual subprocess is less, and the consuming time is shorter, but the pre-stack time migration velocity analysis knot obtained by subprocess
The signal to noise ratio of fruit can decline.When being embodied as, it is long to amount of calculation size, calculating time that operating personnel can weigh actual job
Short, noise compares so as to see who is superior etc. and to require, determines and selects conitnuous forms or alternate form.
In a kind of preferred embodiment, when alternate form is selected, in order to ensure to distribute to the geological data of subprocess
Road truly reflects geological condition as far as possible, and host process sets quantity by first（As survey line, CMP CMP road collection, 1,000
Road, 10,000 roads etc.）Continuous earthquake data before superposition road is defined as one group of earthquake data before superposition road；In whole earthquake data before superposition roads
In, the continuous earthquake data before superposition road group of quantity is set every second, institute is distributed in subsequent group earthquake data before superposition road
Subprocess is stated, such as host process often skips several earthquake data before superposition road groups in whole earthquake data before superposition roads, will be follow-up
An earthquake data before superposition road group distribute to the subprocess.
Step S15, each subprocess is using the pre-stack time migration velocity analysis position and the earthquake data before superposition
Road carries out pre-stack time migration, generates pre-stack time migration result and is sent to the host process.
Step S16, host process obtain pre-stack time migration speed according to the pre-stack time migration result that all subprocess send
Degree analysis result, the pre-stack time migration velocity analysis result include whole pre-stack time migration velocity analysis positions and all
The corresponding common reflection point CRP road collection in earthquake data before superposition road, migrated section section and Migration velocity spectra.
In step S15 and S16, each subprocess executed in parallel is calculated, and each subprocess is distributed to itself to host process and folded
Front time migration velocity analysis position and earthquake data before superposition road carry out pre-stack time migration, by the parallel work of all subprocess
Industry, completes the calculating task in whole pre-stack time migration velocity analysis positions and whole earthquake data before superposition roads on the whole, carries
The high efficiency of pre-stack time migration velocity analysis, shortens duty cycle；Finally, then collect all subprocess by host process to send out
The pre-stack time migration result result that sends, obtains pre-stack time migration velocity analysis result, i.e., all pre-stack time migration is fast
Degree analysis position and the corresponding CRP road collection in whole earthquake data before superposition roads, migrated section section and Migration velocity spectra.
CRP road collection, migrated section section and the Migration velocity spectra obtained using the present invention, it is possible to determine one optimal folded
Front time migration speed, as the present invention is remarkably improved the calculating speed of pre-stack time migration velocity analysis, duty cycle is short,
The aborted risk of calculating process will reduce, so as to for folded during oil-gas resource Seismic Exploration Data Processing
Front time migration provides advantage.
In a kind of preferred embodiment, step S15 is specific as follows：Each subprocess determines pre-stack time migration speed sequence
Row, using the every velocity function in the pre-stack time migration velocity series as pre-stack time migration speed, divide to host process
The pre-stack time migration velocity analysis position that joins and earthquake data before superposition road carry out pre-stack time migration, generate the subprocess and correspond to
CRP road collection, migrated section section and migration velocity modal data.In the embodiment, the calculating process that each subprocess is carried out is as follows：
Every velocity function in pre-stack time migration velocity series is scanned, to scan the velocity function for arriving as prestack time
Migration velocity, based on the pre-stack time migration speed to distributing to pre-stack time migration velocity analysis position and the prestack ground of itself
Shake data track carries out pre-stack time migration calculating, and all speed in complete pre-stack time migration velocity series to be scanned are obtained
The pre-stack time migration velocity analysis position that distributes corresponding to the pre-stack time migration velocity series and host process and prestack
The CRP road collection, migrated section section of seismic-data traces, migration velocity modal data, wherein, migration velocity modal data be for host process
The basic data information needed for Migration velocity spectra is generated in subsequent step.Additionally, the task of host process is concrete in the embodiment
As follows：The CRP road collection, migrated section section of all subprocess calculating of synthesis, obtains whole pre-stack time migration velocity analysis
The corresponding CRP road collection in position and whole earthquake data before superposition roads, migrated section section, and the skew calculated according to all subprocess
The corresponding migration velocity in the whole pre-stack time migration velocity analysis positions of normal-moveout spectrum data genaration and whole earthquake data before superposition roads
Spectrum.
Center trace in view of migrated section section can be obtained according to the direct read group total of CRP road collection, and Migration velocity spectra
If the arterial highway information of data also can be obtained according to the direct read group total of CRP road collection, in order to avoid computing repeatedly, improve and calculate effect
Rate, in another kind of preferred embodiment, can by total to CRP road collection and migrated section section and migration velocity modal data, can
Out, used as co-information, each subprocess is in concrete calculating process for the information flag obtained according to the direct read group total of CRP road collection
In no longer calculate complete migrated section section and migration velocity modal data, but only calculate CRP road collection and migrated section section and
In migration velocity modal data, removing lacks the part beyond co-information；Then, each subprocess by corresponding CRP road collection, lack
The migrated section section and Migration velocity spectra data is activation of co-information is to host process；Directly sued for peace using CRP road collection by host process again
Co-information is calculated, and it is corresponding each subprocess to be obtained using co-information and the migrated section section combination for lacking co-information
Complete migrated section section, and each subprocess pair is obtained using co-information with the Migration velocity spectra data combination for lacking co-information
The complete migration velocity modal data that answers；Finally, host process is also by corresponding for all subprocess CRP road collection, complete migrated section
Section synthesis is corresponding in whole pre-stack time migration velocity analysis positions and whole earthquake data before superposition roads together, is obtained
CRP road collection, migrated section section, and according to the corresponding complete migration velocity modal data of all subprocess, generate whole prestack time
The corresponding Migration velocity spectra in migration velocity analysis position and whole earthquake data before superposition roads.
For example, in one embodiment, pre-stack time migration velocity series, sequence number are become using 61 sweep speed group of functions
Respectively -30, -29 ..., -1,0,1 ..., 29,30.9 CRP road collections and 9 migrated section sections are had, is respectively adopted sequence number
Sweep speed function for -8, -6, -4, -2,0,2,4,6,8 is calculated.11 road of each migrated section Duan You, represents velocity analysis position
Put and before, each 5 road afterwards.Subprocess carry out pre-stack time migration calculate when, from sequence number -30 to 30 cyclically use more than
Sweep speed function, the sweep speed function for serial number -8, -6, -4, -2,0,2,4,6,8 are skipped and are not calculated, and are cutd open in skew
Center trace is skipped in the calculating of face section（Velocity analysis position）Do not calculate.After the completion of the calculating of subprocess, host process is to CRP road collection
Summation, just defines the center trace of 9 migrated section sections, and be further converted into vacancy on normal-moveout spectrum, serial number -8, -
6th, the result of calculation of -4, -2,0,2,4,6,8 sweep speed function.This measure can make sweep speed pre-stack time migration
Amount of calculation reduce about 15%.Each numeral is only merely illustrative above, and treatment people can specify others number according to actual conditions
Word.
Due to significantly improving the calculating speed of pre-stack time migration velocity analysis, shorten duty cycle, therefore this
The bright large-area three-dimensional seismic prospecting work area that can be applied to more velocity analysis position.
Embodiment one
The optimum velocity section for pre-stack time migration calculating side that the present embodiment is provided with the specific embodiment explanation present invention
The beneficial effect of method.In the present embodiment, it is 553-2130 that the survey line scope of input data is 871-900, CMP scope, has
4600323 roads, sampling interval 4ms, 2000 sampling points of per pass, data scale 35115.23MB, this are data in three-dimensional work area
A part.Velocity analysis position 42, is 151 using the original mode job run time:16:11, using 6 nodes, 24 are entered
The parallel mode job run time of journey is 9:43:36, it shorten to the 6.43% of original mode；And each node needs only to 7
The internal memory scale of individual normal-moveout spectrum.1 survey line is input into using per 3 surveys line, 1 CMP is input into per 3 CMP, same using 6 sections
The parallel mode job run time of 24 processes of point is 1:02:14, it shorten to the 0.69% of original mode further.Fig. 3 is for adopting
With the normal-moveout spectrum fully entered at the survey line 890, CMP1500 that data are calculated（Locally）, Fig. 4 is that some sweep speeds are corresponding
CRP road collection（Locally）, Fig. 5 is the corresponding migrated section section of some sweep speeds（Locally）.Fig. 6 is to be input into numbers using 1/9th
According to the normal-moveout spectrum at the survey line 890, CMP1500 for calculating（Locally）, Fig. 7 is the corresponding CRP road collection of some sweep speeds（Locally）,
Fig. 8 is the corresponding migrated section section of some sweep speeds（Locally）.From the contrast of corresponding map as can be seen that adopting 1/9th
The result that input data is calculated is essentially identical, only fine distinction, illustrates that this mode is feasible.
Except using in addition to parallel computation, the pre-stack time migration velocity analysis of the present invention and the handling process for building field and folded
Acceleration analysis with to build field closely similar, as shown in Figure 9.Therefore, the treatment people with common seismic data processing experience
Pre-stack time migration speed parameter is can determine, makes pre-stack time migration easily facilitate large-scale promotion use.
Particular embodiments described above, has been carried out to the purpose of the present invention, technical scheme and beneficial effect further in detail
Describe in detail bright, the specific embodiment that the foregoing is only the present invention is should be understood that, the guarantor being not intended to limit the present invention
Shield scope, all any modification, equivalent substitution and improvement that within the spirit and principles in the present invention, is done etc., should be included in this
Within the protection domain of invention.
Claims (5)
1. a kind of pre-stack time migration velocity analysis method, it is characterised in that include：
Apply for multiple processes, in the plurality of process is defined as host process, remaining is defined as subprocess；
Determine distribution information, the distribution information includes each subprocess and will distribute to the pre-stack time migration of the subprocess
Corresponding relation between velocity analysis position；
For each subprocess storage allocation space；
Host process is each subprocess distribution pre-stack time migration velocity analysis position and pre-stack seismic according to the distribution information
Data track；
Each subprocess executed in parallel is calculated, and each subprocess is using the pre-stack time migration velocity analysis position and the prestack
Seismic-data traces carry out pre-stack time migration velocity analysis, generate pre-stack time migration velocity analysis result and are sent to the master
Process；
It is inclined that the host process obtains total prestack time according to the prestack time velocity analysis migration result that all subprocess send
Velocity analysis result is moved, total pre-stack time migration velocity analysis result includes whole pre-stack time migration velocity analysis positions
Put and the corresponding common reflection point CRP road collection in whole earthquake data before superposition roads, migrated section section and Migration velocity spectra；
Described host process distributes the earthquake data before superposition road for each subprocess, specially：The host process all will be folded
Front seismic-data traces are distributed on the subprocess；Or, part earthquake data before superposition road is distributed to described by the host process
Subprocess；
The subprocess is distributed in part earthquake data before superposition road by wherein described host process, specially：
Continuous for first setting quantity earthquake data before superposition road is defined as one group of earthquake data before superposition road by the host process；
In whole earthquake data before superposition roads, the continuous earthquake data before superposition road group of quantity is set every second, by follow-up one
The subprocess is distributed in group earthquake data before superposition road.
2. method according to claim 1, it is characterised in that described determination distribution information, specially：
According to the subprocess quantity of application, determine and whole pre-stack time migration velocity analysis positions are evenly distributed to enter to each height
Journey.
3. method according to claim 1, it is characterised in that described for each subprocess storage allocation space, specifically
For：
For each subprocess, according to the corresponding pre-stack time migration velocity analysis number of positions of the subprocess, pre-stack seismic number
Memory headroom is calculated according to the sampling point quantity and default pre-stack time migration calculating parameter in road, and distribute to the subprocess.
4. method according to claim 1, it is characterised in that each described subprocess utilizes the pre-stack time migration
Velocity analysis position and the earthquake data before superposition road carry out pre-stack time migration velocity analysis, generate pre-stack time migration speed
Analysis result, specially：
Each subprocess described determines pre-stack time migration velocity series, with the items in the pre-stack time migration velocity series
Velocity function as pre-stack time migration speed, the pre-stack time migration velocity analysis position to host process distribution and
The earthquake data before superposition road carries out pre-stack time migration velocity analysis, generates the corresponding CRP road collection of the subprocess, migrated section
Section and migration velocity modal data.
5. method according to claim 4, it is characterised in that it is corresponding that each described subprocess generates the subprocess
CRP road collection, migrated section section and migration velocity modal data, specially：
Each subprocess described generates the corresponding CRP road collection of the subprocess, the migrated section section for lacking co-information and lacks altogether
The migration velocity modal data of information, wherein, the co-information is common with migrated section section and migration velocity modal data for CRP road collection
Information that is having and being obtained using the direct read group total of CRP road collection；
Then the host process obtains total prestack time according to the pre-stack time migration velocity analysis result that all subprocess send
Migration velocity analysis result, specially：
For each subprocess, the host process is obtained described common using the direct read group total of the corresponding CRP road collection of the subprocess
Information, and using the co-information and the migrated section section for lacking co-information and the migration velocity for lacking co-information
Modal data, obtains the corresponding migrated section section of the subprocess and migration velocity modal data；
The host process comprehensively corresponding CRP road collection of all subprocess, migrated section section and migration velocity modal data, obtain all
The corresponding CRP road collection in pre-stack time migration velocity analysis position and whole earthquake data before superposition roads, migrated section section and skew speed
Degree spectrum.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410003240.7A CN104765067B (en) | 2014-01-03 | 2014-01-03 | A kind of efficient pre-stack time migration velocity analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410003240.7A CN104765067B (en) | 2014-01-03 | 2014-01-03 | A kind of efficient pre-stack time migration velocity analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104765067A CN104765067A (en) | 2015-07-08 |
| CN104765067B true CN104765067B (en) | 2017-03-08 |
Family
ID=53647014
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410003240.7A Active CN104765067B (en) | 2014-01-03 | 2014-01-03 | A kind of efficient pre-stack time migration velocity analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104765067B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109471173B (en) * | 2018-10-08 | 2021-04-30 | 中国石油天然气集团有限公司 | Residual static correction method and device |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7388808B2 (en) * | 2004-09-23 | 2008-06-17 | Pgs Americas, Inc. | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration |
| CN101839999B (en) * | 2009-03-20 | 2011-12-07 | 中国石油集团东方地球物理勘探有限责任公司 | Method for determining optimum velocity section for pre-stack time migration |
| CN102141633B (en) * | 2010-12-10 | 2012-08-29 | 中国科学院地质与地球物理研究所 | Anisotropic three-dimensional prestack time migration method |
-
2014
- 2014-01-03 CN CN201410003240.7A patent/CN104765067B/en active Active
Non-Patent Citations (3)
| Title |
|---|
| 叠前正演模拟与偏移的网络并行计算;胡光岷等;《物探化探计算技术》;20011130;第23卷(第4期);第322页 * |
| 基于MPI+OpenMP的多层次并行偏移算法研究;陈辉 等;《成都理工大学学报（自然科学版）》;20101031;第534页 * |
| 基于共成像点道集评价的偏移速度建模;陈辉 等;《石油地球物理勘探》;20100831;第45卷(第4期);摘要 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104765067A (en) | 2015-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dai et al. | Deep learning for extracting dispersion curves | |
| CN104769459A (en) | High-efficient velocity analysis method for pre-stack time migration | |
| CN104133240B (en) | Large-scale collateral kirchhoff prestack depth migration method and device | |
| CN101598807B (en) | Method for determining attribute difference of different binning in seismological observation system design | |
| CN105005072B (en) | The PML border dimensionally seismic wave propagating mode utilizing CUDA intends method | |
| CN106353811A (en) | Method for identifying fluid of multi-attribute fusion reservoirs on basis of well information attribute optimization | |
| US20210262329A1 (en) | Method for Generating Initial Models For Least Squares Migration Using Deep Neural Networks | |
| Xue et al. | An efficient GPU implementation for locating micro-seismic sources using 3D elastic wave time-reversal imaging | |
| CN107587871A (en) | Determine the method and device of horizontal fracture width | |
| CN106405644B (en) | Crack determines method and apparatus | |
| CN105403913A (en) | Pre-stack depth migration method and device | |
| CN108415073A (en) | Angle domain back scattering offset imaging method and device | |
| CN103698809A (en) | Kirchhoff pre-stack time migration parallel method with no speed-up ratio bottleneck | |
| CN104765067B (en) | A kind of efficient pre-stack time migration velocity analysis method | |
| CN106033126B (en) | A kind of quantitative division methods of oil gas unit degree of prospecting | |
| CN107329174A (en) | The method and device of well logging information is shown on seismic profile | |
| CN106199705B (en) | The detection method of seismic acquisition data | |
| CN102692651A (en) | Preliminary wave residual static correction method with space-variant velocity | |
| CN104570115B (en) | A kind of surface wave attenuation method and device | |
| CN106526667B (en) | Prestack time migration method and device | |
| CN1292263C (en) | Ray traction in earthquake prospection | |
| CN104267434A (en) | Three-dimensional multi-component earthquake observation system geophone offset distribution obtaining method and device | |
| CN105738949B (en) | A kind of nine bin uniformity method for parallel processing for time-lapse seismic | |
| CN103217707B (en) | A kind of method of direct extraction longitudinal wave time domain transformed wave angle gathers | |
| CN106094020A (en) | A kind of seismic inversion method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |