CN104777514A - Geometric spreading compensation method based on uniform horizontal layered medium model - Google Patents
Geometric spreading compensation method based on uniform horizontal layered medium model Download PDFInfo
- Publication number
- CN104777514A CN104777514A CN201510181316.XA CN201510181316A CN104777514A CN 104777514 A CN104777514 A CN 104777514A CN 201510181316 A CN201510181316 A CN 201510181316A CN 104777514 A CN104777514 A CN 104777514A
- Authority
- CN
- China
- Prior art keywords
- layer
- information
- layered medium
- offset
- geometrical attenuation
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a geometric spreading compensation method based on a uniform horizontal layered medium model. The geometric spreading compensation method comprises the following steps: 1) reading root mean square speed information, two-way vertical traveltime information and shot and positioning point distribution information in a computer; 2) obtaining the uniform horizontal layered medium model comprising information about layer speed and layer thickness of each layer according to the root mean square speed information and the two-way vertical traveltime information; 3) performing ray tracing on the uniform horizontal layered medium model from a shot point to obtain ray position and traveltime information; 4) determining relative geometric spreading compensation values of reflecting layers according to a geometric spreading compensation formula; 5) obtaining a compensation value at detection points through linear interpolation, interpolating the relative geometric spreading compensation values to sampling points of corresponding seismic data volumes according to a seismic data format, and outputting data to be used for following data processing. The geometric spreading compensation method has the advantages that the accuracy of a geometric spreading compensation technology is improved, and can be widely applied to the geometric spreading compensation technology in the seismic data processing process.
Description
Technical field
The present invention relates to a kind of geometrical attenuation compensation method, particularly about the 2.5 dimension geometrical attenuation compensation methodes based on even level's layered medium model in a kind of seismic data process.
Background technology
Geometrical attenuation compensation technique is a ring of can not ignore in seismic data process, it is the basis of carrying out the follow-up data treatment technologies such as AVO (Amplitude Versus Offset, the change of amplitude offset distance) attributive analysis, full waveform inversion, skew further.It is whether accurate that geometrical attenuation compensation technique is applied, and directly affects the effect of other data processing and interpretation technology when applying.Therefore the necessary accuracy problems specialized in geometrical attenuation and compensate, and the key problem studying its accuracy is the problem of more accurate computational geometry diffusion.In numerical evaluation problem, the core technology that geometrical attenuation compensation technique realizes is that the account form adopted is different.Prior art and existing software compensate the geometrical attenuation that model or real data carry out, and are to compensate based on 2-D (2 dimension) the geometrical attenuation compensation formula of zero shot-geophone distance.But this compensation way has its defect:
One is compensate with zero shot-geophone distance formula, namely in whole data volume, each road of the difference along with geophone offset offset is constant, but actual conditions are offsets changes along with the change of geophone offset, especially, when geophone offset is larger, the energy of this undercompensation loss is increasing;
Two is diffusions that 2-D geometrical attenuation compensation technique does not consider outside ray plane, and therefore its accuracy is limited.For the image quality that subsequent treatment data particularly offsets, there is obvious impact.
Summary of the invention
For the problems referred to above, the object of this invention is to provide the 2.5 dimension geometrical attenuation compensation methodes based on even level's layered medium model that can improve accuracy, counting yield can be kept again in a kind of seismic data process.
For achieving the above object, the present invention takes following technical scheme: a kind of geometrical attenuation compensation method based on even level's layered medium model, and it comprises the following steps: when the root-mean-square velocity information of seismic data zoning, round trip that 1) reading in computing machine needs to compensate vertically are walked, distributed intelligence examined by information and big gun; 2) according to root-mean-square velocity information vertical with round trip walk time information, draw and comprise every interval velocity of one deck and even level's layered medium model of layer thickness information in seismic data zoning that needs compensate; 3) integrating step 1) in the big gun that obtains examine distributed intelligence, to step 2) even level's layered medium model of building, information when carrying out ray tracing to obtain ray position and to walk from shot point; 4) according to step 3) in the reflection angle in each reflection horizon that obtains of ray tracing result, determined the relative geometrical attenuation offset in each reflection horizon by geometrical attenuation compensation formula; 5) offset that two rays the most contiguous with geophone station calculate is carried out the offset that linear interpolation obtains geophone station place, according to the form of geological data, the relative geometrical attenuation offset obtained is interpolated on the sampled point of sympathetic earthquakes data volume, and exports data in order to the use of follow-up data treatment technology.
Described step 2) in, interval velocity is determined by Dix formula (1):
Wherein, V
ibe the interval velocity of i-th layer, U
ibe the root-mean-square velocity of i-th layer, U
i+1be the root-mean-square velocity of the i-th+1 layer, T
ifor the round trip of corresponding i-th layer vertically walk time, T
i+1for the round trip of corresponding the i-th+1 layer vertically walk time.
Described step 2) in, layer thickness is determined by following formula (2):
d
i＝V
i×(T
i+1-T
i)/2 (2)
Wherein, d
ibe the layer thickness of i-th layer, V
ibe the interval velocity of i-th layer, T
ifor the round trip of corresponding i-th layer vertically walk time, T
i+1for the round trip of corresponding the i-th+1 layer vertically walk time.
Described step 3) in, ray position is determined by Snell's law (3):
V
1/sinα
1＝V
2/sinβ
2(3)
Wherein, V
1for the interval velocity of ground floor, α
1for incident angle when ray incides the second layer from ground floor, V
2for the interval velocity of the second layer, β
2for ray incides the angle of transmission after the second layer from ground floor.
Described step 3) in, when walking, information is determined by following formula (4):
t
n＝d
n/(V
n·cosβ
n) (4)
Wherein, t
nfor information during the walking of n-th layer, d
nfor the layer thickness of n-th layer, V
nfor the interval velocity of n-th layer, β
nfor the angle of transmission of n-th layer.
Described step 4) in, relative geometrical attenuation offset is determined by following formula (5):
Wherein,
for the relative geometrical attenuation offset of n-th layer, V
1for the interval velocity of ground floor, α
1for incident angle when ray incides the second layer from ground floor, x is geophone offset, d
ibe the layer thickness of i-th layer, α
iit is the incident angle of i-th layer of incident ray.
The present invention is owing to taking above technical scheme, and it has the following advantages: 1, the present invention is due to for different geophone offset computational geometry diffusion compensation, improves the accuracy of geometrical attenuation compensation technique.2, the present invention is owing to describing the diffusion of underground ray accurately, further increases the accuracy of geometrical attenuation compensation technique.3, the present invention calculates at 2-D data field due to the main calculated amount of numerical evaluation, maintains counting yield.The present invention can be widely used in the geometrical attenuation compensation technique in seismic data process.
Accompanying drawing explanation
Fig. 1 is schematic flow sheet of the present invention
Fig. 2 is even level's layered medium model information that the interval velocity obtained by Dix inverting builds
Fig. 3 is the comparison diagram of North Sea model the 5th layer of normalization relative compensation value with the calculating of 2.5-D and 2-D geometrical attenuation
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in detail.
As shown in Figure 1, a kind of geometrical attenuation compensation method based on even level's layered medium model provided by the invention, it comprises the following steps:
1) in computing machine, read in information and big gun cautious (shot point and geophone station) distributed intelligence when the root-mean-square velocity information of seismic data zoning, the round trip that need to compensate vertically are walked.
2) according to root-mean-square velocity information vertical with round trip walk time information, the interval velocity drawn by Dix (Otto Dix) formula (1) inverting constructs the general frame of even level's layered medium model:
Wherein, V
ibe the interval velocity of i-th layer, U
ibe the root-mean-square velocity of i-th layer, U
i+1be the root-mean-square velocity of the i-th+1 layer, T
ifor the round trip of corresponding i-th layer vertically walk time, T
i+1for the round trip of corresponding the i-th+1 layer vertically walk time.
By root-mean-square velocity U
iand U
i+1, round trip information T when vertically walking
iand T
i+1be finally inversed by the interval velocity V of every one deck in the seismic data zoning needing to compensate
i, and information T when vertically walking according to round trip
iand T
i+1, the layer thickness d of i-th layer is determined by formula (2)
i, thus obtain comprising concrete every layer thickness of one deck and even level's layered medium model (as shown in Figure 2) of interval velocity information.
d
i＝V
i×(T
i+1-T
i)/2 (2)
3) by step 2) even level's layered medium model of building, information when carrying out ray tracing to obtain ray position and to walk from shot point:
In the even level's layered medium model established, by Snell (Si Naier) law, integrating step 1) in the shot point that obtains and geophone station distributed intelligence, ray tracing is carried out by shot point, the region of wave detector existence is covered by the ray launching sufficient amount, thus information when determining ray position and walk.
Snell law is shown below:
V
1/sinα
1＝V
2/sinβ
2(3)
Wherein, V
1for the interval velocity of ground floor, α
1for incident angle when ray incides the second layer from ground floor, V
2for the interval velocity of the second layer, β
2for ray incides the angle of transmission after the second layer from ground floor.
As shown in Figure 2, wherein, transverse axis represents geophone station position, and the longitudinal axis represents the degree of depth, V
nwith d
nrepresent interval velocity and the layer thickness of n-th layer respectively, in figure, curve represents the raypath schematic diagram this even level's layered medium model being carried out to ray tracing; α
n-1and β
nbe respectively incident angle when ray is transmitted to n-th layer from (n-1)th layer and angle of transmission, can angle of transmission β be obtained by Snell law (3)
n=argsin (sin α
n-1v
n/ V
n-1), wherein, V
nand V
n-1be respectively the interval velocity of n-th layer and (n-1)th layer; α
nand θ
nbe respectively incident angle and reflection angle, then α that ray reflects from n-th layer
n=θ
n, again due to the incident angle α of n-th layer
nwith the angle of transmission β of n-th layer
nequal, i.e. α
n=β
n=θ
n, then the ray position of n-th layer can by the incident angle α of the n-th layer of trying to achieve
ndetermine.
Information t during the walking of n-th layer
ncan by the layer thickness d of n-th layer
n, n-th layer interval velocity V
nwith the angle of transmission β of n-th layer
ndetermine, that is:
t
n＝d
n/(V
n·cosβ
n) (4)
4) the geometrical attenuation offset in each reflection horizon is calculated:
According to step 3) in the reflection angle in each reflection horizon that obtains of ray tracing result, determined the relative geometrical attenuation offset in each reflection horizon by geometrical attenuation compensation formula (5):
Wherein, x is geophone offset,
for the relative geometrical attenuation offset of n-th layer, α
iit is the incident angle of i-th layer of incident ray.
5) according to the needs of data processing technique, the relative geometrical attenuation offset that two rays the most contiguous with geophone station calculate is carried out the relative geometrical attenuation offset that linear interpolation obtains geophone station place.According to the form of geological data, namely by earthquake number of channels, the seismic data volume that sampling number and sampling interval are determined, is interpolated into the relative geometrical attenuation offset calculated on the sampled point of sympathetic earthquakes data volume, and exports data in order to the use of follow-up data treatment technology.
According to delivering data Ursin, B. (1990). " Offset-dependent geometrical spreadingin a layered medium. " GEOPHYSICS 55 (4): 492-496, give North Sea model in literary composition, its concrete speed, parameter is as shown in table 1 when walking etc.Utilize this model, respectively 2-D and 2.5-D geometrical attenuation is carried out to it and compensate.The basic parameter of big gun: track pitch is 50m, minimum geophone offset is 50m, and geophone station number is 80.The relative geometrical attenuation offset of normalization of the 5th layer of model as shown in Figure 3.
Table 1 North Sea model
As shown in Figure 3, give in figure in 2-D and 2.5-D situation, for the geometrical attenuation normalization relative compensation value comparison diagram of the 5th layer, North Sea model simultaneously.Wherein solid line represents the normalization relative compensation value of 2.5-D, the normalization relative compensation value of represented by dotted arrows 2-D.Can be found by contrast: when layer depth is determined, along with the increase of geophone offset, the growth of the relative compensation value of 2.5-D obviously will be greater than the offset in 2-D situation, and when geophone offset far away, its value even differs more than 50%.This illustrates, geophone offset be can not ignore the impact that geometrical attenuation compensates.
The various embodiments described above are only for illustration of the present invention, and certain methods wherein and step etc. all can change to some extent, and every equivalents of carrying out on the basis of technical solution of the present invention and improvement, all should not get rid of outside protection scope of the present invention.
Claims (10)
1., based on a geometrical attenuation compensation method for even level's layered medium model, it comprises the following steps:
1) in computing machine, read in information and big gun when the root-mean-square velocity information of seismic data zoning, the round trip that need to compensate vertically are walked examine distributed intelligence;
2) according to root-mean-square velocity information vertical with round trip walk time information, draw and comprise every interval velocity of one deck and even level's layered medium model of layer thickness information in seismic data zoning that needs compensate;
3) integrating step 1) in the big gun that obtains examine distributed intelligence, to step 2) even level's layered medium model of building, information when carrying out ray tracing to obtain ray position and to walk from shot point;
4) according to step 3) in the reflection angle in each reflection horizon that obtains of ray tracing result, determined the relative geometrical attenuation offset in each reflection horizon by geometrical attenuation compensation formula;
5) offset that two rays the most contiguous with geophone station calculate is carried out the offset that linear interpolation obtains geophone station place, according to the form of geological data, the relative geometrical attenuation offset obtained is interpolated on the sampled point of sympathetic earthquakes data volume, and exports data in order to the use of follow-up data treatment technology.
2. a kind of geometrical attenuation compensation method based on even level's layered medium model as claimed in claim 1, is characterized in that: described step 2) in, interval velocity is determined by Dix formula (1):
Wherein, V
ibe the interval velocity of i-th layer, U
ibe the root-mean-square velocity of i-th layer, U
i+1be the root-mean-square velocity of the i-th+1 layer, T
ifor the round trip of corresponding i-th layer vertically walk time, T
i+1for the round trip of corresponding the i-th+1 layer vertically walk time.
3. a kind of geometrical attenuation compensation method based on even level's layered medium model as claimed in claim 1, is characterized in that: described step 2) in, layer thickness is determined by following formula (2):
d
i＝V
i×(T
i+1-T
i)/2 (2)
Wherein, d
ibe the layer thickness of i-th layer, V
ibe the interval velocity of i-th layer, T
ifor the round trip of corresponding i-th layer vertically walk time, T
i+1for the round trip of corresponding the i-th+1 layer vertically walk time.
4. a kind of geometrical attenuation compensation method based on even level's layered medium model as claimed in claim 2, is characterized in that: described step 2) in, layer thickness is determined by following formula (2):
d
i＝V
i×(T
i+1-T
i)/2 (2)
Wherein, d
iit is the layer thickness of i-th layer.
5. a kind of geometrical attenuation compensation method based on even level's layered medium model as claimed in claim 1 or 2 or 3 or 4, is characterized in that: described step 3) in, ray position is determined by Snell's law (3):
V
1/sinα
1＝V
2/sinβ
2(3)
Wherein, V
1for the interval velocity of ground floor, α
1for incident angle when ray incides the second layer from ground floor, V
2for the interval velocity of the second layer, β
2for ray incides the angle of transmission after the second layer from ground floor.
6. a kind of geometrical attenuation compensation method based on even level's layered medium model as claimed in claim 1 or 2 or 3 or 4, is characterized in that: described step 3) in, when walking, information is determined by following formula (4):
t
n＝d
n/(V
n·cosβ
n) (4)
Wherein, t
nfor information during the walking of n-th layer, d
nfor the layer thickness of n-th layer, V
nfor the interval velocity of n-th layer, β
nfor the angle of transmission of n-th layer.
7. a kind of geometrical attenuation compensation method based on even level's layered medium model as claimed in claim 5, is characterized in that: described step 3) in, when walking, information is determined by following formula (4):
t
n＝d
n/(V
n·cosβ
n) (4)
Wherein, t
nfor information during the walking of n-th layer, d
nfor the layer thickness of n-th layer, V
nfor the interval velocity of n-th layer, β
nfor the angle of transmission of n-th layer.
8. a kind of geometrical attenuation compensation method based on even level's layered medium model as described in claim 1 or 2 or 3 or 4 or 7, is characterized in that: described step 4) in, relative geometrical attenuation offset is determined by following formula (5):
Wherein,
for the relative geometrical attenuation offset of n-th layer, V
1for the interval velocity of ground floor, α
1for incident angle when ray incides the second layer from ground floor, x is geophone offset, d
ibe the layer thickness of i-th layer, α
iit is the incident angle of i-th layer of incident ray.
9. a kind of geometrical attenuation compensation method based on even level's layered medium model as claimed in claim 5, is characterized in that: described step 4) in, relative geometrical attenuation offset is determined by following formula (5):
Wherein,
for the relative geometrical attenuation offset of n-th layer, x is geophone offset, d
ibe the layer thickness of i-th layer, α
iit is the incident angle of i-th layer of incident ray.
10. a kind of geometrical attenuation compensation method based on even level's layered medium model as claimed in claim 6, is characterized in that: described step 4) in, relative geometrical attenuation offset is determined by following formula (5):
Wherein,
for the relative geometrical attenuation offset of n-th layer, V
1for the interval velocity of ground floor, α
1for incident angle when ray incides the second layer from ground floor, x is geophone offset, d
ibe the layer thickness of i-th layer, α
iit is the incident angle of i-th layer of incident ray.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510181316.XA CN104777514A (en) | 2015-04-16 | 2015-04-16 | Geometric spreading compensation method based on uniform horizontal layered medium model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510181316.XA CN104777514A (en) | 2015-04-16 | 2015-04-16 | Geometric spreading compensation method based on uniform horizontal layered medium model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104777514A true CN104777514A (en) | 2015-07-15 |
Family
ID=53619091
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510181316.XA Pending CN104777514A (en) | 2015-04-16 | 2015-04-16 | Geometric spreading compensation method based on uniform horizontal layered medium model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104777514A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106772577A (en) * | 2016-06-29 | 2017-05-31 | 中国石油大学(华东) | Source inversion method based on microseism data and SPSA optimized algorithms |
| CN107229067A (en) * | 2017-05-27 | 2017-10-03 | 中国科学院测量与地球物理研究所 | Stratified model anisotropy ray-tracing procedure and system |
| CN110031897A (en) * | 2019-03-21 | 2019-07-19 | 北京信息科技大学 | Amplitude energy compensates and bearing calibration and system multi-component earthquake data in the seawater |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1239304A1 (en) * | 2001-03-05 | 2002-09-11 | Compagnie Generale De Geophysique | Improvements to processes for tomographic inversion of picked events in migrated seismic data |
| CN101315427A (en) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | Method and system for processing seismic exploration data of complex area |
| CN101630018A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气股份有限公司 | Method for processing seismic exploration data in process of controlling full acoustic wave equation inversion |
| CN102109615A (en) * | 2010-11-19 | 2011-06-29 | 恒泰艾普石油天然气技术服务股份有限公司 | Amplitude compensation method and device for seismic data processing under condition of true ground surface |
| CN102478667A (en) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | Inversion method of frequency dispersion amplitude versus offset (AVO) |
-
2015
- 2015-04-16 CN CN201510181316.XA patent/CN104777514A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1239304A1 (en) * | 2001-03-05 | 2002-09-11 | Compagnie Generale De Geophysique | Improvements to processes for tomographic inversion of picked events in migrated seismic data |
| CN101315427A (en) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | Method and system for processing seismic exploration data of complex area |
| CN101630018A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气股份有限公司 | Method for processing seismic exploration data in process of controlling full acoustic wave equation inversion |
| CN102109615A (en) * | 2010-11-19 | 2011-06-29 | 恒泰艾普石油天然气技术服务股份有限公司 | Amplitude compensation method and device for seismic data processing under condition of true ground surface |
| CN102478667A (en) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | Inversion method of frequency dispersion amplitude versus offset (AVO) |
Non-Patent Citations (4)
| Title |
|---|
| ALIA.NOWROOZI: "计算层速度的Dix公式的一般形式", 《国外地质勘探技术》 * |
| 刘明忱 等: "基于模型的2.5维几何扩散研究", 《中国地球物理2013》 * |
| 叶云飞: "Kirchhoff型真振幅反偏移脉冲响应分析", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 石殿祥 等: "角道集上的振幅补偿", 《石油地球物理勘探》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106772577A (en) * | 2016-06-29 | 2017-05-31 | 中国石油大学(华东) | Source inversion method based on microseism data and SPSA optimized algorithms |
| CN107229067A (en) * | 2017-05-27 | 2017-10-03 | 中国科学院测量与地球物理研究所 | Stratified model anisotropy ray-tracing procedure and system |
| CN107229067B (en) * | 2017-05-27 | 2019-05-31 | 中国科学院测量与地球物理研究所 | Stratified model anisotropy ray-tracing procedure and system |
| CN110031897A (en) * | 2019-03-21 | 2019-07-19 | 北京信息科技大学 | Amplitude energy compensates and bearing calibration and system multi-component earthquake data in the seawater |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102967882B (en) | The modeling method of the interval velocity model on stratum | |
| WO2015196779A1 (en) | Method for obtaining velocity of high-accuracy seismic waves using stacking velocity | |
| CN104570102B (en) | Method for combining near-surface velocity model with middle-deep stratum velocity model | |
| CN102841379B (en) | Method for analyzing pre-stack time migration and speed based on common scatter point channel set | |
| CN107817516B (en) | Near-surface modeling method and system based on first-motion wave information | |
| CN103576197B (en) | A kind of converted wave angle gathers abstracting method | |
| CN102053260B (en) | Method for acquiring azimuth velocity of primary wave and method for processing earthquake data | |
| CN108845350B (en) | Method and device for inverting two-dimensional velocity model | |
| CN102901984A (en) | Method for constructing true earth surface dip angle trace gathers of seismic data | |
| CN104777514A (en) | Geometric spreading compensation method based on uniform horizontal layered medium model | |
| CN104122582A (en) | Method for calculating high-accuracy seismic wave velocity by utilizing stacking velocity | |
| CN109143366A (en) | Near surface first arrival tomographic statics method and computer readable storage medium | |
| Gong et al. | Combined migration velocity model-building and its application in tunnel seismic prediction | |
| Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
| CN102353992B (en) | Irregular topography dynamic correction method based on central point regularization | |
| CN110646840B (en) | Angle gather extraction method and system | |
| CN112684498A (en) | Reservoir fracture prediction method and system based on wide-azimuth seismic data | |
| CN102478665B (en) | Method for determining incidence angle and amplitude of seismic wave | |
| CN104570091B (en) | A kind of method and apparatus for obtaining first arrival wave ray | |
| CN109490964B (en) | Improved high-precision AVO elastic parameter fast inversion method | |
| Baziw et al. | Deriving interval velocities from downhole seismic data | |
| CN109581521B (en) | Local chromatographic method and system for TTI anisotropy | |
| US10401515B2 (en) | Estimation of water properties from seismic data | |
| CN102033244B (en) | High-precision stacking and imaging method suitable for shallow curved earth surface | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150715 |