叠前时间偏移: 基于波动方程的架构、发展和应用

This dissertation starts from the point that the prestack time migration can been considered as an approximation of the prestack depth migration, giving a wave equation based prestack time migration approach. The new approach includes: analytically getting the travel time and amplitude based on the one way wave equation and the stationary-phase theory, using ‘spread’ imaging method and imaging following the prestack depth migration, updating the velocity model with respect to the flats of the events in CRP gathers. Based on this approach, we present a scheme that can image land seismic data without field static correction. We may determine the correct near surface velocities and stack velocities by picking up the residual correction of the events in the CRP gathers. We may get the rational migration section based on the updated velocities and correct the migration section from a floating datum plane to a universal datum plane. We may adaptively determine the migration aperture according to the dips of the imaging structures. This not only speed up the processing, but may suppress the migration noise produce by the extra aperture. We adopt the deconvolution imaging condition of wave equation migration. It may partially compensate the geometric divergence. In this scheme, we use the table-driven technique which may enhance the computational efficiency. If the subsurface is much more complicated, it may be impossible to distinguish the DTS curve. To solve this problem, we proposed a technique to determine the appropriate range of the DTS curve. We synthesize DTS panel in this range using different velocities and depths, and stack the amplitude around the zero time. Determine the correct velocity and location of the considered grid point by comparing the values.

波动方程模拟与偏移的实用化方法研究

Along with the widespread and in-depth applications in petroleum prospecting and development, the seismic modeling and migration technologies are proposed with a higher requirement by oil industrial, and the related practical demand is getting more and more urgent. Based on theories of modeling and migration methods for wave equation, both related with velocity model, I thoroughly research and develop some methods for the goal of highly effective and practical in this dissertation. In the first part, this dissertation probes into the layout designing by wave equations modeling, focusing on the target-oriented layout designing method guided by wave equation modeling in complicated structure areas. It is implemented by using the fourth order staggered grid finite difference (FD) method in velocity-stress 2D acoustic wave equations plus perfectly matched layer (PML) absorbing boundary condition. To design target-oriented layout: (a) match the synthetic record on the surface with events of subsurface structures by analyzing the snapshots of theoretical model; (b) determine the shot-gather distance by tracking the events of target areas and measuring the receiving range when it reaches the surface; (c) restrict the range of valid shot-gather distance by drawing seismic windows in single shot records; (d) choose the best trace distance by comparing the resolution of prospecting targets from the simulated records with different trace distance. Eventually, we obtained the observation system parameters, which achieve the design requirements. In the second part, this dissertation presents the practical method to improve the 3D Fourier Finite Difference (FFD) migration, and carefully analyzes all the factors which influence 3D FFD migration’s efficiency. In which, one of the most important parameters of migration is the extrapolating step. This dissertation presents an efficient 3D FFD migration algorithm, which use FFD propagator to extrapolate wavefields over big layers, and use Born-Kirchhoff interpolator to image wavefields over small layers between the big ones. Finally, I show the effectiveness of this hybrid migration method by comparing migration results from 3D SEG/EAGE model with different methods.

反射地震转换波成像方法研究

Multi-waves and multi-component get more and more attentions from oil industry. On the basis of existent research results, My research focuses on some key steps of OBC 4C datum processing. OBC datum must be preprocessed quite well for getting a good image. We show a flow chart of preprocess including attenuation of noise on multi-component datum、elimination ghost by summing P and Z and rotation of horizontal components. This is a good foundation for the coming steps about OBC processing. How to get exact converted point location and to analyze velocity are key points in processing reflection seismic converted wave data. This paper includes computing converted point location, analyzing velocity and nonhyperbolic moveout about converted waves. Anisotropic affects deeply the location of converted wave and the nonhyperbolic moveout. Supposed VTI, we research anisotropic effect on converted wave location and the moveout. Since Vp/Vs is important, we research the compute method of Vp/Vs from post-stack data and pre-stack data. It is a part of the paper that inversing anisotropic parameter by traveltime. Pre-stack time migration of converted wave is an focus, using common-offset Kirchhoff migration, we research the velocity model updating in anisotropic media. I have achieved the following results: 1) using continued Fractions, we proposed a new converted point approximate equation, when the offset is long enough ,the thomsen’s 2 order equation can’t approximate to the exact location of converted point, our equation is a good approximate for the exact location. 2) our new methods about scanning nonhyperbolic velocity and Vp/Vs can get a high quality energy spectrum. And the new moveout can fit the middle and long offset events. Processing the field data get a good result. 3) a new moveout equation, which have the same form as Alkhalifah’s long offset P wave moveout equation, have the same degree preciseness as thomsen’s moveout equation by testing model data. 4) using c as a function of the ratio offset to depth, we can uniform the Li’s and thomsen’s moveout equation in a same equation, the model test tell us choice the reasonable function C can improve the exact degree of Li’s and thomsen’s equation. 5) using traveltime inversion ,we can get anisotropic parameter, which can help to flat the large offset event and propose a model of anisotropic parameter which will useful for converted wave pre-stack time migration in anisotropic media. 6）using our pre-stack time migration method and flow, we can update the velocity model and anisotropic parameter model then get good image. Key words: OBC, Common converted Point (CCP), Nonhyperbolic moveout equation, Normal moveout correction, Velocity analysis, Anisotropic parameters inversion, Kirchhoff anisotropic pre-stack time migration, migration velocity model updating

