CN104597484A - Three-dimensional transmission time interval (TTI) earthquake anisotropic medium reverse time migration imaging method and device - Google Patents

Abstract

The invention provides a three-dimensional transmission time interval (TTI) earthquake anisotropic medium reverse time migration imaging method and a device, and is applicable to the technical field of reflected wave earthquake data processing. The method comprises the steps of: determining all guns required to be subject to reverse time migration imaging; executing the following steps for each gun: placing a wavelet at a position corresponding to the gun, and using a coupling second-order partial differential equation to simulate a shot point wave field; performing the coupling second-order partial differential equation on gun data corresponding to the gun, and simulating a detection point wave field; and according to a cross-correlation imaging condition, a shot point wave field simulation result and a detection point wave field simulation result, imaging to obtain a single-gun reverse time migration result. By adopting the differential solution of wave equation by suing the stable coupling second-order partial differential equation, a problem of unstable time caused by mutation of a dip angle of a symmetrical axis by using a TTI medium is solved, and a defect of the three-dimensional complex-structure imaging with greatly changed speed can be overcome.

Description

A kind of three-dimensional TTI seismic anisotropy medium reverse-time migration formation method and device

Technical field

The present invention relates to reflection wave seismic data processing technology field, particularly, relate to a kind of three-dimensional TTI seismic anisotropy medium reverse-time migration formation method and device.

Background technology

The seismic anisotropy of ball medium is proved to be ubiquitous, but seismic prospecting usually regards ground spheric approximation as isotropic medium as.In the past, this approximate problem and the formula that can simplify seismic processing explanation, but pursuing today of meticulous oil reservoir detection, ignore anisotropy that some seismic material (the igneous reservoirs medium as the fractured-porous reservoir type reservoir medium in thin interbed sandstone reservoir medium, laminated slate medium, carbonatite and cranny development) exists seismic processing may be made to explain produce larger error, the dislocation that after one of them of this error shows as skew in seism processing, tomography produces.Therefore anisotropy seismic prospecting technology is studied quite necessary and urgent.

At present, three-dimensional TTI(Tilted Transversely Isotropic, inclination transverse isotropy) seismic anisotropy medium reverse-time migration technology is only serving listed by international each major company, one of several special technologies not for sale, and the domestic research for this technology is in the starting stage, be particularly substantially in blank in conjunction with this technology, the research and development of business software that can be applicable to large-scale industrial production.

In addition, there is calculating instability problem at the mutated site of axis of symmetry dip angle parameter section in prior art, causes calculating to occur exceptional value, calculating is reported an error and exits.

Summary of the invention

The fundamental purpose of the embodiment of the present invention is to provide a kind of three-dimensional TTI seismic anisotropy medium reverse-time migration formation method and device, solving in prior art when carrying out shot point wave-field simulation, existing at the mutated site of axis of symmetry dip angle parameter section and calculating unstable problem.

To achieve these goals, the embodiment of the present invention provides a kind of three-dimensional TTI seismic anisotropy medium reverse-time migration formation method, comprising:

Determine all big guns needing to carry out reverse-time migration imaging;

For each big gun, perform following treatment step:

Wavelet is placed in the sp location corresponding at this big gun, and application coupling partial differential equation of second order carries out shot point wave-field simulation to this big gun;

By to the partial differential equation of second order that is coupled described in big gun market demand corresponding to this big gun, geophone station wave-field simulation is carried out to this big gun;

The result of application cross-correlation image-forming condition to the result of described shot point wave-field simulation and described geophone station wave-field simulation carries out imaging, obtains single big gun reverse-time migration result of this big gun;

After above-mentioned treatment step is all performed to all big guns, single big gun reverse-time migration result of all big guns is stacked up and forms reverse-time migration imaging section;

Wherein, described coupling partial differential equation of second order is:

