广义预测控制方法在网络遥操作机器人系统中的应用

Internet网络的时变时延及网络数据丢包严重影响了遥操作机器人系统的操作性能,甚至造成系统不稳定。为了解决这一问题,提出一种新的基于Internet的遥操作机器人系统控制结构。通过在主端对给定信息加入时间标签获得过去的系统回路时延,采用多元线性回归算法,预测下一时刻系统回路时延,然后在从端设计一个广义预测控制器控制远端机器人,从而改善时变时延对系统性能的影响。应用广义预测控制器产生的冗余控制信息,降低了网络数据丢包对系统的影响。最后根据预测控制稳定性定理,推导出系统的稳定性条件。仿真试验结果表明,该方法能有效解决时变时延以及网络数据丢包引起的性能下降问题。

基于介电泳的电极阵列电场仿真研究

介电泳方法被广泛地应用于微纳颗粒的分离和操纵中,实现介电泳操作的关键是设计满足所需电场分布的电极阵列。针对目前在微电极阵列设计中尚缺乏简单有效的电场解析方法的现状,提出一种基于格林公式的电极阵列电场的解析方法。首先介绍了传统介电泳和行波介电泳的概念和计算模型,分析了介电泳过程与电极上所施加的交变电压的频率和幅度的关系,然后在确立电极电势的边界条件的基础上,采用基于格林公式的电场解析方法,建立了非均匀电场的解析模型,得出不同条件下的电极阵列电场分布的仿真结果,最后利用FEMLAB有限元仿真软件对解析模型进行了对比仿真,验证了该解析模型的可行性。基于格林公式的电场解析求解方法能够有效地提高电极阵列设计中的针对性以及缩短电极设计的时间。

基于测角的自主移动机器人定位算法

给出了自主移动机器人定位的两种算法：解析算法和数值算法。解析法公式较以往的简洁。数值算法结合解析法和高斯－牛顿算法，不仅能避免因初值选取不合理而导致求解过程发散的问题，而且能提高运算精度和速度，通过对两种算法的计算机仿真，表明了解析算法具有运算速度快，而数值算法精度高的特点。其结果已用于自主移动机器人的研制中。

一种多移动机器人避碰规划方法

本文采用集中预规划方法 ,通过调整机器人的运动速度实现多机器人避碰 ,所提算法的基本思想为 :将机器人的运动路径分段 ,然后按避碰要求对机器人通过各段的时间进行约束 ,从而将避碰问题转化为高维线性空间的优化问题 ,并进一步将其转化为线性方程的求解 ,使问题具有明确的解析解 .由于该方法的复杂度较高 ,在实现过程中采用了多种方法降低复杂度 ,简化计算 .本文给出了该算法的基本思路 ,有关定理及证明 ,算法的化简方法 ,最后给出了实验结果及分析 .

两代竞争遗传算法及其应用研究

本文在分析简单遗传算法 (Simple Genetic Algorithm,SGA)的基础上 ,提出了一种新型结构的两代竞争遗传算法 ,并给出了算法演进的模式定理 .通过理论分析和对 TSP(TravelSalesman Problem,TSP)问题的应用研究 ,表明了该算法具有搜索效率高、鲁棒性强的特点

20世纪运动稳定性研究的三个重要结果

本文论述 2 0世纪运动稳定性理论研究的三个重要结果 :李雅普诺夫函数、谢聂稳定判据、卡利托洛夫定理 .这三个结果都是对一般的连续系统作出的 ,结论明确 ,简单实用 ,因而具有广泛的应用性

全方位六足步行机器人运动规划的相对运动算法

本文提出了步行机器人运动控制算法。该方法以相对运动学原理为基础,把机体的运动规划问题转化为腿的足端轨迹规划问题,从而使步行机器人运动控制问题得到大大简化.并应用该方法对全方位三角步态算法及稳定性进行分析求解.

基于弹性参数分析的碳酸盐岩储层预测方法研究

