CN109856676A - A method of realizing earthquake common reflection surface stack parameter optimization - Google Patents
A method of realizing earthquake common reflection surface stack parameter optimization Download PDFInfo
- Publication number
- CN109856676A CN109856676A CN201811451480.8A CN201811451480A CN109856676A CN 109856676 A CN109856676 A CN 109856676A CN 201811451480 A CN201811451480 A CN 201811451480A CN 109856676 A CN109856676 A CN 109856676A
- Authority
- CN
- China
- Prior art keywords
- indicate
- wave
- initial
- curvature radius
- radius
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
The present invention provides a kind of method for realizing earthquake common reflection surface stack parameter optimization, comprising the following steps: step S1 obtains seismic data, based on radius of curvature and inclination angle before seismic data acquisition initial reflection point wave-front curvature radius, initial reflection surface wave;Step S2 is based on objective function, obtains the optimal value at radius of curvature and inclination angle before the reflection point wave-front curvature radius, initial reflection surface wave by radius of curvature and inclination angle before the initial reflection point wave-front curvature radius, initial reflection surface wave.The present invention provides a kind of method for realizing earthquake common reflection surface stack parameter optimization, first finds out three initial parameters by the first step, search initial parameter be it is single search for one by one, simple linear search can be used;Second step is to search for optimum combination solution near initial parameter using global optimization approach, to realize that CMAES algorithm solves the problems, such as that CRS imaging parameters quickly scan.
Description
Technical field
The present invention relates to geological exploration technical fields, more particularly, to a kind of realization earthquake common reflection surface stack parameter
The method of optimization.
Background technique
Common reflective surface (common reflection surface, abbreviation CRS) superposition is a kind of special zero-offset
Imaging method has the characteristics that drive independent of macro-velocity model and complete data.Since CRS theory introduces reflecting surface
Concept, and be limited within the scope of the first Fresnel-zone, so as to be cutd open by improving stacking fold to enhance superposition
While the signal-to-noise ratio in face, and certain resolution ratio is maintained, so CRS stacking method is more suitable for low signal-to-noise ratio, low covering
The seismic data of number.Specifically, common reflection surface (CRS) stacking method utilizes Fresnel principle, and spread reflection face element is (right
For two dimension, expand line element) increase degree of covering, it is theoretical using ray holography, different CRP gathers is corrected
It into same trace gather, is then overlapped, is enhanced the energy for reflecting signal.
Under two-dimensional case, three parameters need to be determined, it is slightly complicated for opposite confirmation CMP stack velocity, in order to find out
The optimum combination of three parameters, we first look for three initial parameter values, then find out optimum combination in initial value attachment.Initial value scanning
Using linear scan, to determine the probable ranges of three parameters, need to come using global optimization approach later on this basis again same
When optimized three parameters.
Under two-dimensional condition, ray theory (Schleicher etc., 1993), classical CRS trip when according to paraxial ray travelling
Equation can be expressed as (Jager etc., 2001) hyperbolic type expression formula when row,
(1);
In formula,When travelling for zero shot-geophone distance,For the speed at zero-offset,For the corresponding coordinate of zero-offset,For
Midpoint coordinates between big gun inspection,For half offset distance,For zero-offset ray emergence angle,Indicate the curvature half in IP wave
Diameter,It indicatesThe radius of curvature of point N wave.
Traditional CMP superposition needs a speed parameter, and CRS superposition needs three parameters, at this time parameter search space
Three-dimensional is risen to from one-dimensional, search calculation amount exponentially increases, and calculation amount is huge.The one-dimensional space is searched for using linear search
Can, but three-dimensional space needs more advanced quick Optimizing Search algorithm.At present using it is more be simulated annealing and base
In the algorithm of gradient, relative efficiency is lower, is easily trapped into local minimum, limits the popularization and application of CRS technology.
Summary of the invention
The present invention provides a kind of method for realizing earthquake common reflection surface stack parameter optimization, and three are first found out by the first step
Initial parameter, search initial parameter be it is single search for one by one, simple linear search can be used;Second step is using global optimization
Algorithm searches for optimum combination solution near initial parameter, to realize that CMAES algorithm solves the problems, such as that CRS imaging parameters quickly scan.
According to an aspect of the present invention, a kind of method for realizing earthquake common reflection surface stack parameter optimization is provided, including
Following steps:
Step S1 obtains seismic data, before seismic data acquisition initial reflection point wave-front curvature radius, initial reflection surface wave
Radius of curvature and initial tilt;
Step S2 is based on mesh by radius of curvature and inclination angle before the initial reflection point wave-front curvature radius, initial reflection surface wave
Scalar functions obtain the wave-front curvature radius, reflecting surface wave-front curvature radius and the optimal value at inclination angle.
On the basis of above scheme preferably, the objective function are as follows:
;
And;
;
;
Wherein,
Indicate reflection point wave-front curvature radius to be optimized;
Indicate inclination angle to be optimized;
Indicate wave-front curvature radius to be optimizedInverse, -1≤≤+1;
Indicate reflecting surface wave-front curvature radius to be optimized;
Indicate initial reflection point wave-front curvature radius;
Indicate initial tilt;
The inverse of radius of curvature before expression initial reflection surface wave, -1≤≤+1;
Indicate the compatibility function that seismic signal changes over time;
A indicates earthquake channel amplitude, is the function of time location t;
Indicate the consistency composite function that seismic signal changes over time;
M indicates the road number of current trace gather;
K indicates the seismic signal number of samples in actual time window;
Indicate the L of sampling point spatial offset initial parameter model2Norm;
Indicate the coefficient of regularization, 0 <<+1.
On the basis of above scheme preferably, the calculation expression of the position time t are as follows:
;
Wherein: h indicates offset distance;
Indicate reflection point wave-front curvature radius to be optimized;
Indicate inclination angle to be optimized;
Indicate the vertical travelling two-way time;
Indicate earth's surface speed.
A kind of method for realizing earthquake common reflection surface stack parameter optimization of the invention, by the way that objective function is designed to one
Cause property function and weighted model difference and structure, make it in search parameter, guarantee obtain relatively good condition for consistence
Under, and without departing from given original model parameter (model difference is small), and regularisation parameter, it can be used for adjusting and deviate model journey
The size of degree.
Reflection point wave-front curvature radius, reflecting surface wave-front curvature radius when the present invention is by acquisition objective function minimum value
And the optimal value combination at inclination angle, formula (1) can calculate CRS hourage when being travelled using CRS, so that it may to CMP trace gather
Or channel set moved compared with be superimposed, generate last CRS superposition achievement section, done for Geophysicist and geological personnel
Further processing and explanation.
Detailed description of the invention
Fig. 1 be certain CDP point position of the invention three parameter optimizations before with the curve graph after optimization;
Fig. 2 is that suburb positive time depth and offset distance variation diagram at any time are moved in certain position CDP of the invention;
Fig. 3 is single CDP stacked profile map of the invention;
Fig. 4 is CRS stacked profile map of the invention;
Fig. 5 is a kind of method flow block diagram for realizing earthquake common reflection surface stack parameter optimization of the invention.
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.Implement below
Example is not intended to limit the scope of the invention for illustrating the present invention.
It please refers to shown in Fig. 5, the present invention provides a kind of methods for realizing earthquake common reflection surface stack parameter optimization, including
Following steps:
Step S1 obtains seismic data, is calculated based on seismic data and obtains initial reflection point wave-front curvature radius, initial reflection face
Wave-front curvature radius and initial tilt;
Step S2 is based on mesh by radius of curvature and initial tilt before initial reflection point wave-front curvature radius, initial reflection surface wave
Scalar functions obtain the reflection point wave-front curvature radius, reflecting surface wave-front curvature radius and the optimal value at inclination angle.
For the technical solution of the present invention of further instruction, introduced detailed in step S1 of the invention below,
For obtaining reflection point wave-front curvature radius, reflecting surface wave-front curvature radius and the initial parameter at inclination angle.
In common seismic Data processing, it would be desirable to artificial pickup velocity, and the stack velocity needed in CRS processingIt is exactly conventional stacking speed, but surveying stack velocity can realize according to CMP travel-time equation automatically scanning.Using artificial
Reference of the speed of pickup as automatically scanning, on the basis of artificial pickup velocity, scanning obtains compatibility function maximum
Corresponding speed when value, the stability for the scanning speed that can guarantee in this way, while avoiding scanning to the low of multiple wave again
Fast range.
Under two-dimensional condition, ray theory (Schleicher etc., 1993), classical CRS trip when according to paraxial ray travelling
Equation can be expressed as (Jager etc., 2001) hyperbolic type expression formula when row,
(1);
Wherein, in formula,When travelling for zero shot-geophone distance,For the speed at zero-offset,For the corresponding seat of zero-offset
Mark,For big gun inspection between midpoint coordinates,For half offset distance,For zero-offset ray emergence angle,It indicates in point NIP
The radius of curvature of wave,It indicatesThe radius of curvature of point N wave.
Due to stack velocity in CRS processingWithIt is directly related, for formula (1), it is assumed that, then the public affairs
Formula is variable are as follows:
(2)
Stack velocity formula is compared,
(3)
It enables,
(4)
It is obtained after transformation,
(5)
It can be concluded that, obtained from upper formulaThe inclination angle and, dependent variable is all constant, can be calculated.It obtains in next step
Scan inclination angle。
After the scanning of previous step stack velocity, superposition is done using the speed of scanning and obtains stacked section.Notice that inclination angle has
It is positive and negative, such as -40 ° to+40 ° of scan position are defined before scanning.After obtaining inclination angle, in conjunction with the first step, we scan obtained superposition speed
Degree, formula (5) calculating can be used in we。
Equally, the radius of curvature for scanning N wave still uses result obtained in the previous step as input, i.e. stacked section.Make
With containing after simplificationEquation linear scan N wave radius of curvature.Pay attention toRange be it is infinite to just infinite from bearing, can be with
Its inverse, which is scanned, by linear scan obtains its initial value.
Specific implementation, is parameter automatically scanning speed, it need to be noted that stack velocity is extremely important in CRS processing,
But the rate pattern automatically scanning without any restriction on the parameters may have many problems, before scanning speed, provide with reference to speed
Degree carries out linear scan according to come the range that sets scanning centered on reference velocity, take consistency preferably when corresponding speed
Degree, then will obtain that carry out is smooth, later using this speed do common midpoint gather it is dynamic compared with and superposition, generate earthquake overlap and cut open
Face.After generating section, angle scanning is carried out to section, obtains the inclination value that each sampling point goes out.The scanning at inclination angle uses greatly first
Angle interval scan, such as 2 ° or so use closely-spaced carry out rescan, such as 0.1 ° or so after the general angular range of determination.
After preliminary sweep obtains speed and inclination angle, NIP radius can be calculated according to formula (5);Assuming that the h of formula (1) is 0, i.e.,
It can scan to obtain N wave radius, take its inverse, extremely N curve rate.
Obtaining three initial values, inclination angle, the radius of curvature of N wave, the radius of curvature of NIP wave, next may be used
To use the algorithm of global optimization, best parameter group is found near initial value.
Wherein, CMAES is a kind of algorithm (Hansen, 2011) of random multi-Dimensional parameters search unrelated with gradient.It should
Algorithm is not based on gradient search, does not also have derivative calculations process, is suitble to various complex search problems, and such as more Local Minimums are asked
Topic, non-linear smooth perturbation problem, the inseparable problem of parameter, discontinuously and containing noise problem etc..CMAES algorithm search is one
The process of a iteration updates stochastic variable using the continuous iteration of stochastical sampling method, can be summarized as three steps according to its feature:
(1) spatial sampling is carried out to data according to sampling point mean value and covariance matrix, (2) estimate new sampling point mean value, and (3) are estimated new
Sampling point covariance.Successively continuous loop iteration updates this three step, can search parameter space be mostly efficiently multidimensional solution, and
And fast convergence rate.
This distinctive random parameter update mode can be expressed as with equation,
(6)
(7)
(8)
~: indicate stochastic variable distribution;
: indicate multiple random variables normal distribution, mean value 0, covariance matrix;
: indicate k-th of offspring from g+1 generation;
: indicate the mean value in g generation;
: indicate g for population variance or step-length;
: indicate g for covariance matrix;
: it indicates total sample, will be generally greater than or equal to 2;
: indicate weighting coefficient space;
: representation space sample;
: indicate g+1 for Space Experiments covariance.
The search of CRS superposition parameter is substantially exactly one containing noise, and non-linear, discontinuously, inseparable parameter is excellent
Change problem, thus CMAES is a kind of extraordinary optimization method for being suitable for CRS parameter search.
Searching for the first step first using CMAES is to define minimum objective function, objective function of the invention are as follows:
(9);
And(10);
(11);
(12);
Wherein,
Indicate reflection point wave-front curvature radius to be optimized;
Indicate inclination angle to be optimized;
Indicate wave-front curvature radius to be optimizedInverse, -1≤≤+1;
Indicate reflecting surface wave-front curvature radius to be optimized;
Indicate initial reflection point wave-front curvature radius;
Indicate initial tilt;
The inverse of radius of curvature before expression initial reflection surface wave, -1≤≤+1;
Indicate the compatibility function that seismic signal changes over time;
A indicates earthquake channel amplitude, is the function of time location t;
Indicate the consistency composite function that seismic signal changes over time;
M indicates the road number of current trace gather;
K indicates the seismic signal number of samples in actual time window;
Indicate the L of sampling point spatial offset initial parameter model2Norm;
Indicate the coefficient of regularization, 0 <<+1.
Wherein, the calculation expression of position time t of the invention are as follows:
;
Wherein: h indicates offset distance;
Indicate reflection point wave-front curvature radius to be optimized;
Indicate inclination angle to be optimized;
Indicate the vertical travelling two-way time;
Indicate earth's surface speed.
It is illustrated below with reference to embodiment and how objective function is based on come the initial ginseng of scanning optimization three using CMAES
Number.
After three initial parameters are calculated, using CMAES based on objective function (9) come three ginsengs after scanning optimization
Number, is illustrated here with example.Fig. 1 be certain CDP point position three parameter optimizations before (blue) optimization after (red) curve graph.Phase
Comparatively, red curve bounce is more, it is a kind of random algorithm that this, which is due to CMAES, in order to find best parameter group,
The optimizing near initial parameter of the parameter of search reaches smaller target function value.
After searching optimal three parameter combinations, followed by the hourage for using formula (1) that curved surface is calculated,
Fig. 2 be shown certain position CDP move the positive time graph in suburb, x-axis is offset distance, the time parameter in Fig. 2 be after smooth,
Due to being random solution procedure, three parameters of calculating contain random element, thus the time calculated need to do it is smooth.
Obtained it is smooth after dynamic correction time, can CDP trace gather data be carried out with dynamic correction and superposition, what Fig. 3 was shown
It is single CDP stacked profile map, Fig. 4 is CRS stacked profile map, belongs to newest fruits of the invention, it can be seen that the signal of Fig. 4
Obtained enhancing, noise are suppressed, and Fig. 4 imaging effect is greatly improved.
Finally, the present processes are only preferable embodiment, it is not intended to limit the scope of the present invention.It is all
Within the spirit and principles in the present invention, any modification, equivalent replacement, improvement and so on should be included in protection of the invention
Within the scope of.
Claims (3)
1. a kind of method for realizing earthquake common reflection surface stack parameter optimization, which comprises the following steps:
Step S1 obtains seismic data, is calculated based on seismic data and obtains initial reflection point wave-front curvature radius, initial reflection face
Wave-front curvature radius and initial tilt;
Step S2 passes through radius of curvature and initial tilt before the initial reflection point wave-front curvature radius, initial reflection surface wave, base
In objective function, the wave-front curvature radius, reflecting surface wave-front curvature radius and the optimal value at inclination angle are obtained.
2. a kind of method for realizing earthquake common reflection surface stack parameter optimization as described in claim 1, which is characterized in that described
Objective function are as follows:
;
And;
;
;
Wherein,
Indicate reflection point wave-front curvature radius to be optimized;
Indicate inclination angle to be optimized;
Indicate wave-front curvature radius to be optimizedInverse, -1≤≤+1;
Indicate reflecting surface wave-front curvature radius to be optimized;
Indicate initial reflection point wave-front curvature radius;
Indicate initial tilt;
The inverse of radius of curvature before expression initial reflection surface wave, -1≤≤+1;
Indicate the compatibility function that seismic signal changes over time;
A indicates earthquake channel amplitude, is the function of time location t;
Indicate the consistency composite function that seismic signal changes over time;
M indicates the road number of current trace gather;
K indicates the seismic signal number of samples in actual time window;
Indicate the L of sampling point spatial offset initial parameter model2Norm;
Indicate the coefficient of regularization, 0 <<+1.
3. a kind of method for realizing earthquake common reflection surface stack parameter optimization as claimed in claim 2, which is characterized in that described
The calculation expression of position time t are as follows:
;
Wherein: h indicates offset distance;
Indicate reflection point wave-front curvature radius to be optimized;
Indicate inclination angle to be optimized;
Indicate the vertical travelling two-way time;
Indicate earth's surface speed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811451480.8A CN109856676B (en) | 2018-11-30 | 2018-11-30 | Method for realizing optimization of stacking parameters of seismic common reflection surface |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811451480.8A CN109856676B (en) | 2018-11-30 | 2018-11-30 | Method for realizing optimization of stacking parameters of seismic common reflection surface |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109856676A true CN109856676A (en) | 2019-06-07 |
CN109856676B CN109856676B (en) | 2021-03-12 |
Family
ID=66890441
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811451480.8A Active CN109856676B (en) | 2018-11-30 | 2018-11-30 | Method for realizing optimization of stacking parameters of seismic common reflection surface |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109856676B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110954955A (en) * | 2019-12-23 | 2020-04-03 | 西南石油大学 | Seismic stochastic inversion method based on evolutionary optimization algorithm |
CN112327364A (en) * | 2020-11-02 | 2021-02-05 | 中国石油大学(华东) | Improved one-step method three-dimensional CRS (Cross-reference Signal) superposition method |
CN113960669A (en) * | 2021-10-21 | 2022-01-21 | 中国石油化工股份有限公司 | Reflection information enhancement method and device based on common imaging point gather combination calculation |
CN113960668A (en) * | 2021-10-21 | 2022-01-21 | 中国石油化工股份有限公司 | Method and device for enhancing reflection information based on pre-stack time migration |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4878205A (en) * | 1987-07-23 | 1989-10-31 | Boris Gelchinsky | Method of analyzing the structure of a medium, particularly useful for seismic prospecting |
CN105301648A (en) * | 2014-07-31 | 2016-02-03 | 中国石油化工股份有限公司 | Method of acquiring common reflection surface stacking parameters |
CN106483559A (en) * | 2015-08-24 | 2017-03-08 | 中国石油化工股份有限公司 | A kind of construction method of subsurface velocity model |
KR20170040569A (en) * | 2015-10-05 | 2017-04-13 | 주식회사 에프에스 | Eismic Data Processing Method Using Common Reflection Surface and Common Diffraction Surface |
CN107179543A (en) * | 2016-03-11 | 2017-09-19 | 中国石油化工股份有限公司 | The method and apparatus that regularization is carried out to pre stack data |
CN109557585A (en) * | 2019-01-25 | 2019-04-02 | 中国石油大学(华东) | A kind of CRS dip decomposition method based on differential evolution algorithm |
-
2018
- 2018-11-30 CN CN201811451480.8A patent/CN109856676B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4878205A (en) * | 1987-07-23 | 1989-10-31 | Boris Gelchinsky | Method of analyzing the structure of a medium, particularly useful for seismic prospecting |
CN105301648A (en) * | 2014-07-31 | 2016-02-03 | 中国石油化工股份有限公司 | Method of acquiring common reflection surface stacking parameters |
CN106483559A (en) * | 2015-08-24 | 2017-03-08 | 中国石油化工股份有限公司 | A kind of construction method of subsurface velocity model |
KR20170040569A (en) * | 2015-10-05 | 2017-04-13 | 주식회사 에프에스 | Eismic Data Processing Method Using Common Reflection Surface and Common Diffraction Surface |
CN107179543A (en) * | 2016-03-11 | 2017-09-19 | 中国石油化工股份有限公司 | The method and apparatus that regularization is carried out to pre stack data |
CN109557585A (en) * | 2019-01-25 | 2019-04-02 | 中国石油大学(华东) | A kind of CRS dip decomposition method based on differential evolution algorithm |
Non-Patent Citations (2)
Title |
---|
杨锴等: "共反射面元叠加的应用实践", 《地球物理学报》 * |
黄亚飞 等: "求解全局优化问题的正交协方差矩阵自适应进化策略算法", 《计算机应用》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110954955A (en) * | 2019-12-23 | 2020-04-03 | 西南石油大学 | Seismic stochastic inversion method based on evolutionary optimization algorithm |
CN112327364A (en) * | 2020-11-02 | 2021-02-05 | 中国石油大学(华东) | Improved one-step method three-dimensional CRS (Cross-reference Signal) superposition method |
CN112327364B (en) * | 2020-11-02 | 2022-04-29 | 中国石油大学(华东) | Improved one-step method three-dimensional CRS (Cross-reference Signal) superposition method |
CN113960669A (en) * | 2021-10-21 | 2022-01-21 | 中国石油化工股份有限公司 | Reflection information enhancement method and device based on common imaging point gather combination calculation |
CN113960668A (en) * | 2021-10-21 | 2022-01-21 | 中国石油化工股份有限公司 | Method and device for enhancing reflection information based on pre-stack time migration |
CN113960668B (en) * | 2021-10-21 | 2024-04-16 | 中国石油化工股份有限公司 | Method and device for enhancing reflection information based on prestack time migration |
Also Published As
Publication number | Publication date |
---|---|
CN109856676B (en) | 2021-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109856676A (en) | A method of realizing earthquake common reflection surface stack parameter optimization | |
US10330807B2 (en) | Device and method for displaying full azimuth angle domain image data | |
Gizon et al. | Local helioseismology | |
Wu et al. | Directional illumination analysis using beamlet decomposition and propagation | |
Medvedev et al. | Relation of internal gravity wave anisotropy with neutral wind characteristics in the upper atmosphere | |
CN106597540B (en) | Gaussian beam offset imaging method and device | |
CN106094029A (en) | The method utilizing offset distance vector sheet geological data Predicating Reservoir Fractures | |
CN112883564B (en) | Water body temperature prediction method and prediction system based on random forest | |
CN105607122B (en) | A kind of earthquake texture blending and Enhancement Method based on full variation geological data decomposition model | |
Abbad et al. | Automatic nonhyperbolic velocity analysis | |
Xiang‐Bo et al. | Anisotropic Radon transform and its application to demultiple | |
Liu et al. | Navigability analysis of local gravity map with projection pursuit-based selection method by using gravitation field algorithm | |
CN106257309A (en) | Post-stack seismic data body processing method and processing device | |
Coimbra et al. | Offset-continuation stacking: Theory and proof of concept | |
CN109581488A (en) | One kind being based on particle group optimizing CRS channel set seismic imaging method | |
CN109188516A (en) | The micro-seismic event localization method of the domain Radon energy scan superposition | |
CN109387872A (en) | Surface-related multiple prediction technique | |
CN112904407B (en) | Micro-motion exploration method under complex terrain and interference conditions | |
CN108983290A (en) | It is a kind of three-dimensional Method in Transverse Isotropic Medium in travel when determine method and system | |
CN106569269B (en) | A kind of preferred trace weighting normal-moveout spectrum calculation method | |
CN109655888A (en) | The quantitative choosing method and system of smooth floating datum in seismic data process | |
Tang et al. | A stable and effective strategy for angle gather extraction from reverse time migration using traveltime gradient | |
Wu et al. | Stacking common-image gathers via the flatten-illumination-similarity method | |
Leung et al. | Liouville partial-differential-equation methods for computing 2D complex multivalued eikonals in attenuating media | |
Zang et al. | Calibration method of petroleum underground prospects based on high precision gravity and magnetic exploration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |