CN106353798A - Multi-component joint Gaussian beam pre-stack reverse-time migration imaging method - Google Patents

Abstract

The invention discloses a multi-component joint Gaussian beam pre-stack reverse-time migration imaging method. The multi-component joint Gaussian beam pre-stack reverse-time migration imaging method includes inputting initial velocity fields of longitudinal waves and transverse waves; reading seismic shot records of the longitudinal waves and converted waves and determining parameters of frequency bandwidths, beam widths, beam center intervals and the like; constructing underground elastic vector wave fields; carrying out imaging by the aid of cross-correlation imaging conditions; migrating and superposing all the shot records to obtain ultimate elastic wave imaging results. Imaging achievement can be effectively migrated by the aid of the multi-component joint Gaussian beam pre-stack reverse-time migration imaging method. Compared with the prior art, the multi-component joint Gaussian beam pre-stack reverse-time migration imaging method has the advantages that the vector wave fields with abundant information can be constructed, and accordingly the problem of multiple imaging solutions due to an insufficient quantity of information of single longitudinal waves or transverse waves can be solved; Gaussian beam migration is high in efficiency and is practical, and accordingly the multi-component joint Gaussian beam pre-stack reverse-time migration imaging method is a depth-domain migration imaging method with computational efficiency and imaging precision and is extremely high in practicality for actual data processing at present; processes and the parameters are easy to set by the aid of the multi-component joint Gaussian beam pre-stack reverse-time migration imaging method, and the multi-component joint Gaussian beam pre-stack reverse-time migration imaging method is high in operation speed and suitable to be applied to processing three-dimensional seismic data.

Description

Multi -components joint Gaussian beam pre-Stack Reverse imaging method

Technical field

The invention belongs to oil-gas exploration seismic data migration imaging processing technology field, specifically a kind of multi -components connection Close Gaussian beam pre-Stack Reverse imaging method.

Prior art

Gaussian beam skew is offset method based on ray theory, at present for mainly have traditional ray tracing Conventional ray-tracing procedure under method and high-frequency approximation.Traditional ray-tracing procedure biases toward to ray road Footpath and description when walking, by it is flexible and efficient, do not have inclination angle to limit and easy expansive approach is to relief surface Advantage, it has also become most widely used aborning kirchhoff integrates the important composition of pre-stack depth migration Part.Conventional ray tracing under high-frequency approximation thinks that central ray represents the main energy of seismic wave, only Seimic wave propagation process is described using central ray, such approximate processing can only reflect the motion of seismic wave Learn feature.For complex dielectricss, conventional ray tracing there is likely to be caustic and multipath in numerical computations Problem, therefore application effect be not very good.

So far, the also main compressional wave data processing method using routine of multi component signal skew, this First multi component signal is carried out in method with P-wave And S separate, then with existing compressional wave skew software respectively Process longitudinal wave field and shear wave field. the major defect of this method is: 1) it multi component signal is regarded as several The simple superposition of scalar wave, have ignored the many features of vector wave, thus affect the process of multi component signal with Explanation Accuracy；2) if seismic wave field is not properly separated as compressional wave and shear wave (this phenomenon is widely present), And assume in the migration process of data handled for scalar wave field, this will necessarily make migration result produce Many artefacts, affect processing accuracy.

During traditional multi-component earthquake data is processed, general way is to vertical by wave field separation Component compressional wave routinely is treated and is carried out conventional treatment to it, and horizontal component is then treated by ps converted shear wave It is carried out with the conventional treatment of similar compressional wave；And tend not to the energy of different wave modes during wave field separation Amount is kept completely separate, and non-remaining type wave energy can lead to the much noise in imaging results to disturb, serious shadow Ring imaging effect；Particularly traditional wave equation class offset method due to computational efficiency relatively low under, in mesh The problems such as there is no practicality before the big data dignity of front three-dimensional high-density acquisition.

With respect to longitudinal wave exploration, multi-wave multicomponent exploration can obtain more underground medium elastic information, many The reverse-time migration research of ripple data also makes certain gains, but is primarily directed to simple compressional wave and converted wave Migration technology, the joint reverse-time migration the method for ripple in length and breadth being currently based on Gaussian beam has no corresponding Research Literature.

Content of the invention

The purpose of the present invention is the problem existing for prior art, and a kind of imaging effect of proposition is good, practicality Strong multi -components joint Gaussian beam pre-Stack Reverse imaging method.

The technical scheme of the multi -components joint Gaussian beam pre-Stack Reverse imaging method of the present invention, comprising:

(1) input initial p-and s-wave velocity field；

(2) read in compressional wave and converted wave earthquake big gun record, and determine bandwidth, beam width, beam center interval Parameter；

(3) structure of underground elasticity vector wave field；

(4) cross-correlation image-forming condition is utilized to be imaged；

(5) all big gun records are entered line displacement and are superimposed, obtain final elastic wave imaging result.

Further:

Described step (3) is the forward and reverse continuation entering traveling-wave field using elastodynamics Gaussian beam equation；

Described step (4) is using zero-lag cross-correlation image-forming condition, the wave field of continuation to be imaged；

Described step (5) is to carry out migration imaging by big gun to all of big gun record, and all big gun records are inclined Move into as being overlapped, just obtain imaging results.

Further:

The structure of step (3) underground elasticity vector wave field includes:

Hypothesis focus is x_s, receiving point x_rThe two component elastic wave earthquake records receiving are u_i(x_r；ω), then instead Can integrate by kirchhoff-helmholtz to represent to the elastic wave field of continuation:

$u_m(x;x_r;\mathrm{\ω})={\&integral;}_s{\mathrm{dx}}_r\[t_i(x_r;\mathrm{\ω})g_{im}^{*}(x;x_r;\mathrm{\ω})-u_i(x_r;\mathrm{\ω}){\mathrm{\σ}}_{im}^{*}(x;x_r;\mathrm{\ω})\]---\left(1\right)$

Wherein, * represents complex conjugate；g_lm(x；x_r；It is ω) displacement Green tensor, it represents x_rPlace's l direction unit The component in the m direction of displacement at x caused by muscle power；t_i(x_r) it is x_rPlace's stress；σ_im(x；x_r) it is stress lattice Woods tensor；

Step (4) utilizes cross-correlation image-forming condition to be imaged:

According to claerbout image-forming principle, by asking for connecing of source wavefield and the backward extension of different wave modes Receive cross-correlation between wave field being calculated as picture value, for this reason, first by focus displacement wave fieldLogical Cross elastodynamics Gaussian beam to represent,

$u_m^p(x;x_s;\mathrm{\ω})\≈\frac{i}{4{\mathrm{\πv}}_p^2\left(x_s\right)}\sqrt{\frac{{\mathrm{\ω}}_r{w}_0^2}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right)}{p_2^p\left(x_s\right)}{\hat{u}}_m^p(x;x_s;\mathrm{\ω})---\left(2\right)$

Next the propagating characteristic according to p ripple and s ripple, is defined as follows imaging formula:

$\begin{array}{c}{i}^{pp}\left(x\right)=\&integral;u_2^{p*}(x;x_s;\mathrm{\ω})u_2^p(x;x_r;\mathrm{\ω})d\mathrm{\ω}\\ =\frac{\mathrm{\δ}l\sqrt{{\mathrm{\ω}}_r{w}_0^2}}{16{\mathrm{\π}}^2}\underset{l}{\σ}\&integral;d\mathrm{\ω}\frac{i\mathrm{\ω}}{v_p^2\left(x_s\right)}\sqrt{\frac{\mathrm{\ρ}\left(l\right)}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right){\mathrm{dp}}_1^p\left(l\right)}{p_2^p\left(x_s\right)p_2^p\left(l\right)}\\ \×{\hat{u}}_2^{p*}(x;x_s;\mathrm{\ω}){\hat{u}}_1^{p*}(x;l;\mathrm{\ω})\[w_1^p\left(l\right)d_1^p(l;p_1^p;\mathrm{\ω})+w_2^p\left(l\right)d_2^p(l;p_1^p;\mathrm{\ω})\]\end{array}---\left(3\right)$

$\begin{array}{c}{i}^{ps}\left(x\right)=\&integral;u_2^{p*}(x;x_s;\mathrm{\ω})u_1^s(x;x_r;\mathrm{\ω})d\mathrm{\ω}\\ =\frac{\mathrm{\δ}l\sqrt{{\mathrm{\ω}}_r{w}_0^2}}{16{\mathrm{\π}}^2}\underset{l}{\σ}\&integral;d\mathrm{\ω}\frac{i\mathrm{\ω}}{v_p^2\left(x_s\right)}\sqrt{\frac{\mathrm{\ρ}\left(l\right)}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right){\mathrm{dp}}_1^s\left(l\right)}{p_2^p\left(x_s\right)p_2^s\left(l\right)}\\ \×{\hat{u}}_2^{p*}(x;x_s;\mathrm{\ω}){\hat{u}}_1^{s*}(x;l;\mathrm{\ω})\[w_1^s\left(l\right)d_1^s(l;p_1^s;\mathrm{\ω})+w_2^s\left(l\right)d_2^s(l;p_1^s;\mathrm{\ω})\]\end{array}---\left(4\right)$

Wherein, i^ppX () becomes picture value, i for p-p single-shot^psX () becomes picture value for p-s single-shot,(v table The type of oscillography) it is the adding window local dip superposition to different wave mode multi component seismic records, l represents Gaussian beam The effective half width of the Gaussian beam relevant with frequency, focus is x_s, receiving point x_r, ρ is Media density, v_p(x_s) For the velocity of longitudinal wave of central ray at hypocentral location, w is weight coefficient, and p is kinetics ray-tracing parameter.

Invention effect

The present invention, divides by the use of multi-component seismic data as input not entering traveling-wave field to input data by directly From in the case of, positive and inverse time extrapolation is carried out based on equations for elastic waves, joint multi -components build underground jointly Elastic vector seismic wave field, and obtain pure compressional wave and the imaging results of pure shear wave using cross-correlation image-forming condition, Substantially increase imaging effect.And Gaussian beam reverse-time migration remains the efficient of conventional kirchhoff skew Property and motility it is also possible to easily adapt to complexity surface conditions.Can to many subwaves to being imaged, Imaging precision is better than conventional kirchhoff skew and close to wave equation migration, is that one kind has calculating effect concurrently Rate and the Depth Domain offset imaging method of imaging precision, have very strong reality at present in the process of real data The property used.

Using multi -components joint Gaussian beam pre-Stack Reverse imaging method can preferably migration imaging achievement, should Method has the advantage that other technologies do not possess, and its concrete advantage and feature show the following aspects:

Firstth, construct informative vector wave field.The method application wave field information that ripple enriches in length and breadth carrys out structure Build underground vector wave field, make full use of the respective advantage wave field characteristics of ripple in length and breadth, solve single compressional wave or Solving problem the imaging that shear wave brings because quantity of information is not enough more.

Secondth, Gaussian beam skew highly effective.Gaussian beam skew remains the efficient of conventional kirchhoff skew Property and motility it is also possible to easily adapt to complexity surface conditions.It is that one kind has computational efficiency and imaging concurrently The Depth Domain offset imaging method of precision, has very strong practicality at present in the process of real data.

3rd, easy realization simple to operate.The method flow process and parameter setting are simple, fast operation, and being suitable for should For D seismic data processing.

Brief description

Fig. 1 multi -components joint Gaussian beam pre-Stack Reverse imaging technique flow process；

The oil reservoir section of Fig. 2 realistic model；

Fig. 3 velocity of longitudinal wave model；

The migration imaging result of Fig. 4 compressional wave；

Fig. 5 changes velocity model；

The migration imaging result of Fig. 6 converted wave.

Specific embodiment

Overview embodiment:

Referring to the drawings 1, the techniqueflow of the multi -components joint Gaussian beam pre-Stack Reverse imaging method of the present invention It is:

(1) input initial p-and s-wave velocity field；

(2) read in compressional wave and converted wave earthquake big gun record, and determine bandwidth, beam width, beam center interval Etc. parameter；

(3) structure of underground elasticity vector wave field:

Hypothesis focus is x_s, receiving point x_rThe two component elastic wave earthquake records receiving are u_i(x_r；ω), then instead Can integrate by kirchhoff-helmholtz to represent to the elastic wave field of continuation:

$u_m(x;x_r;\mathrm{\ω})={\&integral;}_s{\mathrm{dx}}_r\[t_i(x_r;\mathrm{\ω})g_{im}^{*}(x;x_r;\mathrm{\ω})-u_i(x_r;\mathrm{\ω}){\σ}_{im}^{*}(x;x_r;\mathrm{\ω})\]---(4-1-4)$

Wherein, * represents complex conjugate；g_lm(x；x_r；It is ω) displacement Green tensor, it represents x_rPlace's l direction unit The component in the m direction of displacement at x caused by muscle power；t_i(x_r) it is x_rPlace's stress；σ_im(x；x_r) it is stress lattice Woods tensor；

(4) cross-correlation image-forming condition is utilized to be imaged:

According to claerbout image-forming principle, by asking for connecing of source wavefield and the backward extension of different wave modes Receive cross-correlation between wave field being calculated as picture value, for this reason, first by focus displacement wave fieldLogical Cross elastodynamics Gaussian beam to represent,

$u_m^p(x;x_s;\mathrm{\ω})\≈\frac{i}{4{\mathrm{\πv}}_p^2\left(x_s\right)}\sqrt{\frac{{\mathrm{\ω}}_2{w}_0^2}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right)}{p_2^p\left(x_s\right)}{\hat{u}}_m^p(x;x_s;\mathrm{\ω})---(4-1-22)$

Wherein, ω_rFor reference frequency, w₀For the original width of Gaussian beam, p is kinetics ray-tracing parameter.

Next the propagating characteristic according to p ripple and s ripple, is defined as follows imaging formula:

$\begin{array}{c}{i}^{pp}\left(x\right)=\&integral;u_2^{p*}(x;x_s;\mathrm{\ω})u_2^p(x;x_r;\mathrm{\ω})d\mathrm{\ω}\\ =\frac{\mathrm{\δ}l\sqrt{{\mathrm{\ω}}_r{w}_0^2}}{16{\mathrm{\π}}^2}\underset{l}{\σ}\&integral;d\mathrm{\ω}\frac{i\mathrm{\ω}}{v_p^2\left(x_s\right)}\sqrt{\frac{\mathrm{\ρ}\left(l\right)}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right){\mathrm{dp}}_1^p\left(l\right)}{p_2^p\left(x_s\right)p_2^p\left(l\right)}\\ \×{\hat{u}}_2^{p*}(x;x_s;\mathrm{\ω}){\hat{u}}_1^{p*}(x;l;\mathrm{\ω})\[w_1^p\left(l\right)d_1^p(l;p_1^p;\mathrm{\ω})+w_2^p\left(l\right)d_2^p(l;p_1^p;\mathrm{\ω})\]\end{array}---(4-1-23)$

$\begin{array}{c}{i}^{ps}\left(x\right)=\&integral;u_2^{p*}(x;x_s;\mathrm{\ω})u_1^s(x;x_r;\mathrm{\ω})d\mathrm{\ω}\\ =\frac{\mathrm{\δ}l\sqrt{{\mathrm{\ω}}_r{w}_0^2}}{16{\mathrm{\π}}^2}\underset{l}{\σ}\&integral;d\mathrm{\ω}\frac{i\mathrm{\ω}}{v_p^2\left(x_s\right)}\sqrt{\frac{\mathrm{\ρ}\left(l\right)}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right){\mathrm{dp}}_1^s\left(l\right)}{p_2^p\left(x_s\right)p_2^s\left(l\right)}\\ \×{\hat{u}}_2^{p*}(x;x_s;\mathrm{\ω}){\hat{u}}_1^{s*}(x;l;\mathrm{\ω})\[w_1^s\left(l\right)d_1^s(l;p_1^s;\mathrm{\ω})+w_2^s\left(l\right)d_2^s(l;p_1^s;\mathrm{\ω})\]\end{array}---(4-1-24)$