$\left\{\begin{array}{c}\frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}p}{{\∂t}^2}=(1+2\mathrm{\ϵ})H_{2}p+H_{1}q+\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_1(p-q)\\ \frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}q}{{\∂t}^2}=(1+2\mathrm{\δ})H_{2}p+H_{1}q-\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_2(p-q)\\ H_1={\sin }^2\mathrm{\θ}{\cos }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂x}^2}+{\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂y}^2}+{\cos }^2\mathrm{\θ}\frac{{\∂}^2}{{\∂z}^2}+\\ {\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{\∂x\∂y}+\sin 2\mathrm{\θ}\sin \mathrm{\φ}\frac{{\∂}^2}{\∂y\∂z}+\sin 2\mathrm{\θ}\cos \mathrm{\φ}\frac{{\∂}^2}{\∂x\∂z}\\ H_2=\frac{{\∂}^2}{{\∂x}^2}+\frac{{\∂}^2}{{\∂y}^2}+\frac{{\∂}^2}{{\∂z}^2}-H_1\end{array}\right.$

In described equation, δ, ε are respectively Thomson anisotropic parameters corresponding to imaging space;

θ, φ are respectively axis of symmetry dip angle parameter corresponding to imaging space and axis of symmetry position angle parameter.

The present invention also provides a kind of three-dimensional TTI seismic anisotropy medium reverse-time migration imaging device, comprising:

Big gun determination module, for determining all big guns needing to carry out reverse-time migration imaging;

Single big gun processing module, for for each big gun, performs following treatment step:

Wavelet is placed in the sp location corresponding at this big gun, and application coupling partial differential equation of second order carries out shot point wave-field simulation;

Obtain the big gun data that this big gun is corresponding, and apply described coupling partial differential equation of second order and carry out geophone station wave-field simulation;

The result of application cross-correlation image-forming condition to the result of described shot point wave-field simulation and described geophone station wave-field simulation carries out imaging, obtains single big gun reverse-time migration result of this big gun;

Stacking image module, after all performing above-mentioned treatment step to all big guns, stacks up single big gun reverse-time migration result of all big guns and forms reverse-time migration imaging section;

Wherein, described coupling partial differential equation of second order is:

$\left\{\begin{array}{c}\frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}p}{{\∂t}^2}=(1+2\mathrm{\ϵ})H_{2}p+H_{1}q+\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_1(p-q)\\ \frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}q}{{\∂t}^2}=(1+2\mathrm{\δ})H_{2}p+H_{1}q-\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_2(p-q)\\ H_1={\sin }^2\mathrm{\θ}{\cos }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂x}^2}+{\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂y}^2}+{\cos }^2\mathrm{\θ}\frac{{\∂}^2}{{\∂z}^2}+\\ {\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{\∂x\∂y}+\sin 2\mathrm{\θ}\sin \mathrm{\φ}\frac{{\∂}^2}{\∂y\∂z}+\sin 2\mathrm{\θ}\cos \mathrm{\φ}\frac{{\∂}^2}{\∂x\∂z}\\ H_2=\frac{{\∂}^2}{{\∂x}^2}+\frac{{\∂}^2}{{\∂y}^2}+\frac{{\∂}^2}{{\∂z}^2}-H_1\end{array}\right.$

In described equation, δ, ε are respectively Thomson anisotropic parameters corresponding to imaging space;

θ, φ are respectively axis of symmetry dip angle parameter corresponding to imaging space and axis of symmetry position angle parameter.

By means of technique scheme, present invention employs the difference that stable coupling partial differential equation of second order realizes wave equation to solve, TTI medium axis of symmetry inclination angle can be solved to suddenly change the calculating instability problem caused, and final solution speed 3 D complex structure imaging jumpy problem.Method of the present invention has that counting yield is high, imaging effect good and is easy to the advantage that realizes, is suitable for the exploitation of reverse-time migration commercial software and the needs of suitability for industrialized production.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, below the accompanying drawing used required in describing embodiment is briefly described, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

Fig. 1 is three-dimensional TTI seismic anisotropy medium reverse-time migration formation method schematic flow sheet provided by the invention;

Fig. 2 is the TTI reverse-time migration shot point wave-field simulation snapshot in certain TTI seismic anisotropy medium area provided by the invention;

Fig. 3 is shot point wave field snapshot provided by the invention;

Fig. 4 is geophone station wave field snapshot provided by the invention;

Fig. 5 is single big gun reverse-time migration result provided by the invention;

Fig. 6 is the reverse-time migration imaging results in certain TTI seismic anisotropy medium area provided by the invention;

Fig. 7 is the VTI reverse-time migration imaging results ignoring anisotropy dip angle parameter provided by the invention;

Fig. 8 is the isotropy reverse-time migration imaging results ignoring all anisotropic parameterses provided by the invention;

Fig. 9 is the partial enlarged drawing of Fig. 6;

Figure 10 is the partial enlarged drawing of Fig. 7;

Figure 11 is the partial enlarged drawing of Fig. 8;

Figure 12 is provided by the invention by the result of the result of TTI reverse-time migration and the superimposed display of rate pattern;

Figure 13 is the earthquake isotropy of constant speed model provided by the invention, compare before seismic event in VTI, TTI medium;

Figure 14 is three-dimensional TTI seismic anisotropy medium reverse-time migration image device structure schematic diagram provided by the invention;

Figure 15 is the three-dimensional TTI seismic anisotropy medium reverse-time migration formation method schematic flow sheet that the embodiment of the present invention provides.

Embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

Realizing in process of the present invention, inventor finds:

Based on the finite difference migration method (be also reverse-time migration method, or RTM method) of full acoustic wave equation by temporal-spatial field display high-order finite difference method algorithm direct solution full acoustic wave partial differential equation, simulate wave propagation phenomenon truly.The method observes wave equation completely, and there is not inclination angle restriction, be applicable to the sharply change of velocity field, in 3 D complex structure imaging, possess clear superiority, imaging precision is high.Because reverse-time migration algorithm is direct modeling earthquake wave propagation in time-space domain, it is made more to be easy to be applied to the VTI(Vertical Transversely Isotropic of complicated change, vertical transverse isotropy) and TTI(Tilted Transversely Isotropic, inclination transverse isotropy) etc. the imaging problem in anisotropic medium.

When there is anisotropic character in underground medium, the circulation way of seismic event in underground medium is different, Figure 13 (a) and (b), (c) are the earthquake isotropy of constant speed model respectively, compare before seismic event in VTI, TTI medium, the impact due to anisotropy factor can be found out, the propagation trajectories of seismic event there occurs obvious change, wherein, consider seismic wave propagation track in the vertical transverse isotropy TTI medium of more anisotropy factors the most complicated compared with the first two, seismic wave propagation track truly can be reflected.Therefore, TTI reverse-time migration is compared with isotropic medium reverse-time migration, and image quality is higher.

Based on above-mentioned discovery, the invention provides a kind of three-dimensional TTI seismic anisotropy medium reverse-time migration formation method, as shown in Figure 1, comprising:

Step S11, determines all big guns needing to carry out reverse-time migration imaging;

Step S12, for each big gun, performs following treatment step:

Step S121, wavelet is placed in the sp location corresponding at this big gun, and application coupling partial differential equation of second order carries out shot point wave-field simulation to this big gun;

Step S122, by the partial differential equation of second order that is coupled described in big gun market demand corresponding to this big gun, carries out geophone station wave-field simulation to this big gun;

Step S123, the result of application cross-correlation image-forming condition to the result of described shot point wave-field simulation and described geophone station wave-field simulation carries out imaging, obtains single big gun reverse-time migration result of this big gun;

Step S13, after performing above-mentioned treatment step, stacks up single big gun reverse-time migration result of all big guns and forms reverse-time migration imaging section all big guns;

Wherein, described coupling partial differential equation of second order is:

$\left\{\begin{array}{c}\frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}p}{{\∂t}^2}=(1+2\mathrm{\ϵ})H_{2}p+H_{1}q+\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_1(p-q)\\ \frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}q}{{\∂t}^2}=(1+2\mathrm{\δ})H_{2}p+H_{1}q-\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_2(p-q)\\ H_1={\sin }^2\mathrm{\θ}{\cos }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂x}^2}+{\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂y}^2}+{\cos }^2\mathrm{\θ}\frac{{\∂}^2}{{\∂z}^2}+\\ {\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{\∂x\∂y}+\sin 2\mathrm{\θ}\sin \mathrm{\φ}\frac{{\∂}^2}{\∂y\∂z}+\sin 2\mathrm{\θ}\cos \mathrm{\φ}\frac{{\∂}^2}{\∂x\∂z}\\ H_2=\frac{{\∂}^2}{{\∂x}^2}+\frac{{\∂}^2}{{\∂y}^2}+\frac{{\∂}^2}{{\∂z}^2}-H_1\end{array}\right.$

In described equation, p, q are coupling wave field, and x, y, z is spatial axes coordinate, and t is time shaft coordinate, and δ, ε are respectively Thomson anisotropic parameters corresponding to imaging space;

θ, φ are respectively axis of symmetry dip angle parameter corresponding to imaging space and axis of symmetry position angle parameter.

Apply when three-dimensional TTI seismic anisotropy medium reverse-time migration formation method provided by the invention can solve the sudden change of axis of symmetry inclination angle and calculate unstable problem.Fig. 2 is the application TTI reverse-time migration shot point wave-field simulation snapshot in certain TTI seismic anisotropy medium area that obtains of the present invention and the result of the superimposed display of axis of symmetry dip angle parameter θ, can see, sudden change is there is on the right of axis of symmetry dip angle parameter θ section, traditional equation, when simulating, exists at this mutated site and calculates unstable problem.And the equation that the present invention proposes can solve this instability problem, as can be seen from Figure 2, seismic event have passed this mutated site smoothly, there is not calculating instability problem.

Figure 3 shows that the shot point wave field snapshot instance that application the present invention obtains.

Figure 4 shows that the geophone station wave field snapshot instance that application the present invention obtains.

Figure 5 shows that single big gun reverse-time migration result example that application the present invention obtains.

Figure 6 shows that the reverse-time migration imaging section result in certain TTI seismic anisotropy medium area that application the present invention obtains.

Fig. 7-12 is depicted as the reverse-time migration imaging results in certain TTI seismic anisotropy medium area.Wherein Fig. 7 is the VTI reverse-time migration result ignoring anisotropy dip angle parameter, and Fig. 8 is the isotropy reverse-time migration result ignoring all anisotropic parameterses.Known by comparison diagram 6-8, TTI reverse-time migration imaging results is obviously better than latter two result, particularly the border imaging of high speed body is correct, and latter two result exists error in the imaging on high speed body border, and visible TTI reverse-time migration is for the importance of engineering construction system.Above-mentioned three results amplified, obtain the result of Fig. 9-11, can see in the imaging on high speed body border, TTI reverse-time migration is the most clearly, the border that particularly inclination angle is larger, and the result of VTI and isotropy reverse-time migration all exists comparatively big error.Figure 12, by the result of TTI reverse-time migration and the superimposed display of rate pattern, can see that this result can reflect true underground structure.

Present invention employs the difference that stable coupling partial differential equation of second order realizes wave equation to solve, TTI medium axis of symmetry inclination angle can be solved and to suddenly change the calculating instability problem caused, and final solution speed 3 D complex structure imaging jumpy problem.Method of the present invention has that counting yield is high, imaging effect good and is easy to the advantage that realizes, is suitable for the exploitation of reverse-time migration commercial software and the needs of suitability for industrialized production.

Accordingly, the invention provides a kind of three-dimensional TTI seismic anisotropy medium reverse-time migration imaging device, as shown in figure 14, this device comprises:

Big gun determination module 1401, for determining all big guns needing to carry out reverse-time migration imaging;

Single big gun processing module 1402, for for each big gun, performs following treatment step:

Wavelet is placed in the sp location corresponding at this big gun, and application coupling partial differential equation of second order carries out shot point wave-field simulation;

Obtain the big gun data that this big gun is corresponding, and apply described coupling partial differential equation of second order and carry out geophone station wave-field simulation;

The result of application cross-correlation image-forming condition to the result of described shot point wave-field simulation and described geophone station wave-field simulation carries out imaging, obtains single big gun reverse-time migration result of this big gun;

Stacking image module 1403, after all performing above-mentioned treatment step to all big guns, stacks up single big gun reverse-time migration result of all big guns and forms reverse-time migration imaging section;

Wherein, described coupling partial differential equation of second order is:

$\left\{\begin{array}{c}\frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}p}{{\∂t}^2}=(1+2\mathrm{\ϵ})H_{2}p+H_{1}q+\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_1(p-q)\\ \frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}q}{{\∂t}^2}=(1+2\mathrm{\δ})H_{2}p+H_{1}q-\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_2(p-q)\\ H_1={\sin }^2\mathrm{\θ}{\cos }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂x}^2}+{\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂y}^2}+{\cos }^2\mathrm{\θ}\frac{{\∂}^2}{{\∂z}^2}+\\ {\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{\∂x\∂y}+\sin 2\mathrm{\θ}\sin \mathrm{\φ}\frac{{\∂}^2}{\∂y\∂z}+\sin 2\mathrm{\θ}\cos \mathrm{\φ}\frac{{\∂}^2}{\∂x\∂z}\\ H_2=\frac{{\∂}^2}{{\∂x}^2}+\frac{{\∂}^2}{{\∂y}^2}+\frac{{\∂}^2}{{\∂z}^2}-H_1\end{array}\right.$

In described equation, δ, ε are respectively Thomson anisotropic parameters corresponding to imaging space;

θ, φ are respectively axis of symmetry dip angle parameter corresponding to imaging space and axis of symmetry position angle parameter.

Embodiment

The present embodiment provides a kind of specific embodiment three-dimensional TTI seismic anisotropy medium reverse-time migration formation method of the present invention being applied to commercial software, as shown in figure 15, specifically comprises the steps:

Steps A 1, stores all big gun data in this locality, and Depth Domain velocity field, Thomson anisotropic parameters, axis of symmetry dip angle parameter and axis of symmetry position angle parameter that imaging space is corresponding;

Steps A 2, current big gun to be processed is determined from task list, read the big gun data that this big gun is corresponding, and Depth Domain velocity field, Thomson anisotropic parameters, axis of symmetry dip angle parameter and axis of symmetry position angle parameter that imaging space is corresponding, and perform following process for this big gun:

Steps A 21, wavelet is placed in the sp location corresponding at this big gun, and application coupling partial differential equation of second order carries out shot point wave-field simulation (applying calculation stability when this equation can make to suddenly change at axis of symmetry inclination angle), and the shot point wave-field simulation result obtained is stored in this locality;

Steps A 22, the big gun market demand coupling partial differential equation of second order corresponding to this big gun, realizes geophone station wave-field simulation, the detection wave-field simulation result obtained is stored in this locality;

Steps A 23, the shot point wave-field simulation result that application cross-correlation image-forming condition obtains steps A 21 and the geophone station wave-field simulation result that steps A 22 obtains carry out imaging, obtain single big gun reverse-time migration result of this big gun, and are stored in this locality;

Steps A 3, judges the big gun whether not carrying out processing in addition in task list, and if so, then circulation execution steps A 2(comprises steps A 21 ~ steps A 23); Otherwise, perform steps A 4;

Steps A 4, stacks up single big gun reverse-time migration result of all big guns and forms reverse-time migration imaging section, Output rusults.

In order to improve the processing speed of steps A 23 in the present embodiment, can when performing steps A 21, the time interval according to setting is compressed shot point wave-field simulation result substep, then the compression result corresponding each time interval is stored in this locality, afterwards when performing steps A 23, can decompress by the compressed package of point thread synchronization to shot point wave-field simulation result, and then apply cross-correlation image-forming condition imaging is carried out to shot point wave-field simulation result and geophone station wave-field simulation result.Owing to have employed substep compression and point thread synchronization decompression, improve the speed of single big gun process, and then improve the treatment effeciency of whole process.

In sum, the three-dimensional TTI seismic anisotropy medium reverse-time migration formation method that provides of the embodiment of the present invention and device have following beneficial effect:

(1) have employed the difference that stable coupling partial differential equation of second order realizes wave equation to solve, TTI medium axis of symmetry inclination angle can be solved and to suddenly change the calculating instability problem caused, and final solution speed 3 D complex structure imaging jumpy problem;

(2) have that counting yield is high, imaging effect good and be easy to the advantage that realizes, be suitable for the exploitation of reverse-time migration (RTM) commercial software and the needs of suitability for industrialized production.

Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (4)

1. a three-dimensional TTI seismic anisotropy medium reverse-time migration formation method, is characterized in that, comprising:

Determine all big guns needing to carry out reverse-time migration imaging;

For each big gun, perform following treatment step:

Wavelet is placed in the sp location corresponding at this big gun, and application coupling partial differential equation of second order carries out shot point wave-field simulation to this big gun;

By to the partial differential equation of second order that is coupled described in big gun market demand corresponding to this big gun, geophone station wave-field simulation is carried out to this big gun;

The result of application cross-correlation image-forming condition to the result of described shot point wave-field simulation and described geophone station wave-field simulation carries out imaging, obtains single big gun reverse-time migration result of this big gun;

After above-mentioned treatment step is all performed to all big guns, single big gun reverse-time migration result of all big guns is stacked up and forms reverse-time migration imaging section;

Wherein, described coupling partial differential equation of second order is:

$\left\{\begin{array}{c}\frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}p}{{\∂t}^2}=(1+2\mathrm{\ϵ})H_{2}p+H_{1}q+\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_1(p-q)\\ \frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}q}{{\∂t}^2}=(1+2\mathrm{\δ})H_{2}p+H_{1}q-\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_2(p-q)\\ H_1={\sin }^2\mathrm{\θ}{\cos }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂x}^2}+{\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂y}^2}+{\cos }^2\mathrm{\θ}\frac{{\∂}^2}{{\∂z}^2}+\\ {\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{\∂x\∂y}+\sin 2\mathrm{\θ}\sin \mathrm{\φ}\frac{{\∂}^2}{\∂y\∂z}+\sin 2\mathrm{\θ}\cos \mathrm{\φ}\frac{{\∂}^2}{\∂x\∂z}\\ H_2=\frac{{\∂}^2}{{\∂x}^2}+\frac{{\∂}^2}{{\∂y}^2}+\frac{{\∂}^2}{{\∂z}^2}-H_1\end{array}\right.$

In described equation, δ, ε are respectively Thomson anisotropic parameters corresponding to imaging space;

θ, φ are respectively axis of symmetry dip angle parameter corresponding to imaging space and axis of symmetry position angle parameter.

2. method according to claim 1, is characterized in that, also comprises: compress the shot point wave-field simulation result at setting-up time interval.

3. a three-dimensional TTI seismic anisotropy medium reverse-time migration imaging device, is characterized in that, comprising:

Big gun determination module, for determining all big guns needing to carry out reverse-time migration imaging;

Single big gun processing module, for for each big gun, performs following treatment step:

Wavelet is placed in the sp location corresponding at this big gun, and application coupling partial differential equation of second order carries out shot point wave-field simulation;

Obtain the big gun data that this big gun is corresponding, and apply described coupling partial differential equation of second order and carry out geophone station wave-field simulation;

The result of application cross-correlation image-forming condition to the result of described shot point wave-field simulation and described geophone station wave-field simulation carries out imaging, obtains single big gun reverse-time migration result of this big gun;

Stacking image module, after all performing above-mentioned treatment step to all big guns, stacks up single big gun reverse-time migration result of all big guns and forms reverse-time migration imaging section;

Wherein, described coupling partial differential equation of second order is:

$\left\{\begin{array}{c}\frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}p}{{\∂t}^2}=(1+2\mathrm{\ϵ})H_{2}p+H_{1}q+\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_1(p-q)\\ \frac{1}{v_{\mathrm{pz}}^2}\frac{{\∂}^{2}q}{{\∂t}^2}=(1+2\mathrm{\δ})H_{2}p+H_{1}q-\frac{\mathrm{\ϵ}-\mathrm{\δ}}{\mathrm{\σ}}{H}_2(p-q)\\ H_1={\sin }^2\mathrm{\θ}{\cos }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂x}^2}+{\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{{\∂y}^2}+{\cos }^2\mathrm{\θ}\frac{{\∂}^2}{{\∂z}^2}+\\ {\sin }^2\mathrm{\θ}{\sin }^2\mathrm{\φ}\frac{{\∂}^2}{\∂x\∂y}+\sin 2\mathrm{\θ}\sin \mathrm{\φ}\frac{{\∂}^2}{\∂y\∂z}+\sin 2\mathrm{\θ}\cos \mathrm{\φ}\frac{{\∂}^2}{\∂x\∂z}\\ H_2=\frac{{\∂}^2}{{\∂x}^2}+\frac{{\∂}^2}{{\∂y}^2}+\frac{{\∂}^2}{{\∂z}^2}-H_1\end{array}\right.$

In described equation, δ, ε are respectively Thomson anisotropic parameters corresponding to imaging space;

θ, φ are respectively axis of symmetry dip angle parameter corresponding to imaging space and axis of symmetry position angle parameter.

4. device according to claim 3, is characterized in that, also comprises:

Compression module, for compressing the shot point wave-field simulation result at setting-up time interval.