表面相关多次波衰减方法与应用研究

The function of seismic data in prospecting and exploring oil and gas has exceeded ascertaining structural configuration early. In order to determine the advantageous target area more exactly, we need exactly image the subsurface media. So prestack migration imaging especially prestack depth migration has been used increasingly widely. Currently, seismic migration imaging methods are mainly based on primary energy and most of migration methods use one-way wave equation. Multiple will mask primary and sometimes will be regarded as primary and interferes with the imaging of primary, so multiple elimination is still a very important research subject. At present there are three different wavefield prediction and subtraction methods: wavefield extrapolation; feedback loop; and inverse-scattering series. I mainly do research on feedback loop method in this paper. Feedback loop method includs prediction and subtraction.Currently this method has some problems as follows. Firstly, feedback loop method requires the seismic data used to predict multiple is full wavefield data, but usually the original seismic data don’t meet this assumption, so seismic data must be regularized. Secondly, Multiple predicted through feedback loop method usually can’t match the real multiple in seismic data and they are different in amplitude, phase and arrrival time. So we need match the predicted multiple and that in seismic data through estimating filtering factors and subtract multiple from seismic data. It is the key for multiple elimination how to select a correct matching filtering method. There are many matching filtering methods and I put emphasis on Least-square adaptive matching filtering and L1-norm minimizing adaptive matching filtering methods. Least-square adaptive matching filtering method is computationally very fast, but it has two assumptions: the signal has minimum energy and is orthogonal to the noise. When seismic data don’t meet the two assumptions, this method can’t get good matching results and then can’t attenuate multiple correctly. L1-norm adaptive matching filtering methods can avoid these two assumptions and then get good matching results, but this method is computationally a little slow. The results of my research are as follows: 1. Proposed a method that interpolates seismic traces based on F-K migration and demigration. The main advantage of this method is that it can interpolate seismic traces in any offsets. It shows this method is valid through a simple model. 2. Comparing different Least-square adaptive matching filtering methods. The results show that equipose multi-channel adaptive matching filtering methods can get better results of multiple elimination than other matcing methods through three model data and two field data. 3. Proposed equipose multi-channel L1-norm adaptive matching filtering method. Because L1-norm is robust to large amplitude differences, there are no assumption on the signal has minimum energy and orthogonality, this method can get better results of multiple elimination. 4. Research on multiple elimination in inverse data space. The method is a new multiple elimination method and it is different from those methods mentioned above.The advantages of this method is that it is simple in theory and no need for the adaptive subtraction and computationally very fast. The disadvantage of this method is that it is not stabilized in its solution. The results show that equipose multi-channel and equipose pesudo-multi-channel least-square matching filtering and equipose multi-channel and equipose pesudo-multi-channel L1-norm matching filtering methods can get better results of multiple elimination than other matcing methods through three model data and many field data.

叠前时间偏移：各向异性、多次波衰减及应用