Wherein, i^ppX () becomes picture value, i for p-p single-shot^psX () becomes picture value for p-s single-shot,(v table The type of oscillography) it is the adding window local dip superposition to different wave mode multi component seismic records.L represents Gaussian beam The effective half width of the Gaussian beam relevant with frequency, focus is x_s, receiving point x_r, ρ is Media density, v_p(x_s) For the velocity of longitudinal wave of central ray at hypocentral location, w is weight coefficient, and p is kinetics ray-tracing parameter.

(5) all big gun records are entered line displacement and are superimposed, just obtain final elastic wave imaging result.

This invention major technique key point is following two: (1), how to build underground elasticity vector wave field；(2)、 The cross-correlation imaging of ripple in length and breadth.

Concrete application embodiment:

The detailed technology operation sketch of this invention is as shown in Figure 1.

This research is applied for target to xx oil field xx block real data model, with the method pair Multi component signal has carried out pre-stack depth migration imaging, and to verify the effect of this method, idiographic flow is shown in Fig. 1.

1) according to step 1, the stack velocity analysis being utilized respectively ripple in length and breadth obtain initial p-and s-wave velocity field.

2) according to step 2, compressional wave and converted waves data are read in.

3) according to step 3, input compressional wave and converted waves data are analyzed, the bandwidth of determination skew, Beam width, beam center interval and initial ray parameter space.

4) according to step 5, using elastodynamics Gaussian beam equation, the wave field in each moment is just carried out To and backward extension.

5) according to step 6, all whether continuation finishes, if do not had migration program automatic decision all moment wave field Continuation completes, and starts to proceed continuation calculating from step 5.

6) according to step 7, using zero-lag cross-correlation image-forming condition, the wave field of continuation is imaged, obtains Imaging results corresponding to a beam center.

7) according to step 8, whether all bundles of migration program automatic decision are all imaged and finish, if it is not complete, Start to proceed to calculate from step 4.

8) according to step 9, by big gun, migration imaging is carried out to all of big gun record, obtain shifting into of all big guns As result.

9) according to step 10, whether all big guns of migration program automatic decision are all imaged and finish, if it is not complete, Start to proceed to calculate from step 3.

10) according to step 11, the migration imaging of all big gun records is overlapped and exports, is just imaged Result.

In order to check the imaging effect of the method, we align the model data drilled and have been multi -components joint Gauss Bundle reverse-time migration.Fig. 2 is the oil reservoir section of realistic model.Fig. 3 velocity of longitudinal wave model, Fig. 4 compressional wave inclined Move imaging results, Fig. 5 changes velocity model, the migration imaging result of Fig. 6 converted wave.Permissible from figure Find out, imaging effect is more satisfactory, corresponding with model layer position preferable.This illustrates this imaging method accurately and reliably.

Claims (3)

1. a kind of multi -components joint Gaussian beam pre-Stack Reverse imaging method, it is characterized in that including:

(1) input initial p-and s-wave velocity field；

(2) read in compressional wave and converted wave earthquake big gun record, and determine bandwidth, beam width, beam center interval Parameter；

(3) structure of underground elasticity vector wave field；

(4) cross-correlation image-forming condition is utilized to be imaged；

(5) all big gun records are entered line displacement and are superimposed, obtain final elastic wave imaging result.

2. multi -components joint Gaussian beam pre-Stack Reverse imaging method according to claim 1, it is special Levying is:

Described step (3) is the forward and reverse continuation entering traveling-wave field using elastodynamics Gaussian beam equation；

Described step (4) is using zero-lag cross-correlation image-forming condition, the wave field of continuation to be imaged；

Described step (5) is to carry out migration imaging by big gun to all of big gun record, and all big gun records are inclined Move into as being overlapped, just obtain imaging results.

3. multi -components joint Gaussian beam pre-Stack Reverse imaging method according to claim 2, it is special Levying is,

The structure of step (3) underground elasticity vector wave field includes:

Hypothesis focus is x_s, receiving point x_rThe two component elastic wave earthquake records receiving are u_i(x_r；ω), then instead Can integrate by kirchhoff-helmholtz to represent to the elastic wave field of continuation:

$u_m(x;x_r;\mathrm{\ω})={\&integral;}_s{\mathrm{dx}}_r\[t_i(x_r;\mathrm{\ω})g_{im}^{*}(x;x_r;\mathrm{\ω})-u_i(x_r;\mathrm{\ω}){\mathrm{\σ}}_{im}^{*}(x;x_r;\mathrm{\ω})\]---\left(1\right)$

