Abstarct: The main aim of Seismic reflection imaging is obtaining detailed information about geological structures to form a Seismic depth imaging. This image can be obtained by applying process depth migration or depth imaging on seismic data processing. Kirchhoff migration is a popular method for depth migration because of its flexibility and efficiency in practice. However, it fails in complex structures with varying velocities, which cause multi-pathing of seismic ray. To overcome this problem, Gaussian beam migration was introduced that is an extension of Kirchhoff migration. Like as Fourier finite-difference migration, it is especially compatible with lateral velocity variations and is capable to image steeply dipping reflectors. In this study, the Gaussian beam migration, Kirchhoff migration and Fourier finite-difference migration methods were used on a CRS stacked section. Results from the Gaussian beam migration had shown its ability in seismic imaging of complex structure. The migrated section from Gaussian beam migration revealed not only the steeply dipping events, but also the existing anticline which was not imaged by other methods. One can conclude to integrate CRS stacking with Gaussian beam migration for complex structures imaging. Generally, producing the wave propagation velocity model in the earth is a very complex process. In this study, stacking velocity model with Kinematic wavefield attributes is used as the initial velocity model. The purpose of this paper is studying different relation to prepare the migration velocity model and comparison of the migration sections with each other. Then the initial model of the best migration sections as the most practical initial velocity model will be introduced. In this study, three real seismic data with complex structures are considered. Gaussian beam migration method was applied to all data.