At present, in order to image complex structures more accurately, the seismic migration methods has been developed from isotropic media to the anisotropic media. This dissertation develops a prestack time migration algorithm and application aspects for complex structures systematically. In transversely isotropic media with a vertical symmetry axis (VTI media), the dissertation starts from the theory that the prestack time migration is an approximation of the prestack depth migration, based on the one way wave equation and VTI time migration dispersion relation, by combining the stationary-phase theory gives a wave equation based VTI prestack time migration algorithm. Based on this algorithm, we can analytically obtain the travel time and amplitude expression in VTI media, as while conclude how the anisotropic parameter influence the time migration, and by analyzing the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, we can update the velocity model and estimate the anisotropic parameter model through the time migration. When anisotropic parameter is zero, this algorithm degenerates to the isotropic time migration algorithm naturally, so we can propose an isotopic processing procedure for imaging. This procedure may keep the main character of time migration such as high computational efficiency and velocity estimation through the migration, and, additionally, partially compensate the geometric divergence by adopting the deconvolution imaging condition of wave equation migration. Application of this algorithm to the complicated synthetic dataset and field data demonstrates the effectiveness of the approach. In the dissertation we also present an approach for estimating the velocity model and anisotropic parameter model. After analyzing the velocity and anisotropic parameter impaction on the time migration, and based on the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, through migration we can update the velocity model and estimate the anisotropic parameter model by combining the advantages of velocity analysis in isotropic media and anisotropic parameter estimation in VTI media. Testing on the synthetic and field data, demonstrates the method is effective and very steady. Massive synthetic dataset、2D sea dataset and 3D field datasets are used for VTI prestack time migration and compared to the stacked section after NMO and prestack isotropic time migration stacked section to demonstrate that VTI prestack time migration method in this paper can obtain better focusing and less positioning errors of complicated dip reflectors. When subsurface is more complex, primaries and multiples could not be separated in the Radon domain because they can no longer be described with simple functions (parabolic). We propose an attenuating multiple method in the image domain to resolve this problem. For a given velocity model，since time migration takes the complex structures wavefield propagation in to account, primaries and multiples have different offset-domain moveout discrepancies, then can be separated using techniques similar to the prior migration with Radon transform. Since every individual offset-domain common-reflection point gather incorporates complex 3D propagation effects, our method has the advantage of working with 3D data and complicated geology. Testing on synthetic and real data, we demonstrate the power of the method in discriminating between primaries and multiples after prestack time migration, and multiples can be attenuated in the image space considerably.

人工源宽角地震资料中的广角效应及其对地壳速度成像的影响

The seismic wide-angle reflection/refraction method is the one of the most effective method for probing the crustal and upper mantle structure. It mainly uses the wide-angle reflection information from the boundary in the crust and the top boundary of the upper mantle to rebuild the crust and upper mantle structure. Through analyzing the reflection and transmission coefficients of various incident waves on the interface, we think relative to the pre-critical angle reflection information the post critical angle reflection information that received by wide-angle seismic data exists a time-shift effect with the offset variation, and then it must cause the error for velocity analysis and structure image. The feature of the wide-angle seismic wave field of the fourteen representative crust columns tell us that the wide-angle effects in the different representative tectonic units for the interface depth and the interval velocity in crust. We studied the features of the wide-angle seismic wave field through building the crust model and inverse its travel time by GA method to know the wide-angle influence on crustal velocity image. At last we finished the data processing of the Tunxi-Wenzhou wide-angle seismic profile. The results are as following: (1) Through building crust model, we labeled the travel time for all the phases by ray tracing method and remove wide-angle effects method, it revealed the wide-angle effect exists in the seismic data. (2) The travel time inversion by GA method can tell us that the depth by traditional ray tracing method is shallower than the result by remove wide-angle effects method, the latter can recover the crust structure model in effect. (3) We applied the two method mentioned before to the fourteen representative crust columns in China. It indicates that the removed wide-angle effect method in travel time inversion is reasonable and effective. (4) The real data processing from Tunxi-Wenzhou wide-angle seismic profile give us the basic structure through the two ways. The main influence exhibits in the difference of the interval velocity of the curst, and the wide-angle effects in shallow interface are stronger than the deep interface.

波动方程偏移速度分析