Wherein, * represents complex conjugate；g_lm(x；x_r；It is ω) displacement Green tensor, it represents x_rPlace's l direction unit The component in the m direction of displacement at x caused by muscle power；t_i(x_r) it is x_rPlace's stress；σ_im(x；x_r) it is stress lattice Woods tensor；

Step (4) utilizes cross-correlation image-forming condition to be imaged:

According to claerbout image-forming principle, by asking for connecing of source wavefield and the backward extension of different wave modes Receive cross-correlation between wave field being calculated as picture value, for this reason, first by focus displacement wave fieldLogical Cross elastodynamics Gaussian beam to represent,

$u_m^p(x;x_s;\mathrm{\ω})\≈\frac{i}{4{\mathrm{\πv}}_p^2\left(x_s\right)}\sqrt{\frac{{\mathrm{\ω}}_r{w}_0^2}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right)}{p_2^p\left(x_s\right)}{\hat{u}}_m^p(x;x_s;\mathrm{\ω})---\left(2\right)$

Next the propagating characteristic according to p ripple and s ripple, is defined as follows imaging formula:

$\begin{array}{c}{i}^{pp}\left(x\right)=\&integral;u_2^{p*}(x;x_s;\mathrm{\ω})u_2^p(x;x_r;\mathrm{\ω})d\mathrm{\ω}\\ =\frac{\mathrm{\δ}l\sqrt{{\mathrm{\ω}}_r{w}_0^2}}{16{\mathrm{\π}}^2}\underset{l}{\σ}\&integral;d\mathrm{\ω}\frac{i\mathrm{\ω}}{v_p^2\left(x_s\right)}\sqrt{\frac{\mathrm{\ρ}\left(l\right)}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right){\mathrm{dp}}_1^p\left(l\right)}{p_2^p\left(x_s\right)p_2^p\left(l\right)}\\ \×{\hat{u}}_2^{p*}(x;x_s;\mathrm{\ω}){\hat{u}}_1^{p*}(x;l;\mathrm{\ω})\[w_1^p\left(l\right)d_1^p(l;p_1^p;\mathrm{\ω})+w_2^p\left(l\right)d_2^p(l;p_1^p;\mathrm{\ω})\]\end{array}---\left(3\right)$

$\begin{array}{c}{i}^{ps}\left(x\right)=\&integral;u_2^{p*}(x;x_s;\mathrm{\ω})u_1^s(x;x_r;\mathrm{\ω})d\mathrm{\ω}\\ =\frac{\mathrm{\δ}l\sqrt{{\mathrm{\ω}}_r{w}_0^2}}{16{\mathrm{\π}}^2}\underset{l}{\σ}\&integral;d\mathrm{\ω}\frac{i\mathrm{\ω}}{v_p^2\left(x_s\right)}\sqrt{\frac{\mathrm{\ρ}\left(l\right)}{\mathrm{\ρ}\left(x_s\right)}}\&integral;\&integral;\frac{{\mathrm{dp}}_1^p\left(x_s\right){\mathrm{dp}}_1^s\left(l\right)}{p_2^p\left(x_s\right)p_2^s\left(l\right)}\\ \×{\hat{u}}_2^{p*}(x;x_s;\mathrm{\ω}){\hat{u}}_1^{s*}(x;l;\mathrm{\ω})\[w_1^s\left(l\right)d_1^s(l;p_1^s;\mathrm{\ω})+w_2^s\left(l\right)d_2^s(l;p_1^s;\mathrm{\ω})\]\end{array}---\left(4\right)$

Wherein, i^ppX () becomes picture value, i for p-p single-shot^psX () becomes picture value for p-s single-shot,(v table The type of oscillography) it is the adding window local dip superposition to different wave mode multi component seismic records, l represents Gaussian beam The effective half width of the Gaussian beam relevant with frequency, focus is x_s, receiving point x_r, ρ is Media density, v_p(x_s) For the velocity of longitudinal wave of central ray at hypocentral location, w is weight coefficient, and p is kinetics ray-tracing parameter.