CN105572728A - Reverse illumination speed inversion method based on least-square objective functional - Google Patents
Reverse illumination speed inversion method based on least-square objective functional Download PDFInfo
- Publication number
- CN105572728A CN105572728A CN201410547340.6A CN201410547340A CN105572728A CN 105572728 A CN105572728 A CN 105572728A CN 201410547340 A CN201410547340 A CN 201410547340A CN 105572728 A CN105572728 A CN 105572728A
- Authority
- CN
- China
- Prior art keywords
- grid
- ray
- depth
- velocity
- degree
- 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.)
- Pending
Links
- 238000005286 illumination Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000011159 matrix material Substances 0.000 claims abstract description 24
- 238000003384 imaging method Methods 0.000 claims abstract description 13
- 230000003459 anti-dromic effect Effects 0.000 claims description 42
- 238000013508 migration Methods 0.000 claims description 13
- 230000005012 migration Effects 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 3
- 238000005452 bending Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
Landscapes
- Image Generation (AREA)
Abstract
A reverse illumination speed inversion method based on a least-square objective functional is disclosed. The method comprises the following steps of calculating a reverse illumination matrix; according to top boundary and bottom boundary positions of a current layer position, using a bending degree of a co-phase shaft on an angle domain common imaging point gather to calculate a residual depth; under least-square sense, through searching a speed model, using a minimum difference of actual observation travel time and forward modeling travel time to construct the least-square objective functional. The method has advantages that the method can be compatible with current processing software and implementation is convenient; and calculating efficiency and calculating precision are high. In the method, the angle domain common imaging point gather without illusion is used to calculate the residual depth so that the calculating precision is high; ray tracing from a bottom interface of a reflecting layer to a top interface is used to calculate the reverse illumination matrix and reflection ray tracking with tedious calculating is not needed so that the calculating efficiency is high. By using the method, a requirement to initial speed model precision is low, an inversed non-linear degree is low and adaptability is high.
Description
Technical field
The invention belongs to oil-gas exploration seismic data process field, is that a kind of earthquake data before superposition that is applicable to is for realizing the method for Depth Domain velocity modeling.
Prior art
Along with the develop rapidly of seismic exploration technique and computer level, pre-stack depth migration has become main force's technology of complex hydrocarbon exploratory area high precision imaging.But, pre-stack depth migration is very strong to the susceptibility of velocity accuracy, if Depth Domain rate pattern is inaccurate, offset method accurate, advanced so more all can not carry out accurate imaging to underground structure, therefore, a good Depth Domain velocity modeling method seems most important.
At present, existing Depth Domain velocity modeling method mainly contains three kinds: one utilizes time domain speed to carry out time and depth transfer, but all there is larger error due to conversion formula, the general method as acquisition initial depth territory rate pattern; The second is full waveform inversion, and its modeling accuracy is high, but it is high to initial velocity accuracy requirement, in refutation process, nonlinear degree is strong, and needs to estimate source wavelet etc., and calculated amount is huge, the rarely seen example that should make good use of, is difficult to popularization and carries out large-scale commercial production; The third is migration velocity analysis, and its advantage is moderate accuracy, and counting yield can accept, and current most of business softwares all adopt the method to carry out Depth Domain velocity modeling.But migration velocity analysis uses offset domain common image gathers, produce in the process of this road collection and there is multipath imaging problem, larger to follow-up velocity modeling process influence.
Summary of the invention
The object of the invention is for the problems referred to above, propose a kind of antidromic illumination velocity inversion method based on least square cost functional, more and more higher requirement seismic processing technique proposed with the increase meeting oil-gas exploration and development difficulty.Pre-stack depth migration is as main force's technology of high precision imaging, and it key of large-scale promotion application can be the accurate foundation of Depth Domain rate pattern.
The present invention is under Least squares inversion framework, set up the cost functional of Depth Domain rate pattern inverting, adopt current reflective layer bottom boundary to interface, top antidromic illumination matrix by the residue depth assignment that produces in communication process to the grid of wave path process to obtain residual velocity, iterative inversion is repeatedly to obtain high accuracy depth domain rate pattern.The angle domain common image gathers that the present invention applies without illusion asks for the residue degree of depth, and reflection horizon bottom boundary asks for antidromic illumination matrix to the ray tracing at interface, top, simple to operate, counting yield and computational accuracy high, be beneficial to large-scale promotion application in actual production.
Technical scheme of the present invention comprises:
(1) antidromic illumination matrix is calculated, interface, top and the bottom boundary of current layer position is determined in pre-stack depth migration result, the starting point of ray tracing fixes on the bottom boundary of current layer, given ray exit direction, given fixed step size carries out ray tracing, and the path addition in each grid is obtained the total length of ray in each grid, total length in each grid is lined up smoothly according to grid, obtain the antidromic illumination vector of a ray, the antidromic illumination vector of each ray is lined up according to ray order, obtains antidromic illumination matrix;
(2) according to interface, top and the bottom boundary position of current layer position, utilize the degree of crook of lineups on angle domain common image gathers to calculate the residue degree of depth;
(3) under least square meaning, by finding a rate pattern, such that the difference of actual observation whilst on tour and forward simulation whilst on tour is minimum builds least square cost functional.
Such scheme comprises further:
The raypath computing formula of step (1) described calculating antidromic illumination matrix is:
Wherein, p
i,jbe the path in a grid, dl is fixed step size, l
1, l
2be respectively in grid and start and terminate path; n
0be ray exit direction vector, v is grid local velocity, and λ is velocity gradient
The spread pattern of antidromic illumination matrix P is:
Wherein, i ∈ [1, M] represents ray sequence number, and M is ray sum; J ∈ [1, N] represents grid numbering, and N is grid sum.
Step (2) described calculating residue degree of depth method: the reference mark of first selected velocity inverting, then utilizes the maximum energy value of common imaging gather road, energy cross-correlation method pick-up angles territory collection lineups, obtain the residue depth delta z of different angles:
Wherein, energy (τ
max) be cross-correlation maximum energy value, τ is depth-sampling point amount of movement, and span is [1, kd
down-kd
up+ 1]; F (z) for the degree of depth be the offset data energy value at z place, the span of z is top, the end degree of depth [d
up, d
down]; Max [] is max function, z
0for the lineups degree of depth that zero degree in angle domain common image gathers is corresponding.
Step (3) builds least square cost functional and is expressed from the next:
Wherein, E (s) is least square cost functional, and s represents grid slowness, is the inverse of Grid Velocity; Δ s represents slowness disturbance, for carrying out velocity inversion renewal; P represents antidromic illumination matrix, and Δ z represents the residue degree of depth; θ
1for ray emergence angle, θ
2for bottom boundary inclination angle, current layer position; || || the L2 mould of representation vector
Its Grid Velocity inverting more new formula is:
Wherein, v
n-1and v
nbe respectively the velocity vector that n-th (n>=1) inverting upgrades front and back.
Antidromic illumination velocity inversion method based on least square cost functional of the present invention, its technical advantage is mainly manifested in following three aspects:
The first, method and current process software can be good at compatibility and implement conveniently.The Depth Domain rate pattern that the time domain rate pattern that the method utilizes current conventional process software to provide is become by time and depth transfer inputs as initial velocity model, the pre-stack depth migration module when pre-processing software is utilized to obtain angle domain common image gathers, compatible good, it is convenient to implement;
The second, counting yield and computational accuracy high.The method application asks for the residue degree of depth without the angle domain common image gathers of illusion, and computational accuracy is high; Application reflection horizon bottom boundary asks for antidromic illumination matrix to the ray tracing at interface, top, and do not need to calculate loaded down with trivial details reflected ray and follow the trail of, counting yield is high;
Three, the requirement of method to initial velocity model precision is low, and the nonlinear degree of inverting is weak, strong adaptability.
Accompanying drawing explanation
Fig. 1 is based on the antidromic illumination velocity inversion method operational flowchart of least square cost functional
The complicated rate pattern of Fig. 2
The input initial velocity model of Fig. 3 complex model
Fig. 4 antidromic illumination inversion speed model
The pre-stack depth migration result that Fig. 5 inversion speed model is corresponding
The input initial velocity model of Fig. 6 real data
The rate pattern that the inverting of Fig. 7 real data antidromic illumination obtains
The pre-stack depth migration result that Fig. 8 Real data inversion rate pattern is corresponding
Embodiment
Below by specific embodiment, technical scheme of the present invention is described further.
The detailed technology operation steps of this invention as shown in Figure 1.This invention major technique key point is following three: (1), calculating antidromic illumination matrix; (2) the residue degree of depth, is calculated; (3), least square cost functional is built.
(1) antidromic illumination matrix is calculated
Interface, top and the bottom boundary of current layer position is determined in pre-stack depth migration result, the starting point of ray tracing fixes on the bottom boundary of current layer, given ray exit direction, given fixed step size carries out ray tracing, and the path addition in each grid is obtained the total length of ray in each grid, total length in each grid is lined up smoothly according to grid, obtain the antidromic illumination vector of a ray, the antidromic illumination vector of each ray is lined up according to ray order, obtains antidromic illumination matrix.Raypath computing formula is:
Wherein, p
i,jbe the path in a grid, dl is fixed step size, l
1, l
2be respectively in grid and start and terminate path; n
0be ray exit direction vector, v is grid local velocity, and λ is velocity gradient.
The spread pattern of antidromic illumination matrix P is:
Wherein, i ∈ [1, M] represents ray sequence number, and M is ray sum; J ∈ [1, N] represents grid numbering, and N is grid sum.
(2) the residue degree of depth is calculated
According to interface, top and the bottom boundary position of current layer position, utilize the degree of crook of lineups on angle domain common image gathers to calculate the residue degree of depth.First the reference mark (according to accuracy requirement, the amount doesn't matter) of selected velocity inverting, then utilizes the maximum energy value of common imaging gather road, energy cross-correlation method pick-up angles territory collection lineups, obtains the residue depth delta z of different angles:
Wherein, energy (τ
max) be cross-correlation maximum energy value, τ is depth-sampling point amount of movement, and span is [1, kd
down-kd
up+ 1]; F (z) for the degree of depth be the offset data energy value at z place, the span of z is top, the end degree of depth [d
up, d
down]; Max [] is max function, z
0for the lineups degree of depth that zero degree in angle domain common image gathers is corresponding.
(3) least square cost functional is built
Under least square meaning, can searching rate pattern be passed through, make the difference of actual observation whilst on tour and forward simulation whilst on tour minimum come establishing target functional:
Wherein, E (s) is least square cost functional, and s represents grid slowness, is the inverse of Grid Velocity; Δ s represents slowness disturbance, for carrying out velocity inversion renewal; P represents antidromic illumination matrix, and Δ z represents the residue degree of depth; θ
1for ray emergence angle, θ
2for bottom boundary inclination angle, current layer position; || || the L2 mould of representation vector.
Its Grid Velocity inverting more new formula is:
Wherein, v
n-1and v
nbe respectively the velocity vector that n-th (n>=1) inverting upgrades front and back.
Below for complex model and certain exploratory area real data processing procedure, the effect of this invention is described.
Fig. 2-Fig. 5 is the test result of this invention to complex model.Fig. 2 is real complicated rate pattern, and discrete grid block is horizontal, longitudinal sampling number is 801*638, and corresponding sampling interval is respectively 10m and 20m, and this model exists the complex structures such as high steep dip, tomography and thin layer.
This invention to the implementation step of complex model is:
First, the initial rate pattern of input shown in Fig. 3, antidromic illumination matrix is calculated: the bottom boundary starting point of ray tracing being fixed on ground floor according to summary of the invention (1), given ray exit direction and fixed step size carry out ray tracing, path in each grid is added and obtains the total length of ray in each grid, total length in each grid is lined up smoothly according to grid, obtain the antidromic illumination vector of a ray, the antidromic illumination vector of each ray is lined up according to ray order, obtains antidromic illumination matrix;
The residue degree of depth is calculated: according to 80 reference mark of data sampling number selected velocity inverting according to summary of the invention (2), then utilize the maximum energy value of common imaging gather road, energy cross-correlation method pick-up angles territory collection lineups, obtain the residue degree of depth of different angles;
Then the least square cost functional as shown in formula (4) is built, formula (5) is utilized to carry out velocity inversion renewal, ground floor iteration ends starts the inverting of the second layer, until all layer positions iteration ends exports final inversion speed model (as shown in Figure 4).Can find out: rate pattern (Fig. 4) same to true model (Fig. 1) goodness of fit that inverting obtains is very high, the display of pre-stack depth migration imaging result in Figure 5, visible inverse model precision is higher, can obtain high-quality imaging section for carrying out skew.
Fig. 6-Fig. 8 is the test result of this invention to real data.Real data discrete grid block is horizontal, longitudinal sampling number is 901*1001, and corresponding sampling interval is respectively 15m and 8m.Implementation step illustrates: the initial velocity model of input shown in Fig. 6, first antidromic illumination matrix is calculated: the bottom boundary starting point of ray tracing being fixed on ground floor, given ray exit direction and fixed step size carry out ray tracing, path in each grid is added and obtains the total length of ray in each grid, total length in each grid is lined up smoothly according to grid, obtain the antidromic illumination vector of a ray, the antidromic illumination vector of each ray is lined up according to ray order, obtains antidromic illumination matrix; Secondly the residue degree of depth is calculated: according to 90 reference mark of data sampling number selected velocity inverting, then utilize the maximum energy value of common imaging gather road, energy cross-correlation method pick-up angles territory collection lineups, obtain the residue degree of depth of different angles; Then the least square cost functional as shown in formula (4) is built, formula (5) is utilized to carry out velocity inversion renewal, ground floor iteration ends starts the inverting of the second layer, until all layer positions iteration ends exports final inversion speed model, display in the figure 7.Can see: from shallow-layer to deep layer, speed details is obtained for well portrays, and the position, pre-stack depth migration result middle level shown in Fig. 8 is clearly demarcated, tomography is clear, and the velocity inversion precision demonstrating this invention is higher.
Claims (4)
1., based on an antidromic illumination velocity inversion method for least square cost functional, it is characterized in that comprising:
(1) antidromic illumination matrix is calculated, interface, top and the bottom boundary of current layer position is determined in pre-stack depth migration result, the starting point of ray tracing fixes on the bottom boundary of current layer, given ray exit direction, given fixed step size carries out ray tracing, and the path addition in each grid is obtained the total length of ray in each grid, total length in each grid is lined up smoothly according to grid, obtain the antidromic illumination vector of a ray, the antidromic illumination vector of each ray is lined up according to ray order, obtains antidromic illumination matrix;
(2) according to interface, top and the bottom boundary position of current layer position, utilize the degree of crook of lineups on angle domain common image gathers to calculate the residue degree of depth;
(3) under least square meaning, by finding a rate pattern, such that the difference of actual observation whilst on tour and forward simulation whilst on tour is minimum builds least square cost functional.
2. the antidromic illumination velocity inversion method based on least square cost functional according to claim 1, is characterized in that, the raypath computing formula of step (1) described calculating antidromic illumination matrix is:
Wherein, p
i,jbe the path in a grid, dl is fixed step size, l
1, l
2be respectively in grid and start and terminate path; n
0be ray exit direction vector, v is grid local velocity, and λ is velocity gradient
The spread pattern of antidromic illumination matrix P is:
Wherein, i ∈ [1, M] represents ray sequence number, and M is ray sum; J ∈ [1, N] represents grid numbering, and N is grid sum.
3. the antidromic illumination velocity inversion method based on least square cost functional according to claim 2, it is characterized in that, step (2) described calculating residue degree of depth method: the reference mark of first selected velocity inverting, then utilize the maximum energy value of common imaging gather road, energy cross-correlation method pick-up angles territory collection lineups, obtain the residue depth delta z of different angles:
Wherein, energy (τ
max) be cross-correlation maximum energy value, τ is depth-sampling point amount of movement, and span is [1, kd
down-kd
up+ 1]; F (z) for the degree of depth be the offset data energy value at z place, the span of z is top, the end degree of depth [d
up, d
down]; Max [] is max function, z
0for the lineups degree of depth that zero degree in angle domain common image gathers is corresponding.
4. the antidromic illumination velocity inversion method based on least square cost functional according to claim 1,2 or 3, is characterized in that, step (3) builds least square cost functional and is expressed from the next:
Wherein, E (s) is least square cost functional, and s represents grid slowness, is the inverse of Grid Velocity; Δ s represents slowness disturbance, for carrying out velocity inversion renewal; P represents antidromic illumination matrix, and Δ z represents the residue degree of depth; θ
1for ray emergence angle, θ
2for bottom boundary inclination angle, current layer position; || || the L2 mould of representation vector
Its Grid Velocity inverting more new formula is:
Wherein, v
n-1and v
nbe respectively the velocity vector that n-th (n>=1) inverting upgrades front and back.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410547340.6A CN105572728A (en) | 2014-10-16 | 2014-10-16 | Reverse illumination speed inversion method based on least-square objective functional |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410547340.6A CN105572728A (en) | 2014-10-16 | 2014-10-16 | Reverse illumination speed inversion method based on least-square objective functional |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105572728A true CN105572728A (en) | 2016-05-11 |
Family
ID=55883069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410547340.6A Pending CN105572728A (en) | 2014-10-16 | 2014-10-16 | Reverse illumination speed inversion method based on least-square objective functional |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105572728A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106959467A (en) * | 2017-03-20 | 2017-07-18 | 中国科学院地质与地球物理研究所 | Seimic wave velocity inversion method and device |
CN110023790A (en) * | 2016-12-02 | 2019-07-16 | Bp北美公司 | Earthquake-capturing geometry full waveform inversion |
CN111190224A (en) * | 2020-01-09 | 2020-05-22 | 中国石油大学(华东) | Dynamic sampling full-waveform inversion system and method based on three-dimensional seismic wave reverse illumination |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070282535A1 (en) * | 2006-05-31 | 2007-12-06 | Bp Corporation North America Inc. | System and method for 3d frequency domain waveform inversion based on 3d time-domain forward modeling |
US20120075954A1 (en) * | 2010-09-24 | 2012-03-29 | CGGVeritas Services (U.S.) Inc. | Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration |
CN102841375A (en) * | 2012-09-06 | 2012-12-26 | 中国石油大学(华东) | Method for tomography velocity inversion based on angle domain common imaging gathers under complicated condition |
CN102866422A (en) * | 2012-09-10 | 2013-01-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Depth domain geological entity model generation method |
-
2014
- 2014-10-16 CN CN201410547340.6A patent/CN105572728A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070282535A1 (en) * | 2006-05-31 | 2007-12-06 | Bp Corporation North America Inc. | System and method for 3d frequency domain waveform inversion based on 3d time-domain forward modeling |
US20120075954A1 (en) * | 2010-09-24 | 2012-03-29 | CGGVeritas Services (U.S.) Inc. | Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration |
CN102841375A (en) * | 2012-09-06 | 2012-12-26 | 中国石油大学(华东) | Method for tomography velocity inversion based on angle domain common imaging gathers under complicated condition |
CN102866422A (en) * | 2012-09-10 | 2013-01-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Depth domain geological entity model generation method |
Non-Patent Citations (3)
Title |
---|
潘兴祥等: "叠前深度偏移层析速度建模及应用", 《地球物理学进展》 * |
秦宁等: "基于角道集的井约束层析速度反演", 《石油地球物理勘探》 * |
秦宁等: "自动拾取的成像空间域走时层析速度反演", 《石油地球物理勘探》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110023790A (en) * | 2016-12-02 | 2019-07-16 | Bp北美公司 | Earthquake-capturing geometry full waveform inversion |
CN110023790B (en) * | 2016-12-02 | 2022-03-08 | Bp北美公司 | Seismic acquisition geometric full-waveform inversion |
CN106959467A (en) * | 2017-03-20 | 2017-07-18 | 中国科学院地质与地球物理研究所 | Seimic wave velocity inversion method and device |
CN106959467B (en) * | 2017-03-20 | 2019-02-19 | 中国石油天然气集团有限公司 | Seismic wave velocity inversion method and device |
CN111190224A (en) * | 2020-01-09 | 2020-05-22 | 中国石油大学(华东) | Dynamic sampling full-waveform inversion system and method based on three-dimensional seismic wave reverse illumination |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102692644B (en) | Depth domain common-image gather generation method | |
Woodham et al. | Propagation characteristics of coastally trapped waves on the Australian Continental Shelf | |
US11079509B1 (en) | Building accurate training images for automatic seismic interpretation | |
Araya-Polo et al. | Fast and accurate seismic tomography via deep learning | |
Christian et al. | Integrating radar stratigraphy with high resolution visible stratigraphy of the north polar layered deposits, Mars | |
CN102901985B (en) | A kind of Depth Domain interval velocity modification method being applicable to relief surface | |
CN104181598A (en) | Method and device for calculating discontinuity attribute value of stratum | |
Yuan et al. | Spatial deep deconvolution u-net for traffic analyses with distributed acoustic sensing | |
CN105954802A (en) | Lithology data volume conversion method and device | |
CN103733089A (en) | System and method for subsurface characterization including uncertainty estimation | |
CN105093281A (en) | Earthquake multi-wave modeling method under inverse framework | |
CN104237937A (en) | Pre-stack seismic inversion method and system thereof | |
Sain et al. | Meta-attributes and artificial networking: A new tool for seismic interpretation | |
Muller et al. | Deep-tomography: iterative velocity model building with deep learning | |
Wu et al. | Variable seismic waveforms representation: Weak-supervised learning based seismic horizon picking | |
Weekley et al. | Tracking lake surface elevations with proportional hypsometric relationships, Landsat imagery, and multiple DEMs | |
CN105572728A (en) | Reverse illumination speed inversion method based on least-square objective functional | |
Bai et al. | Seismic wavefront evolution of multiply reflected, transmitted, and converted phases in 2D/3D triangular cell model | |
CN102866422B (en) | A kind of depth domain geological entity model generation method | |
Liu et al. | Bathymetry inversion using a deep‐learning‐based surrogate for shallow water equations solvers | |
Chen et al. | A total error-based multiquadric method for surface modeling of digital elevation models | |
CN103440358B (en) | Slope fitting method based on DEM (dynamic effect model) data | |
Chen et al. | Surface modeling of DEMs based on a sequential adjustment method | |
CN103852789A (en) | Nonlinear chromatography method and device for seismic data | |
CN107340537A (en) | A kind of method of P-SV converted waves prestack reverse-time depth migration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160511 |
|
WD01 | Invention patent application deemed withdrawn after publication |