On the subject of oil and gas exploration, migration is an efficacious technique for imagining structures underground. Wave-equation migration (WEM) dominates over other migration methods in accuracy, despite of higher computational cost. However, the advantages of WEM will emerge as the progress of computer technology. WEM is sensitive to velocity model more than others. Small velocity perturbations result in grate divergence in the image pad. Currently, Kirrchhoff method is still very popular in the exploration industry for the reason of difficult to provide precise velocity model. It is very urgent to figure out a way to migration velocity modeling. This dissertation is mainly devoted to migration velocity analysis method for WEM: 1. In this dissertation, we cataloged wave equation prestack depth migration. The concept of migration is introduced. Then, the analysis is applied to different kinds of extrapolate operator to demonstrate their accuracy and applicability. We derived the DSR and SSR migration method and apply both to 2D model. 2. The output of prestack WEM is in form of common image gathers (CIGs). Angle domain common image gathers (ADCIGs) gained by wave equation are proved to be free of artifacts. They are also the most potential candidates for migration velocity analysis. We discussed how to get ADCIGs by DSR and SSR, and obtained ADCIGs before and after imaging separately. The quality of post stack image is affected by CIGs, only the focused or flattened CIGs generate the correct image. Based on wave equation migration, image could be enhanced by special measures. In this dissertation we use both prestack depth residual migration and time shift imaging condition to improve the image quality. 3. Inaccurate velocities lead to errors of imaging depth and curvature of coherent events in CIGs. The ultimate goal of migration velocity analysis (MVA) is to focus scattered event to correct depth and flatten curving event by updating velocities. The kinematic figures are implicitly presented by focus depth aberration and kinetic figure by amplitude. The initial model of Wave-equation migration velocity analysis (WEMVA) is the output of RMO velocity analysis. For integrity of MVA, we review RMO method in this dissertation. The dissertation discusses the general ideal of RMO velocity analysis for flat and dipping events and the corresponding velocity update formula. Migration velocity analysis is a very time consuming work. Respect to computational convenience, we discus how RMO works for synthetic source record migration. In some extremely situation, RMO method fails. Especially in the areas of poorly illuminated or steep structure, it is very difficult to obtain enough angle information for RMO. WEMVA based on wave extrapolate theory, which successfully overcome the drawback of ray based methods. WEMVA inverses residual velocities with residual images. Based on migration regression, we studied the linearized scattering operator and linearized residual image. The key to WEMVA is the linearized residual image. Residual image obtained by Prestack residual migration, which based on DSR is very inefficient. In this dissertation, we proposed obtaining residual migration by time shift image condition, so that, WEMVA could be implemented by SSR. It evidently reduce the computational cost for this method.

大别－苏鲁造山带与华北陆块东部岩石层三维速度结构