Prediction of Carbonate Reservoir Based on the Elastic Parameter Analysis Zhang Guangzhi (Solid Geophysics) Directed by Professor Liu Hong Abstract With the exploration and development of Puguang Oilfield, oil-gas exploration of carbonate rock in China has shown good prospects. Research on earthquake prediction methods for carbonate reservoir becomes the key of oil and gas exploration. Starting with analysis of geological characteristics of carbonate rock, prestack AVO inversion method, prestack elastic impedance inversion and parameter calculation method and seismic attribute extraction and optimization method were studied based on the analysis of rock physics in this work. First, variation characteristic and law of carbonate rock reservoir parameters were studied based on experimental data of rock physics, log data, analysis assay data, mud logging data and seismic data, so as to lay a foundation for the further reservoir identification and description. Then, the structure, type and propagation law of seismic wave field were analyzed through seismic forward modeling of the reservoir, and contact between information from log and geology data with elastic parameters, such as compressional wave and shear wave velocity and density were established, so as to provide a standard for reservoir identification and hydrocarbon detection using seismic reflection characteristics of the research area. Starting with the general concept of inverse problem, through analysis of Zoeppritz equation, three kinds of pre-stack inversion methods were derived and analyzed in detail, the AVO 3-parameter inversion based on Bayesian theory, the prestack AVO waveform inversion method and the simultaneous inversion method, based on the statistical hypothesis of inversion parameters and observation data and the Gauss distribution assumption of noise. The three methods were validated by model data and real data. Then, the elastic wave impedance inversion method of carbonate reservoir was investigated and the method of elastic parameter extraction from elastic impedance data was put forward. Based on the analysis of conventional methods of seismic attribute extraction and optimization, the time-frequency attributes and the wavelet attributes with time and amplitude feature were presented, and the prestack seismic attribute calculation method which can characterize the reservoir rock and fluid characteristic was presented. And the optimization of seismic attribute using the nonlinear KPCA method was also put forward. A series of seismic prediction technologies for carbonate reservoir were presented based on analysis of rock physics and seismic forward simulation technology. Practical application of these technologies was implemented in A oil field of Southern China and good effect has been achieved. Key words: carbonate rock; reservoir prediction; rock physics, prestack seismic inversion; seismic attribute

地幔压力下CaSiO3 Perovskite状态方程和弹性常数的第一性原理研究

In the past decade density functional theory (DFT) has made its way from a peripheral position in quantum chemistry to center. Of course the often excellent accuracy of the DFT based methods has provided the primary driving force of this development. This dissertation is devoted to the study of physical and chemical properties of planetary materials by first-principle calculation. The concerned properties include the geometry, elastic constants and anisotropy. In the first chapter, we give a systematic introduction to theoretical background and review its progress. Development of quantum chemistry promotes the establishment of DFT. Theorem of Hohenberg-Kohn is the fundament of DFT and is developed to Kohn-Sham equation, which can be used to perform real calculations. Now, new corrections and extensions, together with developed exchange-correlation, have made DFT more accurate and suitable for larger systems. In the second chapter, we focus on the calculational methods and technical aspects of DFT. Although it is important to develop methods and program, external package are still often used. At the end of this chapter, we briefly some widely used simulation package and the application of DFT. In the third chapter, we begin to focus on properties of real materials by first principles calculation. We study a kind of minerals named Ca perovskite, investigate its possible structure and anisotropy at Earth’s mental condition. By understanding and predicting geo-physically important materials properties at extreme conditions, we can get the most accurate information to interpret seismic data in the context of likely geophysical processes.

有源电磁扩散场克希霍夫积分偏移

