VTI各向异性叠前深度偏移及速度分析

Theoretical research, laboratory test and field observation show that most of sediment rock has anisotropic features. It will produce some notable errors when applying isotropic methods such as prestack depth migration and velocity analysis to dada acquired under anisotropic condition; it also has a bad effect on geologic interpretation. Generally speaking, the vertical transverse isotropic media is a good approximation to geologic structure, thus it has an important realistic meaning for anisotropic prestack depth migration theory researching and precise complex geologic imaging if considering anisotropic effect of seismic wave propagation. There are two indispensable parts in prestack depth migration of realistic records, one is proper prestack depth migration algorithm, and the other is velocity analysis using prestack seismic data. The paper consists of the two aspects. Based on implicit finite difference research proposed by Dietrich Ristow et al （1997） about VTI media prestack depth migration, the paper proposed split-step Fourier prestack depth migration algorithm （VTISSF） and Fourier finite difference algorithm （VTIFFD） based on wave equation for VTI media, program are designed and the depth migration method are tested using synthetic model. The result shows that VTISSF is a stable algorithm, it generally gets a good result if the reflector dip is not very steep, while undermigration phenomena appeared in steep dips case; the VTIFFD algorithm bring us better result in steep dips with lower efficiency and frequency dispersion. For anisotropic prestack depth migration velocity analysis of VTI media, The paper discussed the basic hypothesis of VTI model in velocity analysis algorithm, basis of anisotropic prestack depth migration velocity analysis and travel time table calculation of VTI media in integral prestack depth migration. Then , analyzed the P-wave common imaging gather in the case of homogeneous velocity and vertically variable velocity . studied the residual correction in common imaging gather produced by media parameter error, analyzed the condition of flat event and correct depth in common imaging gather . In this case, the anisotropic model parameter vector is , is vertical velocity of a point at top surface, is vertical velocity gradient, and are anisotropic parameter. We can get vertical velocity gradient from seismic data; then the P-wave common imaging gather of VTI media whose velocity varies in vertical and horizontal direction, the relationship between media parameter and event residual time shift of common image gather are studied. We got the condition of flattening common imaging gather with correct depth. In this case the anisotropic model parameter vector is , is velocity gradient in horizontal direction. As a result, the vertical velocity grads can be decided uniquely, but horizontal velocity grads and anisotropic parameter can’t be distinguished if no priori information available, our method is to supply parameter by velocity scanning; then, as soon as is supplied we can get another four parameters of VTI media from seismic data. Based on above analysis, the paper discussed the feasibility of migration velocity analysis in vertically and horizontally varied VTI media, synthetic record of three models are used to test the velocity analysis method . Firstly, anisotropic velocity analysis test is done using a simple model with one block, then we used a model with multiple blocks, thirdly, we analyzed the anisotropic velocity using a part of Marmousi model. The model results show that this velocity analysis method is feasible and correct.