The 3-D velocity images of the crest and upper mantle beneath the region of 112° -124°E, 28°-39°N including the Dabie-Sulu orogenic belt are reconstructed by using 36405 P-wave arrivals of 3437 regional and 670 distant earthquakes during the period from 1981 to 1996, and gridding the area of 0.5° * 0.5°. The results of tomography demonstrate that: 1. The results of tomographic imaging show a broad heterogeneity in P wave velocity structure for the lithosphere beneath the Dabie-Sulu orogenic belt. 2. In the Dabie orogenic belt, the velocity patterns in the crust are different among various tectonic units. The Dabie and Qinling orogenic belts are remarkable in the tomographic images, and in mm the Hongan and Dabie blocks in the Dabie orogenic belt are also imaged very distinguishably. 3. A velocity (about 5.9～6.0 km/s) layer exists in the Dabie block at depth between 15～25 km, which is coincident with the low-resistance layer at the depth of 12-23 km, being inferred to be the tectonic detachment zone and suggesting that the extension detachment structure was formed in the middle crust. Beneath the southern and northerm Dabie tectonic units, the north-dipping high-velocity (at level of 6.5 ～ 6.6 km/s) block was developed in the crust, which might be correlated with the UHP rockswith low content of the meta-ultramafic rocks. This result is in agreement with the geological observation on the surface. 4. The velocity image at 40 km depth reveals the features at the top of mantle and the configuration of the Moho discontinuity. The depth of the Moho changes slightly along the trend of the orogenic belt. It in Hongan block is less than 40 km, but it is different in the western and eastern parts of the Dabie block, the former is more than 40 km, and the latter less than or equal to 40 km. The remnant of the mountain root exists between the Shangcheng-Macheng fault and the line of Huoshan-Yuexi-Yingshan in the Dabie orogenic belt, and beneath the southern and northern Dabie tectonic units. However, the thickness of the Moho is about 40 km and there is no obvious changes, which suggest that the Dabie orogenic belt has been experienced quite in the gravity equilibration. The Moho's depth in the Sulu is less than 40 km. 5. There is a dipping slab-like high-velocity body in the uppermost mantle. It is sandwiched by slow velocities and exists beneath the Dabie-Sulu orogenic belt in the range of depths between the Moho discontinuity and 110 km at least. This high-velocity body outlines a picture of the slab interpreted as the remnant of the Triassic subducted YZ. 6. The Sulu orogenic belt displays "crocodilian" velocity structure, the upper crust of the Yangtze thrusted over the Huabei crest, and the Huabei crust indented into the Yangtze crust, where the ancient subduction zone of the Yangtze lithosphere located. Based on the previous geological data, this structure is not related with the collision between the Yangtze and Sino-Korean Blocks, but caused by the sinistral offset of the Tan-Lu Fault. Studied on the velocity structure of the eastern Huabei lithosphere indicates: 1. The 'present-day' lithosphere of the eastern Huabei is between 40-100 km thick with greatly thinned lithosphere around the Bohai Sea. Generally, thickness of the lithosphere in this region decreased eastwards. 2. The attenuation of the lithosphere is attributed to the strongly uplift of the asthenosphere. In the area between the Taihang Mountains and the Tan-Lu Fault, there is a 'lever' with red low velocity belt, it is clearly defined, transverse continuity, depth between 100-150 km, local variations visible, and an upwards trend towards the Bohai Sea. Generally, the velocity structure in the mantle beneath the lithosphere displays irregular column-shape consisting of alternating high and low velocities, and when cold high velocity ancient lithosphere connects with the hot low velocity mantle materials forming precipitous compact structure. More heat pathways from the mantle occur towards the Tan-Lu Fault. 3. The strongly irregular characteristics of the contact between the asthenosphere and the lithosphere is induced by the long-term hot, chemical erosion and alteration on the contact. 4. There are still preserved high velocity lithosphedc root beneath Huabei with 'block-shape' distribution and surrounded by hot materials. Results of our studies indicate that the evolution models of the eastern China mantle are characterized by the direct contact between the uplifted lithosphere and the Huabei Craton accompanying the upwelling of the deep mantle materials. At the contact betwen the lithosphere and the asthenosphere, the upwelled mantle materials replaced and altered the lower lithosphere forming the metasome through the hot and chemical modifications impacted on the Craton lithosphere, and changed it into the lithosphere gradually, resulting in the lithospheric thinning. Thus, the lithospheric thinning is the result of the upwelling of the asthenosphere.

改善三维波动方程叠前深度偏移处理效果的方法研究

3D wave equation prestack depth migration is the effective tool for obtaining the exact imaging result of complex geology structures. It's a part of the 3D seismic data processing. 3D seismic data processing belongs to high dimension signal processing, and there are some difficult problems to do with. They are: How to process high dimension operators? How to improve the focusing? and how to construct the deconvolution operator? The realization of 3D wave equation prestack depth migration, not only realized the leap from poststack to prestack, but also provided the important means to solve the difficult problems in high dimension signal processing. In this thesis, I do a series research especially for the solve of the difficult problems around the 3D wave equation prestack depth migration and using it as a mean. So this thesis service for the realization of 3D wave equation prestack depth migration for one side and improve the migration effect for another side. This thesis expatiates in five departs. Summarizes the main contents as the follows: In the first part, I have completed the projection from 3D data point area to low dimension are using de big matrix transfer and trace rearrangement, and realized the liner processing of high dimension signal. Firstly, I present the mathematics expression of 3D seismic data and the mean according to physics, present the basic ideal of big matrix transfer and describe the realization of five transfer models for example. Secondly, I present the basic ideal and rules for the rearrange and parallel calculate of 3D traces, and give a example. In the conventional DMO focusing method, I recall the history of DM0 process firstly, give the fundamental of DMO process and derive the equation of DMO process and it's impulse response. I also prove the equivalence between DMO and prestack time migration, from the kinematic character of DMO. And derive the relationship between DMO base on wave equation and prestack time migration. Finally, I give the example of DMO process flow and synthetic data of theoretical models. In the wave equation prestak depth migration, I firstly recall the history of migration from time to depth, from poststack to prestack and from 2D to 3D. And conclude the main migration methods, point out their merit and shortcoming. Finally, I obtain the common image point sets using the decomposed migration program code.In the residual moveout, I firstly describe the Viterbi algorithm based on Markov process and compound decision theory and how to solve the shortest path problem using Viterbi algorithm. And based on this ideal, I realized the residual moveout of post 3D wave equation prestack depth migration. Finally, I give the example of residual moveout of real 3D seismic data. In the migration Green function, I firstly give the concept of migration Green function and the 2D Green function migration equation for the approximate of far field. Secondly, I prove the equivalence of wave equation depth extrapolation algorithms. And then I derive the equation of Green function migration. Finally, I present the response and migration result of Green function for point resource, analyze the effect of migration aperture to prestack migration result. This research is benefit for people to realize clearly the effect of migration aperture to migration result, and study on the Green function deconvolution to improve the focusing effect of migration.