As active electromagnetic method, field data of CSAMT method follow the equation of diffusion. Propagting in solid earth media, diffusion EM signal has strong attenuation and dispersion, otherwise seismic wave shows weak attenuation and dispersion, therefore the resolution power of CSAMT method is not better than seismic reflection method. However, there is consistence and similarity between EM signal and seismic wave in wave equation, we can apply Kirchhoff integral migration technique, a proven one in seismic method in time domain, to carry out seduo-seismic processing for CSAMT signal in frequency domain so that the attenuation and dispersion could be made compensated in some extent, and the resolution power and interpretation precision of active EM wave could be improved. Satisfying passive homogeneous Helmholtz quation, we proceed with Green theorem and combine the active inhomogenous Helmholtz quation, the Kirchhoff integral formula could be derived. Given practical problems, if we only consider the surface integral value, and assume that the intergral value in other interface is zero, combined with Green theorem in uniform half space, the expression could be simplified, and we can obtain frequency-domain Kirchhoff integral formula in surface, which is also called downward continuation of EM field in frequency domain. With image conditions and energy compensation considered, in order to get image conditions in time domain Fourier inverse transformation in frequency domain can be performed, so we can formulate the active Kirchhoff integral migration expression. At first, we construct relative stratified model, with different frequency series taken into account, then we change the distances between transmitter and reciever, the EM response can be obtained. Analyzing the EM properties, we can clarify near and far zone that can instruct us to carry out transmitter layout in practical application. Combined with field data surveyed in far zone, We perform Kirchhoff integral migration and compare the results with model to interpret. Secondly, with far field EM data, we apply TM mode to get EM response of given 2D model, then apply Kirchhoff integral migration on modelling data and interpret the results.

甘肃临夏盆地龙担剖面晚新生代磁性地层与古环境

Linxia Basin, situated in the northeast belt of the Tibetan Plateau, is a late Cenozoic depression basin bounded by the Tibetan Plateau and the Chinese Loess Plateau. The Cenozoic deposition, spanning over 30Ma, in which very abundant mammal fossils were discovered, is very suitable for study of uplift processes and geo-morphological evolution of the Tibetan Plateau. The Longdan section (35°31′31.6″N，103°29′0.6″E) is famous for the middle Miocene Platybelodon fauna and the late Miocene Hipparion fauna for a long time and is also one of the earliest known places for wooly rhino, which lies on the east slope of Longdan, a small village of township Nalesi in the south of the Dongxiang Autonomous County, Linxia Hui Nationallity Autonomous Prefecture. The Longdan mammal fauna was discovered at the base of the Early Pleistocene loess deposits at Dongxiang, where the lithology is different from the typical Wucheng Loess on the Chinese Loess Plateau. The rich fossils contain many new species and the major two layers of fossils are in the loess beds. Geologically the fossiliferous area is located in the central part of the Linxia Cenozoic sedimentary basin. Tectonically the Linxia Basin is an intermountain fault basin, bordered by the Leijishan major fault in the south and the north Qinling and Qilianshan major faults in the north. The section is 51.6m thick above the gravel layer, including the 1.6m Late Pleistocene Malan Loess on the top and the other loess-paleosol sequences in the middle of the section. The base of the section is the Jishi Formation, consisting of gravel layer of 13 ~ 17m thick. In this study, 972 bulk samples were collected with an interval of 5cm and other 401 orientied samples were taken with a magnetic compass. In the laboratory, the paleomagnetism, medium grain size, susceptibility, color, micromorphology, anisotropy of magnetic susceptibility were analyzed. From the stratigraphic analysis, the Longdan section from the top 0.3m to the bottom 51.6m, containing 5 normal polarities (N1-N5) and 5 reversal polarities (R1-R5). The paleomagnetic results show N3 is the Olduvai subchron in the middle of the Matuyama chron, and then the chronology of the Longdan mammal fauna is constructed along the section. The Matuyama-Gauss boundary is 45m and N5 enters Gauss chron. The Olduvai subchron with the age of 1.77 ~ 1.95Ma is found just in the upper fossiliferous level of Longdan mammal fauna. Taking the deposit rate of the section into account, the geological age of the upper fossiliferous level of Longdan mammal fauna is estimated to be about 1.9Ma. The lower fossiliferous level is just below the Reunion subchron and its age is estimated to be 2.25Ma. In addition, anisotropy of magnetic susceptibility of the loess-paleosol and other climatic indexes were used for discussing the late Cenozoic paleoenvironmental changes at Longdan, from which the Longdan area should have been an area of predominantly steppe the same as the Longdan mammal fauna.

边坡工程的极限分析上限法