地震波场延拓+菲涅尔体叠加方法及成像效果

In the increasingly enlarged exploration target, deep target layer(especially for the reservoir of lava) is a potential exploration area. As well known, the reflective energy becomes weak because the seismic signals of reflection in deep layer are absorbed and attenuate by upper layer. Caustics and multi-values traveltime in wavefield are aroused by the complexity of stratum. The ratio of signal to noise is not high and the fold numbers are finite(no more than 30). All the factors above affect the validity of conventional processing methods. So the high S/N section of stack can't always be got with the conventional stack methods even if the prestack depth migration is used. So it is inevitable to develop another kind of stack method instead. In the last a few years, the differential solution of wave equation was hold up by the condition of computation. Kirchhoff integral method rose in the initial stages of the ninetieth decade of last century. But there exist severe problems in it, which is are too difficult to resolve, so new method of stack is required for the oil and gas exploration. It is natural to think about upgrading the traditionally physic base of seismic exploration methods and improving those widely used techniques of stack. On the other hand, great progress is depended on the improvement in the wave differential equation prestack depth migration. The algorithm of wavefield continuation in it is utilized. In combination with the wavefield extrapolation and the Fresnel zone stack, new stack method is carried out It is well known that the seismic wavefield observed on surface comes from Fresnel zone physically, and doesn't comes from the same reflection points only. As to the more complex reflection in deep layer, it is difficult to describe the relationship between the reflective interface and the travel time. Extrapolation is used to eliminate caustic and simplify the expression of travel time. So the image quality is enhanced by Fresnel zone stack in target. Based on wave equation, high-frequency ray solution and its character are given to clarify theoretical foundation of the method. The hyperbolic and parabolic travel time of the reflection in layer media are presented in expression of matrix with paraxial ray theory. Because the reflective wave field mainly comes from the Fresnel Zone, thereby the conception of Fresnel Zone is explained. The matrix expression of Fresnel zone and projected Fresnel zone are given in sequence. With geometrical optics, the relationship between object point in model and image point in image space is built for the complex subsurface. The travel time formula of reflective point in the nonuniform media is deduced. Also the formula of reflective segment of zero-offset and nonzero offset section is provided. For convenient application, the interface model of subsurface and curve surface derived from conventional stacks DMO stack and prestack depth migration are analyzed, and the problem of these methods was pointed out in aspects of using data. Arc was put forward to describe the subsurface, thereby the amount of data to stack enlarged in Fresnel Zone. Based on the formula of hyperbolic travel time, the steps of implementation and the flow of Fresnel Zone stack were provided. The computation of three model data shows that the method of Fresnel Zone stack can enhance the signal energy and the ratio of signal to noise effectively. Practical data in Xui Jia Wei Zhi, a area in Daqing oilfield, was processed with this method. The processing results showed that the ability in increasing S/N ratio and enhancing the continuity of weak events as well as confirming the deep configuration of volcanic reservoir is better than others. In deeper target layer, there exists caustic caused by the complex media overburden and the great variation of velocity. Travel time of reflection can't be exactly described by the formula of travel time. Extrapolation is bring forward to resolve the questions above. With the combination of the phase operator and differential operator, extrapolating operator adaptable to the variation of lateral velocity is provided. With this method, seismic records were extrapolated from surface to any different deptlis below. Wave aberration and caustic caused by the inhomogenous layer overburden were eliminated and multi-value curve was transformed into the curve.of single value. The computation of Marmousi shows that it is feasible. Wave field continuation extends the Fresnel Zone stack's application.