Theory of limit analysis include upper bound theorem and lower bound theorem. To deal with slope stability analysis by limit analysis is to approximate the real solution from upper limit and lower limit. The most used method of limit analysis is upper bound theorem, therefore it is often applied to slope engineering in many cases. Although upper bound approach of limit analysis can keep away from vague constitutive relation and complex stress analyses, it also can obtain rigorous result. Assuming the critical surface is circular slip surface, two kinematically admissible velocity fields for perpendicular slice method and radial slice method can be established according to the limit analysis of upper bound theorem. By means of virtual work rate equation and strength reduction method, the upper-bound solution of limit analysis for homogeneous soil slope can be obtained. A log-spiral rotational failure mechanism for homogeneous slope is discussed from two different conditions which represent the position of shear crack passing the toe and below the toe. In the dissertition, the author also establishes a rotational failure mechanics with combination of different logarithmic spiral arcs. Furthermore, the calculation formula of upper bound solution for inhomogeneous soil slope stability problem can be deduced based on the upper bound approach of rigid elements. Through calculating the external work rate caused by soil nail, anti-slide pile, geotechnological grid and retaining wall, the upper bound solution of safety factor of soil nail structure slope, slip resistance of anti-slide pile, critical height of reinforced soil slope and active earth pressure of retaining wall can be obtained by upper bound limit analysis method. Taking accumulated body slope as subject investigated, with study on the limit analysis method to calculate slope safety factor, the kinematically admissible velocity fields of perpendicular slice method for slope with broken slip surface is proposed. Through calculating not only the energy dissipation rate produced in the broken slip surfaces and the vertical velocity discontinuity, but also the work rate produced by self-weight and external load, the upper bound solution of slope with broken slip surface is deduced. As a case study, the slope stability of the Sanmashan landslide in the area of the Three Gorges reservoir is analyzed. Based on the theory of limit analysis, the upper bound solution for rock slope with planar failure surface is obtained. By means of virtual work-rate equation, energy dissipation caused by dislocation of thin-layer and terrane can be calculated; furthermore, the formulas of safety factor for upper bound approach of limit analysis can be deduced. In the end, a new computational model of stability analysis for anchored rock slope is presented after taking into consideration the supporting effect of rock-bolts, the action of seismic force and fissure water pressure. By using the model, not only the external woke-rate done by self-weight, seismic force, fissure water pressure and anchorage force but also the internal energy dissipation produced in the slip surface and structural planes can be totally calculated. According to the condition of virtual work rate equation in limit state, the formula of safety factor for upper bound limit analysis can be deduced.

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

指数矩阵近似计算在地震波场延拓中的应用

Seismic wave field numerical modeling and seismic migration imaging based on wave equation have become useful and absolutely necessarily tools for imaging of complex geological objects. An important task for numerical modeling is to deal with the matrix exponential approximation in wave field extrapolation. For small value size matrix exponential, we can approximate the square root operator in exponential using different splitting algorithms. Splitting algorithms are usually used on the order or the dimension of one-way wave equation to reduce the complexity of the question. In this paper, we achieve approximate equation of 2-D Helmholtz operator inversion using multi-way splitting operation. Analysis on Gauss integral and coefficient of optimized partial fraction show that dispersion may accumulate by splitting algorithms for steep dipping imaging. High-order symplectic Pade approximation may deal with this problem, However, approximation of square root operator in exponential using splitting algorithm cannot solve dispersion problem during one-way wave field migration imaging. We try to implement exact approximation through eigenfunction expansion in matrix. Fast Fourier Transformation (FFT) method is selected because of its lowest computation. An 8-order Laplace matrix splitting is performed to achieve a assemblage of small matrixes using FFT method. Along with the introduction of Lie group and symplectic method into seismic wave-field extrapolation, accurate approximation of matrix exponential based on Lie group and symplectic method becomes the hot research field. To solve matrix exponential approximation problem, the Second-kind Coordinates (SKC) method and Generalized Polar Decompositions (GPD) method of Lie group are of choice. SKC method utilizes generalized Strang-splitting algorithm. While GPD method utilizes polar-type splitting and symmetric polar-type splitting algorithm. Comparing to Pade approximation, these two methods are less in computation, but they can both assure the Lie group structure. We think SKC and GPD methods are prospective and attractive in research and practice.

基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.